fortranshape.h

Go to the documentation of this file.

// include guard
#ifndef AFF_FORTRANSHAPE_H_VERSION

#define AFF_FORTRANSHAPE_H_VERSION \
  "AFF_FORTRANSHAPE_H   V1.2"

#include<aff/array.h>
#include<aff/lib/checkedcast.h>

namespace aff {

namespace util {

  class FortranShape {
    public:
      typedef aff::Strided Tshape;
      typedef Tshape::TIndexVec TIndexVec;
      typedef Tshape::TSizeVec TSizeVec;
      explicit FortranShape(const Tshape& shape, const bool& BaseOne=true);
      const Tsubscript& first(const Tdim& i) const { return(Mfirst[i]); }
      const Tsubscript& last(const Tdim& i) const { return(Mlast[i]); }
      const Tsubscript& dimlast(const Tdim& i) const { return(Mdimlast[i]); }
      const Tsubscript& offset() const { return(Moffset); }
      const TIndexVec& first() const { return(Mfirst); }
      const TIndexVec& last() const { return(Mlast); }
      const TIndexVec& dimlast() const { return(Mdimlast); }
    private:
      TIndexVec Mfirst;
      TIndexVec Mlast;
      TIndexVec Mdimlast;
      Tsubscript Moffset;
  }; // class FortranShape
} // namespace util

/*----------------------------------------------------------------------*/

template<class C>
class FortranArray: private aff::util::FortranShape {
  public:
    typedef C Tarray;
    typedef typename Tarray::Tshape Tshape;
    typedef typename Tarray::Trepresentation Trepresentation;
    typedef typename Trepresentation::Tpointer Treppointer;
    typedef aff::util::FortranShape Tbase;
    typedef typename Tarray::Tpointer Tpointer;
    typedef typename Tarray::Tvalue Tvalue;
    typedef typename Trepresentation::Tvalue Trepvalue;
    FortranArray(Tarray array, const bool& BaseOne=true):
      Tbase(array.shape(), BaseOne)
      {
        Trepresentation repr=array.representation();
        Mpointer=&repr[this->Tbase::offset()];
      }

    Tbase::first;
    Tbase::last;
    Tbase::dimlast;
    Treppointer pointer() const { return(Mpointer); }
    template<class TT>
      TT* castedpointer() const
      { return(SizeCheckedCast<Trepvalue,TT>::cast(Mpointer)); }
  private:
    Treppointer Mpointer;
}; // class FortranArray

} // namespace aff

#endif // AFF_FORTRANSHAPE_H_VERSION (includeguard)

/* ----- END OF fortranshape.h ----- */