primaldirichletdata.h

Go to the documentation of this file.

#ifndef _PRIMAL_DIRICHLET_DATA_H_
#define _PRIMAL_DIRICHLET_DATA_H_

#include "function_wrapper.h"
#include "dofhandler_wrapper.h"
#include "fevalues_wrapper.h"

namespace DOpE
{

  template<typename DD, typename VECTOR, int dealdim>
    class PrimalDirichletData : public DOpEWrapper::Function<dealdim>
  {
  public:
  PrimalDirichletData(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(
                             _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