Greetings, I am trying to use TetraInside in a user solver, but it fails if the element is very small. I work with MEMS so element are 1e-7 to 1e-6 across. The line that appears to cause the failure is in TetraInside subroutine in ElementDescription:
Line 4530:
I assume this is to check for if the tetrahedron is degenerate, but it also causes my solver to not work. Can the check be made relative to the element's size?
Attached is my files, including a some mesh (not used) to avoid a segmentation fault.The four node coordinates used with TetraInside (nx, ny, nz) are written into the user solver, but the coordinate to be checked can be changed in the sif file.
Thank you,