controlvector.h

Go to the documentation of this file.

#ifndef _CONTROL_VECTOR_H_
#define _CONTROL_VECTOR_H_

#include "spacetimehandler_base.h"

#include <lac/vector.h>
#include <lac/block_vector_base.h>
#include <lac/block_vector.h>

#include <vector>
#include <iostream>
#include <sstream>


namespace DOpE
{

template<typename VECTOR>
  class ControlVector
  {
  public:
    //TODO: Currently we only consider one fixed control 
    //      for all timesteps, if more is desired one needs to augment the 
    //      Spacetimehandler to have a time discretization for the control,
    //      Then one can update this vector similar to the statevector 
    //      with different meshes for Vectors.
    //      Note that this requires to keep track of the interpolation
    //      between state and control time points...
    ControlVector(const ControlVector& ref);
    ControlVector(const SpaceTimeHandlerBase<VECTOR>* STH, std::string behavior);
    ~ControlVector();

    void SetTime(double t, const TimeIterator& interval) const;
    void SetTimeDoFNumber(unsigned int time_point) const;
    VECTOR& GetSpacialVector();
    const VECTOR& GetSpacialVector() const;
    const dealii::Vector<double>& GetSpacialVectorCopy() const;

    void operator=(double value);
    void operator=(const ControlVector& dq);
    void operator+=(const ControlVector& dq);
    void operator*=(double a);
    double operator*(const ControlVector& dq) const;
    void add(double s, const ControlVector& dq);
    void equ(double s, const ControlVector& dq);
    void max(const ControlVector& dq);
    void min(const ControlVector& dq);

    void comp_mult(const ControlVector& dq);
    
    void comp_invert();
    
    void init_by_sign(double smaller, double larger, double unclear, double TOL = 1.e-10);
    
    void PrintInfos(std::stringstream& out);
    void UnLockCopy() const { _lock = false; }

    std::string GetBehavior() const { return _behavior; }

    const SpaceTimeHandlerBase<VECTOR>* GetSpaceTimeHandler() const { return _STH; }

    void ReInit();
 
    double Norm(std::string name,std::string restriction = "all") const;

  private:
    void ReSizeSpace(unsigned int ndofs, const std::vector<unsigned int>& dofs_per_block);

    std::vector<VECTOR* > _control;
    mutable VECTOR _local_control;
    mutable dealii::Vector<double> _copy_control;

    mutable int _accessor;
    mutable bool _lock;

    const SpaceTimeHandlerBase<VECTOR>* _STH;
    std::string _behavior;
    unsigned int _sfh_ticket;
  };


}
#endif