removed .cpp and .h files

1 files changed, 0 insertions, 1205 deletions

diff --git a/video_input/consistencyChecker/CTensor.h b/video_input/consistencyChecker/CTensor.h
deleted file mode 100644
index 0f5af5c..0000000
--- a/video_input/consistencyChecker/CTensor.h
+++ /dev/null
@@ -1,1205 +0,0 @@
-// CTensor

-// A three-dimensional array

-//

-// Author: Thomas Brox

-

-#ifndef CTENSOR_H

-#define CTENSOR_H

-

-#include <iostream>

-#include <fstream>

-#include <string>

-#include <sstream>

-#include <CMatrix.h>

-#include <NMath.h>

-

-inline int int_min(int x, int& y) { return (x<y)?x:y; }

-inline int int_max(int x, int& y) { return (x<y)?y:x; }

-

-template <class T>

-class CTensor {

-public:

-  // standard constructor

-  inline CTensor();

-  // constructor

-  inline CTensor(const int aXSize, const int aYSize, const int aZSize);

-  // copy constructor

-  CTensor(const CTensor<T>& aCopyFrom);

-  // constructor with implicit filling

-  CTensor(const int aXSize, const int aYSize, const int aZSize, const T aFillValue);

-  // destructor

-  virtual ~CTensor();

-

-  // Changes the size of the tensor, data will be lost

-  void setSize(int aXSize, int aYSize, int aZSize);

-  // Downsamples the tensor

-  void downsample(int aNewXSize, int aNewYSize);

-  void downsample(int aNewXSize, int aNewYSize, CMatrix<float>& aConfidence);

-  void downsample(int aNewXSize, int aNewYSize, CTensor<float>& aConfidence);

-  // Upsamples the tensor

-  void upsample(int aNewXSize, int aNewYSize);

-  void upsampleBilinear(int aNewXSize, int aNewYSize);

-  // Fills the tensor with the value aValue (see also operator =)

-  void fill(const T aValue);

-  // Fills a rectangular area with the value aValue

-  void fillRect(const CVector<T>& aValue, int ax1, int ay1, int ax2, int ay2);

-  // Copies a box from the tensor into aResult, the size of aResult will be adjusted

-  void cut(CTensor<T>& aResult, int x1, int y1, int z1, int x2, int y2, int z2);

-  // Copies aCopyFrom at a certain position of the tensor

-  void paste(CTensor<T>& aCopyFrom, int ax, int ay, int az);

-  // Mirrors the boundaries, aFrom is the distance from the boundaries where the pixels are copied from,

-  // aTo is the distance from the boundaries they are copied to

-  void mirrorLayers(int aFrom, int aTo);

-  // Transforms the values so that they are all between aMin and aMax

-  // aInitialMin/Max are initializations for seeking the minimum and maximum, change if your

-  // data is not in this range or the data type T cannot hold these values

-  void normalizeEach(T aMin, T aMax, T aInitialMin = -30000, T aInitialMax = 30000);

-  void normalize(T aMin, T aMax, int aChannel, T aInitialMin = -30000, T aInitialMax = 30000);

-  void normalize(T aMin, T aMax, T aInitialMin = -30000, T aInitialMax = 30000);

-  // Converts from RGB to CIELab color space and vice-versa

-  void rgbToCielab();

-  void cielabToRGB();

-  // Draws a line into the image (only for mZSize = 3)

-  void drawLine(int dStartX, int dStartY, int dEndX, int dEndY, T aValue1 = 255, T aValue2 = 255, T aValue3 = 255);

-  void drawRect(int dStartX, int dStartY, int dEndX, int dEndY, T aValue1 = 255, T aValue2 = 255, T aValue3 = 255);

-

-  // Applies a similarity transform (translation, rotation, scaling) to the image

-  void applySimilarityTransform(CTensor<T>& aWarped, CMatrix<bool>& aOutside, float tx, float ty, float cx, float cy, float phi, float scale);

-  // Applies a homography (linear projective transformation) to the image

-  void applyHomography(CTensor<T>& aWarped, CMatrix<bool>& aOutside, const CMatrix<float>& H);

-

-  // Reads the tensor from a file in Mathematica format

-  void readFromMathematicaFile(const char* aFilename);

-  // Writes the tensor to a file in Mathematica format

-  void writeToMathematicaFile(const char* aFilename);

-  // Reads the tensor from a movie file in IM format

-  void readFromIMFile(const char* aFilename);

-  // Writes the tensor to a movie file in IM format

-  void writeToIMFile(const char* aFilename);

-  // Reads an image from a PGM file

-  void readFromPGM(const char* aFilename);

-  // Writes the tensor in PGM-Format

-  void writeToPGM(const char* aFilename);

-  // Extends a XxYx1 tensor to a XxYx3 tensor with three identical layers

-  void makeColorTensor();

-  // Reads a color image from a PPM file

-  void readFromPPM(const char* aFilename);

-  // Writes the tensor in PPM-Format

-  void writeToPPM(const char* aFilename);

-  // Reads the tensor from a PDM file

-  void readFromPDM(const char* aFilename);

-  // Writes the tensor in PDM-Format

-  void writeToPDM(const char* aFilename, char aFeatureType);

-

-  // Gives full access to tensor's values

-  inline T& operator()(const int ax, const int ay, const int az) const;

-  // Read access with bilinear interpolation

-  CVector<T> operator()(const float ax, const float ay) const;

-  // Fills the tensor with the value aValue (equivalent to fill())

-  inline CTensor<T>& operator=(const T aValue);

-  // Copies the tensor aCopyFrom to this tensor (size of tensor might change)

-  CTensor<T>& operator=(const CTensor<T>& aCopyFrom);

-  // Adds a tensor of same size

-  CTensor<T>& operator+=(const CTensor<T>& aMatrix);

-  // Adds a constant to the tensor

-  CTensor<T>& operator+=(const T aValue);

-  // Multiplication with a scalar

-  CTensor<T>& operator*=(const T aValue);

-

-  // Returns the minimum value

-  T min() const;

-  // Returns the maximum value

-  T max() const;

-  // Returns the average value

-  T avg() const;

-  // Returns the average value of a specific layer

-  T avg(int az) const;

-  // Gives access to the tensor's size

-  inline int xSize() const;

-  inline int ySize() const;

-  inline int zSize() const;

-  inline int size() const;

-  // Returns the az layer of the tensor as matrix (slow and fast version)

-  CMatrix<T> getMatrix(const int az) const;

-  void getMatrix(CMatrix<T>& aMatrix, const int az) const;

-  // Copies the matrix components of aMatrix into the az layer of the tensor

-  void putMatrix(CMatrix<T>& aMatrix, const int az);

-  // Gives access to the internal data representation (use sparingly)

-  inline T* data() const;

-

-  // Possible interpretations of the third tensor dimension for PDM format

-  static const char cSpacial = 'S';

-  static const char cVector = 'V';

-  static const char cColor = 'C';

-  static const char cSymmetricMatrix = 'Y';

-protected:

-  int mXSize,mYSize,mZSize;

-  T *mData;

-};

