On Thu, 8 Jan 2015, Jed Brown wrote: > Roy Stogner <royst...@ices.utexas.edu> writes: > >> By "neither seems to do much"? >> >> I mean that if I set either "-ksp_pc_side right" or "-ksp_norm_type >> unpreconditioned", then instead of ending up falsely "converged" at >> 10 KSP preconditioned resid norm 1.580464771307e-13 true resid norm >> 3.084021824157e-04 ||r(i)||/||b|| 2.711608086867e-03 >> I end up at >> 10 KSP unpreconditioned resid norm 4.683505266751e-14 true resid norm >> 3.820485023471e-04 ||r(i)||/||b|| 3.359139032108e-03 > > Interesting. Can you include the full unpreconditioned output? And > please compare with and without -ksp_gmres_modifiedgramschmidt.
>From "./main-opt -ksp_monitor_true_residual -ksp_norm_type unpreconditioned": 0 KSP unpreconditioned resid norm 1.137340546775e-01 true resid norm 1.137340546775e-01 ||r(i)||/||b|| 1.000000000000e+00 1 KSP unpreconditioned resid norm 1.851221771429e-02 true resid norm 1.850940946785e-02 ||r(i)||/||b|| 1.627428963148e-01 2 KSP unpreconditioned resid norm 8.058732139350e-04 true resid norm 8.487837216519e-04 ||r(i)||/||b|| 7.462881052281e-03 3 KSP unpreconditioned resid norm 1.410067346760e-04 true resid norm 3.750759001371e-04 ||r(i)||/||b|| 3.297832836441e-03 4 KSP unpreconditioned resid norm 3.036297998627e-06 true resid norm 3.296862054294e-04 ||r(i)||/||b|| 2.898746609924e-03 5 KSP unpreconditioned resid norm 1.628004374290e-07 true resid norm 4.048903454605e-04 ||r(i)||/||b|| 3.559974596954e-03 6 KSP unpreconditioned resid norm 6.839537794328e-09 true resid norm 2.859920502301e-04 ||r(i)||/||b|| 2.514568315013e-03 7 KSP unpreconditioned resid norm 3.580827398197e-10 true resid norm 3.424888565589e-04 ||r(i)||/||b|| 3.011313168514e-03 8 KSP unpreconditioned resid norm 1.371032175862e-11 true resid norm 3.818671006476e-04 ||r(i)||/||b|| 3.357544068312e-03 9 KSP unpreconditioned resid norm 1.130911552385e-12 true resid norm 2.475276237056e-04 ||r(i)||/||b|| 2.176372102511e-03 10 KSP unpreconditioned resid norm 4.683505266751e-14 true resid norm 3.820485023471e-04 ||r(i)||/||b|| 3.359139032108e-03 And after adding "-ksp_gmres_modifiedgramschmidt" too: 0 KSP unpreconditioned resid norm 1.137340546775e-01 true resid norm 1.137340546775e-01 ||r(i)||/||b|| 1.000000000000e+00 1 KSP unpreconditioned resid norm 1.851221771429e-02 true resid norm 1.850940946785e-02 ||r(i)||/||b|| 1.627428963148e-01 2 KSP unpreconditioned resid norm 8.058732139350e-04 true resid norm 8.290754235563e-04 ||r(i)||/||b|| 7.289596998076e-03 3 KSP unpreconditioned resid norm 1.409277416992e-04 true resid norm 3.366737464418e-04 ||r(i)||/||b|| 2.960184154134e-03 4 KSP unpreconditioned resid norm 3.016318237424e-06 true resid norm 3.904608063041e-04 ||r(i)||/||b|| 3.433103720879e-03 5 KSP unpreconditioned resid norm 1.684322222043e-07 true resid norm 1.987083874964e-04 ||r(i)||/||b|| 1.747131833643e-03 6 KSP unpreconditioned resid norm 7.043704865544e-09 true resid norm 3.021936047142e-04 ||r(i)||/||b|| 2.657019531846e-03 7 KSP unpreconditioned resid norm 3.538690849000e-10 true resid norm 3.937474047247e-04 ||r(i)||/||b|| 3.462000944582e-03 8 KSP unpreconditioned resid norm 2.298664686412e-11 true resid norm 2.516665836236e-04 ||r(i)||/||b|| 2.212763664649e-03 9 KSP unpreconditioned resid norm 8.918329988448e-13 true resid norm 3.883206450876e-04 ||r(i)||/||b|| 3.414286479003e-03 10 KSP unpreconditioned resid norm 3.074816936477e-14 true resid norm 2.987382118533e-04 ||r(i)||/||b|| 2.626638192935e-03 > How big is this matrix? Can you write it with "-ksp_view_mat binary" > and send it to me (or post somewhere)? Not even 100x100; it's a very-well-simplified version of the real problem. http://users.ices.utexas.edu/~roystgnr/binaryoutput > Dammit, these lists of methods are always wrong. We give useful error > messages if it's not supported. I'll fix/delete these doc problems. Thanks! In my head this thread has now been upgraded to "useful bug report", a big step up from "wasting Jed's time". --- Roy ------------------------------------------------------------------------------ Dive into the World of Parallel Programming! The Go Parallel Website, sponsored by Intel and developed in partnership with Slashdot Media, is your hub for all things parallel software development, from weekly thought leadership blogs to news, videos, case studies, tutorials and more. Take a look and join the conversation now. http://goparallel.sourceforge.net _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
