We can check this sort of thing by looking at the VTK file that is output initially. I<br>
do not have anything to check split nodes yet, but we can dream it up.<br>
<br>
&nbsp;&nbsp; Matt<br><br><div><span class="gmail_quote">On 8/31/06, <b class="gmail_sendername">Brad Aagaard</b> &lt;<a href="mailto:baagaard@usgs.gov">baagaard@usgs.gov</a>&gt; wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Eric-<br><br>I am guessing the order of the nodes for the elements needs to be<br>reveresed. Swapping nodes 2 &amp; 3 would do it. I tried this on my<br>machine and it is running.<br><br>Brad<br><br><br>On Thursday 31 August 2006 5:16 pm, Eric Andreas Hetland wrote:
<br>&gt; I found several errors with the original input files, and I cleaned<br>&gt; them up, the new version is on the CIG website. But only for the<br>&gt; first mesh. I will do the second and third ones once I know this
<br>&gt; one is rights. I still get an error with this input deck (error is<br>&gt; below). As much as I can see, the mesh is being constructed<br>&gt; correctly, and I think that I am constructing the split nodes right<br>
&gt; (I checked them over and over, but I still might have missed<br>&gt; something). The BC's may be the problem, but I will have to look at<br>&gt; them tomorrow morning. I built pylith input files using Carl's<br>&gt; uniform mesh, and get the same crash error on it... Probably the
<br>&gt; BC's, could still be the split nodes... Whatever it is, it is going<br>&gt; to have to wait till morning since I have a prior commitment this<br>&gt; evening.<br>&gt;<br>&gt; - Eric.<br>&gt;<br>&gt; Hello from pl3dscan.preinitialize
 (begin)!<br>&gt; Reforming the stiffness matrix:<br>&gt;<br>&gt; Traceback (most recent call last):<br>&gt;&nbsp;&nbsp; File &quot;/opt/geomod/pylith-0.8-linux-x86/bin/pylith3dapp.py&quot;, line<br>&gt; 49, in ? app.run()<br>&gt;&nbsp;&nbsp; File
<br>&gt; &quot;/home/leif/monolith/pylith-0.8-linux-x86/lib/python2.3/site-packag<br>&gt;es/pythia -0.8/pyre/applications/Application.py&quot;, line 77, in run<br>&gt;&nbsp;&nbsp; File<br>&gt; &quot;/home/leif/monolith/pylith-0.8-linux-x86
/lib/python2.3/site-packag<br>&gt;es/pythia -0.8/pyre/applications/Stager.py&quot;, line 19, in execute<br>&gt;&nbsp;&nbsp; File<br>&gt; &quot;/home/leif/monolith/pylith-0.8-linux-x86/lib/python2.3/site-packag<br>&gt;es/pylith 3d/Application.py&quot;, line 65, in main
<br>&gt;&nbsp;&nbsp; File<br>&gt; &quot;/home/leif/monolith/pylith-0.8-linux-x86/lib/python2.3/site-packag<br>&gt;es/pylith 3d/Pylith3d_run.py&quot;, line 384, in run<br>&gt; ValueError: getjac:&nbsp;&nbsp;element #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;-6.25000000E+07:
 Zero or<br>&gt; negative jacobian for element.<br>&gt;<br>&gt;<br>&gt; _______________________________________________<br>&gt; CIG-SHORT mailing list<br>&gt; <a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org
</a><br>&gt; <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br>_______________________________________________<br>CIG-SHORT mailing list<br>
<a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br><a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br></blockquote>
</div><br><br clear="all"><br>-- <br>&quot;Failure has a thousand explanations. Success doesn't need one&quot; -- Sir Alec Guiness