-

-// Provides basic output functionality (only appropriate for very small tensors)

-template <class T> std::ostream& operator<<(std::ostream& aStream, const CTensor<T>& aTensor);

-

-// Exceptions thrown by CTensor-------------------------------------------------

-

-// Thrown when one tries to access an element of a tensor which is out of

-// the tensor's bounds

-struct ETensorRangeOverflow {

-  ETensorRangeOverflow(const int ax, const int ay, const int az) {

-    using namespace std;

-    cerr << "Exception ETensorRangeOverflow: x = " << ax << ", y = " << ay << ", z = " << az << endl;

-  }

-};

-

-// Thrown when the size of a tensor does not match the needed size for a certain operation

-struct ETensorIncompatibleSize {

-  ETensorIncompatibleSize(int ax, int ay, int ax2, int ay2) {

-    using namespace std;

-    cerr << "Exception ETensorIncompatibleSize: x = " << ax << ":" << ax2;

-    cerr << ", y = " << ay << ":" << ay2 << endl;

-  }

-  ETensorIncompatibleSize(int ax, int ay, int az) {

-    std::cerr << "Exception ETensorIncompatibleTensorSize: x = " << ax << ", y = " << ay << ", z= " << az << std::endl;

-  }

-};

-

-// I M P L E M E N T A T I O N --------------------------------------------

-//

-// You might wonder why there is implementation code in a header file.

-// The reason is that not all C++ compilers yet manage separate compilation

-// of templates. Inline functions cannot be compiled separately anyway.

-// So in this case the whole implementation code is added to the header

-// file.

-// Users of CTensor should ignore everything that's beyond this line :)

-// ------------------------------------------------------------------------

-

-// P U B L I C ------------------------------------------------------------

-

-// standard constructor

-template <class T>

-inline CTensor<T>::CTensor() {

-  mData = 0;

-  mXSize = mYSize = mZSize = 0;

-}

-

-// constructor

-template <class T>

-inline CTensor<T>::CTensor(const int aXSize, const int aYSize, const int aZSize)

-  : mXSize(aXSize), mYSize(aYSize), mZSize(aZSize) {

-  mData = new T[aXSize*aYSize*aZSize];

-}

-

-// copy constructor

-template <class T>

-CTensor<T>::CTensor(const CTensor<T>& aCopyFrom)

-  : mXSize(aCopyFrom.mXSize), mYSize(aCopyFrom.mYSize), mZSize(aCopyFrom.mZSize) {

-  int wholeSize = mXSize*mYSize*mZSize;

-  mData = new T[wholeSize];

-  for (register int i = 0; i < wholeSize; i++)

-    mData[i] = aCopyFrom.mData[i];

-}

-

-// constructor with implicit filling

-template <class T>

-CTensor<T>::CTensor(const int aXSize, const int aYSize, const int aZSize, const T aFillValue)

-  : mXSize(aXSize), mYSize(aYSize), mZSize(aZSize) {

-  mData = new T[aXSize*aYSize*aZSize];

-  fill(aFillValue);

-}

-

-// destructor

-template <class T>

-CTensor<T>::~CTensor() {

-  delete[] mData;

-}

-

-// setSize

-template <class T>

-void CTensor<T>::setSize(int aXSize, int aYSize, int aZSize) {

-  if (mData != 0) delete[] mData;

-  mData = new T[aXSize*aYSize*aZSize];

-  mXSize = aXSize;

-  mYSize = aYSize;

-  mZSize = aZSize;

-}

-

-//downsample

-template <class T>

-void CTensor<T>::downsample(int aNewXSize, int aNewYSize) {

-  T* mData2 = new T[aNewXSize*aNewYSize*mZSize];

-  int aSize = aNewXSize*aNewYSize;

-  for (int z = 0; z < mZSize; z++) {

-    CMatrix<T> aTemp(mXSize,mYSize);

-    getMatrix(aTemp,z);

-    aTemp.downsample(aNewXSize,aNewYSize);

-    for (int i = 0; i < aSize; i++)

-      mData2[i+z*aSize] = aTemp.data()[i];

-  }

-  delete[] mData;

-  mData = mData2;

-  mXSize = aNewXSize;

-  mYSize = aNewYSize;

-}

-

-template <class T>

-void CTensor<T>::downsample(int aNewXSize, int aNewYSize, CMatrix<float>& aConfidence) {

-  T* mData2 = new T[aNewXSize*aNewYSize*mZSize];

-  int aSize = aNewXSize*aNewYSize;

-  for (int z = 0; z < mZSize; z++) {

-    CMatrix<T> aTemp(mXSize,mYSize);

-    getMatrix(aTemp,z);

-    aTemp.downsample(aNewXSize,aNewYSize,aConfidence);

-    for (int i = 0; i < aSize; i++)

-      mData2[i+z*aSize] = aTemp.data()[i];

-  }

-  delete[] mData;

-  mData = mData2;

-  mXSize = aNewXSize;

-  mYSize = aNewYSize;

-}

-

-template <class T>

-void CTensor<T>::downsample(int aNewXSize, int aNewYSize, CTensor<float>& aConfidence) {

-  T* mData2 = new T[aNewXSize*aNewYSize*mZSize];

-  int aSize = aNewXSize*aNewYSize;

-  CMatrix<float> aConf(mXSize,mYSize);

-  for (int z = 0; z < mZSize; z++) {

-    CMatrix<T> aTemp(mXSize,mYSize);

-    getMatrix(aTemp,z);

-    aConfidence.getMatrix(aConf,z);

-    aTemp.downsample(aNewXSize,aNewYSize,aConf);

-    for (int i = 0; i < aSize; i++)

-      mData2[i+z*aSize] = aTemp.data()[i];

-  }

-  delete[] mData;

-  mData = mData2;

-  mXSize = aNewXSize;

-  mYSize = aNewYSize;

-}

-

-// upsample

-template <class T>

