Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
combustion_continuation.cc
Go to the documentation of this file.
1
25#include "rheolef.h"
26using namespace rheolef;
27using namespace std;
28#include "combustion.h"
29int main(int argc, char**argv) {
30 environment rheolef (argc, argv);
31 cin >> noverbose;
32 geo omega (argv[1]);
33 string approx = (argc > 2) ? argv[2] : "P1";
34 Float eps = numeric_limits<Float>::epsilon();
36 opts.ini_delta_parameter = 0.1;
37 opts.max_delta_parameter = 1;
38 opts.min_delta_parameter = 1e-7;
39 opts.tol = eps;
40 derr << setprecision(numeric_limits<Float>::digits10)
41 << "# continuation in lambda:" << endl
42 << "# geo = " << omega.name() << endl
43 << "# approx = " << approx << endl
44 << "# dlambda_ini = " << opts.ini_delta_parameter << endl
45 << "# dlambda_min = " << opts.min_delta_parameter << endl
46 << "# dlambda_max = " << opts.max_delta_parameter << endl
47 << "# tol = " << opts.tol << endl;
48 combustion F (omega, approx);
49 field uh = F.initial();
50 F.put (dout, uh);
51 continuation (F, uh, &dout, &derr, opts);
52}
see the Float page for the full documentation
see the field page for the full documentation
see the geo page for the full documentation
see the environment page for the full documentation
The combustion problem: class header for the Newton method.
int main()
Definition field2bb.cc:58
This file is part of Rheolef.
void continuation(Problem &F, typename Problem::value_type &uh, odiststream *p_out, odiststream *p_err, const continuation_option &opts=continuation_option())
see the continuation page for the full documentation
STL namespace.
rheolef - reference manual
odiststream & put(odiststream &os, const field &uh) const
field initial(std::string restart="")
see the continuation_option page for the full documentation