Rheolef  7.1
an efficient C++ finite element environment
combustion_newton.cc
Go to the documentation of this file.
1 #include "rheolef.h"
26 using namespace rheolef;
27 using namespace std;
28 #include "combustion.h"
29 int main(int argc, char**argv) {
30  environment rheolef (argc,argv);
31  geo omega (argv[1]);
33  string approx = (argc > 2) ? argv[2] : "P1";
34  Float lambda = (argc > 3) ? atof(argv[3]) : 0.1;
35  Float tol = (argc > 4) ? atof(argv[4]) : eps;
36  size_t max_iter = (argc > 5) ? atoi(argv[5]) : 100;
37  combustion F (omega, approx);
39  field uh = F.initial ();
40  Float residue = tol;
41  size_t n_iter = max_iter;
42  damped_newton (F, uh, residue, n_iter, &derr);
43  F.put (dout, uh);
44  return (residue <= sqrt(tol)) ? 0 : 1;
45 }
field
see the field page for the full documentation
residue
field residue(Float p, const field &uh)
Definition: p_laplacian_post.cc:35
combustion::set_parameter
void set_parameter(float_type lambda1)
Definition: combustion.h:35
combustion.h
The combustion problem: class header for the Newton method.
main
int main(int argc, char **argv)
Definition: combustion_newton.cc:29
combustion
Definition: combustion.h:25
rheolef.h
rheolef - reference manual
rheolef::environment
see the environment page for the full documentation
Definition: environment.h:115
rheolef
This file is part of Rheolef.
Definition: compiler_eigen.h:37
rheolef::derr
odiststream derr(cerr)
see the diststream page for the full documentation
Definition: diststream.h:436
Float
see the Float page for the full documentation
rheolef::damped_newton
int damped_newton(const Problem &P, const Preconditioner &T, Field &u, Real &tol, Size &max_iter, odiststream *p_derr=0)
see the damped_newton page for the full documentation
Definition: damped-newton-generic.h:29
rheolef::dout
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:430
epsilon
Float epsilon
Definition: transmission_error.cc:25
rheolef::std
Definition: vec_expr_v2.h:402
combustion::put
odiststream & put(odiststream &os, const field &uh) const
Definition: combustion1.icc:43
geo
see the geo page for the full documentation
combustion::initial
field initial(std::string restart="")
Definition: combustion1.icc:35
lambda
Definition: yield_slip_circle.h:34