instatoptproblemcontainer.h

Go to the documentation of this file.

#ifndef INSTAT_OPT_PROBLEM_CONTAINER_
#define INSTAT_OPT_PROBLEM_CONTAINER_

#include "optproblemcontainer.h"

namespace DOpE
{
  template<template<typename BASE_PROB, typename SPARSITYPATTERN, typename VECTOR,int dopedim, int dealdim, template<int, int> class FE, template<int, int> class DH> class PRIMALTSPROBLEM,
    template<typename BASE_PROB, typename SPARSITYPATTERN, typename VECTOR,int dopedim, int dealdim, template<int, int> class FE, template<int, int> class DH> class  ADJOINTTSPROBLEM,
    typename FUNCTIONAL_INTERFACE, typename FUNCTIONAL, typename PDE,
    typename DD, typename CONSTRAINTS, typename SPARSITYPATTERN,
    typename VECTOR, int dopedim, int dealdim,
    template<int, int> class FE = FESystem,
    template<int, int> class DH = dealii::DoFHandler>
  class InstatOptProblemContainer : public OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,
  SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>
{
  public:
  InstatOptProblemContainer(FUNCTIONAL& functional,PDE& pde,CONSTRAINTS& constraints,
                            SpaceTimeHandler<FE, DH, SPARSITYPATTERN, VECTOR, dopedim,dealdim>& STH)
  : OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>(
    functional,pde,constraints,STH), ts_state_problem_(NULL), ts_adjoint_problem_(NULL), 
  ts_tangent_problem_(NULL), ts_adjoint_hessian_problem_(NULL)
  {
  }
  
  ~InstatOptProblemContainer()
  {
    if(ts_state_problem_ != NULL)
    {
      delete ts_state_problem_;
    }
    if(ts_adjoint_problem_ != NULL)
    {
      delete ts_adjoint_problem_;
    }
    if(ts_tangent_problem_ != NULL)
    {
      delete ts_tangent_problem_;
    }
    if(ts_adjoint_hessian_problem_ != NULL)
    {
      delete ts_adjoint_hessian_problem_;
    }
  }

  void ReInit(std::string algo_type)
  {  
    if(ts_state_problem_ != NULL)
    {
      delete ts_state_problem_;
      ts_state_problem_ = NULL;
    }
    if(ts_adjoint_problem_ != NULL)
    {
      delete ts_adjoint_problem_;
      ts_adjoint_problem_ = NULL;
    }
    if(ts_tangent_problem_ != NULL)
    {
      delete ts_tangent_problem_;
      ts_tangent_problem_ = NULL;
    }
    if(ts_adjoint_hessian_problem_ != NULL)
    {
      delete ts_adjoint_hessian_problem_;
      ts_adjoint_hessian_problem_ = NULL;
    }
    OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>::ReInit(algo_type);
        
  }

  std::string GetName() const
  {
    return "InstatOptProblemContainer";
  }
  
  /*******************************************************************************************/
  PRIMALTSPROBLEM<StateProblem<
  OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
  PDE, DD, SPARSITYPATTERN, VECTOR, dealdim>,
  SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>& GetStateProblem()
  {
    if(ts_state_problem_ == NULL)
    {
      ts_state_problem_ = new PRIMALTSPROBLEM<StateProblem<
      OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
      PDE, DD, SPARSITYPATTERN, VECTOR, dealdim>,
      SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>(OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,
                                              SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>::GetStateProblem());
    }
    return *ts_state_problem_;
  }
  /*******************************************************************************************/

  //FIXME: This should use the GetAdjointProblem of OptProblemContainer once availiable simillar for all calls
  // to ADJOINTTSPROBLEM ...
  ADJOINTTSPROBLEM<OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
  SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>& GetAdjointProblem()
  {
    if(ts_adjoint_problem_ == NULL)
    {
      ts_adjoint_problem_ = new ADJOINTTSPROBLEM<OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
      SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>(OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,
                                                                 SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>::GetBaseProblem());
    }
    return *ts_adjoint_problem_;
  }

  /*******************************************************************************************/
  //FIXME: This should use the GetTangentProblem of OptProblemContainer once availiable simillar for all 
  // other related calls
  PRIMALTSPROBLEM<OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
  SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>& GetTangentProblem()
  {
    if(ts_tangent_problem_ == NULL)
    {
      ts_tangent_problem_ = new PRIMALTSPROBLEM<OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
      SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>(OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,
                                              SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>::GetBaseProblem());
    }
    return *ts_tangent_problem_;
  }
  /*******************************************************************************************/

  //FIXME: This should use the GetAdjointHessianProblem of OptProblemContainer once availiable simillar for all 
  // other related calls
  ADJOINTTSPROBLEM<OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
  SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>& GetAdjointHessianProblem()
  {
    if(ts_adjoint_hessian_problem_ == NULL)
    {
      ts_adjoint_hessian_problem_ = new ADJOINTTSPROBLEM<OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
      SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>(OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,
                                                                 SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>::GetBaseProblem());
    }
    return *ts_adjoint_hessian_problem_;
  }

private:
  PRIMALTSPROBLEM<StateProblem<
      OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
      PDE, DD, SPARSITYPATTERN, VECTOR, dealdim>,
      SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>* ts_state_problem_;
  ADJOINTTSPROBLEM<OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
      SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>* ts_adjoint_problem_;
  PRIMALTSPROBLEM<OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
      SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>* ts_tangent_problem_;
  ADJOINTTSPROBLEM<OptProblemContainer<FUNCTIONAL_INTERFACE,FUNCTIONAL,PDE,DD,CONSTRAINTS,SPARSITYPATTERN,VECTOR,dopedim,dealdim,FE, DH>,
      SPARSITYPATTERN, VECTOR, dopedim, dealdim, FE, DH>* ts_adjoint_hessian_problem_;
};
}
#endif