---

Grid1D.h

Go to the documentation of this file.

#ifndef Grid1D_H
#define Grid1D_H

template <class T>
class Grid1D {
public:

  typedef T Scalar;

  Grid1D() {}

  Grid1D( Scalar lower, Scalar upper, int nodes) : 
    lower_(lower), upper_(upper), nodes_(nodes) {
    step_ = (upper - lower) / (nodes-1);
  }


  Scalar step() const {return step_;}
  Scalar lower() const {return lower_;}
  Scalar upper() const {return upper_;}
  int nodes() const {return nodes_;}
  int cells() const {return nodes()-1;}

  Scalar node( int i) const { return i*step() + lower();}

  Scalar closestNode( Scalar a) const {
    Scalar b = (a-lower())/step();
    Scalar c = floor(b);
    Scalar tmp = (b-c < 0.5) ? std::max(c,0.) : std::min(c+1.,static_cast<Scalar>(nodes()-1));
    return tmp*step()+lower();
  }

  int index( Scalar a) const {
    int ind = static_cast<int>((a-lower())/step());
    // FIXME: this causes an exception to be thrown later. Should be tested
    // more carefully before release
    //  if (ind < -1 || ind > cells()) {
    //     std::cout << "**** ind = " << ind << " cells: " << cells() << std::endl;
    //    return -1;
    //  }
    return std::max(0, std::min( cells()-1, ind));
  }

private:

  Scalar step_;
  Scalar lower_;
  Scalar upper_;
  int    nodes_;

};

#endif

---