verbosity=0;

load "medit"
load "msh3"
load "mshmet"
load "aniso"
load "tetgen"

mesh3 ThBackup = readmesh3("toto.0.mesh");
medit("Th",ThBackup);
func Pk6 = [P1,P1,P1,P1,P1,P1];
fespace Wh1Backup(ThBackup,P1);
fespace Wh6Backup(ThBackup,Pk6);
Wh6Backup [m11,m21,m22,m31,m32,m33];// for the metric of the adaptation
real[int] hmetric1(6*ThBackup.nv);
Wh1Backup TBackup=x+y+z;
real hmin=.5,hmax=1.,errh=1.,adaptratio=1.5;
hmetric1=mshmet(ThBackup,TBackup,normalization=1,aniso=1,nbregul=1,hmax=hmax,hmin=hmin,err=errh);
boundaniso(6,hmetric1,40);
m11[]=hmetric1;
savemesh(ThBackup,"ThBackup.mesh");
savesol("ThBackup.sol",ThBackup,[m11,m21,m22,m31,m32,m33]);
exec("~/mmg/build/bin/mmg3d_O3 ThBackup.mesh -out ThBackupo.mesh -sol ThBackup.sol -hgrad "+adaptratio+" -m 10000 -hmax "+hmax+" -hmin "+hmin+" -v 0");
ThBackup=readmesh3("ThBackupo.mesh");
exec("rm ThBackup.mesh ThBackupo.sol ThBackupo.mesh ThBackup.sol");
savemesh(ThBackup,"toto.1.mesh");
medit("Th",ThBackup);