[CIG-SHORT] Cubit tet mesh not being read into pylith

[Brad Aagaard]

Oliver-

This error is related to creating cohesive cells. First, try reading the 
mesh without any faults (this error should disappear). Then try with one 
fault. If you have multiple faults, try to isolate which fault is 
causing the error. The restrictions are: (1) no vertex can be on more 
than one fault and (2) more than one of a cell's faces cannot be on a fault.

Brad


Oliver Boyd wrote:
> Having trouble reading a cubit generated tet mesh into pylith. I get the
> following output and error.
> 
> $pylith dislocation.cfg
>  >> 
> /Users/buildbot/install/x86_darwin_single_nosched_binbot/lib/python2.3/site-
> packages/pylith/utils/PetscManager.py:48:initialize
>  -- petsc(info)
>  -- Initializing PETSc.
>  >> 
> /Users/buildbot/install/x86_darwin_single_nosched_binbot/lib/python2.3/site-
> packages/pylith/meshio/MeshIOObj.py:46:read
>  -- meshiocubit(info)
>  -- Reading finite-element mesh
>  >> meshio/MeshIOCubit.cc:151:<unknown>
>  -- meshiocubit(info)
>  -- Reading 37250 vertices.
>  >> meshio/MeshIOCubit.cc:207:<unknown>
>  -- meshiocubit(info)
>  -- Reading 199518 cells in 2 blocks.
>  >> meshio/MeshIOCubit.cc:283:<unknown>
>  -- meshiocubit(info)
>  -- Found 7 node sets.
>  >> meshio/MeshIOCubit.cc:303:<unknown>
>  -- meshiocubit(info)
>  -- Reading node set 10 with 788 nodes.
>  >> meshio/MeshIOCubit.cc:303:<unknown>
>  -- meshiocubit(info)
>  -- Reading node set 11 with 85 nodes.
>  >> meshio/MeshIOCubit.cc:303:<unknown>
>  -- meshiocubit(info)
>  -- Reading node set 12 with 85 nodes.
>  >> meshio/MeshIOCubit.cc:303:<unknown>
>  -- meshiocubit(info)
>  -- Reading node set 13 with 85 nodes.
>  >> meshio/MeshIOCubit.cc:303:<unknown>
>  -- meshiocubit(info)
>  -- Reading node set 14 with 85 nodes.
>  >> meshio/MeshIOCubit.cc:303:<unknown>
>  -- meshiocubit(info)
>  -- Reading node set 15 with 590 nodes.
>  >> meshio/MeshIOCubit.cc:303:<unknown>
>  -- meshiocubit(info)
>  -- Reading node set 16 with 590 nodes.
> faults/TopologyOps.cc:111: failed assertion `classifySize <
> vReplaceCells.size() + vNoReplaceCells.size()'
> [0]0:Return code = 0, signaled with Abort trap
> 
> Thanks for your help,
> Oliver
> 
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short