-void CTensor<T>::upsample(int aNewXSize, int aNewYSize) {

-  T* mData2 = new T[aNewXSize*aNewYSize*mZSize];

-  int aSize = aNewXSize*aNewYSize;

-  for (int z = 0; z < mZSize; z++) {

-    CMatrix<T> aTemp(mXSize,mYSize);

-    getMatrix(aTemp,z);

-    aTemp.upsample(aNewXSize,aNewYSize);

-    for (int i = 0; i < aSize; i++)

-      mData2[i+z*aSize] = aTemp.data()[i];

-  }

-  delete[] mData;

-  mData = mData2;

-  mXSize = aNewXSize;

-  mYSize = aNewYSize;

-}

-

-// upsampleBilinear

-template <class T>

-void CTensor<T>::upsampleBilinear(int aNewXSize, int aNewYSize) {

-  T* mData2 = new T[aNewXSize*aNewYSize*mZSize];

-  int aSize = aNewXSize*aNewYSize;

-  for (int z = 0; z < mZSize; z++) {

-    CMatrix<T> aTemp(mXSize,mYSize);

-    getMatrix(aTemp,z);

-    aTemp.upsampleBilinear(aNewXSize,aNewYSize);

-    for (int i = 0; i < aSize; i++)

-      mData2[i+z*aSize] = aTemp.data()[i];

-  }

-  delete[] mData;

-  mData = mData2;

-  mXSize = aNewXSize;

-  mYSize = aNewYSize;

-}

-

-// fill

-template <class T>

-void CTensor<T>::fill(const T aValue) {

-  int wholeSize = mXSize*mYSize*mZSize;

-  for (register int i = 0; i < wholeSize; i++)

-    mData[i] = aValue;

-}

-

-// fillRect

-template <class T>

-void CTensor<T>::fillRect(const CVector<T>& aValue, int ax1, int ay1, int ax2, int ay2) {

-  for (int z = 0; z < mZSize; z++) {

-    T val = aValue(z);

-    for (int y = int_max(0,ay1); y <= int_min(ySize()-1,ay2); y++)

-      for (register int x = int_max(0,ax1); x <= int_min(xSize()-1,ax2); x++)

-        operator()(x,y,z) = val;

-  }

-}

-

-// cut

-template <class T>

-void CTensor<T>::cut(CTensor<T>& aResult, int x1, int y1, int z1, int x2, int y2, int z2) {

-  aResult.mXSize = x2-x1+1;

-  aResult.mYSize = y2-y1+1;

-  aResult.mZSize = z2-z1+1;

-  delete[] aResult.mData;

-  aResult.mData = new T[aResult.mXSize*aResult.mYSize*aResult.mZSize];

-  for (int z = z1; z <= z2; z++)

-    for (int y = y1; y <= y2; y++)

-      for (int x = x1; x <= x2; x++)

-        aResult(x-x1,y-y1,z-z1) = operator()(x,y,z);

-}

-

-// paste

-template <class T>

-void CTensor<T>::paste(CTensor<T>& aCopyFrom, int ax, int ay, int az) {

-  for (int z = 0; z < aCopyFrom.zSize(); z++)

-    for (int y = 0; y < aCopyFrom.ySize(); y++)

-      for (int x = 0; x < aCopyFrom.xSize(); x++)

-        operator()(ax+x,ay+y,az+z) = aCopyFrom(x,y,z);

-}

-

-// mirrorLayers

-template <class T>

-void CTensor<T>::mirrorLayers(int aFrom, int aTo) {

-  for (int z = 0; z < mZSize; z++) {

-    int aToXIndex = mXSize-aTo-1;

-    int aToYIndex = mYSize-aTo-1;

-    int aFromXIndex = mXSize-aFrom-1;

-    int aFromYIndex = mYSize-aFrom-1;

-    for (int y = aFrom; y <= aFromYIndex; y++) {

-      operator()(aTo,y,z) = operator()(aFrom,y,z);

-      operator()(aToXIndex,y,z) = operator()(aFromXIndex,y,z);

-    }

-    for (int x = aTo; x <= aToXIndex; x++) {

-      operator()(x,aTo,z) = operator()(x,aFrom,z);

-      operator()(x,aToYIndex,z) = operator()(x,aFromYIndex,z);

-    }

-  }

-}

-

-// normalize

-template <class T>

-void CTensor<T>::normalizeEach(T aMin, T aMax, T aInitialMin, T aInitialMax) {

-  for (int k = 0; k < mZSize; k++)

-    normalize(aMin,aMax,k,aInitialMin,aInitialMax);

-}

-

-template <class T>

-void CTensor<T>::normalize(T aMin, T aMax, int aChannel, T aInitialMin, T aInitialMax) {

-  int aChannelSize = mXSize*mYSize;

-  T aCurrentMin = aInitialMax;

-  T aCurrentMax = aInitialMin;

-  int aIndex = aChannelSize*aChannel;

-  for (int i = 0; i < aChannelSize; i++) {

-    if (mData[aIndex] > aCurrentMax) aCurrentMax = mData[aIndex];

-    else if (mData[aIndex] < aCurrentMin) aCurrentMin = mData[aIndex];

-    aIndex++;

-  }

-  T aTemp1 = aCurrentMin - aMin;

-  T aTemp2 = (aCurrentMax-aCurrentMin);

-  if (aTemp2 == 0) aTemp2 = 1;

-  else aTemp2 = (aMax-aMin)/aTemp2;

-  aIndex = aChannelSize*aChannel;

-  for (int i = 0; i < aChannelSize; i++) {

-    mData[aIndex] -= aTemp1;

-    mData[aIndex] *= aTemp2;

-    aIndex++;

-  }

-}

-

-// drawLine

-template <class T>

