int main(
int argc,
char**argv) {
size_t subdivide = (argc > 1) ? atoi(argv[1]) : 3;
Float field_error = (argc > 2) ? atof(argv[2]) : 1e-3;
bool do_verbose = true;
bool do_clean = true;
string cleanlist;
check_macro (communicator().size() == 1, "zalesak_dg_adapt: command may be used as mono-process only");
din >> ievent (t0, phi0_h);
while (
din >> ievent (tf, phi_h));
<< oevent (t0, ch, phi0_h, phi_h);
cleanlist +=
tmp+
"criterion-0.vtk";
string py_name =
tmp+
"adapt.py";
cleanlist += " "+py_name;
bool view_2d = phi_h.
get_geo().map_dimension() < 3;
py << "#!/usr/bin/env pvbatch --script=" << endl
<< "# This is a paraview script automatically generated by rheolef." << endl
<< endl
<< "from paraview.simple import *" << endl
<< "from paraview_rheolef import * # load rheolef specific functions" << endl
<< endl
<< "opt = { \\" << endl
<< " 'mark' : 'c', \\" << endl
<< " 'view_2d' : " << view_2d << ", \\" << endl
<< " 'geo_error' : " << geo_error << ", \\" << endl
<< " 'field_error' : " << field_error << ", \\" << endl
<< " 'subdivide' : " << subdivide << " \\" << endl
<< " }" << endl
<< endl
<<
"adapt_Pk_iso_P1 (paraview, \""<<
tmp<<
"criterion-0.vtk\", \""<<
tmp<<
"adapt.vtk\", opt)" << endl
<< endl
;
string prog = "pvbatch --force-offscreen-rendering ";
cleanlist +=
" "+
tmp+
"adapt.vtk";
idiststream in_vtk (
tmp+
"adapt.vtk");
iorheo::setbasename(in_vtk.is(),"adapt");
in_vtk >> ia_event (t, c_ha, phi0_ha, phi_ha);
<< oa_event (t0, phi0_ha)
<< oa_event (tf, phi_ha);
if (do_clean) {
}
if (do_verbose) {
cerr << "zalesak_dg_adapt: mesh size " << Xh.get_geo().size()
<<
" -> " << phi_ha.
get_geo().size() << endl;
}
}