CartesianGrid< nDim >::azimuthalBbox Struct Reference

#include <cartesiangrid.h>

Public Member Functions
void	init (MFloat angle, MInt periodicDir1, MInt periodicDir2, MInt noLayers)
MInt	azimuthalSide (const MFloat phi)
MFloat	azimuthalCenter ()
MFloat	azimuthalBoundary (const MFloat *coords, MInt side)
MInt	index (const MFloat *coord, MInt dir)

Public Attributes
MFloat	azimuthalAngle = F0
MFloat	minCoord [nDim]
MFloat	boxLength [nDim]
MInt	nBins [nDim]
MFloat *	minPer3 = nullptr
MFloat *	maxPer3 = nullptr
MFloat	center = F0
MInt	axDir = -1
MInt	perDir1 = -1
MInt	perDir2 = -1
MInt	noHaloLayers = -1

Detailed Description

template<MInt nDim>
struct CartesianGrid< nDim >::azimuthalBbox

Definition at line 276 of file cartesiangrid.h.

Member Function Documentation

azimuthalBoundary()

template<MInt nDim>

MFloat CartesianGrid< nDim >::azimuthalBbox::azimuthalBoundary ( const MFloat *  coords, MInt  side  )

inline

Definition at line 310 of file cartesiangrid.h.

                                                              {
      MInt perSize2 = nBins[perDir2];
      MInt ind1 = index(coords, perDir1);
      MInt ind2 = index(coords, perDir2);
      if(side < 0) {
        MFloat minPhi = center - F1B2 * azimuthalAngle; // std::numeric_limits<MFloat>::max();
        for(MInt i = -(noHaloLayers + 1); i < (noHaloLayers + 2); i++) {
          MInt ind1_ = ind1 + i;
          ind1_ = mMax(0, mMin(ind1_, nBins[perDir1] - 1));
          MInt ind2_ = ind2 + i;
          ind2_ = mMax(0, mMin(ind2_, nBins[perDir2] - 1));
          minPhi = mMin(minPhi, minPer3[ind1 * perSize2 + ind2_]);
          minPhi = mMin(minPhi, minPer3[ind1_ * perSize2 + ind2]);
        }
        return minPhi;
      } else if(side > 0) {
        MFloat maxPhi = center + F1B2 * azimuthalAngle; //-std::numeric_limits<MFloat>::max();
        for(MInt i = -(noHaloLayers + 1); i < (noHaloLayers + 2); i++) {
          MInt ind1_ = ind1 + i;
          ind1_ = mMax(0, mMin(ind1_, nBins[perDir1] - 1));
          MInt ind2_ = ind2 + i;
          ind2_ = mMax(0, mMin(ind2_, nBins[perDir2] - 1));
          maxPhi = mMax(maxPhi, maxPer3[ind1 * perSize2 + ind2_]);
          maxPhi = mMax(maxPhi, maxPer3[ind1_ * perSize2 + ind2]);
        }
        return maxPhi;
      } else {
        mTerm(1, AT_, "side = 0?");
      }
    }

azimuthalCenter()

template<MInt nDim>

MFloat CartesianGrid< nDim >::azimuthalBbox::azimuthalCenter ( )

inline

Definition at line 309 of file cartesiangrid.h.

{ return center; }

azimuthalSide()

template<MInt nDim>

MInt CartesianGrid< nDim >::azimuthalBbox::azimuthalSide ( const MFloat  phi )

inline

Definition at line 301 of file cartesiangrid.h.

                                         {
      MFloat mid = center;
      if(mid > phi) {
        return -1;
      } else {
        return 1;
      }
    }

index()

template<MInt nDim>

MInt CartesianGrid< nDim >::azimuthalBbox::index ( const MFloat *  coord, MInt  dir  )

inline

Definition at line 340 of file cartesiangrid.h.

                                              {
      MInt ind = -1;
      ind = (MInt)((coord[dir] - minCoord[dir]) * nBins[dir] / boxLength[dir]);
      if(ind < 0)
        ind = 0;
      else if(ind >= nBins[dir])
        ind = nBins[dir] - 1;
 
      return ind;
    }

init()

template<MInt nDim>

void CartesianGrid< nDim >::azimuthalBbox::init ( MFloat  angle, MInt  periodicDir1, MInt  periodicDir2, MInt  noLayers  )

inline

Definition at line 289 of file cartesiangrid.h.

                                                                                 {
      azimuthalAngle = angle;
      perDir1 = periodicDir1;
      perDir2 = periodicDir2;
      noHaloLayers = noLayers;
      MInt perSize1 = nBins[perDir1];
      MInt perSize2 = nBins[perDir2];
      mDeallocate(minPer3);
      mDeallocate(maxPer3);
      mAlloc(minPer3, perSize1 * perSize2, "minPer3", std::numeric_limits<MFloat>::max(), AT_);
      mAlloc(maxPer3, perSize1 * perSize2, "maxPer3", -std::numeric_limits<MFloat>::max(), AT_);
    }

Member Data Documentation

axDir

template<MInt nDim>

MInt CartesianGrid< nDim >::azimuthalBbox::axDir = -1

Definition at line 284 of file cartesiangrid.h.

azimuthalAngle

template<MInt nDim>

MFloat CartesianGrid< nDim >::azimuthalBbox::azimuthalAngle = F0

Definition at line 277 of file cartesiangrid.h.

boxLength

template<MInt nDim>

MFloat CartesianGrid< nDim >::azimuthalBbox::boxLength[nDim]

Definition at line 279 of file cartesiangrid.h.

center

template<MInt nDim>

MFloat CartesianGrid< nDim >::azimuthalBbox::center = F0

Definition at line 283 of file cartesiangrid.h.

maxPer3

template<MInt nDim>

MFloat* CartesianGrid< nDim >::azimuthalBbox::maxPer3 = nullptr

Definition at line 282 of file cartesiangrid.h.

minCoord

template<MInt nDim>

MFloat CartesianGrid< nDim >::azimuthalBbox::minCoord[nDim]

Definition at line 278 of file cartesiangrid.h.

minPer3

template<MInt nDim>

MFloat* CartesianGrid< nDim >::azimuthalBbox::minPer3 = nullptr

Definition at line 281 of file cartesiangrid.h.

nBins

template<MInt nDim>

MInt CartesianGrid< nDim >::azimuthalBbox::nBins[nDim]

Definition at line 280 of file cartesiangrid.h.

noHaloLayers

template<MInt nDim>

MInt CartesianGrid< nDim >::azimuthalBbox::noHaloLayers = -1

Definition at line 287 of file cartesiangrid.h.

perDir1

template<MInt nDim>

MInt CartesianGrid< nDim >::azimuthalBbox::perDir1 = -1

Definition at line 285 of file cartesiangrid.h.

perDir2

template<MInt nDim>

MInt CartesianGrid< nDim >::azimuthalBbox::perDir2 = -1

Definition at line 286 of file cartesiangrid.h.

---

The documentation for this struct was generated from the following file:

• /home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/GRID/cartesiangrid.h