-void CTensor<T>::drawLine(int dStartX, int dStartY, int dEndX, int dEndY, T aValue1, T aValue2, T aValue3) {

-  int aOffset1 = mXSize*mYSize;

-  int aOffset2 = 2*aOffset1;

-	// vertical line

-	if (dStartX == dEndX) {

-    if (dStartX < 0 || dStartX >= mXSize)	return;

-		int x = dStartX;

-		if (dStartY < dEndY) {

-			for (int y = dStartY; y <= dEndY; y++)

-				if (y >= 0 && y < mYSize) {

-          mData[x+y*mXSize] = aValue1;

-          mData[x+y*mXSize+aOffset1] = aValue2;

-          mData[x+y*mXSize+aOffset2] = aValue3;

-        }

-  	}

-		else {

-			for (int y = dStartY; y >= dEndY; y--)

-				if (y >= 0 && y < mYSize) {

-          mData[x+y*mXSize] = aValue1;

-          mData[x+y*mXSize+aOffset1] = aValue2;

-          mData[x+y*mXSize+aOffset2] = aValue3;

-        }

-    }

-    return;

-  }

-	// horizontal line

-	if (dStartY == dEndY) {

-    if (dStartY < 0 || dStartY >= mYSize) return;

- 		int y = dStartY;

-		if (dStartX < dEndX) {

-			for (int x = dStartX; x <= dEndX; x++)

-				if (x >= 0 && x < mXSize) {

-          mData[x+y*mXSize] = aValue1;

-          mData[x+y*mXSize+aOffset1] = aValue2;

-          mData[x+y*mXSize+aOffset2] = aValue3;

-        }

-  	}

-		else {

-			for (int x = dStartX; x >= dEndX; x--)

-				if (x >= 0 && x < mXSize) {

-          mData[x+y*mXSize] = aValue1;

-          mData[x+y*mXSize+aOffset1] = aValue2;

-          mData[x+y*mXSize+aOffset2] = aValue3;

-        }

-    }

-    return;

-  }

-  float m = float(dStartY - dEndY) / float(dStartX - dEndX);

-  float invm = 1.0/m;

-  if (fabs(m) > 1.0) {

-    if (dEndY > dStartY) {

-      for (int y = dStartY; y <= dEndY; y++) {

-        int x = (int)(0.5+dStartX+(y-dStartY)*invm);

-        if (x >= 0 && x < mXSize &&	y >= 0 && y < mYSize) {

-          mData[x+y*mXSize] = aValue1;

-          mData[x+y*mXSize+aOffset1] = aValue2;

-          mData[x+y*mXSize+aOffset2] = aValue3;

-        }

-      }

-    }

-    else {

-      for (int y = dStartY; y >= dEndY; y--) {

-        int x = (int)(0.5+dStartX+(y-dStartY)*invm);

-        if (x >= 0 && x < mXSize &&	y >= 0 && y < mYSize) {

-          mData[x+y*mXSize] = aValue1;

-          mData[x+y*mXSize+aOffset1] = aValue2;

-          mData[x+y*mXSize+aOffset2] = aValue3;

-        }

-      }

-    }

-  }

-  else {

-    if (dEndX > dStartX) {

-      for (int x = dStartX; x <= dEndX; x++) {

-        int y = (int)(0.5+dStartY+(x-dStartX)*m);

-        if (x >= 0 && x < mXSize &&	y >= 0 && y < mYSize) {

-          mData[x+y*mXSize] = aValue1;

-          mData[x+y*mXSize+aOffset1] = aValue2;

-          mData[x+y*mXSize+aOffset2] = aValue3;

-        }

-      }

-    }

-    else {

-      for (int x = dStartX; x >= dEndX; x--) {

-        int y = (int)(0.5+dStartY+(x-dStartX)*m);

-        if (x >= 0 && x < mXSize &&	y >= 0 && y < mYSize) {

-          mData[x+y*mXSize] = aValue1;

-          mData[x+y*mXSize+aOffset1] = aValue2;

-          mData[x+y*mXSize+aOffset2] = aValue3;

-        }

-      }

-    }

-  }

-}

-

-// drawRect

-template <class T>

-void CTensor<T>::drawRect(int dStartX, int dStartY, int dEndX, int dEndY, T aValue1, T aValue2, T aValue3) {

-  drawLine(dStartX,dStartY,dEndX,dStartY,aValue1,aValue2,aValue3);

-  drawLine(dStartX,dEndY,dEndX,dEndY,aValue1,aValue2,aValue3);

-  drawLine(dStartX,dStartY,dStartX,dEndY,aValue1,aValue2,aValue3);

-  drawLine(dEndX,dStartY,dEndX,dEndY,aValue1,aValue2,aValue3);

-}

-

-template <class T>

-void CTensor<T>::normalize(T aMin, T aMax, T aInitialMin, T aInitialMax) {

-  int aSize = mXSize*mYSize*mZSize;

-  T aCurrentMin = aInitialMax;

-  T aCurrentMax = aInitialMin;

-  for (int i = 0; i < aSize; i++) {

-    if (mData[i] > aCurrentMax) aCurrentMax = mData[i];

-    else if (mData[i] < aCurrentMin) aCurrentMin = mData[i];

-  }

-  T aTemp1 = aCurrentMin - aMin;

-  T aTemp2 = (aCurrentMax-aCurrentMin);

-  if (aTemp2 == 0) aTemp2 = 1;

-  else aTemp2 = (aMax-aMin)/aTemp2;

-  for (int i = 0; i < aSize; i++) {

-    mData[i] -= aTemp1;

-    mData[i] *= aTemp2;

-  }

-}

-

-template <class T>

-void CTensor<T>::rgbToCielab() {

-  for (int y = 0; y < mYSize; y++)

-    for (int x = 0; x < mXSize; x++) {

-      float R = operator()(x,y,0)*0.003921569;

-      float G = operator()(x,y,1)*0.003921569;

-      float B = operator()(x,y,2)*0.003921569;

-      if (R>0.0031308) R = pow((R + 0.055)*0.9478673, 2.4); else R *= 0.077399381;

-      if (G>0.0031308) G = pow((G + 0.055)*0.9478673, 2.4); else G *= 0.077399381;

-      if (B>0.0031308) B = pow((B + 0.055)*0.9478673, 2.4); else B *= 0.077399381;

-      //Observer. = 2?, Illuminant = D65

-      float X = R * 0.4124 + G * 0.3576 + B * 0.1805;

-      float Y = R * 0.2126 + G * 0.7152 + B * 0.0722;

-      float Z = R * 0.0193 + G * 0.1192 + B * 0.9505;

-      X *= 1.052111;

-      Z *= 0.918417;

-      if (X > 0.008856) X = pow(X,0.33333333333); else X = 7.787*X + 0.137931034;

-      if (Y > 0.008856) Y = pow(Y,0.33333333333); else Y = 7.787*Y + 0.137931034;

-      if (Z > 0.008856) Z = pow(Z,0.33333333333); else Z = 7.787*Z + 0.137931034;

-      operator()(x,y,0) = 1000.0*((295.8*Y) - 40.8)/255.0;

-      operator()(x,y,1) = 128.0+637.5*(X-Y);

-      operator()(x,y,2) = 128.0+255.0*(Y-Z);

-    }

-}

