<div dir="ltr"><div><div>Hi,<br><br></div>There is a heat conduction problem. When superlu_dist is used as a preconditioner, we have random results from different runs. Is there a random algorithm in superlu_dist? If we use ASM or MUMPS as the preconditioner, we then don't have this issue. <br><br></div>run 1:<br><div><br> 0 Nonlinear |R| = 9.447423e+03<br>      0 Linear |R| = 9.447423e+03<br>      1 Linear |R| = 1.013384e-02<br>      2 Linear |R| = 4.020995e-08<br> 1 Nonlinear |R| = 1.404678e-02<br>      0 Linear |R| = 1.404678e-02<br>      1 Linear |R| = 5.104757e-08<br>      2 Linear |R| = 7.699637e-14<br> 2 Nonlinear |R| =<span style="color:rgb(255,0,0)"> 5.106418e-08</span><br><br><br></div><div>run 2:<br><br> 0 Nonlinear |R| = 9.447423e+03<br>      0 Linear |R| = 9.447423e+03<br>      1 Linear |R| = 1.013384e-02<br>      2 Linear |R| = 4.020995e-08<br> 1 Nonlinear |R| = 1.404678e-02<br>      0 Linear |R| = 1.404678e-02<br>      1 Linear |R| = 5.109913e-08<br>      2 Linear |R| = 7.189091e-14<br> 2 Nonlinear |R| = <span style="color:rgb(255,0,0)">5.111591e-08</span><br><br></div><div>run 3:<br><br> 0 Nonlinear |R| = 9.447423e+03<br>      0 Linear |R| = 9.447423e+03<br>      1 Linear |R| = 1.013384e-02<br>      2 Linear |R| = 4.020995e-08<br> 1 Nonlinear |R| = 1.404678e-02<br>      0 Linear |R| = 1.404678e-02<br>      1 Linear |R| = 5.104942e-08<br>      2 Linear |R| = 7.465572e-14<br> 2 Nonlinear |R| = <span style="color:rgb(255,0,0)">5.106642e-08</span><br><br></div><div>run 4:<br><br> 0 Nonlinear |R| = 9.447423e+03<br>      0 Linear |R| = 9.447423e+03<br>      1 Linear |R| = 1.013384e-02<br>      2 Linear |R| = 4.020995e-08<br> 1 Nonlinear |R| = 1.404678e-02<br>      0 Linear |R| = 1.404678e-02<br>      1 Linear |R| = 5.102730e-08<br>      2 Linear |R| = 7.132220e-14<br> 2 Nonlinear |R| = <span style="color:rgb(255,0,0)">5.104442e-08</span><br><br></div><div>Solver details:<br><br>SNES Object: 8 MPI processes<br>  type: newtonls<br>  maximum iterations=15, maximum function evaluations=10000<br>  tolerances: relative=1e-08, absolute=1e-11, solution=1e-50<br>  total number of linear solver iterations=4<br>  total number of function evaluations=7<br>  norm schedule ALWAYS<br>  SNESLineSearch Object: 8 MPI processes<br>    type: basic<br>    maxstep=1.000000e+08, minlambda=1.000000e-12<br>    tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08<br>    maximum iterations=40<br>  KSP Object: 8 MPI processes<br>    type: gmres<br>      restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<br>      happy breakdown tolerance 1e-30<br>    maximum iterations=100, initial guess is zero<br>    tolerances:  relative=1e-06, absolute=1e-50, divergence=10000.<br>    right preconditioning<br>    using UNPRECONDITIONED norm type for convergence test<br>  PC Object: 8 MPI processes<br>    type: lu<br>      out-of-place factorization<br>      tolerance for zero pivot 2.22045e-14<br>      matrix ordering: natural<br>      factor fill ratio given 0., needed 0.<br>        Factored matrix follows:<br>          Mat Object: 8 MPI processes<br>            type: superlu_dist<br>            rows=7925, cols=7925<br>            package used to perform factorization: superlu_dist<br>            total: nonzeros=0, allocated nonzeros=0<br>            total number of mallocs used during MatSetValues calls =0<br>              SuperLU_DIST run parameters:<br>                Process grid nprow 4 x npcol 2 <br>                Equilibrate matrix TRUE <br>                Matrix input mode 1 <br>                Replace tiny pivots FALSE <br>                Use iterative refinement TRUE <br>                Processors in row 4 col partition 2 <br>                Row permutation LargeDiag <br>                Column permutation METIS_AT_PLUS_A<br>                Parallel symbolic factorization FALSE <br>                Repeated factorization SamePattern<br>    linear system matrix followed by preconditioner matrix:<br>    Mat Object: 8 MPI processes<br>      type: mffd<br>      rows=7925, cols=7925<br>        Matrix-free approximation:<br>          err=1.49012e-08 (relative error in function evaluation)<br>          Using wp compute h routine<br>              Does not compute normU<br>    Mat Object: () 8 MPI processes<br>      type: mpiaij<br>      rows=7925, cols=7925<br>      total: nonzeros=63587, allocated nonzeros=63865<br>      total number of mallocs used during MatSetValues calls =0<br>        not using I-node (on process 0) routines<br><br><br></div><div>Fande,<br></div><div><br></div><div><br></div></div>