ManufacturedSolution.h

Go to the documentation of this file.

#ifndef MMS_H
#define MMS_H
 
#include <array>
#include <functional>
#include <memory>
#include <string>
#include <valarray>
 
 
namespace AMP {
 
class Database;
 
 
class ManufacturedSolution final
{
public:
    explicit ManufacturedSolution( std::shared_ptr<Database> db );
 
    std::array<double, 10> evaluate( double x, double y, double z ) const;
 
    size_t getNumberOfParameters() const { return d_NumberOfParameters; }
 
    size_t getNumberOfInputs() const { return d_NumberOfInputs; }
 
    void setTricubicParams( const std::valarray<double> &cin, const std::valarray<double> &ain )
    {
        d_c = cin;
        d_a = ain;
    }
 
    const std::string &get_name() const { return d_Name; }
 
public:
    const auto &getc() const { return d_c; }
    const auto &geta() const { return d_a; }
    const auto &geth() const { return d_h; }
    const auto &geths() const { return d_hs; }
 
private:
    enum class FunctionType { POLYNOMIAL, GENERALQUADRATIC };
    enum class Geometry { BRICK, CYLROD, CYLRODRZ, CYLSHELL, QTRCYLSHELL, LASTGeometry };
    enum class Order { QUADRATIC, CUBIC, FOURIER, GAUSSIAN, LASTOrder };
    enum class BCType { NEUMANN, DIRICHLET1, DIRICHLET2, DIRICHLETZ2, NONE, LASTType };
 
    FunctionType d_FunctionType;
 
    Geometry d_geom;
    Order d_order;
    BCType d_bcType;
    size_t d_NumberOfParameters;
    size_t d_NumberOfInputs;
 
    std::function<std::array<double, 10>( double, double, double, const ManufacturedSolution * )>
        d_functionPointer;
 
    std::valarray<double> d_c;
    std::valarray<double> d_a;
 
    double d_MinX, d_MaxX, d_ScaleX;
    double d_MinY, d_MaxY, d_ScaleY;
    double d_MinZ, d_MaxZ, d_ScaleZ;
    double d_MinR, d_MaxR, d_ScaleR;
    double d_MinTh, d_MaxTh, d_ScaleTh;
 
    std::array<std::array<double, 3>, 3> d_h;
    std::array<std::array<double, 3>, 3> d_hs; // symmetrized h
 
    double d_MaximumTheta;
    bool d_CylindricalCoords;
 
    std::string d_Name;
};
} // namespace AMP
#endif