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]);
32
Float
eps =
numeric_limits<Float>::epsilon
();
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);
38
F.
set_parameter
(
lambda
);
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:391
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