The p-Laplacian problem by the damped 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]) : eps;
size_t max_iter = (argc > 5) ? atoi(argv[5]) : 500;
derr <<
"# P-Laplacian problem by damped Newton:" << endl
<< "# geo = " << omega.name() << endl
<< "# approx = " << approx << endl
<<
"# p = " <<
p << endl;
}