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

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";
continuation_option opts;
opts.ini_delta_parameter = 0.1;
opts.max_delta_parameter = 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
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
Float epsilon