all_parameters.h

#ifndef __ALL_PARAMETERS_H__
#define __ALL_PARAMETERS_H__

#include <deal.II/base/conditional_ostream.h>
#include <deal.II/base/parameter_handler.h>
#include <deal.II/base/tensor.h>

#include "parameters.h"
#include "parameters/parameters_ode_solver.h"
#include "parameters/parameters_linear_solver.h"
#include "parameters/parameters_manufactured_convergence_study.h"

#include "parameters/parameters_euler.h"
#include "parameters/parameters_navier_stokes.h"
#include "parameters/parameters_physics_model.h"

#include "parameters/parameters_reduced_order.h"
#include "parameters/parameters_hyper_reduction.h"
#include "parameters/parameters_burgers.h"
#include "parameters/parameters_grid_refinement_study.h"
#include "parameters/parameters_grid_refinement.h"
#include "parameters/parameters_artificial_dissipation.h"
#include "parameters/parameters_limiter.h"
#include "parameters/parameters_flow_solver.h"
#include "parameters/parameters_mesh_adaptation.h"
#include "parameters/parameters_functional.h"
#include "parameters/parameters_time_refinement_study.h"

namespace PHiLiP {
namespace Parameters {

class AllParameters
{
public:
    AllParameters();

    ManufacturedConvergenceStudyParam manufactured_convergence_study_param;
    ODESolverParam ode_solver_param;
    LinearSolverParam linear_solver_param;
    EulerParam euler_param;
    NavierStokesParam navier_stokes_param;
    ReducedOrderModelParam reduced_order_param;
    HyperReductionParam hyper_reduction_param;
    BurgersParam burgers_param;
    PhysicsModelParam physics_model_param;
    GridRefinementStudyParam grid_refinement_study_param;
    ArtificialDissipationParam artificial_dissipation_param;
    LimiterParam limiter_param;
    FlowSolverParam flow_solver_param;
    MeshAdaptationParam mesh_adaptation_param;
    FunctionalParam functional_param;
    TimeRefinementStudyParam time_refinement_study_param;

    unsigned int dimension;

    enum RunType {
        integration_test,
        flow_simulation
        };
    RunType run_type; 

    enum MeshType {
        default_triangulation,
        triangulation,
        parallel_shared_triangulation,
        parallel_distributed_triangulation,
        };
    MeshType mesh_type;
    

    int overintegration;

    bool use_weak_form;

    enum FluxNodes { GL, GLL };
    FluxNodes flux_nodes_type;

    bool use_collocated_nodes;

    bool use_split_form;

    enum TwoPointNumericalFlux { KG, IR, CH, Ra };
    TwoPointNumericalFlux two_point_num_flux_type;

    bool use_curvilinear_split_form;

    bool use_weight_adjusted_mass;

    bool store_residual_cpu_time;


    bool use_periodic_bc;

    bool use_curvilinear_grid;

    bool use_energy;

    bool use_L2_norm;

    //The default ESFR scheme is the Nonlinearly Stable FR where the volume is also reconstructed
    bool use_classical_FR;

    bool store_global_mass_matrix;

    double sipg_penalty_factor;

    bool use_invariant_curl_form;

    bool use_inverse_mass_on_the_fly;

    bool check_valid_metric_Jacobian;

    std::string energy_file;

    int nstate;

    enum TestType {
        run_control,
        grid_refinement_study,
        stability_fr_parameter_range,
        advection_limiter,
        burgers_limiter,
        burgers_energy_stability,
        diffusion_exact_adjoint,
        euler_gaussian_bump,
        euler_gaussian_bump_enthalpy,
        euler_gaussian_bump_adjoint,
        euler_cylinder,
        euler_cylinder_adjoint,
        euler_vortex,
        euler_entropy_waves,
        euler_split_taylor_green,
        taylor_green_scaling,
        burgers_split_form,
        optimization_inverse_manufactured,
        euler_bump_optimization,
        euler_naca_optimization,
        shock_1d,
        euler_naca0012,
        reduced_order,
        unsteady_reduced_order,
        convection_diffusion_periodicity,
        POD_adaptation,
        POD_adaptive_sampling_run,
        adaptive_sampling_testing,
        finite_difference_sensitivity,
        advection_periodicity,
        dual_weighted_residual_mesh_adaptation,
        anisotropic_mesh_adaptation,
        taylor_green_vortex_energy_check,
        taylor_green_vortex_restart_check,
        time_refinement_study,
        time_refinement_study_reference,
        rrk_numerical_entropy_conservation_check,
        euler_entropy_conserving_split_forms_check,
        h_refinement_study_isentropic_vortex,
        khi_robustness,
        naca0012_unsteady_check_quick,
        homogeneous_isotropic_turbulence_initialization_check,
        build_NNLS_problem,
        hyper_reduction_comparison,
        hyper_adaptive_sampling_run,
        hyper_reduction_post_sampling,
        ROM_error_post_sampling,
        HROM_error_post_sampling,
        hyper_adaptive_sampling_new_error,
        halton_sampling_run,
        low_density
    };
    TestType test_type;

    enum PartialDifferentialEquation {
        advection,
        diffusion,
        convection_diffusion,
        advection_vector,
        burgers_inviscid,
        burgers_viscous,
        burgers_rewienski,
        euler,
        mhd,
        navier_stokes,
        physics_model,
    };
    PartialDifferentialEquation pde_type;

    enum ModelType {
        large_eddy_simulation,
        reynolds_averaged_navier_stokes,
    };
    ModelType model_type;

    enum BoundaryType {
        manufactured_dirichlet,
        manufactured_neumann,
        manufactured_inout_flow,
    };

    enum SourceTerm {
        zero,
        manufactured,
    };

    enum ConvectiveNumericalFlux { 
        lax_friedrichs, 
        roe, 
        l2roe, 
        central_flux,
        two_point_flux,
        two_point_flux_with_lax_friedrichs_dissipation,
        two_point_flux_with_roe_dissipation,
        two_point_flux_with_l2roe_dissipation
    };
    ConvectiveNumericalFlux conv_num_flux_type;

    enum DissipativeNumericalFlux { symm_internal_penalty, bassi_rebay_2, central_visc_flux };
    DissipativeNumericalFlux diss_num_flux_type;

    enum Flux_Reconstruction {cDG, cSD, cHU, cNegative, cNegative2, cPlus, c10Thousand, cHULumped, user_specified_value};
    Flux_Reconstruction flux_reconstruction_type;

    double FR_user_specified_correction_parameter_value;

    enum Flux_Reconstruction_Aux {kDG, kSD, kHU, kNegative, kNegative2, kPlus, k10Thousand};
    Flux_Reconstruction_Aux flux_reconstruction_aux_type;

    enum NonPhysicalBehaviorEnum {return_big_number, abort_run, print_warning};
    NonPhysicalBehaviorEnum non_physical_behavior_type;

    std::string solution_vtk_files_directory_name;

    bool output_high_order_grid;

    bool enable_higher_order_vtk_output;

    bool output_face_results_vtk;

    bool do_renumber_dofs;

    enum RenumberDofsType { CuthillMckee };
    RenumberDofsType renumber_dofs_type;

    double matching_surface_jac_det_tolerance;


    static void declare_parameters (dealii::ParameterHandler &prm);


    void parse_parameters (dealii::ParameterHandler &prm);

    //FunctionParser<dim> initial_conditions;
    //BoundaryConditions  boundary_conditions[max_n_boundaries];
protected:
    dealii::ConditionalOStream pcout; 

};  

} // Parameters namespace
} // PHiLiP namespace

#endif