<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_GNHEP);CHKERRQ(ierr);¬<br> ierr = EPSSetWhichEigenpairs(eps,EPS_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><br><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!<br><br></div>Wenbo<br></div>