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.
_______________________________________________
ff3d-users mailing list
[email protected]
http://lists.nongnu.org/mailman/listinfo/ff3d-users
