Array4D_Inline.cxx

// Copyright (C) 2010 Lin Wu
// Copyright (C) 2001-2009 Vivien Mallet
// Copyright (C) 2003-2009 Marc Duruflé
//
// 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/.
 
 
#ifndef SELDON_FILE_ARRAY4D_INLINE_CXX
 
#include "Array4D.hxx"
 
namespace Seldon
{
 
  /**************
   * DESTRUCTOR *
   **************/
 
 
  template <class T, class Allocator>
  inline Array4D<T, Allocator>::~Array4D()
  {
    Clear();
  }
 
 
  /*****************
   * BASIC METHODS *
   *****************/
 
 
 
  template <class T, class Allocator>
  inline int Array4D<T, Allocator>::GetLength1() const
  {
    return length1_;
  }
 
 
 
  template <class T, class Allocator>
  inline int Array4D<T, Allocator>::GetLength2() const
  {
    return length2_;
  }
 
 
 
  template <class T, class Allocator>
  inline int Array4D<T, Allocator>::GetLength3() const
  {
    return length3_;
  }
 
 
 
  template <class T, class Allocator>
  inline int Array4D<T, Allocator>::GetLength4() const
  {
    return length4_;
  }
 
 
 
  template <class T, class Allocator>
  inline long Array4D<T, Allocator>::GetSize() const
  {
    return length1_ * length234_;
  }
 
 
 
  template <class T, class Allocator>
  inline long Array4D<T, Allocator>::GetDataSize() const
  {
    return length1_ * length234_;
  }
 
 
 
  template <class T, class Allocator>
  inline typename Array4D<T, Allocator>::pointer Array4D<T, Allocator>
  ::GetData() const
  {
    return data_;
  }
 
 
 
  template <class T, class Allocator>
  inline typename Array4D<T, Allocator>::pointer Array4D<T, Allocator>
  ::GetDataPointer(int i, int j, int k, int l) const
  {
    return data_ + i * length234_ + j * length34_ + k * length4_ + l;
  }
 
 
  /**********************************
   * ELEMENT ACCESS AND AFFECTATION *
   **********************************/
 
 
 
  template <class T, class Allocator>
  inline typename Array4D<T, Allocator>::reference
  Array4D<T, Allocator>::operator() (int i, int j, int k, int l)
  {
 
#ifdef SELDON_CHECK_BOUNDS
    CheckBounds(i, j, k, l, length1_,
                length2_, length3_, length4_, "Array4D");
#endif
 
    return data_[i * length234_ + j * length34_ + k * length4_ + l];
  }
 
 
 
  template <class T, class Allocator>
  inline typename Array4D<T, Allocator>::const_reference
  Array4D<T, Allocator>::operator() (int i, int j, int k, int l) const
  {
 
#ifdef SELDON_CHECK_BOUNDS
    CheckBounds(i, j, k, l, length1_,
                length2_, length3_, length4_, "Array4D");
#endif
 
    return data_[i*length234_ + j*length34_ + k*length4_ + l];
  }
 
 
  template <class T, class Allocator>
  inline Array4D<T, Allocator>& Array4D<T, Allocator>::operator=
  (const Array4D<T, Allocator>& A)
  {
    Copy(A);
 
    return *this;
  }
 
 
  template <class T, class Allocator>
  inline void Array4D<T, Allocator>::Copy(const Array4D<T, Allocator>& A)
  {
    Reallocate(A.GetLength1(), A.GetLength2(),
               A.GetLength3(), A.GetLength4());
 
    Allocator::memorycpy(data_, A.GetData(), GetDataSize());
  }
  
} // namespace Seldon.
 
#define SELDON_FILE_ARRAY4D_INLINE_CXX
#endif