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;

 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;



 // 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;


};
};
#endif