56 int main(
int argc,
char *argv[])
60 MPI_Init(&argc, &argv);
61 MPI_Comm_size(MPI_COMM_WORLD, &num_procs);
62 MPI_Comm_rank(MPI_COMM_WORLD, &myid);
65 const char *mesh_file =
"../data/star.mesh";
66 int ser_ref_levels = 2;
67 int par_ref_levels = 1;
71 bool slu_solver =
false;
72 bool sp_solver =
false;
73 bool visualization = 1;
76 args.
AddOption(&mesh_file,
"-m",
"--mesh",
78 args.
AddOption(&ser_ref_levels,
"-rs",
"--refine-serial",
79 "Number of times to refine the mesh uniformly in serial.");
80 args.
AddOption(&par_ref_levels,
"-rp",
"--refine-parallel",
81 "Number of times to refine the mesh uniformly in parallel.");
83 "Finite element order (polynomial degree) or -1 for"
84 " isoparametric space.");
86 "Number of desired eigenmodes.");
88 "Random seed used to initialize LOBPCG.");
89 #ifdef MFEM_USE_SUPERLU
90 args.
AddOption(&slu_solver,
"-slu",
"--superlu",
"-no-slu",
91 "--no-superlu",
"Use the SuperLU Solver.");
93 #ifdef MFEM_USE_STRUMPACK
94 args.
AddOption(&sp_solver,
"-sp",
"--strumpack",
"-no-sp",
95 "--no-strumpack",
"Use the STRUMPACK Solver.");
97 args.
AddOption(&visualization,
"-vis",
"--visualization",
"-no-vis",
99 "Enable or disable GLVis visualization.");
101 if (slu_solver && sp_solver)
104 cout <<
"WARNING: Both SuperLU and STRUMPACK have been selected,"
105 <<
" please choose either one." << endl
106 <<
" Defaulting to SuperLU." << endl;
131 Mesh *mesh =
new Mesh(mesh_file, 1, 1);
137 for (
int lev = 0; lev < ser_ref_levels; lev++)
148 for (
int lev = 0; lev < par_ref_levels; lev++)
173 cout <<
"Number of unknowns: " << size << endl;
214 #if defined(MFEM_USE_SUPERLU) || defined(MFEM_USE_STRUMPACK)
216 #ifdef MFEM_USE_SUPERLU
222 #ifdef MFEM_USE_STRUMPACK
237 if (!slu_solver && !sp_solver)
245 #ifdef MFEM_USE_SUPERLU
256 #ifdef MFEM_USE_STRUMPACK
295 ostringstream mesh_name, mode_name;
296 mesh_name <<
"mesh." << setfill(
'0') << setw(6) << myid;
298 ofstream mesh_ofs(mesh_name.str().c_str());
299 mesh_ofs.precision(8);
300 pmesh->
Print(mesh_ofs);
302 for (
int i=0; i<nev; i++)
307 mode_name <<
"mode_" << setfill(
'0') << setw(2) << i <<
"."
308 << setfill(
'0') << setw(6) << myid;
310 ofstream mode_ofs(mode_name.str().c_str());
311 mode_ofs.precision(8);
320 char vishost[] =
"localhost";
323 mode_sock.precision(8);
325 for (
int i=0; i<nev; i++)
329 cout <<
"Eigenmode " << i+1 <<
'/' << nev
330 <<
", Lambda = " << eigenvalues[i] << endl;
336 mode_sock <<
"parallel " << num_procs <<
" " << myid <<
"\n"
337 <<
"solution\n" << *pmesh << x << flush
338 <<
"window_title 'Eigenmode " << i+1 <<
'/' << nev
339 <<
", Lambda = " << eigenvalues[i] <<
"'" << endl;
344 cout <<
"press (q)uit or (c)ontinue --> " << flush;
347 MPI_Bcast(&c, 1, MPI_CHAR, 0, MPI_COMM_WORLD);
362 #if defined(MFEM_USE_SUPERLU) || defined(MFEM_USE_STRUMPACK)
int Size() const
Logical size of the array.
HypreParVector & GetEigenvector(unsigned int i)
Extract a single eigenvector.
Subclass constant coefficient.
void SetFromCommandLine()
virtual void Save(std::ostream &out) const
void SetPreconditioner(Solver &precond)
void SetMassMatrix(Operator &M)
Abstract parallel finite element space.
void SetPrintLevel(int logging)
int main(int argc, char *argv[])
void SetSymmetricPattern(bool sym)
void SetOperator(const Operator &op)
Set/update the solver for the given operator.
The BoomerAMG solver in hypre.
void SetPrintFactorStatistics(bool print_stat)
void SetMaxIter(int max_iter)
void SetColumnPermutation(superlu::ColPerm col_perm)
void UniformRefinement(int i, const DSTable &, int *, int *, int *)
void SetPrintLevel(int print_level)
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 SetPrintSolveStatistics(bool print_stat)
void SetPrintStatistics(bool print_stat)
void PrintUsage(std::ostream &out) const
void SetReorderingStrategy(strumpack::ReorderingStrategy method)
Array< int > bdr_attributes
A list of all unique boundary attributes used by the Mesh.
void GetEigenvalues(Array< double > &eigenvalues)
Collect the converged eigenvalues.
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)
void SetSize(int nsize)
Change logical size of the array, keep existing entries.
void SetRandomSeed(int s)
void SetKrylovSolver(strumpack::KrylovSolver method)
void SetOperator(Operator &A)
void PrintOptions(std::ostream &out) const
void SetOperator(const Operator &op)
Set/update the solver for the given operator.
void GetNodes(Vector &node_coord) const
void SetPrecondUsageMode(int pcg_mode)
Arbitrary order H1-conforming (continuous) finite elements.
Class for parallel grid function.
Wrapper for hypre's ParCSR matrix class.
void SetNumModes(int num_eigs)
Class for parallel meshes.
void Solve()
Solve the eigenproblem.