On Thu, Oct 4, 2012 at 10:51 AM, Hongfeng Yang <span dir="ltr">&lt;<a href="mailto:hyang@whoi.edu" target="_blank">hyang@whoi.edu</a>&gt;</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">

  
    
  
  <div text="#000000" bgcolor="#FFFFFF"><div class="im">
    <div>On 10/04/2012 10:41 AM, Matthew Knepley
      wrote:<br>
    </div>
    </div><div class="im"><blockquote type="cite">On Thu, Oct 4, 2012 at 10:39 AM, Hongfeng Yang <span dir="ltr">&lt;<a href="mailto:hyang@whoi.edu" target="_blank">hyang@whoi.edu</a>&gt;</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&lt;pylith::topology::Mesh&gt;&amp;,
          PylithScalar,<br>
          pylith::topology::SolutionFields*): Assertion
          `fabs(tractionNormal) &lt;<br>
          _zeroTolerance&#39; 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>
    </blockquote></div>
    Hi Matt,<br>
    <br>
    Do you mean the v_fault not only represents the ID of vertices on
    the fault, but any vertex in the volume/full mesh?<br></div></blockquote><div><br></div><div>Yes. That is the way it is currently setup.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div>
 </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div text="#000000" bgcolor="#FFFFFF">
    If that is the case, I know where I go wrong.<span class="HOEnZb"><font color="#888888"><br>
    <br>
    Hongfeng<br>
  </font></span></div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>