ArrayBase.h

Go to the documentation of this file.

// Copyright András Vukics 2006–2014. Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE.txt)
// -*- C++ -*-
#ifndef   CPPQEDCORE_QUANTUMDATA_ARRAYBASE_H_INCLUDED
#define   CPPQEDCORE_QUANTUMDATA_ARRAYBASE_H_INCLUDED

#include "ArrayBaseFwd.h"

#include "BlitzArrayExtensions.h"
#include "CMatrix.h"
#include "ComplexArrayExtensions.h"

#include <boost/utility.hpp>


namespace quantumdata {


template<int RANK>
class ArrayBase : private boost::noncopyable
{
protected:
  typedef CArray<RANK> ArrayLow; 
  typedef linalg::CVector CVector;

  explicit ArrayBase(const ArrayLow& arrayLow) : arrayLow_(arrayLow) {} 

  virtual ~ArrayBase() {}


  ArrayBase& operator=(const ArrayLow& arrayLow ) {arrayLow_=arrayLow; return *this;}
  

  const ArrayLow& getArray() const {return arrayLow_;}
        ArrayLow& getArray()       {return const_cast<ArrayLow&>(static_cast<const ArrayBase*>(this)->getArray());}
  

  void operator+=(const ArrayBase& arrayBase) {arrayLow_+=arrayBase.arrayLow_;}
  void operator-=(const ArrayBase& arrayBase) {arrayLow_-=arrayBase.arrayLow_;}
  

  template<typename OTHER>
  void operator*=(const OTHER& dc) {arrayLow_*=dc;}

  template<typename OTHER>
  void operator/=(const OTHER& dc) {arrayLow_/=dc;}



  const CVector vectorView() const {return blitzplusplus::unaryArray(arrayLow_);}
        CVector vectorView()       {return blitzplusplus::unaryArray(arrayLow_);} 
  // The usual technique of defining the non-const in terms of the const doesn't work here, because `CVector is not a pointer, reference, nor a pointer-to-data-member type'


  double frobeniusNorm() const {return sqrt(sum(blitzplusplus::sqrAbs(arrayLow_)));}

private:
  ArrayLow arrayLow_;

};


} // quantumdata



#endif // CPPQEDCORE_QUANTUMDATA_ARRAYBASE_H_INCLUDED