The elasticity problem for the embankment geometry with adaptive mesh
int main(
int argc,
char**argv) {
string approx = (argc > 2) ? argv[2] : "P1";
options.
err = (argc > 3) ? atof(argv[3]) : 5e-3;
size_t n_adapt = (argc > 4) ? atoi(argv[4]) : 5;
for (size_t i = 0; true; i++) {
if (i == n_adapt) break;
omega =
adapt(ch, options);
og << omega;
}
}