Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037 #include "prefix2d.hpp"
00038 #include "assert.h"
00039
00040 template<typename T, bool DEBUG>
00041 template<typename ar2d>
00042 util::Prefix2D<T,DEBUG>::Prefix2D(int x, int y, const ar2d& a)
00043 :sizeX(x),sizeY(y)
00044 {
00045 if (DEBUG)
00046 assert (! (x<=0 || y<=0));
00047
00048 data = new T[(sizeX+1)*(sizeY+1)];
00049 for (int i=0; i<=sizeX; i++)
00050 {
00051 T sum = 0;
00052 for (int j=0; j<=sizeY; j++)
00053 {
00054 if (i==0)
00055 data[i*(sizeY+1)+j] = 0;
00056 else if (j==0)
00057 data[i*(sizeY+1)+j] = 0;
00058 else
00059 {
00060 sum += a[i-1][j-1];
00061 data[i*(sizeY+1)+j] = data[(i-1)*(sizeY+1)+j] + sum;
00062 }
00063 }
00064 }
00065 }
00066
00067 template <typename T, bool DEBUG>
00068 int util::Prefix2D<T,DEBUG>::prefixsizeX() const
00069 {
00070 return sizeX + 1;
00071 }
00072
00073 template <typename T, bool DEBUG>
00074 int util::Prefix2D<T,DEBUG>::prefixsizeY() const
00075 {
00076 return sizeY + 1;
00077 }
00078
00079 template <typename T, bool DEBUG>
00080 T* util::Prefix2D<T,DEBUG>::operator[] (int x)
00081 {
00082 if (DEBUG)
00083 assert(! (x<0 || x > sizeX));
00084
00085
00086 return data+(sizeY+1)*x;
00087 }
00088
00089 template <typename T, bool DEBUG>
00090 const T* util::Prefix2D<T,DEBUG>::operator[] (int x) const
00091 {
00092 if (DEBUG)
00093 assert(! (x<0 || x > sizeX));
00094
00095 return data+(sizeY+1)*x;
00096 }
00097
00098 template <typename T, bool DEBUG>
00099 util::Prefix2D<T,DEBUG>::~Prefix2D()
00100 {
00101 delete[] data;
00102 }