MED_SliceArray.hxx

Go to the documentation of this file.

// Copyright (C) 2007-2023  CEA, EDF, OPEN CASCADE
//
// Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
//
// This library 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.
//
// This library 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 this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
//
// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
//
 
#ifndef MED_SliceArray_HeaderFile
#define MED_SliceArray_HeaderFile
 
#ifdef WIN32 // for correct compiling of "valarray" in modules, which include this file
  #undef max
  #undef min
#endif
 
#include <valarray>
#include <stdexcept>
 
#  define MED_TCSLICE_CHECK_RANGE
 
namespace MED
{
  //---------------------------------------------------------------
 
  template<class TValueType>
  class TCSlice
  {
    const TValueType* myCValuePtr; 
    size_t mySourceSize; 
    std::slice mySlice; 
 
  protected:
    void
    check_id(size_t theId) const
    {
      long int anId = -1;
      if(theId < mySlice.size()){
        anId = mySlice.start() + theId*mySlice.stride();
        if(anId < (long int)mySourceSize)
          return;
      }
      throw std::out_of_range("TCSlice::check_id");
    }
 
    size_t
    calculate_id(size_t theId) const
    {
      return mySlice.start() + theId*mySlice.stride();
    }
 
    size_t
    get_id(size_t theId) const
    {
#ifdef MED_TCSLICE_CHECK_RANGE
      check_id(theId);
#endif
      return calculate_id(theId);
    }
 
    size_t
    get_id_at(size_t theId) const
    {
      check_id(theId);
      return calculate_id(theId);
    }
 
  public:
    typedef TValueType value_type;
 
    TCSlice(const value_type* theValuePtr,
            size_t theSourceSize,
            const std::slice& theSlice):
      myCValuePtr(theValuePtr),
      mySourceSize(theSourceSize),
      mySlice(theSlice)
    {}
 
    TCSlice(const TVector<value_type>& theContainer,
            const std::slice& theSlice):
      myCValuePtr(&theContainer[0]),
      mySourceSize(theContainer.size()),
      mySlice(theSlice)
    {}
 
    TCSlice():
      myCValuePtr(NULL)
    {}
 
    const value_type&
    operator[](size_t theId) const
    {
      return *(myCValuePtr + get_id(theId));
    }
 
    const value_type&
    at(size_t theId) const
    {
      return *(myCValuePtr + get_id_at(theId));
    }
 
    size_t
    size() const
    {
      return mySlice.size();
    }
  };
 
  //---------------------------------------------------------------
  template<class TValueType>
  class TSlice: public TCSlice<TValueType>
  {
    TValueType* myValuePtr;
 
  public:
    typedef TValueType value_type;
    typedef TCSlice<TValueType> TSupperClass;
 
    TSlice(value_type* theValuePtr,
           size_t theSourceSize,
           const std::slice& theSlice):
      TSupperClass(theValuePtr, theSourceSize, theSlice),
      myValuePtr(theValuePtr)
    {}
 
    TSlice(TVector<value_type>& theContainer,
           const std::slice& theSlice):
      TSupperClass(theContainer, theSlice),
      myValuePtr(&theContainer[0])
    {}
 
    TSlice():
      myValuePtr(NULL)
    {}
 
    value_type&
    operator[](size_t theId)
    {
      return *(myValuePtr + this->get_id(theId));
    }
 
    value_type&
    at(size_t theId)
    {
      return *(myValuePtr + this->get_id_at(theId));
    }
  };
}
 
#undef MED_TCSLICE_CHECK_RANGE
 
#endif // MED_SliceArray_HeaderFile