testArray.cc

Go to the documentation of this file.

//
// tests matrix implementation using vector of vector
// run using:
// g++ -Wall testArray.cxx -o testArray && ./testArray
//


#include <iostream>
#include <complex>

#include "boost/multi_array.hpp"

#include "arrayUtils.hpp"
#include "reportingUtils.hpp"


using namespace std;
using namespace boost;
using namespace rpwa;


int main(int argc, char** argv)
{
  double a[2][3][4];
  for (unsigned int i = 0; i < 2; ++i)
    for (unsigned int j = 0; j < 3; ++j)
      for (unsigned int k = 0; k < 4; ++k)
        a[i][j][k] = i * 100 + j * 10 + k;

  if (0) {
  
    double***          x      = 0;
    const unsigned int dim[3] = {2, 3, 4};
    allocate3DArray(x, dim);
    for (unsigned int i = 0; i < 2; ++i)
      for (unsigned int j = 0; j < 3; ++j)
        for (unsigned int k = 0; k < 4; ++k)
          x[i][j][k] = i * 100 + j * 10 + k;

    for (unsigned int i = 0; i < 2; ++i)
      for (unsigned int j = 0; j < 3; ++j)
        for (unsigned int k = 0; k < 4; ++k) {
          if (a[i][j][k] != a[i][j][k])
            cout << "ERROR!" << endl;
          cout << "x[" << i << "][" << j << "][" << k << "] = " << x[i][j][k] << " vs. "
               << "a[" << i << "][" << j << "][" << k << "] = " << a[i][j][k] << endl;
        }

    delete3DArray(x, dim);
    
  }


  if (0) {

    double*            x      = 0;
    const unsigned int dim[3] = {2, 3, 4};
    const unsigned int nmbDim = 3;
    allocatePseudoNdimArray(x, dim, nmbDim);
    for (unsigned int i = 0; i < 2; ++i)
      for (unsigned int j = 0; j < 3; ++j)
        for (unsigned int k = 0; k < 4; ++k) {
          const unsigned int indices[3] = {i, j, k};
          const unsigned int offset     = indicesToOffset(indices, dim, nmbDim);
          x[offset] = i * 100 + j * 10 + k;
        }

    double* y = (double*)a;
    for (unsigned int i = 0; i < 2; ++i)
      for (unsigned int j = 0; j < 3; ++j)
        for (unsigned int k = 0; k < 4; ++k) {
          const unsigned int indices[3] = {i, j, k};
          const unsigned int offset     = indicesToOffset(indices, dim, nmbDim);
          if (x[offset] != a[i][j][k])
            cout << "ERROR!" << endl;
          cout << "x[" << offset << "] = " << x[offset] << " vs. "
               << "a[" << i << "][" << j << "][" << k << "] = " << a[i][j][k] << endl;
          if (x[offset] != y[offset])
            cout << "ERROR!" << endl;
          cout << "x[" << offset << "] = " << x[offset] << " vs. "
               << "y[" << offset << "] = " << y[offset] << endl;
          unsigned int testIndices[3];
          offsetToIndices(offset, dim, nmbDim, testIndices);
          cout << " indices {";
          for (unsigned int d = 0; d < nmbDim; ++d) {
            if (testIndices[d] != indices[d])
              cout << "ERROR!  ";
            cout << testIndices[d] << "  ";
          }
          cout << "} vs. {" << i << "  " << j << "  " << k << "}" << endl;
        }
    
    delete[] x;
  }

  if (1) {

    const unsigned int dim[3] = {1000, 300, 400};
    // double             a[dim[0]][dim[1]][dim[2]];
    double***          x      = 0;
    allocate3DArray(x, dim);
    double*            y      = 0;
    const unsigned int nmbDim = 3;
    allocatePseudoNdimArray(y, dim, nmbDim);
    multi_array<double, 3> z(extents[dim[0]][dim[1]][dim[2]]);

    for (unsigned int i = 0; i < dim[0]; ++i)
      for (unsigned int j = 0; j < dim[1]; ++j)
        for (unsigned int k = 0; k < dim[2]; ++k) {
          const unsigned int indices[3] = {i, j, k};
          const unsigned int offset     = indicesToOffset(indices, dim, nmbDim);
          // if (a[i][j][k] != 0)
          //   cout << "a ";
          if (x[i][j][k] != 0)
            cout << "x ";
          if (y[offset] != 0)
            cout << "y ";
          if (z[i][j][k] != 0)
            cout << "z ";
        }
    cout << endl;
  }
  
  return 0;
}