BlockStructure.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 BLOCKSTRUCTURE_H
#define BLOCKSTRUCTURE_H

#include <iostream>
#include <memory>
#include <vector>

#include "Matrix.h"
#include "Vector.h"

namespace doci2DM
{

template<class BlockType>
class BlockStructure
{
   template<class MyBlockType>
   friend std::ostream &operator<< (std::ostream &output,const doci2DM::BlockStructure<MyBlockType> &blocks_p);

   public:

      BlockStructure(int);

      //construct with as input a BlockMatrixType
      //   BlockStructure(BlockType &);

      //copy constructor
      BlockStructure(const BlockStructure<BlockType> &);

      BlockStructure(BlockStructure<BlockType> &&);

      //destructor
      virtual ~BlockStructure() = default;

      void setDim(int,int,int);

      //overload equality operator
      BlockStructure &operator=(const BlockStructure<BlockType> &);

      BlockStructure &operator=(BlockStructure<BlockType> &&);

      BlockStructure &operator=(double);

      //overload += operator
      BlockStructure &operator+=(const BlockStructure<BlockType> &);

      //overload -= operator
      BlockStructure &operator-=(const BlockStructure<BlockType> &);

      BlockStructure &daxpy(double alpha,const BlockStructure<BlockType> &);

      BlockStructure &operator*=(double);

      BlockStructure &operator/=(double);

      BlockType &operator[](int block);

      //const version
      const BlockType &operator[](int block) const;

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

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

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

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

//      virtual void diagonalize(BlockStructure<Matrix> &);

      int gnr() const;

      int gdim(int) const;

      int gdeg(int) const;

      double trace() const;

      double ddot(const BlockStructure<BlockType> &) const;

      void dscal(double alpha);

      double min() const;

      double max() const;

      void fill_Random();

      void fill_Random(int);

      virtual void invert();

      virtual void sqrt(int);

      virtual void L_map(const BlockStructure<BlockType> &,const BlockStructure<BlockType> &);

      virtual BlockStructure<BlockType> &mprod(const BlockStructure<BlockType> &, const BlockStructure<BlockType> &);

      void symmetrize();

      void sort();

      virtual void sep_pm(BlockStructure<BlockType> &, BlockStructure<BlockType> &);

   private:

      std::vector< std::unique_ptr<BlockType> > blocks;

      std::vector<int> degen;
};

typedef BlockStructure<Matrix> BlockMatrix;
typedef BlockStructure<Vector> BlockVector;

}

#endif

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