Functions for dealing with geometry information


Functions

void sharp_make_weighted_healpix_geom_info (int nside, int stride, const double *weight, sharp_geom_info **geom_info)
void sharp_make_gauss_geom_info (int nrings, int nphi, double phi0, int stride_lon, int stride_lat, sharp_geom_info **geom_info)
void sharp_make_fejer1_geom_info (int nrings, int nphi, double phi0, int stride_lon, int stride_lat, sharp_geom_info **geom_info)
void sharp_make_cc_geom_info (int nrings, int ppring, double phi0, int stride_lon, int stride_lat, sharp_geom_info **geom_info)
void sharp_make_fejer2_geom_info (int nrings, int ppring, double phi0, int stride_lon, int stride_lat, sharp_geom_info **geom_info)
void sharp_make_mw_geom_info (int nrings, int ppring, double phi0, int stride_lon, int stride_lat, sharp_geom_info **geom_info)
void sharp_make_geom_info (int nrings, const int *nph, const ptrdiff_t *ofs, const int *stride, const double *phi0, const double *theta, const double *wgt, sharp_geom_info **geom_info)
void sharp_destroy_geom_info (sharp_geom_info *info)

Detailed Description


Function Documentation

void sharp_make_weighted_healpix_geom_info ( int  nside,
int  stride,
const double *  weight,
sharp_geom_info **  geom_info 
)

Creates a geometry information describing a HEALPix map with an Nside parameter nside. weight contains the relative ring weights and must have 2*nside entries.

Note:
if weight is a null pointer, all weights are assumed to be 1.

Definition at line 38 of file sharp_geomhelpers.c.

void sharp_make_gauss_geom_info ( int  nrings,
int  nphi,
double  phi0,
int  stride_lon,
int  stride_lat,
sharp_geom_info **  geom_info 
)

Creates a geometry information describing a Gaussian map with nrings iso-latitude rings and nphi pixels per ring. The azimuth of the first pixel in each ring is phi0 (in radians). The index difference between two adjacent pixels in an iso-latitude ring is stride_lon, the index difference between the two start pixels in consecutive iso-latitude rings is stride_lat.

Definition at line 158 of file sharp_geomhelpers.c.

void sharp_make_fejer1_geom_info ( int  nrings,
int  nphi,
double  phi0,
int  stride_lon,
int  stride_lat,
sharp_geom_info **  geom_info 
)

Creates a geometry information describing an ECP map with nrings iso-latitude rings and nphi pixels per ring. The azimuth of the first pixel in each ring is phi0 (in radians). The index difference between two adjacent pixels in an iso-latitude ring is stride_lon, the index difference between the two start pixels in consecutive iso-latitude rings is stride_lat.

Note:
The spacing of pixel centers is equidistant in colatitude and longitude.

The sphere is pixelized in a way that the colatitude of the first ring is 0.5*(pi/nrings) and the colatitude of the last ring is pi-0.5*(pi/nrings). There are no pixel centers at the poles.

This grid corresponds to Fejer's first rule.

Definition at line 193 of file sharp_geomhelpers.c.

void sharp_make_cc_geom_info ( int  nrings,
int  ppring,
double  phi0,
int  stride_lon,
int  stride_lat,
sharp_geom_info **  geom_info 
)

Creates a geometry information describing an ECP map with nrings iso-latitude rings and nphi pixels per ring. The azimuth of the first pixel in each ring is phi0 (in radians). The index difference between two adjacent pixels in an iso-latitude ring is stride_lon, the index difference between the two start pixels in consecutive iso-latitude rings is stride_lat.

Note:
The spacing of pixel centers is equidistant in colatitude and longitude.

The sphere is pixelized in a way that the colatitude of the first ring is 0 and that of the last ring is pi.

This grid corresponds to Clenshaw-Curtis integration.

Definition at line 240 of file sharp_geomhelpers.c.

void sharp_make_fejer2_geom_info ( int  nrings,
int  ppring,
double  phi0,
int  stride_lon,
int  stride_lat,
sharp_geom_info **  geom_info 
)

Creates a geometry information describing an ECP map with nrings iso-latitude rings and nphi pixels per ring. The azimuth of the first pixel in each ring is phi0 (in radians). The index difference between two adjacent pixels in an iso-latitude ring is stride_lon, the index difference between the two start pixels in consecutive iso-latitude rings is stride_lat.

Note:
The spacing of pixel centers is equidistant in colatitude and longitude.

The sphere is pixelized in a way that the colatitude of the first ring is pi/(nrings+1) and that of the last ring is pi-pi/(nrings+1).

This grid corresponds to Fejer's second rule.

Definition at line 289 of file sharp_geomhelpers.c.

void sharp_make_mw_geom_info ( int  nrings,
int  ppring,
double  phi0,
int  stride_lon,
int  stride_lat,
sharp_geom_info **  geom_info 
)

Creates a geometry information describing a map with nrings iso-latitude rings and nphi pixels per ring. The azimuth of the first pixel in each ring is phi0 (in radians). The index difference between two adjacent pixels in an iso-latitude ring is stride_lon, the index difference between the two start pixels in consecutive iso-latitude rings is stride_lat.

Note:
The spacing of pixel centers is equidistant in colatitude and longitude.

The sphere is pixelized in a way that the colatitude of the first ring is pi/(2*nrings-1) and that of the last ring is pi.

This is the grid introduced by McEwen & Wiaux 2011.

This function does not define any quadrature weights.

Definition at line 336 of file sharp_geomhelpers.c.

void sharp_make_geom_info ( int  nrings,
const int *  nph,
const ptrdiff_t *  ofs,
const int *  stride,
const double *  phi0,
const double *  theta,
const double *  wgt,
sharp_geom_info **  geom_info 
)

Creates a geometry information from a set of ring descriptions. All arrays passed to this function must have nrings elements.

Parameters:
nrings the number of rings in the map
nph the number of pixels in each ring
ofs the index of the first pixel in each ring in the map array
stride the stride between consecutive pixels
phi0 the azimuth (in radians) of the first pixel in each ring
theta the colatitude (in radians) of each ring
wgt the pixel weight to be used for the ring in map2alm and adjoint map2alm transforms. Pass NULL to use 1.0 as weight for all rings.
geom_info will hold a pointer to the newly created data structure

Definition at line 179 of file sharp.c.

void sharp_destroy_geom_info ( sharp_geom_info *  info  ) 

Deallocates the geometry information in info.

Definition at line 229 of file sharp.c.


Generated on Thu Oct 8 14:48:51 2015 for LevelS SHT library