DefaultParams.hpp

Go to the documentation of this file.

const std::string default_options = "-order 2 \n"
                                    "-load_history load_history.in \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 basic \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";
 
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();
  }
}