[Getfem-users] Problem with mixed meshes

[Huu Phuoc BUI]

Dear Getfem-users,

I have a problem with computing the outer faces when using a mixed mesh 
(having tetras, hexa, prisms and pyramids).

The error is:

face_left = m.outer_faces_with_direction([0, 0., -1.0], 0.01)
  File "/usr/local/lib/python2.7/dist-packages/getfem/getfem.py", line 
1372, in outer_faces_with_direction
    return self.get("outer_faces_with_direction", v, angle, CVIDs)
  File "/usr/local/lib/python2.7/dist-packages/getfem/getfem.py", line 
1152, in get
    return getfem('mesh_get',self.id, *args)
RuntimeError: (Getfem::InterfaceError) -- Error in 
bgeot_geometric_trans.cc, line 381 const base_matrix& 
bgeot::geotrans_interpolation_context::B() const:
Non invertible matrix


Please let me know if you have any idea. Thank you.

Best regards,

Phuoc