GMGPolar/include/ConfigParser/config_parser.h at 6a9b7334d6ecddfaf5d8cd06950d9c94175e20c2 · SciCompMod/GMGPolar · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
#pragma once

#include <optional>
#include <string>
#include <stdexcept>
#include <memory>

#include "../../include/ConfigParser/cmdline.h"
#include "../../include/Definitions/global_definitions.h"
#include "../../include/PolarGrid/polargrid.h"
#include "../../include/GMGPolar/test_cases.h"
#include "../../include/GMGPolar/igmgpolar.h"
#include "test_selection.h"

class ConfigParser
{
public:
    /* Specifies whether user input is required */
    static constexpr bool OPTIONAL = false;
    static constexpr bool REQUIRED = true;

    ConfigParser();

    bool parse(int argc, char* argv[]);

    // Test Case
    const DomainGeometryVariant& domainGeometry() const;
    const DensityProfileCoefficientsVariant& densityProfileCoefficients() const;
    const BoundaryConditionsVariant& boundaryConditions() const;
    const SourceTerm& sourceTerm() const;
    const ExactSolution& exactSolution() const;
    std::unique_ptr<IGMGPolar> solver() const;

    // Control Parameters
    int verbose() const;
    bool paraview() const;
    int maxOpenMPThreads() const;
    bool DirBC_Interior() const;
    StencilDistributionMethod stencilDistributionMethod() const;
    bool cacheDensityProfileCoefficients() const;
    bool cacheDomainGeometry() const;

    const PolarGrid& grid() const;

    // Full Multigrid Method
    bool FMG() const;
    int FMG_iterations() const;
    MultigridCycleType FMG_cycle() const;

    // Preconditioned Conjugate Gradient Method
    bool PCG() const;
    bool PCG_FMG() const;
    int PCG_FMG_iterations() const;
    MultigridCycleType PCG_FMG_cycle() const;
    int PCG_MG_iterations() const;
    MultigridCycleType PCG_MG_cycle() const;

    // Multigrid Parameters
    ExtrapolationType extrapolation() const;
    int maxLevels() const;
    int preSmoothingSteps() const;
    int postSmoothingSteps() const;
    MultigridCycleType multigridCycle() const;
    int maxIterations() const;
    ResidualNormType residualNormType() const;
    std::optional<double> absoluteTolerance() const;
    std::optional<double> relativeTolerance() const;

private:
    // Parse command-line arguments to extract problem configuration
    cmdline::parser parser_;
    // Input Functions
    std::unique_ptr<const DomainGeometryVariant> domain_geometry_;
    std::unique_ptr<const DensityProfileCoefficientsVariant> density_profile_coefficients_;
    std::unique_ptr<const BoundaryConditionsVariant> boundary_conditions_;
    std::unique_ptr<const SourceTerm> source_term_;
    std::unique_ptr<const ExactSolution> exact_solution_;
    // General solver output and visualization settings
    int verbose_;
    bool paraview_;
    // Parallelization and threading settings
    int max_omp_threads_;
    // Numerical method setup
    bool DirBC_Interior_;
    StencilDistributionMethod stencil_distribution_method_;
    bool cache_density_profile_coefficients_;
    bool cache_domain_geometry_;
    // Grid configuration
    PolarGrid grid_;
    // Multigrid settings
    ExtrapolationType extrapolation_;
    int max_levels_;
    int pre_smoothing_steps_;
    int post_smoothing_steps_;
    MultigridCycleType multigrid_cycle_;
    bool FMG_;
    int FMG_iterations_;
    MultigridCycleType FMG_cycle_;
    bool PCG_;
    bool PCG_FMG_;
    int PCG_FMG_iterations_;
    MultigridCycleType PCG_FMG_cycle_;
    int PCG_MG_iterations_;
    MultigridCycleType PCG_MG_cycle_;
    // Iterative solver controls
    int max_iterations_;
    ResidualNormType residual_norm_type_;
    std::optional<double> absolute_tolerance_;
    std::optional<double> relative_tolerance_;

    void selectTestCase(GeometryType geometry_type, ProblemType problem_type, AlphaCoeff alpha_type,
                        BetaCoeff beta_type, double Rmax, double kappa_eps, double delta_e, double alpha_jump);
};