@@ -310,11 +310,13 @@ def fix_surface_conflicts(dimension, surfaces, pname_map, name_priority):
310310 for lid , lvertexes in priority_vertexes .items ():
311311 if lvertexes .intersection (nset ):
312312 hpname = pname_map [lid ][1 ]
313- errors += "overlapping elements between priority_name %s and boundary of higher priority %s\n " % (n , hpname )
313+ #errors += "overlapping elements between priority_name %s and boundary of higher priority %s\n" % (n, hpname)
314+ errors += 'WARNING: boundaries "%s" and "%s" are touching\n ' % (n , hpname )
314315 priority_vertexes [nid ] = nset
315316 all_vertexes |= nset
316317 if errors :
317- raise RuntimeError (errors )
318+ errors += "WARNING: this may cause issues when the boundaries are solving the same equations on the same regions\n "
319+ print (errors )
318320
319321 new_surfaces = []
320322 removed_surfaces = set ([])
@@ -338,9 +340,9 @@ def fix_surface_conflicts(dimension, surfaces, pname_map, name_priority):
338340 new_surfaces .extend (local_new_elements )
339341 all_vertexes |= local_vertexes
340342 if removed > 0 :
341- print (" removed %d/%d elements from %s for overlap with other boundary" % (removed , removed + kept , pname_map [phys_id ][1 ]))
343+ print ('INFO: removed %d/%d elements from generated surface "%s" for overlap with other boundary' % (removed , removed + kept , pname_map [phys_id ][1 ]))
342344 if kept == 0 :
343- print (" generated surface %s removed for 0 elements" % pname_map [phys_id ][1 ])
345+ print ('INFO: generated surface "%s" removed for 0 elements' % pname_map [phys_id ][1 ])
344346 removed_surfaces .add (phys_id )
345347 return new_surfaces , removed_surfaces
346348
0 commit comments