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