stokes-slip/parameters.hh at master · jstebel/stokes-slip · 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
#ifndef PARAMETERS_HH
#define PARAMETERS_HH

#include <deal.II/base/parameter_handler.h>
#include <deal.II/base/point.h>


using namespace dealii;


namespace Parameters {

//     struct Materials {
//         double Young_modulus_matrix;
//         double Poisson_ratio_matrix;
//         double Young_modulus_fiber;
//         double Poisson_ratio_fiber;
//         double Fiber_volume_ratio;
//         unsigned int Reinforcement_material_id;
//         bool use_1d_fibers;
//
//         static void declare_parameters(ParameterHandler &prm);
//         void parse_parameters(ParameterHandler &prm);
//     };

    struct BoundaryConditions {
        std::map<int, std::string > velocity;
        std::map<int, std::string > velocity_tangent;
        std::map<int, std::string > velocity_normal;
        std::map<int, std::string > velocity_slip;
        std::map<int, std::string > traction;

        static void declare_parameters(ParameterHandler &prm);
        void parse_parameters(ParameterHandler &prm);
    };

    struct BulkParameters {
        std::map<int, std::string> force;
        bool sym_grad;

        static void declare_parameters(ParameterHandler &prm);
        void parse_parameters(ParameterHandler &prm);
    };

    struct Numerics {
      double epsilon;
      double a_tol;
      double r_tol;
      unsigned int max_iter;

      static void declare_parameters(ParameterHandler &prm);
      void parse_parameters(ParameterHandler &prm);
    };

    struct IO {
        unsigned int dim;
        std::string mesh_file;
        std::string output_file_base;

        static void declare_parameters(ParameterHandler &prm);
        void parse_parameters(ParameterHandler &prm);
    };

    struct CostFunction {
      std::map<int, std::string> boundary_integral;
      std::map<int, std::string> volume_integral;

      static void declare_parameters(ParameterHandler &prm);
      void parse_parameters(ParameterHandler &prm);
    };

    struct ShapeOptimization {
      Point<3> p1;
      Point<3> p2;
      Point<3> ptop;
      double height;
      unsigned int np;
      unsigned int maxit;
      std::vector<double> init_guess;
      std::string f_max;
      std::string f_min;
      double g_max;
      double h_max;
      std::string logfile;

      static void declare_parameters(ParameterHandler &prm);
      void parse_parameters(ParameterHandler &prm);
    };


    struct AllParameters :
//             public Materials,
            public BoundaryConditions,
            public BulkParameters,
            public CostFunction,
            public Numerics,
            public IO,
            ShapeOptimization
    {
        AllParameters(const std::string &input_file);
        static void declare_parameters(ParameterHandler &prm);
        void parse_parameters(ParameterHandler &prm);
        void print(std::ostream &out);

        ParameterHandler prm;
    };

}


#endif