<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">As further input, when I ran Birendra's problem on my Linux machine, I didn't get any output either, except for output_info.vtk and output_t0000000.vtk, even though I had created the output directory. Also, on my Mac the linear solution converged in 134 iterations, while on the Linux machine there was an initial residual of 0 (no iterations). I'm attaching the pylithinfo output using test.cfg (run on my Linux machine). I'm not sure what the problem might be at present.<div><br></div><div>Charles</div><div><br></div><div></div><span style="white-space:pre-wrap">[see attached file: pylith_parameters.txt.gz]</span><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div></div><div><br><div><div>On 14/03/2012, at 11:47 AM, Brad Aagaard wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Birenda,<br><br>I was able to run your test.cfg file after creating the output directory <br>(as Charles noted), but initializing the integrators is taking a rather <br>long time. Your spatial databases that describe the variation in <br>physical properties are quite large. For every quadrature point in every <br>cell PyLith has to find which points to use for the interpolation. Using <br>a cloud of points with a much finer spacing than the size of the cell <br>will cause PyLith to run much slower than necessary.<br><br>The solution did converge, but when I looked at your mesh I see some <br>cells that have very poor aspect ratios so the solution may not be very <br>accurate. See the PyLith Wiki <br>(<a href="http://www.geodynamics.org/cig/community/workinggroups/short/workarea/pylith-wiki/index_html">http://www.geodynamics.org/cig/community/workinggroups/short/workarea/pylith-wiki/index_html</a>) <br>for hints on checking the aspect ratio. Using tetehedral cells will give <br>you much greater flexibility in varying the discretization size while <br>keeping the aspect ratio of cells close to 1. See examples/2d/subduction <br>for an example of using the bias scheme.<br><br>Regards,<br>Brad<br><br><br><br>On 03/13/2012 03:08 PM, Birendra jha wrote:<br><blockquote type="cite">Dear developers,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I am trying to run a simple case like Step16 in the tutorial except<br></blockquote><blockquote type="cite">that density now varies with depth. I set the initial geostatic<br></blockquote><blockquote type="cite">stress from this density profile. I expect to see a zero<br></blockquote><blockquote type="cite">displacement. I have attached the files as zipped.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Why I am not getting any output (It must be some small mistake on my<br></blockquote><blockquote type="cite">part. Just can't figure it out)?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Thanks and regards Birendra<br></blockquote><blockquote type="cite"><br></blockquote><br>_______________________________________________<br>CIG-SHORT mailing list<br><a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br>http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short<br></div></blockquote></div><br><div>
<span class="Apple-style-span" style="font-size: 12px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="3" style="font: normal normal normal 12px/normal Helvetica; ">Charles A. Williams</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Scientist<br>GNS Science<br>1 Fairway Drive, Avalon<br>PO Box 30368<br>Lower Hutt 5040<br>New Zealand<br>ph (office): 0064-4570-4566<br>fax (office): 0064-4570-4600</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><a href="mailto:C.Williams@gns.cri.nz">C.Williams@gns.cri.nz</a></div></div></div></span>
</div>
<br></div></body></html>