fix issue with IBM

Author: arciyer123

src/common/m_constants.fpp

@@ -39,6 +39,6 @@ module m_constants
     real(wp), parameter :: threshold_vector_zero = 1e-10 !< Threshold to treat the component of a vector to be zero
     real(wp), parameter :: threshold_edge_zero = 1e-10 !< Threshold to treat two edges to be overlapped
     real(wp), parameter :: threshold_bary = 1e-1 !< Threshold to interpolate a barycentric facet
-    real(wp), parameter :: initial_distance_buffer = 1e-12 !< Initialized levelset distance for the shortest path pair algorithm
+    real(wp), parameter :: initial_distance_buffer = 1e12 !< Initialized levelset distance for the shortest path pair algorithm
 
-end module m_constants
+end module m_constants

src/pre_process/m_model.fpp

@@ -1055,8 +1055,8 @@ contains
         t_vec3 :: dist_buffer !< Distance between the cell center and the vertices
         integer :: i, j, tri_idx !< Iterator
 
-        dist_min = 1e-12
-        dist_min_normal = 1e-12
+        dist_min = 1e12
+        dist_min_normal = 1e12
         distance = 0._wp
 
         tri_idx = 0
@@ -1221,8 +1221,8 @@ contains
         cross(1) = AB(2)*AC(3) - AB(3)*AC(2)
         cross(2) = AB(3)*AC(1) - AB(1)*AC(3)
         cross(3) = AB(1)*AC(2) - AB(2)*AC(1)
-        tri_area = 0.5d0*sqrt(cross(1)**2 + cross(2)**2 + cross(3)**2)
+        tri_area = 0.5_wp*sqrt(cross(1)**2 + cross(2)**2 + cross(3)**2)
 
     end function f_tri_area
 
-end module m_model
+end module m_model

toolchain/mfc/test/test.py

@@ -53,14 +53,14 @@ def __filter(cases_) -> typing.List[TestCase]:
                 cases.remove(case)
 
     for case in cases[:]:
-        if case.ppn > 1 and not ARG("mpi") or 'IBM' in case.trace:
+        if case.ppn > 1 and not ARG("mpi"):
             cases.remove(case)
             skipped_cases.append(case)
 
     for case in cases[:]:
         if ARG("single"):
             skip = ['low_Mach', 'Hypoelasticity', 'teno', 'Chemistry', 'Phase Change model 6'
-            ,'Axisymmetric', 'Transducer', 'Transducer Array', 'Cylindrical']
+            ,'Axisymmetric', 'Transducer', 'Transducer Array', 'Cylindrical', 'Example']
             if any(label in case.trace for label in skip):
                 cases.remove(case)
 
@@ -248,8 +248,8 @@ def handle_case(case: TestCase, devices: typing.Set[int]):
     global errors
 
     nAttempts = 0
-    if ARG('single') and ARG('test_all'):
-        max_attempts = max(ARG('max_attempts'), 20)
+    if ARG('single'):
+        max_attempts = max(ARG('max_attempts'), 3)
     else:
         max_attempts = ARG('max_attempts')