-

-template <class T>

-void CTensor<T>::cielabToRGB() {

-  for (int y = 0; y < mYSize; y++)

-    for (int x = 0; x < mXSize; x++) {

-      float L = operator()(x,y,0)*0.255;

-      float A = operator()(x,y,1);

-      float B = operator()(x,y,2);

-      float Y = (L+40.8)*0.00338066;

-      float X = (A-128.0+637.5*Y)*0.0015686;

-      float Z = (128.0+255.0*Y-B)*0.00392157;

-      float temp = Y*Y*Y;

-      if (temp > 0.008856) Y = temp;

-      else Y = (Y-0.137931034)*0.12842;

-      temp = X*X*X;

-      if (temp > 0.008856) X = temp;

-      else X = (X-0.137931034)*0.12842;

-      temp = Z*Z*Z;

-      if (temp > 0.008856) Z = temp;

-      else Z = (Z-0.137931034)*0.12842;

-      X *= 0.95047;

-      Y *= 1.0;

-      Z *= 1.08883;

-      float r = 3.2406*X-1.5372*Y-0.4986*Z;

-      float g = -0.9689*X+1.8758*Y+0.0415*Z;

-      float b = 0.0557*X-0.204*Y+1.057*Z;

-      if (r < 0) r = 0;

-      temp = 1.055*pow(r,0.41667)-0.055;

-      if (temp > 0.0031308) r = temp;

-      else r *= 12.92;

-      if (g < 0) g = 0;

-      temp = 1.055*pow(g,0.41667)-0.055;

-      if (temp > 0.0031308) g = temp;

-      else g *= 12.92;

-      if (b < 0) b = 0;

-      temp = 1.055*pow(b,0.41667)-0.055;

-      if (temp > 0.0031308) b = temp;

-      else b *= 12.92;

-      operator()(x,y,0) = 255.0*r;

-      operator()(x,y,1) = 255.0*g;

-      operator()(x,y,2) = 255.0*b;

-    }

-}

-

-// applySimilarityTransform

-template <class T>

-void CTensor<T>::applySimilarityTransform(CTensor<T>& aWarped, CMatrix<bool>& aOutside, float tx, float ty, float cx, float cy, float phi, float scale) {

-  float cosphi = scale*cos(phi);

-  float sinphi = scale*sin(phi);

-  int aSize = mXSize*mYSize;

-  int aWarpedSize = aWarped.xSize()*aWarped.ySize();

-  float ctx = cx+tx-cx*cosphi+cy*sinphi;

-  float cty = cy+ty-cy*cosphi-cx*sinphi;

-  aOutside = false;

-  int i = 0;

-  for (int y = 0; y < aWarped.ySize(); y++)

-    for (int x = 0; x < aWarped.xSize(); x++,i++) {

-      float xf = x; float yf = y;

-      float ax = xf*cosphi-yf*sinphi+ctx;

-      float ay = yf*cosphi+xf*sinphi+cty;

-      int x1 = (int)ax; int y1 = (int)ay;

-      float alphaX = ax-x1; float alphaY = ay-y1;

-      float betaX = 1.0-alphaX; float betaY = 1.0-alphaY;

-      if (x1 < 0 || y1 < 0 || x1+1 >= mXSize || y1+1 >= mYSize) aOutside.data()[i] = true;

-      else {

-        int j = y1*mXSize+x1;

-        for (int k = 0; k < mZSize; k++) {

-          float a = betaX*mData[j]       +alphaX*mData[j+1];

-          float b = betaX*mData[j+mXSize]+alphaX*mData[j+1+mXSize];

-          aWarped.data()[i+k*aWarpedSize] = betaY*a+alphaY*b;

-          j += aSize;

-        }

-      }

-    }

-}

-

-// applyHomography

-template <class T>

-void CTensor<T>::applyHomography(CTensor<T>& aWarped, CMatrix<bool>& aOutside, const CMatrix<float>& H) {

-  int aSize = mXSize*mYSize;

-  int aWarpedSize = aWarped.xSize()*aWarped.ySize();

-  aOutside = false;

-  int i = 0;

-  for (int y = 0; y < aWarped.ySize(); y++)

-    for (int x = 0; x < aWarped.xSize(); x++,i++) {

-      float xf = x; float yf = y;

-      float ax = H.data()[0]*xf+H.data()[1]*yf+H.data()[2];

-      float ay = H.data()[3]*xf+H.data()[4]*yf+H.data()[5];

-      float az = H.data()[6]*xf+H.data()[7]*yf+H.data()[8];

-      float invaz = 1.0/az;

-      ax *= invaz; ay *= invaz;

-      int x1 = (int)ax; int y1 = (int)ay;

-      float alphaX = ax-x1; float alphaY = ay-y1;

-      float betaX = 1.0-alphaX; float betaY = 1.0-alphaY;

-      if (x1 < 0 || y1 < 0 || x1+1 >= mXSize || y1+1 >= mYSize) aOutside.data()[i] = true;

-      else {

-        int j = y1*mXSize+x1;

-        for (int k = 0; k < mZSize; k++) {

-          float a = betaX*mData[j]       +alphaX*mData[j+1];

-          float b = betaX*mData[j+mXSize]+alphaX*mData[j+1+mXSize];

-          aWarped.data()[i+k*aWarpedSize] = betaY*a+alphaY*b;

-          j += aSize;

-        }

-      }

-    }

-}

-

-// -----------------------------------------------------------------------------

-// File I/O

-// -----------------------------------------------------------------------------

-

-// readFromMathematicaFile

-template <class T>

