<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Sep 9, 2015 at 5:29 PM, Gideon Simpson <span dir="ltr"><<a href="mailto:gideon.simpson@gmail.com" target="_blank">gideon.simpson@gmail.com</a>></span> wrote:<br><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 a fresh install of 3.6.1 against gcc 5.2 with —download-mpi and —download-f2cblas with quad math support. I then ran my code with the solver flags:<div><span class=""><div>-snes_max_it 100 -snes_converged_reason -snes_monitor \</div><div><span style="white-space:pre-wrap"> </span>-pc_type lu \</div></span><div><span style="white-space:pre-wrap"> </span>-snes_max_funcs 500000 -snes_fd -snes_monitor -ksp_monitor_true_residual -binary_read_double</div><div><br></div><div>I get similar results as before, but, clearly, at higher precision:</div></div></div></blockquote><div><br></div><div>Great! The Newton direction is just plain bad. We need to precondition it :) So the first thing to do is find</div><div>the "problem" variables (or the possible coarse variables). We do this by looking at the residual. This should</div><div>be fun.</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><div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 SNES Function norm 2.528207888951e-12 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.518112956043e-12 true resid norm 2.528207888951e-12 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 0.000000000000e+00 true resid norm 1.904732438038e-41 ||r(i)||/||b|| 7.533923323167e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 SNES Function norm 2.299665925647e-23 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 1</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 SNES Function norm 5.066222213176e+03 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.135087050102e+01 true resid norm 5.066222213176e+03 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 2.700962799224e-30 true resid norm 5.578720034265e-27 ||r(i)||/||b|| 1.101159759585e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 SNES Function norm 8.484697184455e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 2.668083368285e+02 true resid norm 8.484697184455e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 8.794084418808e-29 true resid norm 2.437641286218e-26 ||r(i)||/||b|| 2.872985603639e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 2 SNES Function norm 6.549559723473e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.552154272986e+02 true resid norm 6.549559723473e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 2.507070573051e-29 true resid norm 4.361978496405e-27 ||r(i)||/||b|| 6.659956822398e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 3 SNES Function norm 5.770723278332e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.348596652974e+02 true resid norm 5.770723278332e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 3.578833635012e-28 true resid norm 1.182025127733e-26 ||r(i)||/||b|| 2.048313652764e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 4 SNES Function norm 5.237702240779e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.524401223281e+02 true resid norm 5.237702240779e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.348074578423e-29 true resid norm 3.851281508331e-27 ||r(i)||/||b|| 7.352998187537e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 5 SNES Function norm 4.753909020039e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.864261890236e+02 true resid norm 4.753909020039e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 7.878414198958e-29 true resid norm 5.419935080688e-26 ||r(i)||/||b|| 1.140100716661e-28</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 6 SNES Function norm 4.221784590918e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 2.550992139992e+02 true resid norm 4.221784590918e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 4.075971594113e-29 true resid norm 7.094679472303e-27 ||r(i)||/||b|| 1.680493004680e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 7 SNES Function norm 3.806525080588e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 5.183605289859e+02 true resid norm 3.806525080588e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.103573513147e-28 true resid norm 2.144465704039e-26 ||r(i)||/||b|| 5.633657098375e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 8 SNES Function norm 3.762054656159e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 5.320776223634e+03 true resid norm 3.762054656159e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 2.585030728289e-27 true resid norm 5.530627313578e-25 ||r(i)||/||b|| 1.470108177329e-27</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 9 SNES Function norm 3.758975227011e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 6.252711433929e+03 true resid norm 3.758975227011e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 6.414696869631e-28 true resid norm 2.875693139366e-25 ||r(i)||/||b|| 7.650205084359e-28</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 10 SNES Function norm 3.757032042685e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 5.890497158679e+02 true resid norm 3.757032042685e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.912199122959e-28 true resid norm 4.765625322483e-26 ||r(i)||/||b|| 1.268454798452e-28</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 11 SNES Function norm 3.728798164350e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.081356529363e+03 true resid norm 3.728798164350e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.465322278036e-27 true resid norm 8.948779168104e-26 ||r(i)||/||b|| 2.399909776201e-28</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 12 SNES Function norm 3.723078740346e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.745127665225e+03 true resid norm 3.723078740346e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 6.682508949012e-28 true resid norm 2.402288617952e-26 ||r(i)||/||b|| 6.452424956579e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 13 SNES Function norm 3.721848058457e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.495477557738e+03 true resid norm 3.721848058457e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 2.389031769454e-28 true resid norm 4.589073952100e-26 ||r(i)||/||b|| 1.233009483467e-28</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 14 SNES Function norm 3.720227574913e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.265999878583e+04 true resid norm 3.720227574913e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 3.936728612764e-27 true resid norm 6.774363900864e-26 ||r(i)||/||b|| 1.820954165962e-28</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 15 SNES Function norm 3.720051997704e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 2.912354770053e+03 true resid norm 3.720051997704e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.193330450218e-26 true resid norm 9.007685186918e-25 ||r(i)||/||b|| 2.421386903322e-27</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 16 SNES Function norm 3.718945465838e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.082057529463e+03 true resid norm 3.718945465838e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.143696537658e-28 true resid norm 1.173159164030e-26 ||r(i)||/||b|| 3.154547908288e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 17 SNES Function norm 3.700412588801e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.676198400599e+01 true resid norm 3.700412588801e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 3.884191202285e-29 true resid norm 8.141362804702e-28 ||r(i)||/||b|| 2.200122988810e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 18 SNES Function norm 3.351964193518e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.071494896915e+02 true resid norm 3.351964193518e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 5.159408051956e-29 true resid norm 3.785844889574e-27 ||r(i)||/||b|| 1.129440731168e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 19 SNES Function norm 3.096020694819e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.361110777562e+02 true resid norm 3.096020694819e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 3.067163043170e-28 true resid norm 3.097548865490e-26 ||r(i)||/||b|| 1.000493591879e-28</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 20 SNES Function norm 3.008463400554e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 7.151989169521e+01 true resid norm 3.008463400554e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 2.963021196462e-29 true resid norm 3.267842047327e-27 ||r(i)||/||b|| 1.086216321171e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 21 SNES Function norm 2.752372240183e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.737818994441e+02 true resid norm 2.752372240183e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.355051821355e-29 true resid norm 1.159582173983e-27 ||r(i)||/||b|| 4.213028154600e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 22 SNES Function norm 2.707638991136e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 4.095365990261e+02 true resid norm 2.707638991136e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 3.313462531085e-29 true resid norm 2.326415447553e-27 ||r(i)||/||b|| 8.592044416441e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 23 SNES Function norm 2.698513655244e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 7.571732975899e+01 true resid norm 2.698513655244e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 4.706642213758e-29 true resid norm 5.275659070101e-28 ||r(i)||/||b|| 1.955024040679e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 24 SNES Function norm 2.617973287453e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.781133707863e+02 true resid norm 2.617973287453e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 7.624863846438e-29 true resid norm 1.041367251217e-27 ||r(i)||/||b|| 3.977761179642e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 25 SNES Function norm 2.600135095660e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.377746106386e+02 true resid norm 2.600135095660e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.148291819547e-28 true resid norm 1.521069513901e-27 ||r(i)||/||b|| 5.849963397822e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 26 SNES Function norm 2.584547548464e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.778914250872e+02 true resid norm 2.584547548464e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.300336366195e-28 true resid norm 2.069166109298e-27 ||r(i)||/||b|| 8.005912332807e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 27 SNES Function norm 2.570427136855e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 9.377626460831e+01 true resid norm 2.570427136855e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 5.766637486688e-29 true resid norm 5.436682842489e-28 ||r(i)||/||b|| 2.115089264557e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 28 SNES Function norm 2.553675835085e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.196813637571e+02 true resid norm 2.553675835085e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 5.412882946976e-28 true resid norm 4.755028260087e-26 ||r(i)||/||b|| 1.862032837041e-28</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 29 SNES Function norm 2.543107878046e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.373192219793e+02 true resid norm 2.543107878046e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 3.620724126195e-29 true resid norm 2.234654091065e-27 ||r(i)||/||b|| 8.787099085950e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 30 SNES Function norm 2.539436050245e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.351117240495e+02 true resid norm 2.539436050245e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.229207196291e-29 true resid norm 3.225569644851e-27 ||r(i)||/||b|| 1.270191326354e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 31 SNES Function norm 2.536651625665e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.687817755461e+02 true resid norm 2.536651625665e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 3.233376086959e-29 true resid norm 2.522505835235e-27 ||r(i)||/||b|| 9.944234398264e-30</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 32 SNES Function norm 2.536237446853e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 6.074223711615e+02 true resid norm 2.536237446853e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 2.047078682230e-28 true resid norm 5.038976546196e-27 ||r(i)||/||b|| 1.986792109094e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 33 SNES Function norm 2.535658772687e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 8.774992658126e+02 true resid norm 2.535658772687e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.036350427238e-28 true resid norm 4.058678859635e-27 ||r(i)||/||b|| 1.600640789428e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 34 SNES Function norm 2.535415069989e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.619728779960e+03 true resid norm 2.535415069989e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.173925947581e-28 true resid norm 1.460901190689e-26 ||r(i)||/||b|| 5.761980387282e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 35 SNES Function norm 2.535362549500e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 4.120252504215e+03 true resid norm 2.535362549500e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.482004320762e-27 true resid norm 2.515988923022e-26 ||r(i)||/||b|| 9.923586366447e-29</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 36 SNES Function norm 2.535350372131e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 7.065615798134e+03 true resid norm 2.535350372131e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.267231805260e-27 true resid norm 9.127341716497e-26 ||r(i)||/||b|| 3.600031702452e-28</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 37 SNES Function norm 2.535346214155e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.200630987401e+04 true resid norm 2.535346214155e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 6.331744826877e-27 true resid norm 3.434648364475e-25 ||r(i)||/||b|| 1.354705856462e-27</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 38 SNES Function norm 2.535344771151e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 2.045653956385e+04 true resid norm 2.535344771151e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 2.729474459326e-27 true resid norm 4.809771115033e-25 ||r(i)||/||b|| 1.897087595250e-27</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 39 SNES Function norm 2.535344273428e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.496422852684e+04 true resid norm 2.535344273428e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.413460430546e-26 true resid norm 2.800309113539e-25 ||r(i)||/||b|| 1.104508426287e-27</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 40 SNES Function norm 2.535344102918e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 5.986286390122e+04 true resid norm 2.535344102918e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.071099626460e-25 true resid norm 2.269549087060e-24 ||r(i)||/||b|| 8.951641256299e-27</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 41 SNES Function norm 2.535344046223e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.158888344932e+05 true resid norm 2.535344046223e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 9.362165764120e-26 true resid norm 5.723929359721e-24 ||r(i)||/||b|| 2.257653894448e-26</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 42 SNES Function norm 2.535344044297e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 6.015343562155e+05 true resid norm 2.535344044297e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 3.045273774078e-24 true resid norm 5.965631634103e-23 ||r(i)||/||b|| 2.352987022618e-25</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 43 SNES Function norm 2.535344043720e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.032398569969e+06 true resid norm 2.535344043720e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 6.133229797301e-24 true resid norm 2.132052198244e-22 ||r(i)||/||b|| 8.409321028933e-25</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 44 SNES Function norm 2.535344043525e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.771990642532e+06 true resid norm 2.535344043525e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 1.399736530559e-23 true resid norm 2.816123472957e-22 ||r(i)||/||b|| 1.110746085980e-24</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 45 SNES Function norm 2.535344043458e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 3.041524303674e+06 true resid norm 2.535344043458e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 6.700523240131e-24 true resid norm 1.214940402127e-21 ||r(i)||/||b|| 4.792013948806e-24</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 46 SNES Function norm 2.535344043436e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 5.220719804889e+06 true resid norm 2.535344043436e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 7.749472721038e-24 true resid norm 9.462448677416e-22 ||r(i)||/||b|| 3.732214845522e-24</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 47 SNES Function norm 2.535344043428e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 8.961379294572e+06 true resid norm 2.535344043428e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 2.735576070962e-23 true resid norm 3.261102152175e-21 ||r(i)||/||b|| 1.286256261997e-23</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 48 SNES Function norm 2.535344043425e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 1.538234198056e+07 true resid norm 2.535344043425e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 5.924165850746e-23 true resid norm 3.890747245316e-21 ||r(i)||/||b|| 1.534603264360e-23</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 49 SNES Function norm 2.535344043424e+02 </span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 0 KSP preconditioned resid norm 2.640413180607e+07 true resid norm 2.535344043424e+02 ||r(i)||/||b|| 1.000000000000e+00</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px"> 1 KSP preconditioned resid norm 2.163013019454e-23 true resid norm 5.034396892612e-21 ||r(i)||/||b|| 1.985685889719e-23</span></font></div><div style="margin:0px"><font face="Menlo"><span style="font-size:12px">Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 49</span></font></div><span class="HOEnZb"><font color="#888888"><div>
<span style="border-collapse:separate;border-spacing:0px">-gideon</span>
</div></font></span><div><div class="h5">
<br><div><blockquote type="cite"><div>On Sep 9, 2015, at 3:15 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>> wrote:</div><br><div><br> Thanks. I think you should try running this with __float128 (quad precision) to see what happens.<br><br> Do you have access to a relatively recent (last couple of years) gcc (and gfortran) compiler?<br><br> Select a new PETSC_ARCH name say arch-quad<br><br> and run configure with the additional option -with-precision=__float128 <br><br> not use --download-mumps or other such packages since they don't work with __float128. If you get everything build then run<br>with the same options except do not use pc_factor_mat_solver_package mumps just use the default PETSc solver LU<br><br> Note that if your program is reading in binary files this will be slightly tricky since you will need to convert the result to __float128 format. If you use PETSc VecLoad() or PetscBinaryRead() etc you can just use the flag -binary_read_double and will read from a double binary file and automatically put the results into the PETSc __float128 matrices.<br><br> Interpreting results. You may get convergence which might likely indicate the scaling and conditioning of your problem are very nasty but can be overcome with more precision or you may get the same lack of convergence you get now and I don't have a good explanation for that :-)<br><br> Note: it could also be that your SNES function evaluation has a problem.<br><br><br><br> Barry<br><br><blockquote type="cite">On Sep 9, 2015, at 1:48 PM, Gideon Simpson <<a href="mailto:gideon.simpson@gmail.com" target="_blank">gideon.simpson@gmail.com</a>> wrote:<br><br>So here’s one grid sequence with these solver flags:<br>-snes_max_it 100 -snes_converged_reason -snes_monitor \<br>-pc_type lu \<br>-pc_factor_mat_solver_package mumps -mat_mumps_icntl_28 1 -mat_mumps_icntl_7 6 \<br>-snes_max_funcs 500000 -snes_fd -snes_monitor -ksp_monitor_true_residual<br><br>By the way, the problem, I believe is approaching something that is singular, so it is not unreasonable to think there is some kind of inherent ill conditioning in the problem.<br><br> 0 SNES Function norm 2.977789375068e-12 <br> 0 KSP preconditioned resid norm 9.487609584353e-12 true resid norm 2.977789375068e-12 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 3.713064216887e-24 true resid norm 9.247299640119e-23 ||r(i)||/||b|| 3.105424351885e-11<br> 1 SNES Function norm 2.511921299174e-12 <br> Nonlinear solve converged due to CONVERGED_SNORM_RELATIVE iterations 1<br> 0 SNES Function norm 5.066222213175e+03 <br> 0 KSP preconditioned resid norm 3.135073899423e+01 true resid norm 5.066222213175e+03 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 4.667575811795e-11 true resid norm 2.582648659172e-07 ||r(i)||/||b|| 5.097780062738e-11<br> 1 SNES Function norm 8.484664416715e+02 <br> 0 KSP preconditioned resid norm 2.668033649331e+02 true resid norm 8.484664416715e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 3.096350794450e-10 true resid norm 1.126402252137e-08 ||r(i)||/||b|| 1.327574311505e-11<br> 2 SNES Function norm 6.549557200006e+02 <br> 0 KSP preconditioned resid norm 1.552337314505e+02 true resid norm 6.549557200006e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.781790941843e-09 true resid norm 5.596253874028e-08 ||r(i)||/||b|| 8.544476677024e-11<br> 3 SNES Function norm 5.770704424482e+02 <br> 0 KSP preconditioned resid norm 1.348464485234e+02 true resid norm 5.770704424482e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.909040505447e-10 true resid norm 6.992127254821e-09 ||r(i)||/||b|| 1.211659225719e-11<br> 4 SNES Function norm 5.237674957045e+02 <br> 0 KSP preconditioned resid norm 1.524393824592e+02 true resid norm 5.237674957045e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.664845931493e-10 true resid norm 5.713325641481e-09 ||r(i)||/||b|| 1.090813326206e-11<br> 5 SNES Function norm 4.753878819749e+02 <br> 0 KSP preconditioned resid norm 1.864281909532e+02 true resid norm 4.753878819749e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.409973828386e-10 true resid norm 4.289369979836e-09 ||r(i)||/||b|| 9.022884559061e-12<br> 6 SNES Function norm 4.221755368773e+02 <br> 0 KSP preconditioned resid norm 2.551073925659e+02 true resid norm 4.221755368773e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 5.540011959843e-10 true resid norm 1.441047729187e-08 ||r(i)||/||b|| 3.413385199545e-11<br> 7 SNES Function norm 3.806522663388e+02 <br> 0 KSP preconditioned resid norm 5.183358055470e+02 true resid norm 3.806522663388e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 2.063740790885e-09 true resid norm 3.530342872323e-08 ||r(i)||/||b|| 9.274456464634e-11<br> 8 SNES Function norm 3.762033603907e+02 <br> 0 KSP preconditioned resid norm 5.324397409227e+03 true resid norm 3.762033603907e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 6.352530213970e-09 true resid norm 9.590500220593e-08 ||r(i)||/||b|| 2.549286165502e-10<br> 9 SNES Function norm 3.758953286234e+02 <br> 0 KSP preconditioned resid norm 6.246847657999e+03 true resid norm 3.758953286234e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 5.367419000872e-08 true resid norm 8.213904234688e-07 ||r(i)||/||b|| 2.185157305564e-09<br> 10 SNES Function norm 3.756870388529e+02 <br> 0 KSP preconditioned resid norm 5.901626862375e+02 true resid norm 3.756870388529e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.330043229228e-09 true resid norm 1.418363475093e-08 ||r(i)||/||b|| 3.775385702481e-11<br> 11 SNES Function norm 3.728715884883e+02 <br> 0 KSP preconditioned resid norm 1.069749399855e+03 true resid norm 3.728715884883e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 2.911440821785e-10 true resid norm 9.552819316154e-09 ||r(i)||/||b|| 2.561959562241e-11<br> 12 SNES Function norm 3.722328158041e+02 <br> 0 KSP preconditioned resid norm 3.064673278570e+03 true resid norm 3.722328158041e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.105090853307e-08 true resid norm 1.085271466536e-07 ||r(i)||/||b|| 2.915571707969e-10<br> 13 SNES Function norm 3.720989356426e+02 <br> 0 KSP preconditioned resid norm 9.068377570178e+03 true resid norm 3.720989356426e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 3.784084466576e-08 true resid norm 3.777243066339e-07 ||r(i)||/||b|| 1.015117944322e-09<br> 14 SNES Function norm 3.720707403757e+02 <br> 0 KSP preconditioned resid norm 7.198056009453e+04 true resid norm 3.720707403757e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 7.120525808630e-08 true resid norm 1.054867000686e-06 ||r(i)||/||b|| 2.835124846476e-09<br> 15 SNES Function norm 3.720707293983e+02 <br> 0 KSP preconditioned resid norm 4.134314655599e+03 true resid norm 3.720707293983e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 2.003141582850e-08 true resid norm 2.046203596859e-07 ||r(i)||/||b|| 5.499501667784e-10<br> 16 SNES Function norm 3.719027845818e+02 <br> 0 KSP preconditioned resid norm 2.545115269075e+03 true resid norm 3.719027845818e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 4.965948884734e-08 true resid norm 5.071610656715e-07 ||r(i)||/||b|| 1.363692574235e-09<br> 17 SNES Function norm 3.714400762925e+02 <br> 0 KSP preconditioned resid norm 1.109601826815e+03 true resid norm 3.714400762925e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 6.169052076803e-09 true resid norm 6.563603648295e-08 ||r(i)||/||b|| 1.767069324832e-10<br> 18 SNES Function norm 3.701626928766e+02 <br> 0 KSP preconditioned resid norm 2.447547395077e+02 true resid norm 3.701626928766e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.250579098159e-09 true resid norm 1.218652092692e-08 ||r(i)||/||b|| 3.292206686799e-11<br> 19 SNES Function norm 3.556442798407e+02 <br> 0 KSP preconditioned resid norm 1.418496111605e+02 true resid norm 3.556442798407e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 3.627252384359e-10 true resid norm 1.951334080585e-09 ||r(i)||/||b|| 5.486757952241e-12<br> 20 SNES Function norm 3.496125058841e+02 <br> 0 KSP preconditioned resid norm 6.543042326746e+01 true resid norm 3.496125058841e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.377371559702e-10 true resid norm 1.295256556885e-09 ||r(i)||/||b|| 3.704834738706e-12<br> 21 SNES Function norm 3.146776071310e+02 <br> 0 KSP preconditioned resid norm 4.667517635820e+01 true resid norm 3.146776071310e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 2.743789267634e-09 true resid norm 1.644758285998e-08 ||r(i)||/||b|| 5.226804350630e-11<br> 22 SNES Function norm 2.855422195269e+02 <br> 0 KSP preconditioned resid norm 2.243742065865e+02 true resid norm 2.855422195269e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.047816533727e-10 true resid norm 2.299887594213e-09 ||r(i)||/||b|| 8.054457228859e-12<br> 23 SNES Function norm 2.806777486780e+02 <br> 0 KSP preconditioned resid norm 1.561899503179e+03 true resid norm 2.806777486780e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 3.972680147986e-09 true resid norm 7.239772649903e-08 ||r(i)||/||b|| 2.579389596790e-10<br> 24 SNES Function norm 2.805893466480e+02 <br> 0 KSP preconditioned resid norm 1.722111652119e+03 true resid norm 2.805893466480e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 5.170553469009e-09 true resid norm 6.273773919075e-08 ||r(i)||/||b|| 2.235927341513e-10<br> 25 SNES Function norm 2.805171735939e+02 <br> 0 KSP preconditioned resid norm 8.024869901472e+03 true resid norm 2.805171735939e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 2.741187759948e-08 true resid norm 3.211723859501e-07 ||r(i)||/||b|| 1.144929495173e-09<br> 26 SNES Function norm 2.805136061483e+02 <br> 0 KSP preconditioned resid norm 3.239737262200e+04 true resid norm 2.805136061483e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 2.317708795749e-08 true resid norm 3.672739850558e-07 ||r(i)||/||b|| 1.309291160949e-09<br> 27 SNES Function norm 2.805134008969e+02 <br> 0 KSP preconditioned resid norm 6.920567009699e+04 true resid norm 2.805134008969e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.101785240910e-07 true resid norm 1.523556762078e-06 ||r(i)||/||b|| 5.431315428092e-09<br> 28 SNES Function norm 2.805133530790e+02 <br> 0 KSP preconditioned resid norm 3.265182256725e+05 true resid norm 2.805133530790e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 3.903754744633e-07 true resid norm 6.654431764515e-06 ||r(i)||/||b|| 2.372233510980e-08<br> 29 SNES Function norm 2.805133509342e+02 <br> 0 KSP preconditioned resid norm 1.312587966212e+06 true resid norm 2.805133509342e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 1.990958223875e-07 true resid norm 1.311133369898e-05 ||r(i)||/||b|| 4.674049793106e-08<br> 30 SNES Function norm 2.805133508006e+02 <br> 0 KSP preconditioned resid norm 2.049236485623e+07 true resid norm 2.805133508006e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 3.432837714993e-05 true resid norm 4.348842776781e-04 ||r(i)||/||b|| 1.550315792232e-06<br> 31 SNES Function norm 2.805133508001e+02 <br> 0 KSP preconditioned resid norm 7.694877106483e+06 true resid norm 2.805133508001e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 4.401952388358e-06 true resid norm 1.008598135500e-04 ||r(i)||/||b|| 3.595544142990e-07<br> 32 SNES Function norm 2.805133507963e+02 <br> 0 KSP preconditioned resid norm 1.116553807875e+08 true resid norm 2.805133507963e+02 ||r(i)||/||b|| 1.000000000000e+00<br> 1 KSP preconditioned resid norm 7.718443445086e-06 true resid norm 1.327187421086e-03 ||r(i)||/||b|| 4.731280765492e-06<br>Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 32<br><br>-gideon<br><br><blockquote type="cite">On Sep 9, 2015, at 1:52 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>> wrote:<br><br><br> Please run with the -snes_fd and one level of grid sequencing (where Newton was having trouble converging on the second level) with -snes_monitor -ksp_monitor_true_residual and send the output. Too much speculation and not enough data.<br><br></blockquote><br></blockquote><br></div></blockquote></div><br></div></div></div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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</div>
</div></div>