On Thu, Aug 30, 2012 at 4:11 PM, Romain Jolivet <span dir="ltr">&lt;<a href="mailto:rjolivet@caltech.edu" target="_blank">rjolivet@caltech.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 style="word-wrap:break-word">I did just run:<div><br></div><div>pylith step01.cfg --petsc.fp_trap=0</div><div><br></div><div>Is that what you meant? If it is, I still have the same problem :-) and it does not change what I see on the screen. Should it change something?\</div>
</div></blockquote><div><br></div><div>Then it is not underflow.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>
However, the norm computation problem does not appear when I set the zero_tolerance to 1.0 (it does not converge though, but I don&#39;t have this Floating point exception...)</div></div></blockquote><div><br></div><div>As Brad says, change the relative tolerance, but leave the absolute tolerance intact.</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 style="word-wrap:break-word"><div>Romain</div><div><br></div><div><br>
<div>
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px">Romain Jolivet<br>
Postdoctoral Scholar<br>Geological and Planetary Sciences<br>California Institute of Technology, Pasadena, CA <br><a href="mailto:rjolivet@caltech.edu" target="_blank">rjolivet@caltech.edu</a><br><a href="tel:%2B1%20%28626%29%20560%206356" value="+16265606356" target="_blank">+1 (626) 560 6356</a></span>
</div>
<br><div><div>On Aug 30, 2012, at 2:00 PM, Matthew Knepley wrote:</div><br><blockquote type="cite">On Thu, Aug 30, 2012 at 3:38 PM, Romain Jolivet <span dir="ltr">&lt;<a href="mailto:rjolivet@caltech.edu" target="_blank">rjolivet@caltech.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">

Thanks for the advice. I changed the preconditioner settings to what you suggested. I don&#39;t see that jump in the linear norm anymore.<br>
However, I still have the floating point exception problem. This arise after the first round of ksp.<br>
<br>
In the log file (attached), I see to messages saying the ksp solver for friction does converge (twice), then the SNES residual norm is printed, then KSP starts, and when it has to compute the new SNES residual norm, it says that VecNorm() encoutered a Inf or NAN value:<br>


</blockquote><div><br></div><div>It could still be underflow. Try running with <span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">--petsc.fp_trap=0</span></div>
<div><span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif"><br></span></div><div><span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">  Thanks,</span></div>
<div><span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif"><br></span></div><div><span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">      Matt</span></div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  141 KSP Residual norm 1.185626397065e-09<br>
  142 KSP Residual norm 8.558244459183e-10<br>
  143 KSP Residual norm 4.537120666370e-10<br>
  144 KSP Residual norm 2.695590025796e-10<br>
  145 KSP Residual norm 2.263153111449e-10<br>
  146 KSP Residual norm 2.156010427246e-10<br>
  Linear solve converged due to CONVERGED_RTOL iterations 146<br>
<div>[0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[0]PETSC ERROR: Floating point exception!<br>
[0]PETSC ERROR: Infinite or not-a-number generated in norm!<br>
</div>[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Development HG revision: 3bede28f1195e406d4257c9fbeebefe827793801  HG Date: Mon Jun 11 22:14:12 2012 -0400<br>
[0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
[0]PETSC ERROR: See docs/index.html for manual pages.<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: /Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/bin/mpinemesis on a pylith-da named <a href="http://dhcp-74-229.caltech.edu/" target="_blank">dhcp-74-229.caltech.edu</a> by jolivetinsar Thu Aug 30 13:34:55 2012<br>



[0]PETSC ERROR: Libraries linked from /Users/buildbot/install/pylith_darwin_10.6_binbot/lib<br>
[0]PETSC ERROR: Configure run at Tue Jun 12 06:09:00 2012<br>
[0]PETSC ERROR: Configure options --prefix=/Users/buildbot/install/pylith_darwin_10.6_binbot --download-f-blas-lapack=1 --with-sieve=1 --with-c2html=0 --CXXFLAGS=-DMPICH_IGNORE_CXX_SEEK --with-clanguage=C++ --with-mpicompilers=1 --with-debugging=0 --download-chaco=1 --download-boost=1 --download-ml=1 --with-hdf5=1 --with-hdf5-dir=/Users/buildbot/install/pylith_darwin_10.6_binbot --with-x=0 --with-shared-libraries=1<br>



[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: VecNorm() line 169 in src/vec/vec/interface/rvector.c<br>
[0]PETSC ERROR: VecNormalize() line 263 in src/vec/vec/interface/rvector.c<br>
[0]PETSC ERROR: KSPGMRESCycle() line 127 in src/ksp/ksp/impls/gmres/gmres.c<br>
[0]PETSC ERROR: KSPSolve_GMRES() line 231 in src/ksp/ksp/impls/gmres/gmres.c<br>
[0]PETSC ERROR: KSPSolve() line 446 in src/ksp/ksp/interface/itfunc.c<br>
[0]PETSC ERROR: _sensitivitySolve() line 2172 in &quot;unknowndirectory/&quot;faults/FaultCohesiveDyn.cc<br>
Traceback (most recent call last):<br>
  File &quot;/Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/lib/python2.7/site-packages/pylith/apps/PetscApplication.py&quot;, line 65, in onComputeNodes<br>
Fatal error. Calling MPI_Abort() to abort PyLith application.<br>
    self.main(*args, **kwds)<br>
  File &quot;/Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/lib/python2.7/site-packages/pylith/apps/PyLithApp.py&quot;, line 119, in main<br>
    self.problem.run(self)<br>
  File &quot;/Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/lib/python2.7/site-packages/pylith/problems/TimeDependent.py&quot;, line 154, in run<br>
    self.formulation.step(t, dt)<br>
  File &quot;/Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/lib/python2.7/site-packages/pylith/problems/Implicit.py&quot;, line 219, in step<br>
    self.solver.solve(dispIncr, self.jacobian, residual)<br>
  File &quot;/Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/lib/python2.7/site-packages/pylith/problems/problems.py&quot;, line 180, in solve<br>
    def solve(self, *args): return _problems.SolverNonlinear_solve(self, *args)<br>
RuntimeError: std::exception<br>
application called MPI_Abort(MPI_COMM_WORLD, -1) - process 0<br>
/Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/bin/nemesis: mpirun: exit 255<br>
/Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/bin/pylith: /Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/bin/nemesis: exit 1<br>
<br>
Any idea?<br>
<div>Romain<br>
<br>
<br>
<br>
Romain Jolivet<br>
Postdoctoral Scholar<br>
Geological and Planetary Sciences<br>
California Institute of Technology, Pasadena, CA<br>
<a href="mailto:rjolivet@caltech.edu" target="_blank">rjolivet@caltech.edu</a><br>
<a href="tel:%2B1%20%28626%29%20560%206356" value="+16265606356" target="_blank">+1 (626) 560 6356</a><br>
<br>
</div>On Aug 30, 2012, at 12:00 PM, <a href="mailto:cig-short-request@geodynamics.org" target="_blank">cig-short-request@geodynamics.org</a> wrote:<br>
<br>
&gt; Send CIG-SHORT mailing list submissions to<br>
&gt;       <a href="mailto:cig-short@geodynamics.org" target="_blank">cig-short@geodynamics.org</a><br>
&gt;<br>
&gt; To subscribe or unsubscribe via the World Wide Web, visit<br>
&gt;       <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br>
&gt; or, via email, send a message with subject or body &#39;help&#39; to<br>
&gt;       <a href="mailto:cig-short-request@geodynamics.org" target="_blank">cig-short-request@geodynamics.org</a><br>
&gt;<br>
&gt; You can reach the person managing the list at<br>
&gt;       <a href="mailto:cig-short-owner@geodynamics.org" target="_blank">cig-short-owner@geodynamics.org</a><br>
&gt;<br>
&gt; When replying, please edit your Subject line so it is more specific<br>
&gt; than &quot;Re: Contents of CIG-SHORT digest...&quot;<br>
&gt;<br>
&gt;<br>
&gt; Today&#39;s Topics:<br>
&gt;<br>
&gt;   1. Re: Rate and State and Initial Traction (Brad Aagaard)<br>
&gt;   2. Fwd:  Rate and State and Initial Traction (Matthew Knepley)<br>
&gt;<br>
&gt;<br>
&gt; ----------------------------------------------------------------------<br>
&gt;<br>
&gt; Message: 1<br>
&gt; Date: Wed, 29 Aug 2012 12:36:08 -0700<br>
&gt; From: Brad Aagaard &lt;<a href="mailto:baagaard@usgs.gov" target="_blank">baagaard@usgs.gov</a>&gt;<br>
&gt; Subject: Re: [CIG-SHORT] Rate and State and Initial Traction<br>
&gt; To: <a href="mailto:cig-short@geodynamics.org" target="_blank">cig-short@geodynamics.org</a><br>
&gt; Message-ID: &lt;<a href="mailto:503E6F28.6050001@usgs.gov" target="_blank">503E6F28.6050001@usgs.gov</a>&gt;<br>
&gt; Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<div>&gt;<br>
&gt; Romain,<br>
&gt;<br>
&gt; It looks like you are using Schur complement preconditioner settings.<br>
&gt; The jump in the KSP residual at iteration 10 is odd and suggests<br>
&gt; something is not quite right in the linear solve (probably related to<br>
</div><div><div>&gt; the preconditioner). Try using the AMG preconditioner. It may not be<br>
&gt; quite as fast but it is more robust. Matt can probably help you tune the<br>
&gt; Schur complement solver conditions some more if necessary.<br>
&gt;<br>
&gt; Your SNES tolerances are quite large. I suggest using tolerances similar<br>
&gt; to those in examples/3d/hex8/step14.cfg.<br>
&gt;<br>
&gt; SOLVER TOLERANCES<br>
&gt;<br>
&gt; ksp_rtol = 1.0e-16<br>
&gt; ksp_atol = 1.0e-12<br>
&gt;<br>
&gt; snes_rtol = 1.0e-14<br>
&gt; snes_atol = 1.0e-10<br>
&gt;<br>
&gt;<br>
&gt; AMG PRECONDITIONER SETTINGS<br>
&gt;<br>
&gt; [pylithapp.timedependent.formulation]<br>
&gt; split_fields = True<br>
&gt; matrix_type = aij<br>
&gt;<br>
&gt; [pylithapp.petsc]<br>
&gt; ksp_gmres_restart = 100<br>
&gt; fs_pc_type = fieldsplit<br>
&gt; fs_pc_fieldsplit_real_diagonal = true<br>
&gt; fs_pc_fieldsplit_type = multiplicative<br>
&gt; fs_fieldsplit_0_pc_type = ml<br>
&gt; fs_fieldsplit_1_pc_type = jacobi<br>
&gt; fs_fieldsplit_0_ksp_type = preonly<br>
&gt; fs_fieldsplit_1_ksp_type = preonly<br>
&gt;<br>
&gt; Brad<br>
&gt;<br>
&gt;<br>
&gt; On 08/29/2012 11:42 AM, Romain Jolivet wrote:<br>
&gt;&gt; Hello Brad, Matt, Charles and the pylith team,<br>
&gt;&gt;<br>
&gt;&gt; Following your advice from the other day, I am experimenting  how to create earthquakes on a frictional interface by applying a traction change history. It works fine with the SlipWeakening law in 2D. I&#39;m trying to move on to the RateStateAgeing behavior and I have problems applying the initial traction conditions:<br>



&gt;&gt;<br>
&gt;&gt; I use a SimpleDB to apply the frictional properties on the fault: one domain is rate-strengthening, the other is rate-weakening with a linear transition between the two domains. If I do not try to turn on any traction_perturbation option, the code runs great, and I have a slow slip all the along the fault at any time.<br>



&gt;&gt;<br>
&gt;&gt; As I want the rate-weakening domain to be stuck sometimes, I want to turn on traction_perturbation.db_initial. I use a SimpleDB to put a negative traction all along the fault, and a shear traction compatible with the friction coefficients I impose. However, when I run this, I have the following error message from PETSc (I attach the log):<br>



&gt;&gt;<br>
&gt;&gt;&gt;&gt; /Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/lib/python2.7/site-packages/pylith/problems/Implicit.py:215:step<br>
&gt;&gt;  -- implicit(info)<br>
&gt;&gt;  -- Solving equations.<br>
&gt;&gt;   0 SNES Function norm 2.679589636063e+01<br>
&gt;&gt;     0 KSP Residual norm 9.637565795900e-01<br>
&gt;&gt;     1 KSP Residual norm 3.564937135471e-08<br>
&gt;&gt;     2 KSP Residual norm 4.313963401911e-16<br>
&gt;&gt;     3 KSP Residual norm 3.050427683646e-16<br>
&gt;&gt;     4 KSP Residual norm 2.490663081181e-16<br>
&gt;&gt;     5 KSP Residual norm 2.156976950523e-16<br>
&gt;&gt;     6 KSP Residual norm 1.929258577738e-16<br>
&gt;&gt;     7 KSP Residual norm 1.761163907914e-16<br>
&gt;&gt;     8 KSP Residual norm 1.630520838230e-16<br>
&gt;&gt;     9 KSP Residual norm 1.525212509706e-16<br>
&gt;&gt;    10 KSP Residual norm 1.833325155789e-11<br>
&gt;&gt;    11 KSP Residual norm 7.039784754868e-23<br>
&gt;&gt;   Linear solve converged due to CONVERGED_ATOL iterations 11<br>
&gt;&gt; [0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
&gt;&gt; [0]PETSC ERROR: Floating point exception!<br>
&gt;&gt; [0]PETSC ERROR: Infinite or not-a-number generated in norm!<br>
&gt;&gt;<br>
&gt;&gt; If anyone has an idea about that error message....<br>
&gt;&gt;<br>
&gt;&gt; I think it has to do with some slip vector. If I set friction.min_slip_rate=1e-1, this problem does not happen (the nonlinear solver does not converge, but at least, no crash). If I set the normal traction to be positive, it does not crash as well (and it does not converge as well). I join to this message the cfg file I am using for that experiment.<br>



&gt;&gt;<br>
&gt;&gt; Thank you<br>
&gt;&gt; Romain<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; Romain Jolivet<br>
&gt;&gt; Postdoctoral Scholar<br>
&gt;&gt; Geological and Planetary Sciences<br>
&gt;&gt; California Institute of Technology, Pasadena, CA<br>
&gt;&gt; <a href="mailto:rjolivet@caltech.edu" target="_blank">rjolivet@caltech.edu</a><br>
&gt;&gt; <a href="tel:%2B1%20%28626%29%20560%206356" value="+16265606356" target="_blank">+1 (626) 560 6356</a><br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; _______________________________________________<br>
&gt;&gt; CIG-SHORT mailing list<br>
&gt;&gt; <a href="mailto:CIG-SHORT@geodynamics.org" target="_blank">CIG-SHORT@geodynamics.org</a><br>
&gt;&gt; <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br>
&gt;&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
</div></div>&gt; ------------------------------<br>
&gt;<br>
&gt; Message: 2<br>
&gt; Date: Wed, 29 Aug 2012 22:09:25 -0500<br>
&gt; From: Matthew Knepley &lt;<a href="mailto:knepley@mcs.anl.gov" target="_blank">knepley@mcs.anl.gov</a>&gt;<br>
&gt; Subject: [CIG-SHORT] Fwd:  Rate and State and Initial Traction<br>
&gt; To: <a href="mailto:cig-short@geodynamics.org" target="_blank">cig-short@geodynamics.org</a><br>
&gt; Message-ID:<br>
&gt;       &lt;CAMYG4GkMcUqRjhMaYZ8Ooj-x=<a href="mailto:06QKFrc7fiwW1dStmDSDO1qWA@mail.gmail.com" target="_blank">06QKFrc7fiwW1dStmDSDO1qWA@mail.gmail.com</a>&gt;<br>
&gt; Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<div><div>&gt;<br>
&gt; On Wed, Aug 29, 2012 at 2:36 PM, Brad Aagaard &lt;<a href="mailto:baagaard@usgs.gov" target="_blank">baagaard@usgs.gov</a>&gt; wrote:<br>
&gt;<br>
&gt;&gt; Romain,<br>
&gt;&gt;<br>
&gt;&gt; It looks like you are using Schur complement preconditioner settings.<br>
&gt;&gt; The jump in the KSP residual at iteration 10 is odd and suggests<br>
&gt;&gt; something is not quite right in the linear solve (probably related to<br>
&gt;&gt;<br>
&gt;<br>
&gt; Brad is right that the tolerances here will produce some garbage. Going<br>
&gt; below<br>
&gt; machine epsilon (about 10^-16) means you get randomness and ultimately<br>
&gt; underflow which I think you got in a norm calculation since things get<br>
&gt; squared.<br>
&gt; The solver should have converged on iterate 2, and will with the tolerances<br>
&gt; that<br>
&gt; Brad suggested.<br>
&gt;<br>
&gt;  Thanks,<br>
&gt;<br>
&gt;     Matt<br>
&gt;<br>
&gt;<br>
&gt;&gt; the preconditioner). Try using the AMG preconditioner. It may not be<br>
&gt;&gt; quite as fast but it is more robust. Matt can probably help you tune the<br>
&gt;&gt; Schur complement solver conditions some more if necessary.<br>
&gt;&gt;<br>
&gt;&gt; Your SNES tolerances are quite large. I suggest using tolerances similar<br>
&gt;&gt; to those in examples/3d/hex8/step14.cfg.<br>
&gt;&gt;<br>
&gt;&gt; SOLVER TOLERANCES<br>
&gt;&gt;<br>
&gt;&gt; ksp_rtol = 1.0e-16<br>
&gt;&gt; ksp_atol = 1.0e-12<br>
&gt;&gt;<br>
&gt;&gt; snes_rtol = 1.0e-14<br>
&gt;&gt; snes_atol = 1.0e-10<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; AMG PRECONDITIONER SETTINGS<br>
&gt;&gt;<br>
&gt;&gt; [pylithapp.timedependent.formulation]<br>
&gt;&gt; split_fields = True<br>
&gt;&gt; matrix_type = aij<br>
&gt;&gt;<br>
&gt;&gt; [pylithapp.petsc]<br>
&gt;&gt; ksp_gmres_restart = 100<br>
&gt;&gt; fs_pc_type = fieldsplit<br>
&gt;&gt; fs_pc_fieldsplit_real_diagonal = true<br>
&gt;&gt; fs_pc_fieldsplit_type = multiplicative<br>
&gt;&gt; fs_fieldsplit_0_pc_type = ml<br>
&gt;&gt; fs_fieldsplit_1_pc_type = jacobi<br>
&gt;&gt; fs_fieldsplit_0_ksp_type = preonly<br>
&gt;&gt; fs_fieldsplit_1_ksp_type = preonly<br>
&gt;&gt;<br>
&gt;&gt; Brad<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; On 08/29/2012 11:42 AM, Romain Jolivet wrote:<br>
&gt;&gt;&gt; Hello Brad, Matt, Charles and the pylith team,<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Following your advice from the other day, I am experimenting  how to<br>
&gt;&gt; create earthquakes on a frictional interface by applying a traction change<br>
&gt;&gt; history. It works fine with the SlipWeakening law in 2D. I&#39;m trying to move<br>
&gt;&gt; on to the RateStateAgeing behavior and I have problems applying the initial<br>
&gt;&gt; traction conditions:<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; I use a SimpleDB to apply the frictional properties on the fault: one<br>
&gt;&gt; domain is rate-strengthening, the other is rate-weakening with a linear<br>
&gt;&gt; transition between the two domains. If I do not try to turn on any<br>
&gt;&gt; traction_perturbation option, the code runs great, and I have a slow slip<br>
&gt;&gt; all the along the fault at any time.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; As I want the rate-weakening domain to be stuck sometimes, I want to<br>
&gt;&gt; turn on traction_perturbation.db_initial. I use a SimpleDB to put a<br>
&gt;&gt; negative traction all along the fault, and a shear traction compatible with<br>
&gt;&gt; the friction coefficients I impose. However, when I run this, I have the<br>
&gt;&gt; following error message from PETSc (I attach the log):<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt; /Users/jolivetinsar/MYBIN/pylith/pylith-1.7.1-darwin-10.6.8/lib/python2.7/site-packages/pylith/problems/Implicit.py:215:step<br>
&gt;&gt;&gt;  -- implicit(info)<br>
&gt;&gt;&gt;  -- Solving equations.<br>
&gt;&gt;&gt;   0 SNES Function norm 2.679589636063e+01<br>
&gt;&gt;&gt;     0 KSP Residual norm 9.637565795900e-01<br>
&gt;&gt;&gt;     1 KSP Residual norm 3.564937135471e-08<br>
&gt;&gt;&gt;     2 KSP Residual norm 4.313963401911e-16<br>
&gt;&gt;&gt;     3 KSP Residual norm 3.050427683646e-16<br>
&gt;&gt;&gt;     4 KSP Residual norm 2.490663081181e-16<br>
&gt;&gt;&gt;     5 KSP Residual norm 2.156976950523e-16<br>
&gt;&gt;&gt;     6 KSP Residual norm 1.929258577738e-16<br>
&gt;&gt;&gt;     7 KSP Residual norm 1.761163907914e-16<br>
&gt;&gt;&gt;     8 KSP Residual norm 1.630520838230e-16<br>
&gt;&gt;&gt;     9 KSP Residual norm 1.525212509706e-16<br>
&gt;&gt;&gt;    10 KSP Residual norm 1.833325155789e-11<br>
&gt;&gt;&gt;    11 KSP Residual norm 7.039784754868e-23<br>
&gt;&gt;&gt;   Linear solve converged due to CONVERGED_ATOL iterations 11<br>
&gt;&gt;&gt; [0]PETSC ERROR: --------------------- Error Message<br>
&gt;&gt; ------------------------------------<br>
&gt;&gt;&gt; [0]PETSC ERROR: Floating point exception!<br>
&gt;&gt;&gt; [0]PETSC ERROR: Infinite or not-a-number generated in norm!<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; If anyone has an idea about that error message....<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; I think it has to do with some slip vector. If I set<br>
&gt;&gt; friction.min_slip_rate=1e-1, this problem does not happen (the nonlinear<br>
&gt;&gt; solver does not converge, but at least, no crash). If I set the normal<br>
&gt;&gt; traction to be positive, it does not crash as well (and it does not<br>
&gt;&gt; converge as well). I join to this message the cfg file I am using for that<br>
&gt;&gt; experiment.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Thank you<br>
&gt;&gt;&gt; Romain<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Romain Jolivet<br>
&gt;&gt;&gt; Postdoctoral Scholar<br>
&gt;&gt;&gt; Geological and Planetary Sciences<br>
&gt;&gt;&gt; California Institute of Technology, Pasadena, CA<br>
&gt;&gt;&gt; <a href="mailto:rjolivet@caltech.edu" target="_blank">rjolivet@caltech.edu</a><br>
&gt;&gt;&gt; <a href="tel:%2B1%20%28626%29%20560%206356" value="+16265606356" target="_blank">+1 (626) 560 6356</a><br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; _______________________________________________<br>
&gt;&gt;&gt; CIG-SHORT mailing list<br>
&gt;&gt;&gt; <a href="mailto:CIG-SHORT@geodynamics.org" target="_blank">CIG-SHORT@geodynamics.org</a><br>
&gt;&gt;&gt; <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br>
&gt;&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; _______________________________________________<br>
&gt;&gt; CIG-SHORT mailing list<br>
&gt;&gt; <a href="mailto:CIG-SHORT@geodynamics.org" target="_blank">CIG-SHORT@geodynamics.org</a><br>
&gt;&gt; <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br>
&gt;&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; What most experimenters take for granted before they begin their<br>
&gt; experiments is infinitely more interesting than any results to which their<br>
&gt; experiments lead.<br>
&gt; -- Norbert Wiener<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; What most experimenters take for granted before they begin their<br>
&gt; experiments is infinitely more interesting than any results to which their<br>
&gt; experiments lead.<br>
&gt; -- Norbert Wiener<br>
</div></div>&gt; -------------- next part --------------<br>
&gt; An HTML attachment was scrubbed...<br>
&gt; URL: <a href="http://geodynamics.org/pipermail/cig-short/attachments/20120829/9e91899c/attachment.html" target="_blank">http://geodynamics.org/pipermail/cig-short/attachments/20120829/9e91899c/attachment.html</a><br>



&gt;<br>
&gt; ------------------------------<br>
<div>&gt;<br>
&gt; _______________________________________________<br>
&gt; CIG-SHORT mailing list<br>
&gt; <a href="mailto:CIG-SHORT@geodynamics.org" target="_blank">CIG-SHORT@geodynamics.org</a><br>
&gt; <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br>
&gt;<br>
&gt;<br>
</div>&gt; End of CIG-SHORT Digest, Vol 72, Issue 21<br>
&gt; *****************************************<br>
<div><div><br>
_______________________________________________<br>
CIG-SHORT mailing list<br>
<a href="mailto:CIG-SHORT@geodynamics.org" target="_blank">CIG-SHORT@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br>
</div></div></blockquote></div><br><br clear="all"><span class="HOEnZb"><font color="#888888"><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>
</font></span></blockquote></div><br></div></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>