Commits

Peter Brune  committed 30454c4

Added test examples for classical AMG

  • Participants
  • Parent commits e5a0faa

Comments (0)

Files changed (3)

File src/ksp/ksp/examples/tutorials/makefile

          ${DIFF} output/ex54_1.out ex54_sa.tmp || echo ${PWD} "\nPossible problem with with ex54_1, diffs above \n======================================"; \
         ${RM} -f ex54_sa.tmp
 
+runex54_Classical:
+	-@${MPIEXEC} -n 4 ./ex54 -ne 109 -alpha 1.e-3 -ksp_monitor_short -ksp_type cg -pc_gamg_type classical -ksp_converged_reason -pc_gamg_coarse_eq_limit 80 -mg_levels_ksp_type chebyshev -mg_levels_ksp_chebyshev_estimate_eigenvalues 0,0.05,0,1.05 -mg_levels_pc_type jacobi > ex54_classical.tmp 2>&1; \
+         ${DIFF} output/ex54_classical.out ex54_classical.tmp || echo ${PWD} "\nPossible problem with with ex54_Classical, diffs above \n======================================"; \
+        #${RM} -f ex54_classical.tmp
+
 runex54f:
 	-@${MPIEXEC} -n 4 ./ex54f -ne 59 -theta 30.0 -epsilon 1.e-1 -ksp_monitor_short -ksp_type cg -pc_type gamg -pc_gamg_type agg -pc_gamg_agg_nsmooths 1 -ksp_converged_reason  -pc_gamg_coarse_eq_limit 80 -blob_center 0.,0. -mat_coarsen_type hem -pc_gamg_square_graph false > ex54f.tmp 2>&1; \
          ${DIFF} output/ex54f.out ex54f.tmp || echo ${PWD} "\nPossible problem with with ex54f, diffs above \n======================================"; \
          ${DIFF} output/ex55_sa.out ex55_sa.tmp || echo ${PWD} "\nPossible problem with with ex55_SA, diffs above \n========================================="; \
         ${RM} -f ex55_sa.tmp
 
+runex55_Classical:
+	-@${MPIEXEC} -n 4 ./ex55 -ne 29 -alpha 1.e-3 -ksp_monitor_short -ksp_type gmres -pc_gamg_type classical -ksp_converged_reason -pc_gamg_coarse_eq_limit 80 -mg_levels_ksp_type chebyshev -mg_levels_pc_type sor -mg_levels_ksp_max_it 5 > ex55_classical.tmp 2>&1; \
+         ${DIFF} output/ex55_classical.out ex55_classical.tmp || echo ${PWD} "\nPossible problem with with ex55_Classical, diffs above \n========================================="; \
+        #${RM} -f ex55_classical.tmp
+
 runex55_NC:
 	-@${MPIEXEC} -n 4 ./ex55 -ne 29 -alpha 1.e-3 -ksp_monitor_short -ksp_type cg -pc_gamg_type agg -pc_gamg_agg_nsmooths 1  -ksp_converged_reason -pc_gamg_coarse_eq_limit 80 -mg_levels_ksp_type chebyshev -mg_levels_ksp_chebyshev_estimate_eigenvalues 0,0.05,0,1.05 -mg_levels_pc_type jacobi > ex55_nc.tmp 2>&1; \
          ${DIFF} output/ex55_NC.out ex55_nc.tmp || echo ${PWD} "\nPossible problem with with ex55_NC, diffs above \n======================================"; \
                                  ex31.PETSc ex31.rm ex32.PETSc runex32 ex32.rm ex34.PETSc runex34 ex34.rm ex38.PETSc runex38 ex38.rm \
                                  ex43.PETSc runex43 runex43_2 runex43_3 ex43.rm \
                                  ex45.PETSc runex45 runex45_2 ex45.rm \
-                                 ex49.PETSc runex49 runex49_2 runex49_3 ex49.rm ex53.PETSc runex53 ex53.rm ex55.PETSc runex55_SA ex55.rm\
+                                 ex49.PETSc runex49 runex49_2 runex49_3 ex49.rm ex53.PETSc runex53 ex53.rm ex55.PETSc runex55_SA runex55_Classical ex55.rm\
                                  ex58.PETSc runex58 runex58_baij runex58_sbaij ex58.rm
 TESTEXAMPLES_C_X	       = ex2.PETSc runex2_5 ex2.rm ex5.PETSc runex5_5 ex5.rm ex8.PETSc ex8.rm ex28.PETSc runex28 ex28.rm
 TESTEXAMPLES_FORTRAN	       = ex1f.PETSc runex1f ex1f.rm ex2f.PETSc runex2f ex2f.rm ex6f.PETSc ex6f.rm \

File src/ksp/ksp/examples/tutorials/output/ex54_classical.out

+  0 KSP Residual norm 394.784 
+  1 KSP Residual norm 30.4511 
+  2 KSP Residual norm 11.5016 
+  3 KSP Residual norm 0.855039 
+  4 KSP Residual norm 0.107642 
+  5 KSP Residual norm 0.0269386 
+  6 KSP Residual norm 0.00322622 
+Linear solve converged due to CONVERGED_RTOL iterations 6

File src/ksp/ksp/examples/tutorials/output/ex55_classical.out

+  0 KSP Residual norm 138.804 
+  1 KSP Residual norm 1.6505 
+  2 KSP Residual norm 1.44786 
+  3 KSP Residual norm 1.34685 
+  4 KSP Residual norm 1.22207 
+  5 KSP Residual norm 0.838791 
+  6 KSP Residual norm 0.533181 
+  7 KSP Residual norm 0.41106 
+  8 KSP Residual norm 0.144899 
+  9 KSP Residual norm 0.0930366 
+ 10 KSP Residual norm 0.0733036 
+ 11 KSP Residual norm 0.0329768 
+ 12 KSP Residual norm 0.018408 
+ 13 KSP Residual norm 0.00581743 
+ 14 KSP Residual norm 0.00366803 
+ 15 KSP Residual norm 0.00137723 
+Linear solve converged due to CONVERGED_RTOL iterations 15