#include "IMP/MultiArray.h"

template<>
std::ostream& operator<<(std::ostream& ostrm,
                         const MultiArray< std::complex<double> >& v)
{
  if (v.xdim_ == 0) {
    ostrm << "NULL Array\n";
  }
  else {
    ostrm << std::endl;
  }

  for (int k = STARTINGZ(v); k <= FINISHINGZ(v); k++) {
    ostrm << "Slice No. " << k << std::endl;
    for (int i = STARTINGY(v); i <= FINISHINGY(v); i++) {
      for (int j = STARTINGX(v); j <= FINISHINGX(v); j++) {
        ostrm << VOL_ELEM(v, k, i, j) << " ";
      }
      ostrm << std::endl;
    }
  }
  return ostrm;
}