transposedgradientdirichletdata.h

Go to the documentation of this file.

#ifndef TRANSPOSED_GRADIENT_DIRICHLET_DATA_H_
#define TRANSPOSED_GRADIENT_DIRICHLET_DATA_H_

#include "transposeddirichletdatainterface.h"

namespace DOpE
{

  template<typename DD, typename VECTOR, int dealdim>
    class TransposedGradientDirichletData : public TransposedDirichletDataInterface<dealdim>
  {
  public:
  TransposedGradientDirichletData(const DD& data) : TransposedDirichletDataInterface<dealdim>(), 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;
    }


    void value (const dealii::Point<dealdim>   &p,
                const unsigned int  component,
                const unsigned int  dof_number,
                dealii::Vector<double>& local_vector) const
    {
      dirichlet_data_.Data_QT(
                              param_values_,
                              domain_values_,
                              color_,
                              p,
                              component,
                              dof_number,
                              local_vector);
    }

    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