# Face with no cross edge

**URL:** <https://community.geodynamics.org/t/face-with-no-cross-edge/2745>\
**Category:** PyLith\
**Created:** [December 21, 2022, 11:14am UTC](https://community.geodynamics.org/t/face-with-no-cross-edge/2745 "2022-12-21T11:14:31Z")\
**Posts on this page:** 4\
**Page:** 1

<div class="post-metadata">

**Author:** ![VivenSharma](https://yyz2.discourse-cdn.com/flex036/user_avatar/community.geodynamics.org/vivensharma/32/933_2.png) [@VivenSharma](https://community.geodynamics.org/u/VivenSharma)\
**Post date:** [December 21, 2022, 11:14am UTC](https://community.geodynamics.org/t/face-with-no-cross-edge/2745/1 "2022-12-21T11:14:31Z")

</div>

HI Brad,

I am getting an error that says: _Face has no cross edge_. I have embedded my fault in the center of the domain (pics below).  
**Fault** has all the three sides marked as _fault\_edge_ to impose no-slip bc.

 ![image](https://canada1.discourse-cdn.com/flex036/uploads/geodynamics/original/2X/b/b48c34289117658daaa76fe03e6b4e86c6f7ebe2.png)  
**Domain**  
 ![image](https://canada1.discourse-cdn.com/flex036/uploads/geodynamics/original/2X/6/64bd34f75142d5c5e6a84d1d115e8ce5da588edb.png)

How to address this in V2.2?

Thank you  
Viven

---

<div class="post-metadata">

**Author:** ![baagaard](https://yyz2.discourse-cdn.com/flex036/user_avatar/community.geodynamics.org/baagaard/32/615_2.png) [@baagaard](https://community.geodynamics.org/u/baagaard)\
**Post date:** [December 23, 2022, 5:14pm UTC](https://community.geodynamics.org/t/face-with-no-cross-edge/2745/2 "2022-12-23T17:14:46Z")

</div>

I don’t see anything obvious that would be causing this error.

Have you run mesh quality checks? Is the maximum condition number metric (available in Cubit and ParaView) less than about 3.0?

---

<div class="post-metadata">

**Author:** ![VivenSharma](https://yyz2.discourse-cdn.com/flex036/user_avatar/community.geodynamics.org/vivensharma/32/933_2.png) [@VivenSharma](https://community.geodynamics.org/u/VivenSharma)\
**Post date:** [December 31, 2022, 11:50am UTC](https://community.geodynamics.org/t/face-with-no-cross-edge/2745/3 "2022-12-31T11:50:19Z")

</div>

Hi Brad

Yes, I had done the condition number checks beforehand.

 ![image](https://canada1.discourse-cdn.com/flex036/uploads/geodynamics/original/2X/2/20ca7f91f7b4bdfdc50a17b28bdb351400a6e8b3.png)  
Here is the error log

```auto
 -- Counting vertices for fault 'fault'.
 >> /home/wadai0/sharm760/PylithResearch/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/topology/MeshGenerator.py:125:_adjustTopology
 -- meshimporter(info)
 -- Adjusting topology for fault 'fault' with 6949 vertices.
Fatal error. Calling MPI_Abort() to abort PyLith application.
Traceback (most recent call last):
  File "/home/wadai0/sharm760/PylithResearch/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/apps/PetscApplication.py", line 74, in onComputeNodes
    self.main(*args, **kwds)
  File "/home/wadai0/sharm760/PylithResearch/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/apps/PyLithApp.py", line 115, in main
    mesh = self.mesher.create(self.problem.normalizer, interfaces)
  File "/home/wadai0/sharm760/PylithResearch/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/topology/MeshImporter.py", line 120, in create
    self._adjustTopology(mesh, faults)
  File "/home/wadai0/sharm760/PylithResearch/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/topology/MeshGenerator.py", line 127, in _adjustTopology
    interface.adjustTopology(mesh, firstFaultVertex, firstLagrangeVertex, firstFaultCell)
  File "/home/wadai0/sharm760/PylithResearch/pylith-2.2.2-linux-x86_64/lib/python2.7/site-packages/pylith/faults/faults.py", line 344, in adjustTopology
    def adjustTopology(self, *args): return _faults.FaultCohesive_adjustTopology(self, *args)
RuntimeError: Error occurred while adjusting topology to create cohesive cells for fault 'fault'.
Internal error while creating fault mesh. Face 10311061 has no cross edge.
application called MPI_Abort(MPI_COMM_WORLD, -1) - process 0
/home/wadai0/sharm760/PylithResearch/pylith-2.2.2-linux-x86_64/bin/nemesis: mpirun: exit 255
/home/wadai0/sharm760/PylithResearch/pylith-2.2.2-linux-x86_64/bin/pylith: /home/wadai0/sharm760/PylithResearch/pylith-2.2.2-linux-x86_64/bin/nemesis: exit 1

```

Thank you  
Viven

---

<div class="post-metadata">

**Author:** ![baagaard](https://yyz2.discourse-cdn.com/flex036/user_avatar/community.geodynamics.org/baagaard/32/615_2.png) [@baagaard](https://community.geodynamics.org/u/baagaard)\
**Post date:** [January 4, 2023, 2:46am UTC](https://community.geodynamics.org/t/face-with-no-cross-edge/2745/4 "2023-01-04T02:46:43Z")

</div>

Do you get this same error with a coarser mesh? The easiest way to debug these types of issues is to coarsen the mesh and then examine in detail a run with the coarsest mesh that generates the error.