Matrix_Pointers.hxx

// Copyright (C) 2001-2009 Vivien Mallet
//
// This file is part of the linear-algebra library Seldon,
// http://seldon.sourceforge.net/.
//
// Seldon is free software; you can redistribute it and/or modify it under the
// terms of the GNU Lesser General Public License as published by the Free
// Software Foundation; either version 2.1 of the License, or (at your option)
// any later version.
//
// Seldon 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 Lesser General Public License for
// more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with Seldon. If not, see http://www.gnu.org/licenses/.
 
 
// To be included by Seldon.hxx
 
#ifndef SELDON_FILE_MATRIX_POINTERS_HXX
 
#include "../share/Common.hxx"
#include "../share/Properties.hxx"
#include "../share/Storage.hxx"
#include "../share/Errors.hxx"
#include "../share/Allocator.hxx"
 
namespace Seldon
{
 
 
  template <class T, class Prop, class Storage, class Allocator
            = typename SeldonDefaultAllocator<Storage, T>::allocator>
  class Matrix_Pointers: public Matrix_Base<T, Allocator>
  {
    // typedef declaration.
  public:
    typedef typename Allocator::value_type value_type;
    typedef typename Allocator::pointer pointer;
    typedef typename Allocator::const_pointer const_pointer;
    typedef typename Allocator::reference reference;
    typedef typename Allocator::const_reference const_reference;
    typedef typename Allocator::value_type entry_type;
    typedef typename Allocator::reference access_type;
    typedef typename Allocator::const_reference const_access_type;
 
    // Attributes.
  protected:
    pointer* me_;
 
    // Methods.
  public:
    // Constructor.
    Matrix_Pointers();
    explicit Matrix_Pointers(int i, int j);
    Matrix_Pointers(const Matrix_Pointers<T, Prop, Storage, Allocator>& A);
 
    // Destructor.
    ~Matrix_Pointers();
    void Clear();
 
    // Basic methods.
    long GetDataSize() const;
    size_t GetMemorySize() const;
    pointer* GetMe() const;
 
    // Memory management.
    void Reallocate(int i, int j);
    void SetData(int i, int j, pointer data);
    void Nullify();
    void Resize(int i, int j);
 
    // Element access and affectation.
    pointer GetDataPointer(int i, int j) const;
    reference operator() (int i, int j);
#ifndef SWIG
    const_reference operator() (int i, int j) const;
#endif
    reference Val(int i, int j);
    reference Get(int i, int j);
#ifndef SWIG
    const_reference Val(int i, int j) const;
    const_reference Get(int i, int j) const;
    reference operator[] (int i);
    const_reference operator[] (int i) const;
    
    Matrix_Pointers<T, Prop, Storage, Allocator>&
    operator= (const Matrix_Pointers<T, Prop, Storage, Allocator>& A);
#endif
 
    void Set(int i, int j, const T& val);
    void Copy(const Matrix_Pointers<T, Prop, Storage, Allocator>& A);
 
    // Convenient functions.
    int GetLD() const;
    void Zero();
    void SetIdentity();
    void Fill();
    template <class T0>
    void Fill(const T0& x);
#ifndef SWIG
    template <class T0>
    Matrix_Pointers<T, Prop, Storage, Allocator>& operator= (const T0& x);
#endif
    void FillRand();
    void Print() const;
    void Print(int a, int b, int m, int n) const;
    void Print(int l) const;
 
    // Input/output functions.
    void Append(string FileName) const;
    void Write(string FileName, bool with_size = true) const;
#ifdef SELDON_WITH_HDF5
    void WriteHDF5(string FileName, string group_name, string dataset_name)
      const;
#endif
    void Write(ostream& FileStream, bool with_size = true) const;
    void WriteText(string FileName) const;
    void WriteText(ostream& FileStream) const;
    void Read(string FileName, bool with_size = true);
    void Read(istream& FileStream, bool with_size = true);
    void ReadText(string FileName);
    void ReadText(istream& FileStream);
 
#ifdef SELDON_WITH_VIRTUAL
    // ApplySor, MltVector and MltAddVector are not defined for this matrix
    // since this class is used with differents types T (not only real
    // or complex numbers)
    virtual bool IsSymmetric() const;
 
