LayeredGrid.h

Go to the documentation of this file.

#ifndef _LayeredGrid_
#define _LayeredGrid_
#include <vapor/common.h>
#include "RegularGrid.h"

//
//

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

 LayeredGrid(
    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,
    float ** coords,
    int varying_dim
);

 LayeredGrid(
    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,
    float ** coords,
    int varying_dim,
    float missing_value
 );

 virtual ~LayeredGrid();


 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]
 );

 bool InsideGrid(double x, double y, double z) const;


 //
 virtual void SetPeriodic(const bool periodic[3]);

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

 float **GetCoordBlks() const { return(_coords); };


private:
 float **_coords;
 int _varying_dim;
 double _extents[6];

 double _GetValueLinear(double x, double y, double z) const;
 void _GetUserExtents(double extents[6]) const;
 void _GetBoundingBox(
    const size_t min[3],
    const size_t max[3],
    double extents[6]
 ) const;
 double _GetVaryingCoord(size_t i, size_t j, size_t k) const;

 //double verticalLinearInterpolation(double x, double y, double z) const;
 //double quadraticInterpolation(double x, double y, double z) const;

 double _interpolateVaryingCoord(
    size_t i0, size_t j0, size_t k0,
    double x, double y, double z 
 ) const;

};
};
#endif