Re: [ff3d-users] Inner wall definitions

[Stephane Del Pino]

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.