dump.cc

Go to the documentation of this file.

#define AFF_DUMP_CC_VERSION \
  "AFF_DUMP_CC   V1.1"

#include <aff/lib/dump_strided.h>
#include <aff/lib/error.h>

namespace aff {

void dump(const Strided& shape, std::ostream& os)
{
  os << "dump of a Strided object:" << endl;
  os << "  index ranges: ";
  for (Tsize i=0; i<Strided::Mmax_dimen; i++)
  {
    os << "[" << shape.first(i) << ":" << shape.last(i) << "]";
  }
  os << endl;
  os << "  strides: ";
  for (Tsize i=0; i<Strided::Mmax_dimen; i++)
  {
    if (i>0) os << ", ";
    os << shape.stride(i);
  }
  os << endl;
  os << "  total number of mapped elements: "
    << shape.size() << endl;
  os << "  total address range in memory: ["
    << shape.first_offset() << ":" << shape.last_offset() << "] = "
    << shape.memory_size() << " positions" << endl;
}

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

namespace {

  void dump_map_helper(const Strided& shape, const Tdim& i, const Tdim& j,
                  Strided::TIndexVec& index, std::ostream& os)
  {
    if (i>1)
    {
      for (index[i]=shape.first(i); index[i]<=shape.last(i); index[i]++)
      {
        os << "    dimension "; 
        os << i << ": [" << index[i] << "]";
        dump_map_helper(shape, i-1, j, index, os);
      }
    }
    else if (i>0)
    {
      if (i<j) os << endl;
      os.width(5); os << " ";
      for (int k=shape.first(1); k<=shape.last(1); k++)
      {
        os.width(6); os << k;
      }
      os << endl;
      for (index[0]=shape.first(0); index[0]<=shape.last(0); index[0]++)
      {
        os.width(5); os << index[0];
        for (index[1]=shape.first(1); index[1]<=shape.last(1); index[1]++)
        {
          os.width(5); 
          if (j==0)
          { os << shape.offset(index[0]); }
          else if (j==1)
          { os << shape.offset(index[0], index[1]); }
          else if (j==2)
          { os << shape.offset(index[0], index[1], index[2]); }
          else if (j==3)
          { os << shape.offset(index[0], index[1], index[2], index[3]); }
          else 
          { os << shape.offset(index); }
          os << "#";
        }
        os << endl;
      }
    }
    else 
    {
      for (index[0]=shape.first(0); index[0]<=shape.last(0); index[0]++)
      {
        os.width(5); os << index[0];
        os.width(5); 
        if (j==0)
        { os << shape.offset(index[0]); }
        else if (j==1)
        { os << shape.offset(index[0], index[1]); }
        else if (j==2)
        { os << shape.offset(index[0], index[1], index[2]); }
        else if (j==3)
        { os << shape.offset(index[0], index[1], index[2], index[3]); }
        else 
        { os << shape.offset(index); }
        os << "#" << endl;
      }
    }
  }
}

void dump_map(const Strided& shape, const Tdim& i, std::ostream& os)
{
  AFF_assert((i<Strided::Mmax_dimen),
             "ERROR (dump_map): illegal dimensionality");
  dump(shape);
  os << "  index mapping for " << i+1 << " dimensional access:" << endl;
  Strided::TIndexVec index(shape.first());
  dump_map_helper(shape, i, i, index, os); 
}

} // namespace aff

/* ----- END OF dump.cc ----- */