Re: [ff3d-users] Povray geometry

[Stephane Del Pino]

Dear Denis.

It is true that ff3d's POVRay interpreter is far to be complete, but in your 
case, you can define an analytic domain. Here is an example :

--------------------------
vector a = (-1,-0.5,-1);
vector b = (1,1,1);
vector n = (20,20,20);

mesh m = structured(n,a,b);

double pi = acos(-1);
domain D = domain(y<0.1*sin(4*pi*x));

mesh s = surface(D,m);
solve(u) in D by m
{
 pde(u)
   -div(grad(u))=1;
   u = 0 on s;
   u = 1 on m xmin;
}

mesh t = tetrahedrize(D,m);
save(vtk,"u",u,t);
--------------------------

Best regards,
Stéphane.