On Sun, Sep 16, 2012 at 1:36 PM, Birendra jha <span dir="ltr"><<a href="mailto:bjha7333@yahoo.com" target="_blank">bjha7333@yahoo.com</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">
<table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td valign="top" style="font:inherit">Hi Matt<br><br>Attached is the gdb backtrace of the crash when I use LU preconditioner.<br></td></tr></tbody></table></blockquote>
<div><br></div><div>This crashes in the Nested Dissection routine. I have only ever seen this in matrices missing</div><div>some rows. Are you sure that this problem has no faults turned on? If so, please send the</div><div>
cfg files and meshes so I can reproduce the problem here.</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">
<table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td valign="top" style="font:inherit">Thanks <br>Birendra<br><br>--- On <b>Fri, 9/14/12, Matthew Knepley <i><<a href="mailto:knepley@mcs.anl.gov" target="_blank">knepley@mcs.anl.gov</a>></i></b> wrote:<br>
<blockquote style="border-left:2px solid rgb(16,16,255);margin-left:5px;padding-left:5px"><br>From: Matthew Knepley <<a href="mailto:knepley@mcs.anl.gov" target="_blank">knepley@mcs.anl.gov</a>><br>Subject: Re: [CIG-SHORT] Change in dt means new Jacobian<br>
To: "Birendra jha" <<a href="mailto:bjha7333@yahoo.com" target="_blank">bjha7333@yahoo.com</a>><br>Cc: "Brad Aagaard" <<a href="mailto:baagaard@usgs.gov" target="_blank">baagaard@usgs.gov</a>>, <a href="mailto:cig-short@geodynamics.org" target="_blank">cig-short@geodynamics.org</a><br>
Date: Friday, September 14, 2012, 6:40 PM<div><div class="h5"><br><br><div>On Fri, Sep 14, 2012 at 1:53 AM, Birendra jha <span dir="ltr"><<a rel="nofollow" href="http://mc/compose?to=bjha7333@yahoo.com" target="_blank">bjha7333@yahoo.com</a>></span> wrote:<br>
<div><blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear developers<br>
<br>
I simplified my problem setup to see if no. of KSP iterations goes down. There are 46552 nodes, 42525 hex elements in a box. In this very simple setup, both ASM and ML preconditioner takes ~90 iterations for convergence.<br>
Attached are the cfg files and stdout from LU and ML preconditioners. ASM stdout is not attached.<br></blockquote><div><br></div><div>Something in your setup is causing the problem to be nearly or actually singular. Its either bad elements in</div>
<div>the mesh, or an inconsistent boundary condition I think.</div><div> </div><blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Why does the run with LU crash? I just used following 2 lines to use LU:<br>
[pylithapp.petsc]<br>
pc_type = lu<br>
Is it incorrect?<br></blockquote><div><br></div><div>This is interesting. Can you run in the debugger --petsc.start_in_debugger and get a stack trace?</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div>
<div> </div><blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
For both ASM and ML, the solution looks correct in Paraview. It's just the time taken in integrating residual that I would like to reduce.<br>
<br>
Another question: when Pylith is running on a single multicore machine, can it use OpenMPI instead of MPI (I am assuming it's using MPI by default)?<br>
<br>
Thanks and regards<br>
Birendra<br>
<br>
--- On Tue, 9/4/12, Brad Aagaard <<a rel="nofollow" href="http://mc/compose?to=baagaard@usgs.gov" target="_blank">baagaard@usgs.gov</a>> wrote:<br>
<br>
> From: Brad Aagaard <<a rel="nofollow" href="http://mc/compose?to=baagaard@usgs.gov" target="_blank">baagaard@usgs.gov</a>><br>
> Subject: Re: Change in dt means new Jacobian<br>
> To: "Birendra jha" <<a rel="nofollow" href="http://mc/compose?to=bjha7333@yahoo.com" target="_blank">bjha7333@yahoo.com</a>><br>
> Cc: <a rel="nofollow" href="http://mc/compose?to=cig-short@geodynamics.org" target="_blank">cig-short@geodynamics.org</a><br>
> Date: Tuesday, September 4, 2012, 7:59 PM<br>
> On 09/03/2012 08:54 PM, Birendra jha<br>
> wrote:<br>
> > 3. How do I utilize multiple cores on a single machine?<br>
> Sorry if this<br>
> > question is already answered somewhere. I found<br>
> instructions on<br>
> > running parallel using a batch system but don't know<br>
> how to run on 1<br>
> > node and use multiple cores. Does this require<br>
> re-building Pylith?<br>
><br>
> To run in parallel all you have to do is add the command<br>
> line argument --nodes=NUM_CORES. You do not need to<br>
> rebuild.<br>
><br>
> > 4. For HDF5 output, is it correct that the I have to<br>
> wait until<br>
> > finish to get the xmf file which is needed to load h5<br>
> files in<br>
> > Paraview? Manual says "This file is written when PyLith<br>
> closes the<br>
> > HDF5 file at the end of the simulation".<br>
><br>
> The Xdmf file is written when PyLith closes the HDF5 file at<br>
> the end of a simulation. There is a python script included<br>
> with PyLith called pylith_genxdmf that will generate the<br>
> Xdmf file given an HDF5 file. This is handy if a simulation<br>
> did not finish or you add additional fields to an existing<br>
> HDF5 file. Generating the Xdmf file in the middle of a<br>
> simulation can work but it is risky because the HDF5 file<br>
> can get corrupted if multiple processes try to access the<br>
> file simultaneously.<br>
><br>
> pylith_genxdmf --file=HDF5_FILE<br>
><br>
> Brad<br>
><br>_______________________________________________<br>
CIG-SHORT mailing list<br>
<a rel="nofollow" href="http://mc/compose?to=CIG-SHORT@geodynamics.org" target="_blank">CIG-SHORT@geodynamics.org</a><br>
<a rel="nofollow" href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br>
<br></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>
</div></div></div></blockquote></td></tr></tbody></table></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>