[petsc-users] Fwd: same petsc problem
Mark Adams
mfadams at lbl.gov
Mon Jul 13 13:54:25 CDT 2015
Greg, I am forwarding this to the PETSc mailing list.
Please send the entire output from this run. As I recall you were getting
a message that all values were not the same on all processors in GMRES. I
have seen this when I get NaNs in the system.
While you are doing this you should use a simple solver like change:
-pressure_pc_type gamg
to
-pressure_pc_type jacobi
And add:
*-*pressure_*ksp_monitor_true_residual*
Mark
---------- Forwarded message ----------
From: Greg Miller <grgmiller at ucdavis.edu>
Date: Mon, Jul 13, 2015 at 2:08 PM
Subject: same petsc problem
To: Mark Adams <mfadams at lbl.gov>
Cc: David Trebotich <treb at hpcrd.lbl.gov>
Hi Mark. I'm still stuck on the same petsc problem. Would you please try
the attached minimal example and advise me?
I'm running this without MPI:
make DIM=2 DEBUG=TRUE MPI=FALSE USE_PETSC=TRUE test
There is no input file.
Thank you,
Greg
--
Greg Miller
Department of Chemical Engineering and Materials Science
University of California, Davis
One Shields Avenue
Davis, CA 95616
grgmiller at ucdavis.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150713/9ddeca45/attachment-0001.html>
-------------- next part --------------
makefiles+=test
ebase = test
CHOMBO_HOME = ../Chombo/lib
LibNames := EBAMRTimeDependent EBAMRElliptic EBAMRTools Workshop EBTools AMRTimeDependent AMRElliptic AMRTools BoxTools
base_dir = .
src_dirs =
#$(CHOMBO_HOME)/../releasedExamples/EBAMRINS/srcLib
#
include $(CHOMBO_HOME)/mk/Make.example
# valgrind --leak-check=full
# mpiexec -n 4 valgrind -v --leak-check=full --track-origins=yes --log-file=valgrind.foo
-------------- next part --------------
#-help
#-proj_mac_pc_mg_type full
#-proj_mac_pc_mg_cycle_type w
#-proj_mac_pc_type hypre
-proj_mac_pc_type gamg
-proj_mac_pc_gamg_sym_graph true
-proj_mac_pc_gamg_square_graph true
-proj_mac_pc_gamg_type agg
-proj_mac_pc_gamg_repartition false
-proj_mac_pc_gamg_process_eq_limit 100
-proj_mac_pc_gamg_coarse_eq_limit 1000
-proj_mac_pc_gamg_agg_nsmooths 1
-proj_mac_pc_gamg_threshold -.01
#-proj_mac_pc_gamg_verbose 2
#-proj_mac_pc_mg_log
#-proj_mac_pc_gamg_use_agg_gasm true
#-gamg_est_ksp_view
#-gamg_est_ksp_max_it 10
#-gamg_est_ksp_type gmres
#-gamg_est_ksp_gmres_restart 50
-proj_mac_pc_hypre_type boomeramg
-proj_mac_pc_hypre_boomeramg_no_CF
-proj_mac_pc_hypre_boomeramg_agg_nl 1
-proj_mac_pc_hypre_boomeramg_coarsen_type HMIS
-proj_mac_pc_hypre_boomeramg_interp_type ext+i
#-proj_mac_pc_hypre_boomeramg_P_max 4
#-proj_mac_pc_hypre_boomeramg_agg_num_paths 2
#-proj_mac_pc_hypre_boomeramg_print_statistics
#-proj_mac_pc_hypre_boomeramg_relax_type_coarse symmetric-SOR/Jacobi
#-proj_mac_pc_hypre_boomeramg_grid_sweeps_coarse 4
-proj_mac_pc_ml_Threshold 0.01
-proj_mac_pc_ml_EnergyMinimization 2
-proj_mac_pc_ml_PrintLevel 1
-proj_mac_ksp_type gmres
#-proj_mac_ksp_gmres_restart 50
-proj_mac_ksp_max_it 50
#-proj_mac_mg_levels_ksp_type chebyshev
-proj_mac_mg_levels_ksp_type richardson
#-proj_mac_mg_levels_ksp_chebyshev_estimate_eigenvalues 0,0.1,0,1.1
-proj_mac_mg_levels_pc_type sor
-proj_mac_ksp_rtol 1.e-6
-proj_mac_ksp_atol 1.e-30
-proj_mac_mg_levels_ksp_max_it 2
#-mat_ascii_output_large
-options_left
#-mat_view_matlab
#-log_summary
#-info
-visc_ksp_rtol 1.e-6
-visc_pc_type jacobi
-visc_pc_hypre_type boomeramg
-visc_pc_hypre_boomeramg_no CF
-visc_pc_hypre_boomeramg_agg_nl 1
-visc_pc_hypre_boomeramg_coarsen_type HMIS
-visc_pc_hypre_boomeramg_interp_type ext+i
-visc_ksp_type gmres
-visc_ksp_max_it 50
-diff_ksp_rtol 1.e-6
-diff_pc_type jacobi
-diff_ksp_max_it 50
#-ksp_view
#-visc_ksp_view
#-diff_ksp_view
#-proj_mac_ksp_view
#-proj_mac_ksp_monitor
#-proj_mac_ksp_monitor_max
#-visc_ksp_monitor_max
#-proj_mac_ksp_converged_reason
#-visc_ksp_converged_reason
#-diff_ksp_converged_reason
-proj_mac_ksp_norm_type unpreconditioned
-diff_ksp_norm_type unpreconditioned
-visc_ksp_norm_type unpreconditioned
#-proj_mac_ksp_converged_use_initial_residual_norm
-pressure_pc_type gamg
-pressure_pc_gamg_sym_graph true
-pressure_pc_gamg_square_graph true
-pressure_pc_gamg_type agg
-pressure_pc_gamg_repartition false
-pressure_pc_gamg_process_eq_limit 100
-pressure_pc_gamg_coarse_eq_limit 1000
-pressure_pc_gamg_agg_nsmooths 1
-pressure_pc_gamg_threshold -.01
-pressure_ksp_type gmres
-pressure_ksp_gmres_restart 50
-pressure_ksp_max_it 50
-pressure_ksp_rtol 1.e-12
-pressure_ksp_atol 1.e-30
-pressure_ksp_norm_type unpreconditioned
-pressure_mg_levels_pc_type sor
-pressure_mg_levels_ksp_max_it 4
-pressure_ksp_converged_reason
-pressure_ksp_monitor
-pressure_backup_pc_type gamg
-pressure_backup_pc_gamg_sym_graph true
-pressure_backup_pc_gamg_square_graph true
-pressure_backup_pc_gamg_type agg
-pressure_backup_pc_gamg_repartition false
-pressure_backup_pc_gamg_process_eq_limit 100
-pressure_backup_pc_gamg_coarse_eq_limit 1000
-pressure_backup_pc_gamg_agg_nsmooths 1
-pressure_backup_pc_gamg_threshold 0.05
-pressure_backup_ksp_type gmres
-pressure_backup_ksp_gmres_restart 50
-pressure_backup_ksp_max_it 200
-pressure_backup_ksp_rtol 1.e-12
-pressure_backup_ksp_atol 1.e-30
-pressure_backup_ksp_norm_type unpreconditioned
-pressure_backup_mg_levels_pc_type sor
-pressure_backup_mg_levels_ksp_max_it 4
-pressure_backup_ksp_converged_reason
-pressure_backup_ksp_monitor
-viscous_pc_type gamg
-viscous_pc_gamg_sym_graph true
-viscous_pc_gamg_square_graph true
-viscous_pc_gamg_type agg
-viscous_pc_gamg_repartition false
-viscous_pc_gamg_process_eq_limit 100
-viscous_pc_gamg_coarse_eq_limit 1000
-viscous_pc_gamg_agg_nsmooths 1
-viscous_pc_gamg_threshold 0.05
-viscous_ksp_type gmres
-viscous_ksp_gmres_restart 50
-viscous_ksp_max_it 50
-viscous_ksp_rtol 1.e-12
-viscous_ksp_atol 1.e-30
-viscous_ksp_norm_type unpreconditioned
-viscous_mg_levels_pc_type sor
-viscous_mg_levels_ksp_max_it 4
-viscous_ksp_converged_reason
-viscous_ksp_monitor
################################################
-viscous_backup_pc_type sor
#-viscous_backup_pc_gamg_verbose 2
#-viscous_backup_mg_levels_ksp_type chebyshev
#-viscous_backup_pc_type gamg
#-viscous_backup_pc_gamg_sym_graph true
#-viscous_backup_pc_gamg_square_graph true
#-viscous_backup_pc_gamg_type agg
#-viscous_backup_pc_gamg_repartition false
#-viscous_backup_pc_gamg_process_eq_limit 100
#-viscous_backup_pc_gamg_coarse_eq_limit 1000
#-viscous_backup_pc_gamg_agg_nsmooths 1
#-viscous_backup_pc_gamg_threshold 0.05
-viscous_backup_ksp_view
-viscous_backup_ksp_type gmres
-viscous_backup_ksp_gmres_restart 100
-viscous_backup_ksp_max_it 1000
-viscous_backup_ksp_rtol 1.e-12
-viscous_backup_ksp_atol 1.e-30
#-viscous_backup_ksp_norm_type unpreconditioned
#-viscous_backup_mg_levels_pc_type sor
#-viscous_backup_mg_levels_ksp_max_it 4
-viscous_backup_ksp_converged_reason
-viscous_backup_ksp_monitor
-log_summary
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.cpp
Type: text/x-c++src
Size: 11059 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150713/9ddeca45/attachment-0001.cpp>
More information about the petsc-users
mailing list