[CIG-SHORT] pylith error : Assertion failed

[Brad Aagaard]

The error you are getting arises when PyLith has trouble inserting the 
cohesive cells. It looks like you are using PyLith v1.9.0. Is this correct?

How much space do you have between the faults? One cell? Two cells?

Brad


On 6/26/14, 4:40 PM, Marcelo Contreras wrote:
> 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 <http://MeshIOCubit.cc>:142:<unknown>
>   -- meshiocubit(info)
>   -- Reading 52095 vertices.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:197:<unknown>
>   -- meshiocubit(info)
>   -- Reading 272965 cells in 5 blocks.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:257:<unknown>
>   -- meshiocubit(info)
>   -- Found 12 node sets.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'ALLNODES' with id 1 containing 52095 nodes.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'BOTTOM_MANTLE_BORDER' with id 2 containing 558 nodes.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'BOTTOM_SLAB_BORDER' with id 3 containing 655 nodes.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'EAST_BORDER' with id 4 containing 621 nodes.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'FAULT_SEISMO' with id 5 containing 725 nodes.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'FAULT_SLABBOT' with id 6 containing 6696 nodes.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'FAULT_SLABTOP' with id 7 containing 4167 nodes.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'GROUNDSURF' with id 8 containing 11037 nodes.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'NORTH_BORDER' with id 9 containing 857 nodes.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'SOUTH_BORDER' with id 10 containing 859 nodes.
>   >> meshio/MeshIOCubit.cc <http://MeshIOCubit.cc>:285:<unknown>
>   -- meshiocubit(info)
>   -- Reading node set 'WEST_MANTLE_BORDER' with id 11 containing 110 nodes.
>   >> meshio/MeshIOCubit.cc <http://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 <http://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!!  :-(
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>