The combustion problem: class body for the Newton method
:
lambda(0), Xh(),
m(), pm(), a1(), pa1(), event(
"lambda",
"u") {
if (approx !=
"")
reset (omega, approx);
}
}
if (restart ==
"")
return field (
Xh, 0);
idiststream in (restart);
derr <<
"# restart from lambda=" <<
lambda << endl;
return xh0;
}
}
if (!is) return is;
if (
Xh.name() ==
"")
reset (uh.get_geo(), uh.get_approx());
if (uh.b().dis_size() == 0) {
std::copy (uh.begin_dof(), uh.end_dof(),
tmp.begin_dof());
}
return is;
}