Interpolation/src/Grid3D.cc

0001 #include "Grid3D.h"
0002 #include <iostream>
0003
0004 /*
0005 Grid3D::Grid3D( const Grid1D& ga, const Grid1D& gb, const Grid1D& gc,
0006         std::vector<ValueType> const & data) :
0007   grida_(ga), gridb_(gb), gridc_(gc) {
0008   data_.reserve(data.size());
0009   //FIXME use a std algo
0010   for (size_t i=0; i<=data.size(); ++i)
0011     data_.push_back(ValueType(data[i].x(),data[i].y(),data[i].z()));
0012   stride1_ = gridb_.nodes() * gridc_.nodes();
0013   stride2_ = gridc_.nodes();
0014 }
0015 */
0016
0017 void Grid3D::dump() const {
0018   for (int j = 0; j < gridb().nodes(); ++j) {
0019     for (int k = 0; k < gridc().nodes(); ++k) {
0020       for (int i = 0; i < grida().nodes(); ++i) {
0021         std::cout << grida().node(i) << " " << gridb().node(j) << " " << gridc().node(k) << " " << operator()(i, j, k)
0022                   << std::endl;
0023       }
0024     }
0025   }
0026 }