    virtual void AddInteraction(int, int, const T&);
    virtual void AddInteractionRow(int, int, const Vector<int>&,
                                   const Vector<T>& val, bool sorted = false);
 
    virtual void AddInteractionColumn(int, int, const Vector<int>&,
                                      const Vector<T>& val, bool sorted = false);
    
    virtual void ClearRow(int i);
 
    typedef typename ClassComplexType<T>::Treal Treal;
    typedef typename ClassComplexType<T>::Tcplx Tcplx;
        
    virtual void MltAddVector(const Treal& alpha, const Vector<Treal>& x,
                              const Treal& beta, Vector<Treal>& y) const;
 
    virtual void MltAddVector(const Tcplx& alpha, const Vector<Tcplx>& x,
                              const Tcplx& beta, Vector<Tcplx>& y) const;
 
    virtual void MltAddVector(const Treal& alpha, const SeldonTranspose&,
                              const Vector<Treal>& x,
                              const Treal& beta, Vector<Treal>& y) const;
 
    virtual void MltAddVector(const Tcplx& alpha, const SeldonTranspose&,
                              const Vector<Tcplx>& x,
                              const Tcplx& beta, Vector<Tcplx>& y) const;
    
    virtual void MltVector(const Vector<Treal>& x, Vector<Treal>& y) const;
    virtual void MltVector(const Vector<Tcplx>& x, Vector<Tcplx>& y) const;
    
    virtual void MltVector(const SeldonTranspose&,
                           const Vector<Treal>& x, Vector<Treal>& y) const;
 
    virtual void MltVector(const SeldonTranspose&,
                           const Vector<Tcplx>& x, Vector<Tcplx>& y) const;
 
#endif
 
  };
 
 
  template <class T, class Prop, class Allocator>
  class Matrix<T, Prop, ColMajor, Allocator>:
    public Matrix_Pointers<T, Prop, ColMajor, Allocator>
  {
    // typedef declaration.
  public:
    typedef typename Allocator::value_type value_type;
    typedef Prop property;
    typedef ColMajor storage;
    typedef Allocator allocator;
 
  public:
    Matrix();
    explicit Matrix(int i, int j);
    Matrix(const Matrix<T, Prop, ColMajor, Allocator>& A);
 
    void WriteColumn(string FileName, int col) const;
    void WriteColumn(ostream& FileStream, int col) const;
 
#ifndef SWIG
    template <class T0>
    Matrix<T, Prop, ColMajor, Allocator>& operator= (const T0& x);
    Matrix<T, Prop, ColMajor, Allocator>& operator=(const Matrix<T, Prop,
                                                    ColMajor, Allocator>& A);
#endif
    template<class T0>
    Matrix<T, Prop, ColMajor, Allocator>& operator*= (const T0& x);
  };
 
 
  template <class T, class Prop, class Allocator>
  class Matrix<T, Prop, RowMajor, Allocator>:
    public Matrix_Pointers<T, Prop, RowMajor, Allocator>
  {
    // typedef declaration.
  public:
    typedef typename Allocator::value_type value_type;
    typedef Prop property;
    typedef RowMajor storage;
    typedef Allocator allocator;
 
  public:
    Matrix();
    explicit Matrix(int i, int j);
    Matrix(const Matrix<T, Prop, RowMajor, Allocator>& A);
 
    void WriteRow(string FileName, int row) const;
    void WriteRow(ostream& FileStream, int row) const;
 
#ifndef SWIG
 
    template <class T0>
    Matrix<T, Prop, RowMajor, Allocator>& operator= (const T0& x);
    Matrix<T, Prop, RowMajor, Allocator>& operator=(const Matrix<T, Prop,
                                                    RowMajor, Allocator>& A);
#endif
    template<class T0>
    Matrix<T, Prop, RowMajor, Allocator>& operator*= (const T0& x);
  };
 
 
} // namespace Seldon.
 
#define SELDON_FILE_MATRIX_POINTERS_HXX
#endif