RegularGrid.h

Go to the documentation of this file.

#ifndef _RegularGrid_
#define _RegularGrid_
 
#include <ostream>
#include <vector>
#include <vapor/common.h>
#include <vapor/StructuredGrid.h>
 
namespace VAPoR {
 
//
 
class VDF_API RegularGrid : public StructuredGrid {
public:
    //
    RegularGrid(const DimsType &dims, const DimsType &bs, const std::vector<float *> &blks, const CoordType &minu, const CoordType &maxu);
    RegularGrid(const std::vector<size_t> &dims, const std::vector<size_t> &bs, const std::vector<float *> &blks, const std::vector<double> &minu, const std::vector<double> &maxu);
 
    RegularGrid() = default;
    virtual ~RegularGrid() = default;
 
    virtual size_t GetGeometryDim() const override { return (_geometryDim); }
 
    virtual DimsType GetCoordDimensions(size_t dim) const override;
 
    static std::string GetClassType() { return ("Regular"); }
    std::string        GetType() const override { return (GetClassType()); }
 
    //
    virtual void GetBoundingBox(const DimsType &min, const DimsType &max, CoordType &minu, CoordType &maxu) const override;
 
    //
    virtual void GetUserCoordinates(const DimsType &indices, CoordType &coords) const override;
 
    // For grandparent inheritance of
    // Grid::GetUserCoordinates(const size_t indices[], double coords[])
    //
    using Grid::GetUserCoordinates;
 
    virtual bool GetIndicesCell(const CoordType &coords, DimsType &indices) const override;
 
    //
    virtual bool InsideGrid(const CoordType &coords) const override;
 
    class ConstCoordItrRG : public Grid::ConstCoordItrAbstract {
    public:
        ConstCoordItrRG(const RegularGrid *rg, bool begin);
        ConstCoordItrRG(const ConstCoordItrRG &rhs);
 
        ConstCoordItrRG();
        virtual ~ConstCoordItrRG() {}
 
        virtual void            next();
        virtual void            next(const long &offset);
        virtual ConstCoordType &deref() const { return (_coords); }
        virtual const void *    address() const { return this; };
 
        virtual bool equal(const void *rhs) const
        {
            const ConstCoordItrRG *itrptr = static_cast<const ConstCoordItrRG *>(rhs);
 
            return (_index == itrptr->_index);
        }
 
        virtual std::unique_ptr<ConstCoordItrAbstract> clone() const { return std::unique_ptr<ConstCoordItrAbstract>(new ConstCoordItrRG(*this)); };
 
    private:
        DimsType  _index;
        DimsType  _dims;
        CoordType _minu;
        CoordType _delta;
        CoordType _coords;
    };
 
    virtual ConstCoordItr ConstCoordBegin() const override { return ConstCoordItr(std::unique_ptr<ConstCoordItrAbstract>(new ConstCoordItrRG(this, true))); }
    virtual ConstCoordItr ConstCoordEnd() const override { return ConstCoordItr(std::unique_ptr<ConstCoordItrAbstract>(new ConstCoordItrRG(this, false))); }
 
    VDF_API friend std::ostream &operator<<(std::ostream &o, const RegularGrid &rg);
 
protected:
    virtual float GetValueNearestNeighbor(const CoordType &coords) const override;
 
    virtual float GetValueLinear(const CoordType &coords) const override;
 
    //
    virtual void GetUserExtentsHelper(CoordType &minu, CoordType &maxu) const override;
 
private:
    void _regularGrid(const CoordType &minu, const CoordType &maxu);
 
    CoordType _minu = {{0.0, 0.0, 0.0}};
    CoordType _maxu = {{0.0, 0.0, 0.0}};
    size_t    _geometryDim;
    CoordType _delta;    // increment between grid points in user coords
};
};    // namespace VAPoR
#endif