transposedhessiandirichletdata.h

Go to the documentation of this file.

#ifndef _TRANSPOSED_HESSIAN_DIRICHLET_DATA_H_
#define _TRANSPOSED_HESSIAN_DIRICHLET_DATA_H_

#include "transposeddirichletdatainterface.h"

namespace DOpE
{

  template<typename DD, typename VECTOR, int dealdim>
    class TransposedHessianDirichletData : public TransposedDirichletDataInterface<dealdim>
  {
  public:
  TransposedHessianDirichletData(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);
      _dirichlet_data.Data_QQT(
                               _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