Permeability.h

Go to the documentation of this file.

 
#pragma once
 
#include <Eigen/Core>
 
#include "BaseLib/Error.h"
#include "ParameterLib/Parameter.h"
 
namespace MaterialLib
{
namespace PorousMedium
{
class Permeability
{
public:
    explicit Permeability(
        ParameterLib::Parameter<double> const& permeability_parameter,
        int const dimension)
        : _permeability_parameter(permeability_parameter), _dimension(dimension)
    {
        if (permeability_parameter.getNumberOfGlobalComponents() !=
            _dimension * _dimension)
        {
            OGS_FATAL(
                "The given parameter has {:d} components, but the permeability "
                "tensor is defined for a {:d} dimensional problem.",
                permeability_parameter.getNumberOfGlobalComponents(),
                _dimension);
        }
    }
 
    virtual ~Permeability() = default;
 
    virtual Eigen::MatrixXd getValue(const double t,
                                     ParameterLib::SpatialPosition const& pos,
                                     const double variable,
                                     const double temperature) const
    {
        (void)variable;
        (void)temperature;
 
        return Eigen::Map<Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic,
                                        Eigen::RowMajor> const>(
            _permeability_parameter(t, pos).data(), _dimension, _dimension);
    }
 
private:
    ParameterLib::Parameter<double> const& _permeability_parameter;
    int const _dimension;
};
 
}  // namespace PorousMedium
}  // namespace MaterialLib