Vector.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 VECTOR_H
#define VECTOR_H

#include <iostream>

#include "Matrix.h"


namespace doci2DM
{

class Vector
{

   friend std::ostream &operator<<(std::ostream &output,const doci2DM::Vector &vector_p);

   public:

      //construct with as input a dimension
      Vector(int);

      //construct with as input a Matrix
      Vector(Matrix &);

      //copy constructor
      Vector(const Vector &);

      Vector(Vector &&);

      //destructor
      virtual ~Vector();

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

      Vector &operator=(Vector &&);

      Vector &operator=(double);

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

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

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

      Vector &operator/=(double);

      Vector &operator*=(double);

      void diagonalize(Matrix &);

      //easy to change the numbers
      double &operator[](int i);

      //easy to access the numbers
      double operator[](int i) const;

      void fill_Random();

      void fill_Random(int seed);

      //get the pointer to the vector
      double *gVector();

      //const version
      const double *gVector() const;

      int gn() const;

      double sum() const;

      double trace() const;

      double log_product() const;

      double ddot(const Vector &) const;

      void dscal(double alpha);

      double min() const;

      double max() const;

      void invert();

      void sqrt(int);

      void symmetrize();

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

      Vector &mprod(const Vector &,const Vector &);

      void sort();

      void sep_pm(Vector &, Vector &);

   private:

      std::unique_ptr<double []> vector;

      int n;
};

}

#endif

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