Skip to content
Snippets Groups Projects
Commit 759fd4cb authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

*** empty log message ***

parent fffbd52f
No related branches found
No related tags found
No related merge requests found
...@@ -77,10 +77,10 @@ class gmshLinearSystemGmm : public gmshLinearSystem<scalar> { ...@@ -77,10 +77,10 @@ class gmshLinearSystemGmm : public gmshLinearSystem<scalar> {
{ {
// gmm::ilutp_precond<gmm::row_matrix<gmm::rsvector<scalar> > > P(*_a, 10,0.); // gmm::ilutp_precond<gmm::row_matrix<gmm::rsvector<scalar> > > P(*_a, 10,0.);
gmm::ildltt_precond<gmm::row_matrix<gmm::wsvector<scalar> > > P(*_a, 2, 1.e-10); gmm::ildltt_precond<gmm::row_matrix<gmm::wsvector<scalar> > > P(*_a, 2, 1.e-10);
gmm::iteration iter(_prec); // defines an iteration object, with a max residu of 1E-8 gmm::iteration iter(_prec);
iter.set_noisy(_noisy); iter.set_noisy(_noisy);
if(_gmres) gmm::gmres(*_a, *_x, *_b, P, 100, iter); // execute the GMRES algorithm if(_gmres) gmm::gmres(*_a, *_x, *_b, P, 100, iter);
else gmm::cg(*_a, *_x, *_b, P, iter); // execute the CG algorithm else gmm::cg(*_a, *_x, *_b, P, iter);
return 1; return 1;
} }
}; };
...@@ -107,4 +107,5 @@ public : ...@@ -107,4 +107,5 @@ public :
}; };
#endif #endif
#endif #endif
...@@ -48,6 +48,7 @@ void GMSH_FiniteElementPlugin::getInfos(char *author, char *copyright, ...@@ -48,6 +48,7 @@ void GMSH_FiniteElementPlugin::getInfos(char *author, char *copyright,
"Plugin(FiniteElement) solves simple PDEs\n" "Plugin(FiniteElement) solves simple PDEs\n"
"using the finite element method. This is only\n" "using the finite element method. This is only\n"
"intended as a demonstration!\n" "intended as a demonstration!\n"
"\n"
"Plugin(FiniteElement) is creates a new view.\n"); "Plugin(FiniteElement) is creates a new view.\n");
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment