/usr/share/doc/freefem++/examples/examples++-3d/Laplace-Adapt-3d.edp is in freefem++-doc 3.19.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 | load "msh3" load "tetgen" load "mshmet" load "medit"
//build initial mesh
int nn = 6;
int[int] l1111=[1,1,1,1],l01=[0,1],l11=[1,1];// label numbering
mesh3 Th3=buildlayers(square(nn,nn,region=0,label=l1111),
nn, zbound=[0,1], labelmid=l11, labelup = l01, labeldown = l01);
Th3 = trunc(Th3,(x<0.5) | (y < 0.5) | (z < 0.5) ,label=1);// remove the $]0.5,1[^3 cube$
//end of build initial mesh
fespace Vh(Th3,P1);
Vh u,v,usol,h;
macro Grad(u) [dx(u),dy(u),dz(u)] // EOM
problem Poisson(u,v,solver=CG) = int3d(Th3)( Grad(u)'*Grad(v) ) // ') for emacs
-int3d(Th3)( 1*v ) + on(1,u=0);
real errm=1e-2;// level of error
for(int ii=0; ii<5; ii++)
{
Poisson;
cout <<" u min, max = " << u[].min << " "<< u[].max << endl;
h=0. ;// for resizing h
h[]=mshmet(Th3,u,normalization=1,aniso=0,nbregul=1,hmin=1e-3,hmax=0.3,err=errm);//loptions=MSHloptions,doptions=MSHdoptions);
cout <<" h min, max = " << h[].min << " "<< h[].max << " " << h[].n << " " << Th3.nv << endl;
plot(u,wait=1);
errm*= 0.8;// change the level of error
cout << " Th3" << Th3.nv < " " << Th3.nt << endl;
Th3=tetgreconstruction(Th3,switch="raAQ",sizeofvolume=h*h*h/6.);
medit("U-adap-iso-"+ii,Th3,u,wait=1);
}
|