[CIG-SHORT] pylith error : Assertion failed

[Marcelo Contreras]

Hello,

In my simulation I have there faults for my slab:  fault_slabtop and fault_slabbot with a constant slip rate (in subduction) and I would like a third fault  (fault_seismo) with stick-slip behavior like example step14.cfg.  I have solution from pylith with each  isolated fault, but when I use two faults that are one beside the other (fault_seismo+fault_slabtop) Pylith show me the following error.

pylith inter_3D.cfg 
 >> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/utils/PetscManager.py:64:initialize
 -- petsc(info)
 -- Initialized PETSc.
 >> /Users/marcont/pylith190/lib/python2.7/site-packages/pylith/meshio/MeshIOObj.py:55:read
 -- meshiocubit(info)
 -- Reading finite-element mesh
 >> meshio/MeshIOCubit.cc:142:<unknown>
 -- meshiocubit(info)
 -- Reading 52095 vertices.
 >> meshio/MeshIOCubit.cc:197:<unknown>
 -- meshiocubit(info)
 -- Reading 272965 cells in 5 blocks.
 >> meshio/MeshIOCubit.cc:257:<unknown>
 -- meshiocubit(info)
 -- Found 12 node sets.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'ALLNODES' with id 1 containing 52095 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'BOTTOM_MANTLE_BORDER' with id 2 containing 558 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'BOTTOM_SLAB_BORDER' with id 3 containing 655 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'EAST_BORDER' with id 4 containing 621 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'FAULT_SEISMO' with id 5 containing 725 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'FAULT_SLABBOT' with id 6 containing 6696 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'FAULT_SLABTOP' with id 7 containing 4167 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'GROUNDSURF' with id 8 containing 11037 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'NORTH_BORDER' with id 9 containing 857 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'SOUTH_BORDER' with id 10 containing 859 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'WEST_MANTLE_BORDER' with id 11 containing 110 nodes.
 >> meshio/MeshIOCubit.cc:285:<unknown>
 -- meshiocubit(info)
 -- Reading node set 'WEST_SLAB_BORDER' with id 12 containing 265 nodes.
Assertion failed: (classifySize < vReplaceCells.size() + vNoReplaceCells.size()), function classifyCells, file faults/TopologyOps.cc, line 117.
[0]0:Return code = 0, signaled with Abort trap: 6


In the image attached you see there aren’t nodes shared between fault_seismo and fault_slabtop. Please can anyone help me solve the problem ?

best regards
Marcelo.

testing :
only fault_slabtop —> pylith problem working
only fault_slabbot —> pylith problem working
only fault_seismo —> pylith problem working (constant slip rate, coseismic)
only fault_slabbot+fault_slabtop —> pylith problem working (inter-seismic problem)
only fault_slabbot+fault_seismo —> pylith problem working (constant slip rate + static friction)
only fault_slabtop+fault_seismo —> NOT WORKING!!  :-(

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20140626/eb92009d/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Captura de pantalla 2014-06-26 a la(s) 19.01.37.png
Type: image/png
Size: 336517 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20140626/eb92009d/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Captura de pantalla 2014-06-26 a la(s) 19.02.26.png
Type: image/png
Size: 242523 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20140626/eb92009d/attachment-0003.png>