Rheolef  7.2
an efficient C++ finite element environment
combustion_continuation.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  cin >> noverbose;
32  geo omega (argv[1]);
33  string approx = (argc > 2) ? argv[2] : "P1";
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
Definition: environment.h:121
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:467
odiststream derr(cerr)
see the diststream page for the full documentation
Definition: diststream.h:473
The combustion problem: class header for the Newton method.
int main(int argc, char **argv)
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
Definition: continuation.h:303
rheolef - reference manual
odiststream & put(odiststream &os, const field &uh) const
Definition: combustion1.icc:43
field initial(std::string restart="")
Definition: combustion1.icc:35
see the continuation_option page for the full documentation
Float epsilon