Thanks Stephane, this helps a lot.
Is there any way to get a smoother tetrahedrization of the sinusoidal
surface?
I get many sharp angles instead of a smooth surface, see
http://www.cims.nyu.edu/~ddg5/wave.png
Thanks,
Denis Gueyffier
Courant Institute
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
End of ff3d-users Digest, Vol 34, Issue 2
*****************************************
_______________________________________________
ff3d-users mailing list
[email protected]
http://lists.nongnu.org/mailman/listinfo/ff3d-users
