Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
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";
Float eps = numeric_limits<Float>::epsilon();
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
The combustion problem: class header for the Newton method.
int main()
Definition field2bb.cc:58
This file is part of Rheolef.
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