using namespace std;
int main(
int argc,
char**argv) {
string approx = (argc > 2) ? argv[2] : "P1";
string metric = (argc > 3) ? argv[3] : "orthogonal";
continuation_option opts;
opts.ini_delta_parameter = 0.1;
opts.max_delta_parameter = 0.5;
opts.min_delta_parameter = 1e-10;
opts.tol = eps;
derr << setprecision(numeric_limits<Float>::digits10)
<< "# continuation in s:" << endl
<< "# geo = " << omega.name() << endl
<< "# approx = " << approx << endl
<< "# metric = " << metric << endl
<< "# ds_init = " << opts.ini_delta_parameter << endl
<< "# ds_min = " << opts.min_delta_parameter << endl
<< "# ds_max = " << opts.max_delta_parameter << endl
<< "# tol = " << opts.tol << endl;
dout << catchmark(
"metric") << metric << endl;
keller<combustion> F (
combustion(omega,approx), metric);
keller<combustion>::value_type xh = F.initial();
}
see the Float page for the full documentation
see the geo page for the full documentation
idiststream din(cin)
see the diststream page for the full documentation
odiststream dout(cout)
see the diststream page for the full documentation
odiststream derr(cerr)
see the diststream page for the full documentation
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
rheolef - reference manual