62 int main(
int argc,
char *argv[])
66 MPI_Init(&argc, &argv);
67 MPI_Comm_size(MPI_COMM_WORLD, &num_procs);
68 MPI_Comm_rank(MPI_COMM_WORLD, &myid);
71 const char *mesh_file =
"../data/star.mesh";
73 bool static_cond =
false;
75 const char *device_config =
"cpu";
76 bool visualization =
true;
79 args.
AddOption(&mesh_file,
"-m",
"--mesh",
82 "Finite element order (polynomial degree) or -1 for"
83 " isoparametric space.");
84 args.
AddOption(&static_cond,
"-sc",
"--static-condensation",
"-no-sc",
85 "--no-static-condensation",
"Enable static condensation.");
86 args.
AddOption(&pa,
"-pa",
"--partial-assembly",
"-no-pa",
87 "--no-partial-assembly",
"Enable Partial Assembly.");
88 args.
AddOption(&device_config,
"-d",
"--device",
89 "Device configuration string, see Device::Configure().");
90 args.
AddOption(&visualization,
"-vis",
"--visualization",
"-no-vis",
92 "Enable or disable GLVis visualization.");
110 Device device(device_config);
111 if (myid == 0) { device.
Print(); }
116 Mesh mesh(mesh_file, 1, 1);
125 (int)floor(log(10000./mesh.
GetNE())/log(2.)/
dim);
126 for (
int l = 0; l < ref_levels; l++)
135 ParMesh pmesh(MPI_COMM_WORLD, mesh);
138 int par_ref_levels = 2;
139 for (
int l = 0; l < par_ref_levels; l++)
161 cout <<
"Using isoparametric FEs: " << fec->
Name() << endl;
173 cout <<
"Number of finite element unknowns: " << size << endl;
251 ostringstream mesh_name, sol_name;
252 mesh_name <<
"mesh." << setfill(
'0') << setw(6) << myid;
253 sol_name <<
"sol." << setfill(
'0') << setw(6) << myid;
255 ofstream mesh_ofs(mesh_name.str().c_str());
256 mesh_ofs.precision(8);
257 pmesh.
Print(mesh_ofs);
259 ofstream sol_ofs(sol_name.str().c_str());
260 sol_ofs.precision(8);
270 sol_sock <<
"parallel " << num_procs <<
" " << myid <<
"\n";
271 sol_sock.precision(8);
272 sol_sock <<
"solution\n" << pmesh << x << flush;
int Size() const
Return the logical size of the array.
Class for domain integration L(v) := (f, v)
virtual void GetEssentialTrueDofs(const Array< int > &bdr_attr_is_ess, Array< int > &ess_tdof_list, int component=-1)
Conjugate gradient method.
A coefficient that is constant across space and time.
virtual void Mult(const Vector &b, Vector &x) const
Operator application: y=A(x).
Pointer to an Operator of a specified type.
int GetNE() const
Returns number of elements.
virtual void Save(std::ostream &out) const
void Print(std::ostream &out=mfem::out)
Print the configuration of the MFEM virtual device object.
Abstract parallel finite element space.
bool UsesTensorBasis(const FiniteElementSpace &fes)
int main(int argc, char *argv[])
The BoomerAMG solver in hypre.
void SetPrintLevel(int print_lvl)
void Parse()
Parse the command-line options. Note that this function expects all the options provided through the ...
Jacobi smoothing for a given bilinear form (no matrix necessary).
void UniformRefinement(int i, const DSTable &, int *, int *, int *)
void SetMaxIter(int max_it)
T Max() const
Find the maximal element in the array, using the comparison operator < for class T.
HYPRE_Int GlobalTrueVSize() const
virtual void Print(std::ostream &out=mfem::out) const
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.
void SetRelTol(double rtol)
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 Clear()
Clear the contents of the Mesh.
virtual void SetOperator(const Operator &op)
Also calls SetOperator for the preconditioner if there is one.
void GetNodes(Vector &node_coord) const
virtual void SetPreconditioner(Solver &pr)
This should be called before SetOperator.
Arbitrary order H1-conforming (continuous) finite elements.
Class for parallel grid function.
The MFEM Device class abstracts hardware devices such as GPUs, as well as programming models such as ...
Class for parallel meshes.
bool Good() const
Return true if the command line options were parsed successfully.