[CIG-SHORT] How does Pylith number vertices on fault

[Matthew Knepley]

On Thu, Oct 4, 2012 at 10:39 AM, Hongfeng Yang <hyang at whoi.edu> wrote:

> Hi Developers,
>
> I am tracking time history of variables on particular vertices on the
> fault. However, I found the vertex ID in Pylith (v_fault) is larger than
> the total number of vertices on the fault. How can I convert the value
> of v_fault back to the vertex location on the fault surface?
>
> For example, I got the following fault opening error message
> WARNING! Fault opening with nonzero traction., v_fault: 4154, opening:
> 1.03716e-10, normal traction: -0.0042952
> mpinemesis: faults/FaultCohesiveDyn.cc:382: virtual void
> pylith::faults::FaultCohesiveDyn::integrateResidual(const
> pylith::topology::Field<pylith::topology::Mesh>&, PylithScalar,
> pylith::topology::SolutionFields*): Assertion `fabs(tractionNormal) <
> _zeroTolerance' failed.
> [0]0:Return code = 0, signaled with Aborted
>
> v_fault points to vertex ID 4154. The total number of vertices on the
> fault is 1472, however.
>

The fault vertices are just vertices in the mesh, so this refers to the full
mesh vertex numbering.

  Thanks,

     Matt


> Thanks,
>
> Hongfeng
>
> --
> Postdoc Investigator
> Woods Hole Oceanographic Institution
> Dept. Geology and Geophysics
> 360 Woods Hole Rd, MS 24
> Woods Hole, MA 02543
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-short/attachments/20121004/b9747aa8/attachment.htm