SphericalGrid.h

Go to the documentation of this file.

#ifndef _SphericalGrid_
#define _SphericalGrid_
#include <vector>
#include <vapor/common.h>
#include "RegularGrid.h"
#ifdef _WINDOWS
    #pragma warning(disable : 4251 4100)
#endif
 
namespace VAPoR {
 
//
class VDF_API SphericalGrid : public RegularGrid {
public:
    SphericalGrid(const size_t bs[3], const size_t min[3], const size_t max[3], const double extents[6], const size_t permutation[3], const bool periodic[3], float **blks);
 
    SphericalGrid(const size_t bs[3], const size_t min[3], const size_t max[3], const double extents[6], const size_t permutation[3], const bool periodic[3], float **blks, 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] = _extentsC[i];
    };
 
    virtual void GetBoundingBox(const size_t min[3], const size_t max[3], double extents[6]) 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;
 
    bool InsideGrid(double x, double y, double z) const;
 
    // phi is in range -180 to 180, theta is in range -180/2 to 180/2
    //
    static inline void CartToSph(double x, double y, double z, double *phi, double *theta, double *r);
 
    static inline void SphToCart(double phi, double theta, double r, double *x, double *y, double *z);
 
private:
    double            _extentsC[6];    // extents in Cartesian coordinates, ordered x,y,z
    std::vector<long> _permutation;    // permutation vector for coordinate ordering
 
    void _GetUserExtents(double extents[6]) const;
 
    void _permute(const std::vector<long> &permutation, double result[3], double x, double y, double z) const;
};
};    // namespace VAPoR
#endif