Rheolef  7.2
an efficient C++ finite element environment
combustion_continuation.cc

The combustion problem by continuation.

The combustion problem by continuation

#include "rheolef.h"
using namespace rheolef;
using namespace std;
#include "combustion.h"
int main(int argc, char**argv) {
environment rheolef (argc, argv);
cin >> noverbose;
geo omega (argv[1]);
string approx = (argc > 2) ? argv[2] : "P1";
opts.ini_delta_parameter = 0.1;
opts.min_delta_parameter = 1e-7;
opts.tol = eps;
derr << setprecision(numeric_limits<Float>::digits10)
<< "# continuation in lambda:" << endl
<< "# geo = " << omega.name() << endl
<< "# approx = " << approx << endl
<< "# dlambda_ini = " << opts.ini_delta_parameter << endl
<< "# dlambda_min = " << opts.min_delta_parameter << endl
<< "# dlambda_max = " << opts.max_delta_parameter << endl
<< "# tol = " << opts.tol << endl;
combustion F (omega, approx);
field uh = F.initial();
F.put (dout, uh);
continuation (F, uh, &dout, &derr, opts);
}
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
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
STL namespace.
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