47 #ifdef MFEM_USE_SIMMETRIX
54 #include <apfConvert.h>
59 #error This example requires that MFEM is built with MFEM_USE_PUMI=YES
65 int main(
int argc,
char *argv[])
69 MPI_Init(&argc, &argv);
70 MPI_Comm_size(MPI_COMM_WORLD, &num_procs);
71 MPI_Comm_rank(MPI_COMM_WORLD, &myid);
74 const char *mesh_file =
"../../data/pumi/serial/Kova.smb";
75 #ifdef MFEM_USE_SIMMETRIX
76 const char *model_file =
"../../data/pumi/geom/Kova.x_t";
78 const char *model_file =
"../../data/pumi/geom/Kova.dmg";
81 bool static_cond =
false;
82 bool visualization = 1;
86 args.
AddOption(&mesh_file,
"-m",
"--mesh",
89 "Finite element order (polynomial degree) or -1 for"
90 " isoparametric space.");
91 args.
AddOption(&static_cond,
"-sc",
"--static-condensation",
"-no-sc",
92 "--no-static-condensation",
"Enable static condensation.");
93 args.
AddOption(&visualization,
"-vis",
"--visualization",
"-no-vis",
95 "Enable or disable GLVis visualization.");
96 args.
AddOption(&model_file,
"-p",
"--parasolid",
97 "Parasolid model to use.");
98 args.
AddOption(&geom_order,
"-go",
"--geometry_order",
99 "Geometric order of the model");
117 #ifdef MFEM_USE_SIMMETRIX
118 Sim_readLicenseFile(0);
124 apf::Mesh2* pumi_mesh;
125 pumi_mesh = apf::loadMdsMesh(model_file, mesh_file);
130 crv::BezierCurver bc(pumi_mesh, geom_order, 2);
148 (int)floor(log(50000./mesh->
GetNE())/log(2.)/
dim);
149 for (
int l = 0; l < ref_levels; l++)
166 cout <<
"Using isoparametric FEs: " << fec->
Name() << endl;
173 cout <<
"Number of finite element unknowns: "
219 cout <<
"Size of linear system: " << A.
Height() << endl;
221 #ifndef MFEM_USE_SUITESPARSE
225 PCG(A, M, B, X, 1, 200, 1e-12, 0.0);
229 umf_solver.
Control[UMFPACK_ORDERING] = UMFPACK_ORDERING_METIS;
231 umf_solver.
Mult(B, X);
239 ofstream mesh_ofs(
"refined.mesh");
240 mesh_ofs.precision(8);
241 mesh->
Print(mesh_ofs);
242 ofstream sol_ofs(
"sol.gf");
243 sol_ofs.precision(8);
252 sol_sock.precision(8);
253 sol_sock <<
"solution\n" << *mesh << x << flush;
260 if (order > 0) {
delete fec; }
263 pumi_mesh->destroyNative();
264 apf::destroyMesh(pumi_mesh);
266 #ifdef MFEM_USE_SIMMETRIX
268 Sim_unregisterAllKeys();
int Size() const
Return the logical size of the array.
Class for domain integration L(v) := (f, v)
virtual void Print(std::ostream &out=mfem::out) const
Class for grid function - Vector with associated FE space.
A coefficient that is constant across space and time.
virtual void GetEssentialTrueDofs(const Array< int > &bdr_attr_is_ess, Array< int > &ess_tdof_list, int component=-1)
Get a list of essential true dofs, ess_tdof_list, corresponding to the boundary attributes marked in ...
int GetNE() const
Returns number of elements.
Base class for PUMI meshes.
int main(int argc, char *argv[])
Data type for Gauss-Seidel smoother of sparse matrix.
Direct sparse solver using UMFPACK.
virtual void Save(std::ostream &out) const
Save the GridFunction to an output stream.
void Parse()
Parse the command-line options. Note that this function expects all the options provided through the ...
int Height() const
Get the height (size of output) of the Operator. Synonym with NumRows().
void UniformRefinement(int i, const DSTable &, int *, int *, int *)
T Max() const
Find the maximal element in the array, using the comparison operator < for class T.
void PCG(const Operator &A, Solver &B, const Vector &b, Vector &x, int print_iter, int max_num_iter, double RTOLERANCE, double ATOLERANCE)
Preconditioned conjugate gradient method. (tolerances are squared)
virtual int GetTrueVSize() const
Return the number of vector true (conforming) dofs.
void PrintUsage(std::ostream &out) const
Print the usage message.
Array< int > bdr_attributes
A list of all unique boundary attributes used by the Mesh.
double Control[UMFPACK_CONTROL]
Class FiniteElementSpace - responsible for providing FEM view of the mesh, mainly managing the set of...
Collection of finite elements from the same family in multiple dimensions. This class is used to matc...
void AddOption(bool *var, const char *enable_short_name, const char *enable_long_name, const char *disable_short_name, const char *disable_long_name, const char *description, bool required=false)
Add a boolean option and set 'var' to receive the value. Enable/disable tags are used to set the bool...
virtual const char * Name() const
void PrintOptions(std::ostream &out) const
Print the options.
void GetNodes(Vector &node_coord) const
Arbitrary order H1-conforming (continuous) finite elements.
virtual void Mult(const Vector &b, Vector &x) const
Operator application: y=A(x).
virtual void SetOperator(const Operator &op)
Factorize the given Operator op which must be a SparseMatrix.
bool Good() const
Return true if the command line options were parsed successfully.