DefaultParams.hpp

Go to the documentation of this file.

const std::string default_options =
    "-order 2 \n"
    " \n"
    "-pc_type lu \n"
    "-pc_factor_mat_solver_type mumps \n"
    "#-ksp_monitor \n"
    " \n"
    "-snes_atol 1e-8 \n"
    "-snes_rtol 1e-8 \n"
    "-snes_max_linear_solve_fail -1 \n"
    "-snes_linesearch_type bt \n"
    "-snes_linesearch_max_it 3 \n"
    "#-snes_monitor \n"
    "-snes_max_it 25 \n"
    " \n"
    "-ts_max_time 2 \n"
    "-ts_max_snes_failures 2 \n"
    "-ts_type beuler \n"
    "-ts_exact_final_time matchstep \n"
    "-ts_dt 0.05 \n"
    "-ts_error_if_step_fails false \n"
    "-ts_monitor \n"
    " \n"
    "#-ts_adapt_monitor \n"
    "#-ts_adapt_type basic \n"
    "#-ts_adapt_basic_safety 0.8 \n"
    "#-ts_adapt_dt_max 5.0 \n"
    "#-ts_adapt_dt_min 0.005 \n"
    " \n"
    "-ts_adapt_always_accept 1 \n"
    "-ts_theta_initial_guess_extrapolate 1 \n"
    "-ts_theta_theta 1 \n"
    "-ts_adapt_reject_safety 0.9 \n"
    " \n"
    " \n"
    "-snes_ksp_ew 0 \n"
    "# plastic strains \n"
    "-pc_type fieldsplit  \n"
    "-pc_fieldsplit_type schur \n"
    " \n"
    "-fieldsplit_0_ksp_type preonly  \n"
    "# -fieldsplit_0_ksp_monitor  \n"
    "-fieldsplit_0_pc_type lu  \n"
    "-fieldsplit_0_pc_factor_mat_solver_type mumps \n"
    "-fieldsplit_0_mat_mumps_icntl_14 1200 \n"
    "-fieldsplit_0_mat_mumps_icntl_24 1 \n"
    "-fieldsplit_0_mat_mumps_icntl_13 1 \n"
    " \n"
    "-fieldsplit_1_pc_type fieldsplit  \n"
    "-fieldsplit_1_pc_fieldsplit_type schur \n"
    "-fieldsplit_1_ksp_type preonly  \n"
    "# -fieldsplit_1_ksp_monitor  \n"
    " \n"
    "# plastic multiplier \n"
    "-fieldsplit_1_fieldsplit_0_pc_type lu \n"
    "-fieldsplit_1_fieldsplit_0_ksp_type preonly  \n"
    "-fieldsplit_1_fieldsplit_0_pc_factor_mat_solver_type mumps \n"
    "-fieldsplit_1_fieldsplit_0_mat_mumps_icntl_14 1200 \n"
    "-fieldsplit_1_fieldsplit_0_mat_mumps_icntl_24 1 \n"
    "-fieldsplit_1_fieldsplit_0_mat_mumps_icntl_13 1 \n"
    "-fieldsplit_0_mat_mumps_icntl_14 1200 \n"
    "-fieldsplit_0_mat_mumps_icntl_24 1 \n"
    "-fieldsplit_0_mat_mumps_icntl_13 1 \n"
    "# -fieldsplit_1_fieldsplit_0_ksp_monitor  \n"
    " \n"
    "# boundary \n"
    "# -fieldsplit_1_fieldsplit_1_pc_type lu  \n"
    "-fieldsplit_1_fieldsplit_1_pc_type fieldsplit  \n"
    "-fieldsplit_1_fieldsplit_1_pc_fieldsplit_type multiplicative \n"
    "-fieldsplit_1_fieldsplit_1_ksp_type fgmres  \n"
    "-fieldsplit_1_fieldsplit_1_pc_factor_mat_solver_type mumps \n"
    "-fieldsplit_1_fieldsplit_1_mat_mumps_icntl_14 1200 \n"
    "-fieldsplit_1_fieldsplit_1_mat_mumps_icntl_24 1 \n"
    "-fieldsplit_1_fieldsplit_1_mat_mumps_icntl_13 1 \n"
    "# -fieldsplit_1_fieldsplit_1_ksp_monitor  \n"
    " \n"
    " \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_0_pc_type lu \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_0_ksp_type fgmres  \n"
    "# -fieldsplit_1_fieldsplit_1_fieldsplit_0_ksp_monitor  \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_0_pc_factor_mat_solver_type mumps \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_0_mat_mumps_icntl_14 1200 \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_0_mat_mumps_icntl_24 1 \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_0_mat_mumps_icntl_13 1 \n"
    " \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_1_pc_type lu \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_1_ksp_type fgmres  \n"
    "# -fieldsplit_1_fieldsplit_1_fieldsplit_1_ksp_monitor  \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_1_pc_factor_mat_solver_type mumps \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_1_mat_mumps_icntl_14 1200 \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_1_mat_mumps_icntl_24 1 \n"
    "-fieldsplit_1_fieldsplit_1_fieldsplit_1_mat_mumps_icntl_13 1 \n";
 
std::string param_file = "param_file.petsc";
 
if (!static_cast<bool>(ifstream(param_file))) {
  std::ofstream file(param_file.c_str(), std::ios::ate);
  if (file.is_open()) {
    file << default_options;
    file.close();
  }
}