diff --git a/NonLinearEVP/NonLinearEVP.pro b/NonLinearEVP/NonLinearEVP.pro index 95cbe18b0c78650d953b100b2fb1411a53604c91..137232a0221550cd0a4becf544dfb6379d590443 100644 --- a/NonLinearEVP/NonLinearEVP.pro +++ b/NonLinearEVP/NonLinearEVP.pro @@ -498,15 +498,15 @@ PostOperation { } eig_tol = 1e-6; -slepc_options_nep_pol = " -nep_max_it 1000 -nep_type nleigs -nep_target_magnitude "; -// slepc_options_nep_rat = Sprintf(" -nep_max_it 1000 -nep_type nleigs -nep_rational -nep_tol %.10f -nep_view ",eig_tol); -slepc_options_nep_rat = Sprintf(" -nep_type nleigs -nep_rational "); - +// slepc_options_nep_rat = Sprintf(" -nep_max_it 100 -nep_type nleigs -nep_rational -nep_tol %.10f -nep_view ",eig_tol); +// slepc_options_nep_rat = Sprintf(" -nep_type nleigs -nep_rational -nep_view "); +slepc_options_nep_rat = Sprintf(" -nep_view "); // // Warning petsc > 3.9 st_pc_factor_mat_solver_package => st_pc_factor_mat_solver_type -slepc_options_st = " -solve -st_type sinvert -st_ksp_type preonly -st_pc_type lu -st_pc_factor_mat_solver_type mumps -petsc_prealloc 200 -pos"; +// TODO all this should now be default +slepc_options_st = " -st_type sinvert -st_ksp_type preonly -st_pc_type lu -st_pc_factor_mat_solver_type mumps "; -slepc_options_pep = Sprintf(" -pep_max_it 400 -pep_target_magnitude -pep_tol %.10f ",eig_tol); +slepc_options_pep = Sprintf(" -pep_max_it 400 -pep_target_magnitude -pep_tol %.10f -pep_view ",eig_tol); If (flag_res==0) which_res="Aux_E"; @@ -543,5 +543,5 @@ EndIf DefineConstant[ R_ = {StrCat("res_",which_res), Name "GetDP/1ResolutionChoices", Visible 1}, - C_ = {Sprintf(StrCat("-solve -pos -slepc ",str_slepc_options,slepc_options_rg)), Name "GetDP/9ComputeCommand", Visible 1}, + C_ = {Sprintf(StrCat("-solve -pos -slepc -petsc_prealloc 200 ",str_slepc_options,slepc_options_rg)), Name "GetDP/9ComputeCommand", Visible 1}, P_ = {StrCat("postop_",which_res), Name "GetDP/2PostOperationChoices", Visible 1}];