-void CTensor<T>::readFromMathematicaFile(const char* aFilename) {

-  using namespace std;

-  // Read the whole file and store data in aData

-  // Ignore blanks, tabs and lines

-  // Also ignore Mathematica comments (* ... *)

-  ifstream aStream(aFilename);

-  string aData;

-  char aChar;

-  bool aBracketFound = false;

-  bool aStarFound = false;

-  bool aCommentFound = false;

-  while (aStream.get(aChar))

-    if (aChar != ' ' && aChar != '\t' && aChar != '\n') {

-      if (aCommentFound) {

-        if (!aStarFound && aChar == '*') aStarFound = true;

-        else {

-          if (aStarFound && aChar == ')') aCommentFound = false;

-          aStarFound = false;

-        }

-      }

-      else {

-        if (!aBracketFound && aChar == '(') aBracketFound = true;

-        else {

-          if (aBracketFound && aChar == '*') aCommentFound = true;

-          else aData += aChar;

-          aBracketFound = false;

-        }

-      }

-    }

-  // Count the number of braces and double braces to figure out z- and y-Size of tensor

-  int aDoubleBraceCount = 0;

-  int aBraceCount = 0;

-  int aPos = 0;

-  while ((aPos = aData.find_first_of('{',aPos)+1) > 0) {

-    aBraceCount++;

-    if (aData[aPos] == '{' && aData[aPos+1] != '{') aDoubleBraceCount++;

-  }

-  // Count the number of commas in the first section to figure out xSize of tensor

-  int aCommaCount = 0;

-  aPos = 0;

-  while (aData[aPos] != '}') {

-    if (aData[aPos] == ',') aCommaCount++;

-    aPos++;

-  }

-  // Adapt size of tensor

-  if (mData != 0) delete[] mData;

-  mXSize = aCommaCount+1;

-  mYSize = (aBraceCount-1-aDoubleBraceCount) / aDoubleBraceCount;

-  mZSize = aDoubleBraceCount;

-  mData = new T[mXSize*mYSize*mZSize];

-  // Analyse file ---------------

-  aPos = 0;

-  if (aData[aPos] != '{') throw EInvalidFileFormat("Mathematica");

-  aPos++;

-  for (int z = 0; z < mZSize; z++) {

-    if (aData[aPos] != '{') throw EInvalidFileFormat("Mathematica");

-    aPos++;

-    for (int y = 0; y < mYSize; y++) {

-      if (aData[aPos] != '{') throw EInvalidFileFormat("Mathematica");

-      aPos++;

-      for (int x = 0; x < mXSize; x++) {

-        int oldPos = aPos;

-        if (x+1 < mXSize) aPos = aData.find_first_of(',',aPos);

-        else aPos = aData.find_first_of('}',aPos);

-        #ifdef GNU_COMPILER

-        string s = aData.substr(oldPos,aPos-oldPos);

-        istrstream is(s.c_str());

-        #else

-        string s = aData.substr(oldPos,aPos-oldPos);

-        istringstream is(s);

-        #endif

-        T aItem;

-        is >> aItem;

-        operator()(x,y,z) = aItem;

-        aPos++;

-      }

-      if (y+1 < mYSize) {

-        if (aData[aPos] != ',') throw EInvalidFileFormat("Mathematica");

-        aPos++;

-        while (aData[aPos] != '{')

-          aPos++;

-      }

-    }

-    aPos++;

-    if (z+1 < mZSize) {

-      if (aData[aPos] != ',') throw EInvalidFileFormat("Mathematica");

-      aPos++;

-      while (aData[aPos] != '{')

-        aPos++;

-    }

-  }

-}

-

-// writeToMathematicaFile

-template <class T>

-void CTensor<T>::writeToMathematicaFile(const char* aFilename) {

-  using namespace std;

-  ofstream aStream(aFilename);

-  aStream << '{';

-  for (int z = 0; z < mZSize; z++) {

-    aStream << '{';

-    for (int y = 0; y < mYSize; y++) {

-      aStream << '{';

-      for (int x = 0; x < mXSize; x++) {

-        aStream << operator()(x,y,z);

-        if (x+1 < mXSize) aStream << ',';

-      }

-      aStream << '}';

-      if (y+1 < mYSize) aStream << ",\n";

-    }

-    aStream << '}';

-    if (z+1 < mZSize) aStream << ",\n";

-  }

-  aStream << '}';

-}

-

-// readFromIMFile

-template <class T>

-void CTensor<T>::readFromIMFile(const char* aFilename) {

-  FILE *aStream;

-  aStream = fopen(aFilename,"rb");

-  // Read image data

-  for (int i = 0; i < mXSize*mYSize*mZSize; i++)

-    mData[i] = getc(aStream);

-  fclose(aStream);

-}

-

-// writeToIMFile

-template <class T>

-void CTensor<T>::writeToIMFile(const char *aFilename) {

-  FILE *aStream;

-  aStream = fopen(aFilename,"wb");

-  // write data

-  for (int i = 0; i < mXSize*mYSize*mZSize; i++) {

-    char dummy = (char)mData[i];

-    fwrite(&dummy,1,1,aStream);

-  }

-  fclose(aStream);

-}

-

-// readFromPGM

-template <class T>

-void CTensor<T>::readFromPGM(const char* aFilename) {

-  FILE *aStream;

-  aStream = fopen(aFilename,"rb");

-  if (aStream == 0) std::cerr << "File not found: " << aFilename << std::endl;

-  int dummy;

-  // Find beginning of file (P5)

-  while (getc(aStream) != 'P');

-  if (getc(aStream) != '5') throw EInvalidFileFormat("PGM");

-  do

-    dummy = getc(aStream);

-  while (dummy != '\n' && dummy != ' ');

-  // Remove comments and empty lines

-  dummy = getc(aStream);

-  while (dummy == '#') {

-    while (getc(aStream) != '\n');

-    dummy = getc(aStream);

-  }

-  while (dummy == '\n')

-    dummy = getc(aStream);

-  // Read image size

-  mXSize = dummy-48;

-  while ((dummy = getc(aStream)) >= 48 && dummy < 58)

-    mXSize = 10*mXSize+dummy-48;

-  while ((dummy = getc(aStream)) < 48 || dummy >= 58);

-  mYSize = dummy-48;

-  while ((dummy = getc(aStream)) >= 48 && dummy < 58)

-    mYSize = 10*mYSize+dummy-48;

-  mZSize = 1;

-  while (dummy != '\n' && dummy != ' ')

-    dummy = getc(aStream);

-  while (dummy != '\n' && dummy != ' ')

-    dummy = getc(aStream);

-  // Adjust size of data structure

-  delete[] mData;

-  mData = new T[mXSize*mYSize];

-  // Read image data

-  for (int i = 0; i < mXSize*mYSize; i++)

-    mData[i] = getc(aStream);

-  fclose(aStream);

-}

-

-// writeToPGM

-template <class T>

-void CTensor<T>::writeToPGM(const char* aFilename) {

-  int rows = (int)floor(sqrt(mZSize));

-  int cols = (int)ceil(mZSize*1.0/rows);

-  FILE* outimage = fopen(aFilename, "wb");

-  fprintf(outimage, "P5 \n");

-  fprintf(outimage, "%ld %ld \n255\n", cols*mXSize,rows*mYSize);

-  for (int r = 0; r < rows; r++)

-    for (int y = 0; y < mYSize; y++)

-      for (int c = 0; c < cols; c++)

-        for (int x = 0; x < mXSize; x++) {

-          unsigned char aHelp;

-          if (r*cols+c >= mZSize) aHelp = 0;

-          else aHelp = (unsigned char)operator()(x,y,r*cols+c);

-          fwrite (&aHelp, sizeof(unsigned char), 1, outimage);

-        }

-  fclose(outimage);

-}

