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

fix

Browse source
This commit is contained in:
Martin Reinecke 2018-12-19 15:07:17 +01:00
parent 57da18d154
commit d295e4c609
1 changed files with 3 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

9
libsharp/sharp_core.c
View file

@ -261,10 +261,8 @@ NOINLINE static void calc_alm2map (sharp_job * restrict job,
d->p1i[i] += d->lam2[i]*d->corfac[i]*ai1; d->p1i[i] += d->lam2[i]*d->corfac[i]*ai1;
d->lam2[i] = f20*d->cth[i]*d->lam1[i] - f21*d->lam2[i]; d->lam2[i] = f20*d->cth[i]*d->lam1[i] - f21*d->lam2[i];
if (rescale(&d->lam1[i], &d->lam2[i], &d->scale[i], vload(sharp_ftol))) if (rescale(&d->lam1[i], &d->lam2[i], &d->scale[i], vload(sharp_ftol)))
{
getCorfac(d->scale[i], &d->corfac[i], gen->cf); getCorfac(d->scale[i], &d->corfac[i], gen->cf);
full_ieee &= vallTrue(vge(d->scale[i],vload(sharp_minscale))); full_ieee &= vallTrue(vge(d->scale[i],vload(sharp_minscale)));
}
d->p2r[i] += d->lam1[i]*d->corfac[i]*ar2; d->p2r[i] += d->lam1[i]*d->corfac[i]*ar2;
d->p2i[i] += d->lam1[i]*d->corfac[i]*ai2; d->p2i[i] += d->lam1[i]*d->corfac[i]*ai2;
} }
@ -336,16 +334,15 @@ NOINLINE static void calc_map2alm(sharp_job * restrict job,
atmp[1] += d->lam2[i]*d->corfac[i]*d->p1i[i]; atmp[1] += d->lam2[i]*d->corfac[i]*d->p1i[i];
d->lam2[i] = f20*d->cth[i]*d->lam1[i] - f21*d->lam2[i]; d->lam2[i] = f20*d->cth[i]*d->lam1[i] - f21*d->lam2[i];
if (rescale(&d->lam1[i], &d->lam2[i], &d->scale[i], vload(sharp_ftol))) if (rescale(&d->lam1[i], &d->lam2[i], &d->scale[i], vload(sharp_ftol)))
{
getCorfac(d->scale[i], &d->corfac[i], gen->cf); getCorfac(d->scale[i], &d->corfac[i], gen->cf);
full_ieee &= vallTrue(vge(d->scale[i],vload(sharp_minscale))); full_ieee &= vallTrue(vge(d->scale[i],vload(sharp_minscale)));
}
atmp[2] += d->lam1[i]*d->corfac[i]*d->p2r[i]; atmp[2] += d->lam1[i]*d->corfac[i]*d->p2r[i];
atmp[3] += d->lam1[i]*d->corfac[i]*d->p2i[i]; atmp[3] += d->lam1[i]*d->corfac[i]*d->p2i[i];
} }
vhsum_cmplx_special (atmp[0], atmp[1], atmp[2], atmp[3], &alm[l]); vhsum_cmplx_special (atmp[0], atmp[1], atmp[2], atmp[3], &alm[l]);
l+=2; l+=2;
} }
if (l>lmax) return;
for (int i=0; i<nv2; ++i) for (int i=0; i<nv2; ++i)
{ {