Container.h

Go to the documentation of this file.

/* 
 * @BEGIN LICENSE
 *
 * Copyright (C) 2014-2015  Ward Poelmans
 *
 * This file is part of v2DM-DOCI.
 * 
 * v2DM-DOCI is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * Foobar is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with Foobar.  If not, see <http://www.gnu.org/licenses/>.
 *
 * @END LICENSE
 */

#ifndef CONTAINER_H
#define CONTAINER_H

#include <iostream>
#include <cstdlib>
#include <vector>

#include "BlockStructure.h"

namespace doci2DM
{

class Container
{
   friend std::ostream &operator<<(std::ostream &output,const doci2DM::Container &container);

   public:

      //constructor
      Container(int n, int m);

      //copy constructor
      Container(const Container &);

      Container(Container &&);

      //destructor
      virtual ~Container() = default;

      void setMatrixDim(int,int,int);

      void setVectorDim(int,int,int);

      //overload equality operator
      Container &operator=(const Container &);

      Container &operator=(Container &&);

      Container &operator=(double);

      //overload += operator
      Container &operator+=(const Container &);

      //overload -= operator
      Container &operator-=(const Container &);

      Container &daxpy(double alpha,const Container &);

      Container &operator*=(double);

      Container &operator/=(double);

      Container &mprod(const Container &,const Container &);

      //easy to change the numbers
      double &operator()(int block,int i,int j);

      //easy to access the numbers
      double operator()(int block,int i,int j) const;

      double &operator()(int block,int i);

      //easy to access the numbers
      double operator()(int block,int i) const;

      Matrix & getMatrix(int);
      const Matrix & getMatrix(int) const;

      Vector & getVector(int);
      const Vector & getVector(int) const;

      BlockMatrix & getMatrices();
      const BlockMatrix & getMatrices() const;

      BlockVector & getVectors();
      const BlockVector & getVectors() const;

      int gnr() const;

      int gnMatrix() const;

      int gnVector() const;

      int gdimMatrix(int) const;

      int gdimVector(int) const;

      int gdegMatrix(int) const;

      int gdegVector(int) const;

      double trace() const;

      double ddot(const Container &) const;

      void invert();

      void dscal(double);

      void fill_Random();

      void fill_Random(int);

      //positieve of negatieve vierkantswortel uit de matrix
      void sqrt(int);

      void L_map(const Container &,const Container &);

      void symmetrize();

      void sep_pm(Container &,Container &);

   private:

      std::unique_ptr<BlockMatrix> matrix;
      std::unique_ptr<BlockVector> vector;
};

}

#endif

/*  vim: set ts=3 sw=3 expandtab :*/