-

-// makeColorTensor

-template <class T>

-void CTensor<T>::makeColorTensor() {

-  if (mZSize != 1) return;

-  int aSize = mXSize*mYSize;

-  int a2Size = 2*aSize;

-  T* aNewData = new T[aSize*3];

-  for (int i = 0; i < aSize; i++)

-    aNewData[i] = aNewData[i+aSize] = aNewData[i+a2Size] = mData[i];

-  mZSize = 3;

-  delete[] mData;

-  mData = aNewData;

-}

-

-// readFromPPM

-template <class T>

-void CTensor<T>::readFromPPM(const char* aFilename) {

-  FILE *aStream;

-  aStream = fopen(aFilename,"rb");

-  if (aStream == 0)

-    std::cerr << "File not found: " << aFilename << std::endl;

-  int dummy;

-  // Find beginning of file (P6)

-  while (getc(aStream) != 'P');

-  dummy = getc(aStream);

-  if (dummy == '5') mZSize = 1;

-  else if (dummy == '6') mZSize = 3;

-  else throw EInvalidFileFormat("PPM");

-  do dummy = getc(aStream); while (dummy != '\n' && dummy != ' ');

-  // Remove comments and empty lines

-  dummy = getc(aStream);

-  while (dummy == '#') {

-    while (getc(aStream) != '\n');

-    dummy = getc(aStream);

-  }

-  while (dummy == '\n')

-    dummy = getc(aStream);

-  // Read image size

-  mXSize = dummy-48;

-  while ((dummy = getc(aStream)) >= 48 && dummy < 58)

-    mXSize = 10*mXSize+dummy-48;

-  while ((dummy = getc(aStream)) < 48 || dummy >= 58);

-  mYSize = dummy-48;

-  while ((dummy = getc(aStream)) >= 48 && dummy < 58)

-    mYSize = 10*mYSize+dummy-48;

-  while (dummy != '\n' && dummy != ' ')

-    dummy = getc(aStream);

-  while (dummy < 48 || dummy >= 58) dummy = getc(aStream);

-  while ((dummy = getc(aStream)) >= 48 && dummy < 58);

-  if (dummy != '\n') while (getc(aStream) != '\n');

-  // Adjust size of data structure

-  delete[] mData;

-  mData = new T[mXSize*mYSize*mZSize];

-  // Read image data

-  int aSize = mXSize*mYSize;

-  if (mZSize == 1)

-    for (int i = 0; i < aSize; i++)

-      mData[i] = getc(aStream);

-  else {

-    int aSizeTwice = aSize+aSize;

-    for (int i = 0; i < aSize; i++) {

-      mData[i] = getc(aStream);

-      mData[i+aSize] = getc(aStream);

-      mData[i+aSizeTwice] = getc(aStream);

-    }

-  }

-  fclose(aStream);

-}

-

-// writeToPPM

-template <class T>

-void CTensor<T>::writeToPPM(const char* aFilename) {

-  FILE* outimage = fopen(aFilename, "wb");

-  fprintf(outimage, "P6 \n");

-  fprintf(outimage, "%d %d \n255\n", mXSize,mYSize);

-  for (int y = 0; y < mYSize; y++)

-    for (int x = 0; x < mXSize; x++) {

-      unsigned char aHelp = (unsigned char)operator()(x,y,0);

-      fwrite (&aHelp, sizeof(unsigned char), 1, outimage);

-      aHelp = (unsigned char)operator()(x,y,1);

-      fwrite (&aHelp, sizeof(unsigned char), 1, outimage);

-      aHelp = (unsigned char)operator()(x,y,2);

-      fwrite (&aHelp, sizeof(unsigned char), 1, outimage);

-    }

-  fclose(outimage);

-}

-

-// readFromPDM

-template <class T>

-void CTensor<T>::readFromPDM(const char* aFilename) {

-  std::ifstream aStream(aFilename);

-  std::string s;

-  // Read header

-  aStream >> s;

-  if (s != "P9") throw EInvalidFileFormat("PDM");

-  char aFeatureType;

-  aStream >> aFeatureType;

-  aStream >> s;

-  aStream >> mXSize;

-  aStream >> mYSize;

-  aStream >> mZSize;

-  aStream >> s;

-  // Adjust size of data structure

-  delete[] mData;

-  mData = new T[mXSize*mYSize*mZSize];

-  // Read data

-  for (int i = 0; i < mXSize*mYSize*mZSize; i++)

-    aStream >> mData[i];

-}

-

-// writeToPDM

-template <class T>

-void CTensor<T>::writeToPDM(const char* aFilename, char aFeatureType) {

-  std::ofstream aStream(aFilename);

-  // write header

-  aStream << "P9" << std::endl;

-  aStream << aFeatureType << "SS" << std::endl;

-  aStream << mZSize << ' ' << mYSize << ' ' << mXSize << std::endl;

-  aStream << "F" << std::endl;

-  // write data

-  for (int i = 0; i < mXSize*mYSize*mZSize; i++) {

-    aStream << mData[i];

-    if (i % 8 == 0) aStream << std::endl;

-    else aStream << ' ';

-  }

-}

-

-// operator ()

-template <class T>

-inline T& CTensor<T>::operator()(const int ax, const int ay, const int az) const {

-  #ifdef _DEBUG

-    if (ax >= mXSize || ay >= mYSize || az >= mZSize || ax < 0 || ay < 0 || az < 0)

-      throw ETensorRangeOverflow(ax,ay,az);

-  #endif

-  return mData[mXSize*(mYSize*az+ay)+ax];

-}

-

-template <class T>

-CVector<T> CTensor<T>::operator()(const float ax, const float ay) const {

-  CVector<T> aResult(mZSize);

-  int x1 = (int)ax;

-  int y1 = (int)ay;

-  int x2 = x1+1;

-  int y2 = y1+1;

-  #ifdef _DEBUG

-  if (x2 >= mXSize || y2 >= mYSize || x1 < 0 || y1 < 0) throw ETensorRangeOverflow(ax,ay,0);

-  #endif

-  float alphaX = ax-x1; float alphaXTrans = 1.0-alphaX;

-  float alphaY = ay-y1; float alphaYTrans = 1.0-alphaY;

-  for (int k = 0; k < mZSize; k++) {

-    float a = alphaXTrans*operator()(x1,y1,k)+alphaX*operator()(x2,y1,k);

-    float b = alphaXTrans*operator()(x1,y2,k)+alphaX*operator()(x2,y2,k);

-    aResult(k) = alphaYTrans*a+alphaY*b;

-  }

-  return aResult;

-}

