StretchedGrid.h

Go to the documentation of this file.

#ifndef _StretchedGrid_
#define _StretchedGrid_

#include <vector>
#include <vapor/common.h>
#include "RegularGrid.h"
#ifdef _WINDOWS
#pragma warning(disable : 4251 4100)
#endif
//
//

namespace VAPoR {
class VDF_API StretchedGrid : public RegularGrid {
public:

 StretchedGrid(
    const size_t bs[3],
    const size_t min[3],
    const size_t max[3],
    const double extents[6],
    const bool periodic[3],
    float ** blks,
    const std::vector <double> &xcoords,
    const std::vector <double> &ycoords,
    const std::vector <double> &zcoords
);

 StretchedGrid(
    const size_t bs[3],
    const size_t min[3],
    const size_t max[3],
    const double extents[6],
    const bool periodic[3],
    float ** blks,
    const std::vector <double> &xcoords,
    const std::vector <double> &ycoords,
    const std::vector <double> &zcoords,
    float missing_value
 );


 float GetValue(double x, double y, double z) const;

 virtual void GetUserExtents(double extents[6]) const {
    for (int i=0; i<6; i++) extents[i] = _extents[i];
 }


 virtual void GetBoundingBox(
    const size_t min[3],
    const size_t max[3],
    double extents[6]
 ) const;

 virtual void    GetEnclosingRegion(
    const double minu[3], const double maxu[3], size_t min[3], size_t max[3]
 ) const;


 int GetUserCoordinates(
    size_t i, size_t j, size_t k, 
    double *x, double *y, double *z
 ) const;

 void GetIJKIndex(
    double x, double y, double z,
    size_t *i, size_t *j, size_t *k
 ) const;

 void GetIJKIndexFloor(
    double x, double y, double z,
    size_t *i, size_t *j, size_t *k
 ) const;

 int Reshape(
    const size_t min[3],
    const size_t max[3],
    const bool periodic[3]
 );

 virtual void GetMinCellExtents(double *x, double *y, double *z) const;

 void GetUserCoordinateMaps(
    std::vector <double> &xcoords,
    std::vector <double> &ycoords,
    std::vector <double> &zcoords
 ) const {
    xcoords = _xcoords; ycoords = _ycoords; zcoords = _zcoords;
 };



private:
 size_t _min[3];
 size_t _max[3];
 double _delta[3];
 double _extents[6];
 std::vector <double> _xcoords;
 std::vector <double> _ycoords;
 std::vector <double> _zcoords;

 float _GetValueNearestNeighbor(double x, double y, double z) const;
 float _GetValueLinear(double x, double y, double z) const;
 float _GetValueQuadratic(double x, double y, double z) const;


};
};
#endif