primal_ts_base.h

Go to the documentation of this file.

#ifndef PRIMAL_TS_BASE_H_
#define PRIMAL_TS_BASE_H_

#include "ts_base.h"
namespace DOpE
{

  template<typename OPTPROBLEM, typename SPARSITYPATTERN, typename VECTOR,
      int dopedim, int dealdim,
      template <int, int> class FE = dealii::FESystem,
      template <int, int> class DH = dealii::DoFHandler>
    class PrimalTSBase : public TSBase<OPTPROBLEM, SPARSITYPATTERN, VECTOR,
        dopedim, dealdim, FE, DH>
    {
      public:
        PrimalTSBase(OPTPROBLEM& OP) :
            TSBase<OPTPROBLEM, SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE,
                DH>(OP)
        {
        }

        ~PrimalTSBase()
        {
        }

        /******************************************************/
        /****For the initial values ***************/
        template<typename EDC>
          void
          Init_ElementEquation(const EDC& edc,
              dealii::Vector<double> &local_vector, double scale,
              double scale_ico)
          {
            this->GetProblem().Init_ElementEquation(edc, local_vector, scale,
                scale_ico);
          }

        template<typename EDC>
          void
          Init_ElementRhs(const EDC& edc,
              dealii::Vector<double> &local_vector, double scale)
          {
            this->GetProblem().Init_ElementRhs(edc, local_vector, scale);
          }

        template<typename EDC>
          void
          Init_ElementMatrix(const EDC& edc,
              dealii::FullMatrix<double> &local_matrix, double scale, double scale_ico)
          {
            this->GetProblem().Init_ElementMatrix(edc, local_matrix, scale,
                scale_ico);
          }

        void
        Init_PointRhs(
            const std::map<std::string, const dealii::Vector<double>*> &/*param_values*/,
            const std::map<std::string, const VECTOR*> &/*domain_values*/,
            VECTOR& /*rhs_vector*/, double /*scale*/)
        {
        }

       template<typename FDC>
          void
          Init_FaceEquation(const FDC& /*fdc*/,
              dealii::Vector<double> &/*local_vector*/, double /*scale*/, double /*scale_ico*/)
          {
          }

         template<typename FDC>
          void
          Init_InterfaceEquation(const FDC& /*fdc*/,
              dealii::Vector<double> &/*local_vector*/, double /*scale*/, double /*scale_ico*/)
          {
          }

        template<typename FDC>
          void
          Init_BoundaryEquation(const FDC& /*fdc*/,
              dealii::Vector<double> &/*local_vector*/, double /*scale*/, double /*scale_ico*/)
          {
          }

        template<typename FDC>
          void
          Init_FaceMatrix(const FDC& /*fdc*/,
              FullMatrix<double> &/*local_matrix*/, double /*scale*/, double /*scale_ico*/)
          {
          }

        template<typename FDC>
          void
          Init_InterfaceMatrix(const FDC& /*fdc*/,
              FullMatrix<double> &/*local_matrix*/, double /*scale*/, double /*scale_ico*/)
          {
          }

        template<typename FDC>
          void
          Init_BoundaryMatrix(const FDC& /*fdc*/,
              FullMatrix<double> &/*local_matrix*/, double /*scale*/, double /*scale_ico*/)
          {
          }

        /****End the initial values ***************/
    };
}

#endif /* TS_BASE_PRIMAL_H_ */