uniformlyRefinedMesh2.hh

Go to the documentation of this file.

 
#ifndef uniformlyRefinedMesh2_hh
#define uniformlyRefinedMesh2_hh
 
#include "geometry/mesh.hh"
#include <vector>
 
namespace concepts {
 
  class UniformlyRefinedMesh2 : public Mesh2 {
  public:
    UniformlyRefinedMesh2(Mesh2& coarseMesh, uint level);
    virtual ~UniformlyRefinedMesh2() {}
    virtual uint ncell() const { return cell_.size(); }
    virtual Scan2* scan() { return new S(cell_); }
  private:
    class S : public Scan<Cell2> {
    public:
      inline S(std::vector<Cell2*>& cell) :
        idx_(cell.begin()), cell_(cell) {}
      inline S(const S& scan) : idx_(scan.idx_), cell_(scan.cell_) {}
      inline bool eos() const { return idx_ == cell_.end(); }
      inline Cell2& operator++(int) { return *(*idx_++); }
      inline Scan2* clone() const { return new S(*this); }
    private:
      std::vector<Cell2*>::iterator idx_;
      std::vector<Cell2*>& cell_;
    };
    void addCell(Cell2* cell, uint level);
    std::vector<Cell2*> cell_;
  };
 
} // namespace concepts
 
#endif //uniformlyRefinedMesh2_hh