37 int main(
int argc,
char *argv[])
40 const char *mesh_file =
"../data/star.mesh";
45 bool visualization = 1;
48 args.
AddOption(&mesh_file,
"-m",
"--mesh",
50 args.
AddOption(&ref_levels,
"-r",
"--refine",
51 "Number of times to refine the mesh uniformly, -1 for auto.");
53 "Finite element order (polynomial degree) >= 0.");
55 "One of the two DG penalty parameters, typically +1/-1."
56 " See the documentation of class DGDiffusionIntegrator.");
58 "One of the two DG penalty parameters, should be positive."
59 " Negative values are replaced with (order+1)^2.");
60 args.
AddOption(&visualization,
"-vis",
"--visualization",
"-no-vis",
62 "Enable or disable GLVis visualization.");
71 kappa = (order+1)*(order+1);
78 Mesh *mesh =
new Mesh(mesh_file, 1, 1);
88 ref_levels = (int)floor(log(50000./mesh->
GetNE())/log(2.)/
dim);
90 for (
int l = 0; l < ref_levels; l++)
104 cout <<
"Number of unknowns: " << fespace->
GetVSize() << endl;
135 #ifndef MFEM_USE_SUITESPARSE
142 PCG(A, M, *b, x, 1, 500, 1e-12, 0.0);
146 GMRES(A, M, *b, x, 1, 500, 10, 1e-12, 0.0);
151 umf_solver.
Control[UMFPACK_ORDERING] = UMFPACK_ORDERING_METIS;
153 umf_solver.
Mult(*b, x);
158 ofstream mesh_ofs(
"refined.mesh");
159 mesh_ofs.precision(8);
160 mesh->
Print(mesh_ofs);
161 ofstream sol_ofs(
"sol.gf");
162 sol_ofs.precision(8);
168 char vishost[] =
"localhost";
171 sol_sock.precision(8);
172 sol_sock <<
"solution\n" << *mesh << x << flush;
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.
Subclass constant coefficient.
int GetNE() const
Returns number of elements.
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 UniformRefinement(int i, const DSTable &, int *, int *, int *)
void SetCurvature(int order, bool discont=false, int space_dim=-1, int ordering=1)
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)
void PrintUsage(std::ostream &out) const
L2_FECollection DG_FECollection
double Control[UMFPACK_CONTROL]
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)
NURBSExtension * NURBSext
Optional NURBS mesh extension.
void PrintOptions(std::ostream &out) const
virtual void Mult(const Vector &b, Vector &x) const
Operator application: y=A(x).
int GMRES(const Operator &A, Vector &x, const Vector &b, Solver &M, int &max_iter, int m, double &tol, double atol, int printit)
GMRES method. (tolerances are squared)
virtual void SetOperator(const Operator &op)
Factorize the given Operator op which must be a SparseMatrix.