MED_SliceArray.hxx

Go to the documentation of this file.

// Copyright (C) 2007-2016  CEA/DEN, EDF R&D, 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 correctly compiling "valarray" in modules, which includes this file
  #undef max
  #undef min
#endif

#include <valarray>
#include <stdexcept>

//#if defined(_DEBUG_)
#  define MED_TCSLICE_CHECK_RANGE
//#endif

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