<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">On 10/04/2012 10:41 AM, Matthew Knepley
wrote:<br>
</div>
<blockquote
cite="mid:CAMYG4GkE2zUS2WdTN=pHciBE=Fe6gM=QuhDwovTDVbqVtxh9gg@mail.gmail.com"
type="cite">On Thu, Oct 4, 2012 at 10:39 AM, Hongfeng Yang <span
dir="ltr"><<a moz-do-not-send="true"
href="mailto:hyang@whoi.edu" target="_blank">hyang@whoi.edu</a>></span>
wrote:<br>
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Developers,<br>
<br>
I am tracking time history of variables on particular vertices
on the<br>
fault. However, I found the vertex ID in Pylith (v_fault) is
larger than<br>
the total number of vertices on the fault. How can I convert
the value<br>
of v_fault back to the vertex location on the fault surface?<br>
<br>
For example, I got the following fault opening error message<br>
WARNING! Fault opening with nonzero traction., v_fault: 4154,
opening:<br>
1.03716e-10, normal traction: -0.0042952<br>
mpinemesis: faults/FaultCohesiveDyn.cc:382: virtual void<br>
pylith::faults::FaultCohesiveDyn::integrateResidual(const<br>
pylith::topology::Field<pylith::topology::Mesh>&,
PylithScalar,<br>
pylith::topology::SolutionFields*): Assertion
`fabs(tractionNormal) <<br>
_zeroTolerance' failed.<br>
[0]0:Return code = 0, signaled with Aborted<br>
<br>
v_fault points to vertex ID 4154. The total number of vertices
on the<br>
fault is 1472, however.<br>
</blockquote>
<div><br>
</div>
<div>The fault vertices are just vertices in the mesh, so this
refers to the full</div>
<div>mesh vertex numbering.</div>
<div><br>
</div>
</div>
</blockquote>
Does Pylith output the vertex and topology in a re-numbering way? I
read the output value of each cell on the fault and could not find
the value of v_fault in any cell.<br>
<br>
Thanks,<br>
<br>
Hongfeng <br>
<br>
</body>
</html>