<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Oct 2, 2017 at 11:21 AM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Yea, it fails in the eigen estimator, but the Cheby eigen estimator works in the solve that works:<div><br></div><div><div>        eigenvalue estimates used:  min = 0.100004, max = 1.10004</div><div>        eigenvalues estimate via gmres min 0.0118548, max 1.00004</div></div><div><br></div><div>Why would it just give "KSPSolve has not converged". It is not supposed to converge ...</div></div></blockquote><div><br></div><div>This sounds like a mistake with</div><div><br></div><div>  <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/KSP/KSPSetErrorIfNotConverged.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/KSP/KSPSetErrorIfNotConverged.html</a></div><div><br></div><div>somewhere.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="gmail_extra"><div class="gmail_quote"><span class="gmail-">On Mon, Oct 2, 2017 at 11:11 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br></span><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="gmail-"><span>On Mon, Oct 2, 2017 at 11:06 AM, Wenbo Zhao <span dir="ltr"><<a href="mailto:zhaowenbo.npic@gmail.com" target="_blank">zhaowenbo.npic@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Matt,<br></div></div></div></blockquote><div><br></div></span></span><div>Thanks Wenbo.</div><span class="gmail-"><span><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Test 1 nonsmooth<br>zhaowenbo@ubuntu:~/test_slepc/<wbr>SPARK/spark$ make NCORE=1 runkr_nonsmooth <br>mpirun -n 1 ./step-41 \<br>   -st_ksp_type gmres  -st_ksp_view -st_ksp_monitor  \<br>   -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 0 \<br>   -mata AMAT.dat -matb BMAT.dat \<br>   -eps_nev 1 -eps_ncv 10  -eps_monitor  -log_view > log_nonsmooth 2>&1<br><br>Test 2 smooth<br>zhaowenbo@ubuntu:~/test_slepc/<wbr>SPARK/spark$ make NCORE=1 runkr_smooth <br>mpirun -n 1 ./step-41 \<br>   -st_ksp_type gmres -st_ksp_view -st_ksp_monitor  \<br>   -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 1 \<br>   -mata AMAT.dat -matb BMAT.dat \<br>   -eps_nev 1 -eps_ncv 10  -eps_monitor -log_view > log_smooth 2>&1<br>makefile:43: recipe for target 'runkr_smooth' failed<br>make: *** [runkr_smooth] Error 91<br><br></div></div></div></blockquote><div><br></div></span></span><div>Mark, the solve is not failing, its the construction of the interpolator I think. Check out this stack</div><div><br></div><div><div>[0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--</div><div>[0]PETSC ERROR:   </div><div>[0]PETSC ERROR: KSPSolve has not converged</div><div>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/d<wbr>ocumentation/faq.html</a> for trouble shooting.</div><div>[0]PETSC ERROR: Petsc Release Version 3.8.0, unknown </div><div>[0]PETSC ERROR: ./step-41 on a arch-linux2-c-debug named ubuntu by zhaowenbo Mon Oct  2 08:00:58 2017</div><div>[0]PETSC ERROR: Configure options --with-mpi=1 --with-shared-libraries=1 --with-64-bit-indices=1 --with-debugging=1</div><div>[0]PETSC ERROR: #1 KSPSolve() line 855 in /home/zhaowenbo/research/petsc<wbr>/petsc_git/src/ksp/ksp/interfa<wbr>ce/itfunc.c</div><div>[0]PETSC ERROR: #2 PCGAMGOptProlongator_AGG() line 1186 in /home/zhaowenbo/research/petsc<wbr>/petsc_git/src/ksp/pc/impls/<wbr>gamg/agg.c</div><div>[0]PETSC ERROR: #3 PCSetUp_GAMG() line 528 in /home/zhaowenbo/research/petsc<wbr>/petsc_git/src/ksp/pc/impls/<wbr>gamg/gamg.c</div><div>[0]PETSC ERROR: #4 PCSetUp() line 924 in /home/zhaowenbo/research/petsc<wbr>/petsc_git/src/ksp/pc/interfac<wbr>e/precon.c</div><div>[0]PETSC ERROR: #5 KSPSetUp() line 378 in /home/zhaowenbo/research/petsc<wbr>/petsc_git/src/ksp/ksp/interfa<wbr>ce/itfunc.c</div><div>[0]PETSC ERROR: #6 STSetUp_Shift() line 129 in /home/zhaowenbo/research/slepc<wbr>/slepc_git/src/sys/classes/st/<wbr>impls/shift/shift.c</div><div>[0]PETSC ERROR: #7 STSetUp() line 281 in /home/zhaowenbo/research/slepc<wbr>/slepc_git/src/sys/classes/st/<wbr>interface/stsolve.c</div><div>[0]PETSC ERROR: #8 EPSSetUp() line 273 in /home/zhaowenbo/research/slepc<wbr>/slepc_git/src/eps/interface/<wbr>epssetup.c</div><div>[0]PETSC ERROR: #9 solve_diffusion_3d() line 1029 in src/diffu.c</div><div>[0]PETSC ERROR: #10 main() line 25 in src/main.c</div><div>[0]PETSC ERROR: PETSc Option Table entries:</div><div>[0]PETSC ERROR: -eps_monitor</div><div>[0]PETSC ERROR: -eps_ncv 10</div><div>[0]PETSC ERROR: -eps_nev 1</div><div>[0]PETSC ERROR: -log_view</div><div>[0]PETSC ERROR: -mata AMAT.dat</div><div>[0]PETSC ERROR: -matb BMAT.dat</div><div>[0]PETSC ERROR: -st_ksp_monitor</div><div>[0]PETSC ERROR: -st_ksp_type gmres</div><div>[0]PETSC ERROR: -st_ksp_view</div><div>[0]PETSC ERROR: -st_pc_gamg_agg_nsmooths 1</div><div>[0]PETSC ERROR: -st_pc_gamg_type agg</div><div>[0]PETSC ERROR: -st_pc_type gamg</div><div>[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov-------<wbr>---</div><div>------------------------------<wbr>------------------------------<wbr>--------------</div><div>MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD</div><div>with errorcode 91.</div></div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div><div class="gmail-h5"><div><div class="gmail-m_4862697422351733712h5"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div></div>Thanks,</div><div><br></div>Wenbo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 2, 2017 at 10:48 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095h5">On Mon, Oct 2, 2017 at 10:43 AM, Wenbo Zhao <span dir="ltr"><<a href="mailto:zhaowenbo.npic@gmail.com" target="_blank">zhaowenbo.npic@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Mark,<br></div><br>Thanks for your reply.<br><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 2, 2017 at 9:51 PM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Please send the output with -st_ksp_view and -st_ksp_monitor and we can start to debug it.<div><br></div></div></blockquote><div>Test 1 with nonsmooth and preonly is OK<br>zhaowenbo@ubuntu:~/test_slepc/<wbr>SPARK/spark$ make NCORE=1 runkr_nonsmooth <br>mpirun -n 1 ./step-41 \<br>   -st_ksp_type gmres  -st_ksp_view -st_ksp_monitor  \<br>   -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 0 \<br>   -st_ksp_view  -mata AMAT.dat -matb BMAT.dat \<br>   -st_mg_coarse_ksp_type preonly   -st_mg_coarse_ksp_monitor  \<br>   -eps_nev 1 -eps_ncv 10  -eps_monitor  -log_view > log_nonsmooth 2>&1<br><br>Test 2 smooth and preonly is not OK<br>zhaowenbo@ubuntu:~/test_slepc/<wbr>SPARK/spark$ make NCORE=1 runkr_smooth <br>mpirun -n 1 ./step-41 \<br>   -st_ksp_type gmres -st_ksp_view -st_ksp_monitor  \<br>   -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 1 \<br>   -st_ksp_view  -mata AMAT.dat -matb BMAT.dat \<br>   -st_mg_coarse_ksp_type preonly   -st_mg_coarse_ksp_monitor  \<br>   -eps_nev 1 -eps_ncv 10  -eps_monitor -log_view > log_smooth 2>&1<br>makefile:43: recipe for target 'runkr_smooth' failed<br>make: *** [runkr_smooth] Error 91<br><br>Test 3 nonsmooth and gmres is not OK<br>zhaowenbo@ubuntu:~/test_slepc/<wbr>SPARK/spark$ make NCORE=1 runkr_gmres <br>mpirun -n 1 ./step-41 \<br>   -st_ksp_type gmres  -st_ksp_view -st_ksp_monitor  \<br>   -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 0 \<br>   -st_ksp_view  -mata AMAT.dat -matb BMAT.dat \<br>   -st_mg_coarse_ksp_type gmres  -st_mg_coarse_ksp_monitor -st_mg_coarse_ksp_rtol 1.0e-6 \<br></div></div></div></div></blockquote><div><br></div></div></div><div>DO NOT DO THIS. Please send the output where you do NOTHING to the coarse solver.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095h5"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>   -eps_nev 1 -eps_ncv 10  -eps_monitor  -log_view > log_gmres 2>&1<br>makefile:59: recipe for target 'runkr_gmres' failed<br>make: *** [runkr_gmres] Error 91</div><div><br></div><div>log-files is attached.<br> </div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div></div><div>You mentioned that B is not symmetric. I assume it is elliptic (diffusion). Where does the asymmetry come from?</div></div><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-HOEnZb"><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-h5"><div class="gmail_extra"><br></div></div></div></blockquote><div><br></div><div>It is a two-group diffusion equations, where group denotes neutron enegry discretisation. <br>Matrix B consists of neutron diffusion/leakage term, removal term and minus neutron scatter source term between different energies, when matrix A denotes neutron fission source.<br><br>Diffusion term(Laplace operator) is elliptic and symmetric. Removal term is diagonal only. However scatter term is asymmetry since scatter term from high energy to low energy is far greater than the term from low to high. <br></div><div><br></div><div><br></div><div>Wenbo<br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-HOEnZb"><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-h5"><div class="gmail_extra"><div class="gmail_quote">On Mon, Oct 2, 2017 at 9:39 AM, Wenbo Zhao <span dir="ltr"><<a href="mailto:zhaowenbo.npic@gmail.com" target="_blank">zhaowenbo.npic@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><span><div>Matt,<br></div>Thanks for your reply.<br></span><div class="gmail_extra">For the defalt option doesnt work firstly( -st_ksp_type gmres -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 1), I tried to test those options.</div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815HOEnZb"><font color="#888888"><div class="gmail_extra"><br></div><div class="gmail_extra">Wenbo<br></div></font></span><div><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 2, 2017 at 9:08 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-">On Mon, Oct 2, 2017 at 8:30 AM, Wenbo Zhao <span dir="ltr"><<a href="mailto:zhaowenbo.npic@gmail.com" target="_blank">zhaowenbo.npic@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Matt</div><div><br></div><div>Because I am not clear about what will happen using 'preonly' for large scale problem. <br></div></div></div></blockquote><div><br></div></span><div>The size of the problem has nothing to do with 'preonly'. All it means is to apply a preconditioner without a Krylov solver.</div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>It seems to use a direct solver from below, <br><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/KSP/KSPPREONLY.html" target="_blank">http://www.mcs.anl.gov/petsc/p<wbr>etsc-current/docs/manualpages/<wbr>KSP/KSPPREONLY.html</a></div></div></div></blockquote><div><br></div></span><div>However, I still cannot understand why you would change the default?</div><div><br></div><div>  Matt</div><div><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-h5"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div><br></div>Thanks!<br></div>Wenbo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 2, 2017 at 5:09 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span>On Sun, Oct 1, 2017 at 9:53 PM, Wenbo Zhao <span dir="ltr"><<a href="mailto:zhaowenbo.npic@gmail.com" target="_blank">zhaowenbo.npic@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Matt,<br>Thanks for your reply.<br>It DOES make no sense for this problem.<br>But I am not clear about the 'preonly' option. Which solver is used in preonly? I wonder if 'preonly' is suitable for large scale problem such as 400,000,000 unknowns.<br>So I tried 'gmres' option and found these error messages.<br></div></blockquote><div><br></div></span><div>I mean, why are you setting this at all. Just do not set the coarse solver. The default should work fine.</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278m_-3221509347772759455h5"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Could you give me some suggestions?<br><br>Thanks.<br><br>Wenbo<br><br><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 2, 2017 at 12:34 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span>On Sun, Oct 1, 2017 at 6:49 AM, Wenbo Zhao <span dir="ltr"><<a href="mailto:zhaowenbo.npic@gmail.com" target="_blank">zhaowenbo.npic@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi,<br><br>I met some questions when I use PETSC/SLEPC to solve two-group neutron diffusion equations with finite difference method. The grid is 3*3*3, when DOF on each points is 2. So the matrix size is 54*54.<br>It is generalized eigenvalue problem Ax=\lamda Bx, where B is diagonally dominant matrix but not symmetry.<br>EPS is set as below,<br> ierr = EPSSetProblemType(eps,EPS_GNHE<wbr>P);CHKERRQ(ierr);¬<br> ierr = EPSSetWhichEigenpairs(eps,EPS_<wbr>LARGEST_REAL);CHKERRQ(ierr);¬<br><br>Krylovschur is used as eps sovler. GAMG is used as PC.<br>I tried agg_nsmooths and mg_coarse_ksp_type. Only non-smooths and preonly is OK.<br></div></div></blockquote><div><br></div></span><div>Why are you setting the coarse solver. This makes no sense.</div><div><br></div><div>   Thanks,</div><div><br></div><div>    Matt</div><div><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278m_-3221509347772759455m_-8305666638070075194m_3153904134462254503h5"><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br>Test 1<br>$ make NCORE=1 runkr_nonsmooth <br>mpirun -n 1 ./step-41 \<br>   -st_ksp_type gmres  \<br>   -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 0 \<br>   -st_ksp_view  -mata AMAT.dat -matb BMAT.dat \<br>   -st_mg_coarse_ksp_type preonly   -st_mg_coarse_ksp_monitor  \<br>   -eps_nev 1 -eps_ncv 10  -eps_monitor  -log_view > log_nonsmooth 2>&1<br> <br>Test 2<br>$ make NCORE=1 runkr_smooth <br>mpirun -n 1 ./step-41 \<br>   -st_ksp_type gmres  \<br>   -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 1 \<br>   -st_ksp_view  -mata AMAT.dat -matb BMAT.dat \<br>   -st_mg_coarse_ksp_type preonly   -st_mg_coarse_ksp_monitor  \<br>   -eps_nev 1 -eps_ncv 10  -eps_monitor -log_view > log_smooth 2>&1<br>makefile:43: recipe for target 'runkr_smooth' failed<br>make: *** [runkr_smooth] Error 91<br><br>Test 3<br>$ make NCORE=1 runkr_gmres <br>mpirun -n 1 ./step-41 \<br>   -st_ksp_type gmres  \<br>   -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 0 \<br>   -st_ksp_view  -mata AMAT.dat -matb BMAT.dat \<br>   -st_mg_coarse_ksp_type gmres  -st_mg_coarse_ksp_monitor -st_mg_coarse_ksp_rtol 1.0e-6 \<br>   -eps_nev 1 -eps_ncv 10  -eps_monitor  -log_view > log_gmres 2>&1<br>makefile:59: recipe for target 'runkr_gmres' failed<br>make: *** [runkr_gmres] Error 91<br><br>Log files were attched.<br>The matrix file were also attched as AMAT.dat and BMAT.dat.<br><br>Is it correct? Or something wrong with my code or commad-line?<br><br>Thanks!<span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278m_-3221509347772759455m_-8305666638070075194m_3153904134462254503m_5044481091402600956HOEnZb"><font color="#888888"><br><br></font></span></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278m_-3221509347772759455m_-8305666638070075194m_3153904134462254503m_5044481091402600956HOEnZb"><font color="#888888">Wenbo<br></font></span></div>
</blockquote></div></div></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278m_-3221509347772759455m_-8305666638070075194m_3153904134462254503HOEnZb"><font color="#888888"><br><br clear="all"><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278HOEnZb"><font color="#888888"><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278m_-3221509347772759455m_-8305666638070075194HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278m_-3221509347772759455m_-8305666638070075194m_3153904134462254503m_5044481091402600956gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/%7Emk51/" target="_blank">https://www.cse.buffalo.edu/~k<wbr>nepley/</a><br></div></div></div></div></div>
</font></span></font></span></font></span></div></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278HOEnZb"><font color="#888888">
</font></span></blockquote></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278HOEnZb"><font color="#888888"><br></font></span></div></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278HOEnZb"><font color="#888888">
</font></span></blockquote></div></div></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278HOEnZb"><font color="#888888"><div><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278m_-3221509347772759455h5"><br><br clear="all"><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278m_-3221509347772759455m_-8305666638070075194gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/%7Emk51/" target="_blank">https://www.cse.buffalo.edu/~k<wbr>nepley/</a><br></div></div></div></div></div>
</font></span></div></div></font></span></div></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237HOEnZb"><font color="#888888">
</font></span></blockquote></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237HOEnZb"><font color="#888888"><br></font></span></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237HOEnZb"><font color="#888888">
</font></span></blockquote></div></div></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237HOEnZb"><font color="#888888"><div><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-h5"><br><br clear="all"><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237m_51802801014489142gmail-m_-5744122032765884815m_1885042015433017079gmail-m_-7838439549008716278gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/%7Emk51/" target="_blank">https://www.cse.buffalo.edu/~k<wbr>nepley/</a><br></div></div></div></div></div>
</font></span></div></div></font></span></div></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-HOEnZb"><font color="#888888"><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237HOEnZb"><font color="#888888">
</font></span></font></span></blockquote></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-HOEnZb"><font color="#888888"><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237HOEnZb"><font color="#888888"><br></font></span></font></span></div></div></div></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-HOEnZb"><font color="#888888"><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237HOEnZb"><font color="#888888">
</font></span></font></span></blockquote></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-HOEnZb"><font color="#888888"><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237HOEnZb"><font color="#888888"><br></font></span></font></span></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-HOEnZb"><font color="#888888">
</font></span></div></div></blockquote></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-HOEnZb"><font color="#888888"><br></font></span></div></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-HOEnZb"><font color="#888888">
</font></span></blockquote></div></div></div><span class="gmail-m_4862697422351733712m_1581769487803968623gmail-HOEnZb"><font color="#888888"><div><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095h5"><br><br clear="all"><div><br></div>-- <br><div class="gmail-m_4862697422351733712m_1581769487803968623gmail-m_6855455092711105095m_-4426497621381592237gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~k<wbr>nepley/</a><br></div></div></div></div></div>
</div></div></font></span></div></div>
</blockquote></div><br></div>
</blockquote></div></div></div></div></div><div><div class="gmail-h5"><div><div class="gmail-m_4862697422351733712h5"><br><br clear="all"><div><br></div>-- <br><div class="gmail-m_4862697422351733712m_1581769487803968623gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~k<wbr>nepley/</a><br></div></div></div></div></div>
</div></div></div></div></div></div>
</blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>