+ dXFrac * dYFrac * m_ppMatrix[iFloorX+1][iFloorY+1]);
+
+ return result;
+ }
+};
+
+
+template<class T>
+class BicubicPolyInterpolator {
+private:
+ T** const m_ppMatrix;
+ const unsigned int m_nx;
+ const unsigned int m_ny;
+
+public:
+ BicubicPolyInterpolator (T** ppMatrix, unsigned int nx, unsigned int ny)
+ : m_ppMatrix(ppMatrix), m_nx(nx), m_ny(ny)
+ {}
+
+ T interpolate (double dXPos, double dYPos)
+ {
+ // int iFloorX = static_cast<int>(floor (dXPos));
+ // int iFloorY = static_cast<int>(floor (dYPos));
+ // double dXFrac = dXPos - iFloorX;
+ // double dYFrac = dYPos - iFloorY;
+
+ T result = 0;
+
+ // Need to add code