-

-// operator =

-template <class T>

-inline CTensor<T>& CTensor<T>::operator=(const T aValue) {

-  fill(aValue);

-  return *this;

-}

-

-template <class T>

-CTensor<T>& CTensor<T>::operator=(const CTensor<T>& aCopyFrom) {

-  if (this != &aCopyFrom) {

-    delete[] mData;

-    if (aCopyFrom.mData == 0) {

-      mData = 0; mXSize = 0; mYSize = 0; mZSize = 0;

-    }

-    else {

-      mXSize = aCopyFrom.mXSize;

-      mYSize = aCopyFrom.mYSize;

-      mZSize = aCopyFrom.mZSize;

-      int wholeSize = mXSize*mYSize*mZSize;

-      mData = new T[wholeSize];

-      for (register int i = 0; i < wholeSize; i++)

-        mData[i] = aCopyFrom.mData[i];

-    }

-  }

-  return *this;

-}

-

-// operator +=

-template <class T>

-CTensor<T>& CTensor<T>::operator+=(const CTensor<T>& aTensor) {

-  #ifdef _DEBUG

-  if (mXSize != aTensor.mXSize || mYSize != aTensor.mYSize || mZSize != aTensor.mZSize)

-    throw ETensorIncompatibleSize(mXSize,mYSize,mZSize);

-  #endif

-  int wholeSize = size();

-  for (int i = 0; i < wholeSize; i++)

-    mData[i] += aTensor.mData[i];

-  return *this;

-}

-

-// operator +=

-template <class T>

-CTensor<T>& CTensor<T>::operator+=(const T aValue) {

-  int wholeSize = mXSize*mYSize*mZSize;

-  for (int i = 0; i < wholeSize; i++)

-    mData[i] += aValue;

-  return *this;

-}

-

-// operator *=

-template <class T>

-CTensor<T>& CTensor<T>::operator*=(const T aValue) {

-  int wholeSize = mXSize*mYSize*mZSize;

-  for (int i = 0; i < wholeSize; i++)

-    mData[i] *= aValue;

-  return *this;

-}

-

-// min

-template <class T>

-T CTensor<T>::min() const {

-  T aMin = mData[0];

-  int aSize = mXSize*mYSize*mZSize;

-  for (int i = 1; i < aSize; i++)

-    if (mData[i] < aMin) aMin = mData[i];

-  return aMin;

-}

-

-// max

-template <class T>

-T CTensor<T>::max() const {

-  T aMax = mData[0];

-  int aSize = mXSize*mYSize*mZSize;

-  for (int i = 1; i < aSize; i++)

-    if (mData[i] > aMax) aMax = mData[i];

-  return aMax;

-}

-

-// avg

-template <class T>

-T CTensor<T>::avg() const {

-  T aAvg = 0;

-  for (int z = 0; z < mZSize; z++)

-    aAvg += avg(z);

-  return aAvg/mZSize;

-}

-

-template <class T>

-T CTensor<T>::avg(int az) const {

-  T aAvg = 0;

-  int aSize = mXSize*mYSize;

-  int aTemp = (az+1)*aSize;

-  for (int i = az*aSize; i < aTemp; i++) 

-    aAvg += mData[i];

-  return aAvg/aSize;

-}

-

-// xSize

-template <class T>

-inline int CTensor<T>::xSize() const {

-  return mXSize;

-}

-

-// ySize

-template <class T>

-inline int CTensor<T>::ySize() const {

-  return mYSize;

-}

-

-// zSize

-template <class T>

-inline int CTensor<T>::zSize() const {

-  return mZSize;

-}

-

-// size

-template <class T>

-inline int CTensor<T>::size() const {

-  return mXSize*mYSize*mZSize;

-}

-

-// getMatrix

-template <class T>

-CMatrix<T> CTensor<T>::getMatrix(const int az) const {

-  CMatrix<T> aTemp(mXSize,mYSize);

-  int aMatrixSize = mXSize*mYSize;

-  int aOffset = az*aMatrixSize;

-  for (int i = 0; i < aMatrixSize; i++)

-    aTemp.data()[i] = mData[i+aOffset];

-  return aTemp;

-}

-

-// getMatrix

-template <class T>

-void CTensor<T>::getMatrix(CMatrix<T>& aMatrix, const int az) const {

-  if (aMatrix.xSize() != mXSize || aMatrix.ySize() != mYSize)

-    throw ETensorIncompatibleSize(aMatrix.xSize(),aMatrix.ySize(),mXSize,mYSize);

-  int aMatrixSize = mXSize*mYSize;

-  int aOffset = az*aMatrixSize;

-  for (int i = 0; i < aMatrixSize; i++)

-    aMatrix.data()[i] = mData[i+aOffset];

-}

-

-// putMatrix

-template <class T>

-void CTensor<T>::putMatrix(CMatrix<T>& aMatrix, const int az) {

-  if (aMatrix.xSize() != mXSize || aMatrix.ySize() != mYSize)

-    throw ETensorIncompatibleSize(aMatrix.xSize(),aMatrix.ySize(),mXSize,mYSize);

-  int aMatrixSize = mXSize*mYSize;

-  int aOffset = az*aMatrixSize;

-  for (int i = 0; i < aMatrixSize; i++)

-    mData[i+aOffset] = aMatrix.data()[i];

-}

-

-// data()

-template <class T>

-inline T* CTensor<T>::data() const {

-  return mData;

-}

-

-// N O N - M E M B E R  F U N C T I O N S --------------------------------------

-

-// operator <<

-template <class T>

-std::ostream& operator<<(std::ostream& aStream, const CTensor<T>& aTensor) {

-  for (int z = 0; z < aTensor.zSize(); z++) {

-    for (int y = 0; y < aTensor.ySize(); y++) {

-      for (int x = 0; x < aTensor.xSize(); x++)

-        aStream << aTensor(x,y,z) << ' ';

-      aStream << std::endl;

-    }

-    aStream << std::endl;

-  }

-  return aStream;

-}

-

-#endif