The p-Laplacian problem by the Newton method
int main(
int argc,
char**argv) {
string approx = (argc > 2) ? argv[2] : "P1";
Float p = (argc > 3) ? atof(argv[3]) : 1.5;
Float tol = (argc > 4) ? atof(argv[4]) : 1e5*eps;
size_t max_iter = (argc > 5) ? atoi(argv[5]) : 500;
derr <<
"# P-Laplacian problem by Newton:" << endl
<< "# geo = " << omega.name() << endl
<< "# approx = " << approx << endl
<< "# tol = " << tol << endl
<< "# max_iter = " << max_iter << endl;
dout << setprecision(numeric_limits<Float>::digits10)
}