[CIG-SHORT] cohesive cells topology error

Fri Nov 6 07:37:25 PST 2015

I have attached the grid file. Below is the same error on pylith-2.1.0:

 -- Reading node set 'fault1' with id 25 containing 31 nodes.
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Invalid argument
[0]PETSC ERROR: Split face 2369 has 1 != 2 supports
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Development GIT revision: v3.5.3-2251-gb193602  GIT Date: 2015-02-19 19:20:25 -0600
[0]PETSC ERROR: /home/pylith-2.1.0/bin/mpinemesis on a arch-pylith named eofe5 by user Fri Nov  6 10:32:18 2015
[0]PETSC ERROR: Configure options --prefix=/home/pylith-2.1.0 --with-c2html=0 --with-x=0 --with-clanguage=C --with-mpicompilers=1 --with-debugging=1 --with-shared-libraries=1 --download-chaco=1 --download-ml=1 --download-fblaslapack=1 --with-hdf5=1 --with-hdf5-include=/home/pylith-2.1.0/include --with-hdf5-lib=/home/pylith-2.1.0/lib/libhdf5.dylib --LIBS=-lz CPPFLAGS="-I/home/pylith-2.1.0/include " LDFLAGS="-L/home/pylith-2.1.0/lib " CFLAGS="-g -O2" CXXFLAGS="-g -O2 -DMPICH_IGNORE_CXX_SEEK" FCFLAGS="-g -O2" PETSC_DIR=/home/build/pylith-2.1.0/petsc-pylith PETSC_ARCH=arch-pylith
[0]PETSC ERROR: #1 DMPlexLabelCohesiveComplete() line 1431 in /home/build/pylith-2.1.0/petsc-pylith/src/dm/impls/plex/plexsubmesh.c
[0]PETSC ERROR: #2 static void pylith::faults::CohesiveTopology::create(pylith::topology::Mesh*, const pylith::topology::Mesh&, PetscDMLabel, int, int&, int&, int&, bool)() line 163 in ../../../pylith-2.1.0/libsrc/pylith/faults/CohesiveTopology.cc
Fatal error. Calling MPI_Abort() to abort PyLith application.
Traceback (most recent call last):
  File "/home/pylith-2.1.0/lib/python2.7/site-packages/pylith/apps/PetscApplication.py", line 64, in onComputeNodes
    self.main(*args, **kwds)
  File "/home/pylith-2.1.0/lib/python2.7/site-packages/pylith/apps/PyLithApp.py", line 112, in main
    mesh = self.mesher.create(self.problem.normalizer, interfaces)
  File "/home/pylith-2.1.0/lib/python2.7/site-packages/pylith/topology/MeshImporter.py", line 120, in create
    self._adjustTopology(mesh, faults)
  File "/home/pylith-2.1.0/lib/python2.7/site-packages/pylith/topology/MeshGenerator.py", line 127, in _adjustTopology
    interface.adjustTopology(mesh, firstFaultVertex, firstLagrangeVertex, firstFaultCell)
  File "/home/pylith-2.1.0/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 'fault1'.
Error detected while in PETSc function.

--------------------------------------------
On Thu, 11/5/15, Birendra jha <bjha7333 at yahoo.com> wrote:

 Subject: Re: cohesive cells topology error
 To: cig-short at geodynamics.org
 Date: Thursday, November 5, 2015, 7:21 PM

 Here is complete error:

  -- Reading node set 'fault1' with id 25 containing 31
 nodes.
  >>
 ../../../pylith-1.8.0/libsrc/pylith/meshio/MeshIOCubit.cc:285:<unknown>
  -- meshiocubit(info)
  -- Reading node set 'fault2' with id 26 containing 32
 nodes.
 Fatal error. Calling MPI_Abort() to abort PyLith
 application.
 Traceback (most recent call last):
   File
 "/data/pylith/lib/python2.7/site-packages/pylith/apps/PetscApplication.py",
 line 65, in onComputeNodes
     self.main(*args, **kwds)
   File
 "/data/pylith/lib/python2.7/site-packages/pylith/apps/PyLithApp.py",
 line 103, in main
     mesh =
 self.mesher.create(self.problem.normalizer, interfaces)
   File
 "/data/pylith/lib/python2.7/site-packages/pylith/topology/MeshImporter.py",
 line 120, in create
     self._adjustTopology(mesh, faults)
   File
 "/data/pylith/lib/python2.7/site-packages/pylith/topology/MeshGenerator.py",
 line 127, in _adjustTopology
     interface.adjustTopology(mesh,
 firstFaultVertex, firstLagrangeVertex, firstFaultCell)
   File
 "/data/pylith/lib/python2.7/site-packages/pylith/faults/faults.py",
 line 370, in adjustTopology
     def adjustTopology(self, *args): return
 _faults.FaultCohesive_adjustTopology(self, *args)
 RuntimeError: Error occurred while adjusting topology to
 create cohesive cells for fault 'fault1'.
 Error classifying cells during creation of cohesive cells.
   classifySize: 8, classifyTotal: 7

 Is the problem in fault2?

 --------------------------------------------
 On Thu, 11/5/15, Birendra jha <bjha7333 at yahoo.com>
 wrote:

  Subject: cohesive cells topology error
  To: cig-short at geodynamics.org
  Date: Thursday, November 5, 2015, 7:17 PM

  Hello developers

  Can you please help with the following error that shows up
  immediately after running pylith?

    File
  "/data/pylith/lib/python2.7/site-packages/pylith/topology/MeshGenerator.py",
  line 127, in _adjustTopology
      interface.adjustTopology(mesh,
  firstFaultVertex, firstLagrangeVertex, firstFaultCell)
    File
  "/data/pylith/lib/python2.7/site-packages/pylith/faults/faults.py",
  line 370, in adjustTopology
      def adjustTopology(self, *args): return
  _faults.FaultCohesive_adjustTopology(self, *args)
  RuntimeError: Error occurred while adjusting topology to
  create cohesive cells for fault 'fault1'.
  Error classifying cells during creation of cohesive cells.
    classifySize: 8, classifyTotal: 7

  It's a 2D triangle grid with two faults in the grid and I
  only activate one of the two faults in the cfg file.

  The grid quality is good (max condition number <2).

  Many thanks 
  Birendra
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mesh_tri_fault1.exo
Type: application/octet-stream
Size: 63704 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20151106/4cfd0be7/attachment-0001.obj>