tangentdirichletdata.h

Go to the documentation of this file.

#ifndef _TANGENT_DIRICHLET_DATA_H_
#define _TANGENT_DIRICHLET_DATA_H_

#include "function_wrapper.h"

namespace DOpE
{

  template<typename DD, typename VECTOR, int dealdim>
    class TangentDirichletData : public DOpEWrapper::Function<dealdim>
  {
  public:
  TangentDirichletData(const DD& data) : DOpEWrapper::Function<dealdim>(data.n_components(), data.InitialTime()), _dirichlet_data(data)
    {
      _param_values = NULL;
      _domain_values = NULL;
      _color = 0;
    }

    void ReInit(
                const std::map<std::string, const dealii::Vector<double>* > &param_values,
                const std::map<std::string, const VECTOR* > &domain_values,
                unsigned int color)
    {
      _param_values = &param_values;
      _domain_values = &domain_values;
      _color = color;
    }


    double value (const dealii::Point<dealdim>   &p,
                  const unsigned int  component) const
    {
      return _dirichlet_data.Data_Q(
        _param_values,
        _domain_values,
        _color,
        p,
        component);
    }

    void SetTime(double time) const
    {
      _dirichlet_data.SetTime(time);
    }
  private:
    const DD& _dirichlet_data;
    const std::map<std::string, const dealii::Vector<double>* >* _param_values;
    const std::map<std::string, const VECTOR* >* _domain_values;
    unsigned int _color;
  };
}
#endif