Hello.
If I understand your problem: you think that the domain geometry is not
correctly described.
I think that the description is correct: you are using a fictitious domain
approximation. This means that the mesh does not fit the geometry. In fact
the tetrahedrization show all tetrahedra that *intersect* the domain. If you
want to see the effective boundary of the discretization use the instructions
mesh smesh = surface(Omega,M);
save(medit,"surface",smesh);
you will then be able to ensure that the geometry is correct.
I am sorry but I just noticed that there is a bug concerning the surface mesh
generation from the 1.0-pre9 release, so use your 1.0-pre7 to build this
mesh.
Concerning the size of the matrix, you can use the option
solve(salt) in Omega by M
memory(matrix=none)
{
which is a bit slower but consumes much less memory.
Best regards,
Stéphane.
_______________________________________________
ff3d-users mailing list
[email protected]
http://lists.nongnu.org/mailman/listinfo/ff3d-users
