The combustion problem by Keller continuation – post-treatment
int main(
int argc,
char**argv) {
string metric;
<< setprecision(numeric_limits<Float>::digits10)
<< "# metric " << metric << endl
<< "# s lambda umax det(mantissa,base,exp) |u| |grad(u)| |residue|"
<< endl;
for (
size_t n = 0; F.get(
din,xh); ++
n) {
problem::determinant_type det;
if (
n > 0 || metric !=
"spherical") det = F.update_derivative(xh);
const space& Xh = xh.second.get_space();
dout << F.parameter() <<
" " << xh.first
<< " " << xh.second.max_abs()
<< " " << det.mantissa
<< " " << det.base
<< " " << det.exponant
<<
" " << sqrt(
m(xh.second,xh.second))
<<
" " << sqrt(
a(xh.second,xh.second))
<< endl;
dot_xh = F.direction (xh);
F.refresh (F.parameter(), xh, dot_xh);
}
}