Aquila-Consortium/libsharp2
0
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

sharp_make_geom_info: Allow passing NULL as weights for all-one weights

Browse source
This commit is contained in:
Dag Sverre Seljebotn 2012-11-06 14:28:03 +01:00
parent 93a01f769c
commit 0383da73b4
2 changed files with 4 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

2
libsharp/sharp.c
View file

@ -176,7 +176,7 @@ void sharp_make_geom_info (int nrings, const int *nph, const ptrdiff_t *ofs,
infos[m].theta = theta[m]; infos[m].theta = theta[m];
infos[m].cth = cos(theta[m]); infos[m].cth = cos(theta[m]);
infos[m].sth = sin(theta[m]); infos[m].sth = sin(theta[m]);
infos[m].weight = weight[m]; infos[m].weight = (weight != NULL) ? weight[m] : 1.;
infos[m].phi0 = phi0[m]; infos[m].phi0 = phi0[m];
infos[m].ofs = ofs[m]; infos[m].ofs = ofs[m];
infos[m].stride = stride[m]; infos[m].stride = stride[m];

4
libsharp/sharp_lowlevel.h
View file

@ -123,7 +123,9 @@ void sharp_destroy_alm_info (sharp_alm_info *info);
\param stride the stride between consecutive pixels \param stride the stride between consecutive pixels
\param phi0 the azimuth (in radians) of the first pixel in each ring \param phi0 the azimuth (in radians) of the first pixel in each ring
\param theta the colatitude (in radians) of each ring \param theta the colatitude (in radians) of each ring
\param weight the pixel weight to be used for the ring \param weight the pixel weight to be used for the ring. Pass NULL to use
1.0 as weight for all rings. By default weights are used for alm2map
but not map2alm, but the execution flags can override this.
\param geom_info will hold a pointer to the newly created data structure \param geom_info will hold a pointer to the newly created data structure
*/ */
void sharp_make_geom_info (int nrings, const int *nph, const ptrdiff_t *ofs, void sharp_make_geom_info (int nrings, const int *nph, const ptrdiff_t *ofs,