[petsc-users] Fwd: same petsc problem

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>