<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">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:0 0 0 .8ex;border-left:1px #ccc solid;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_<wbr>GNHEP);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><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> <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;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="HOEnZb"><font color="#888888"><br><br></font></span></div><span class="HOEnZb"><font color="#888888">Wenbo<br></font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="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>