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

temporary commit for future reference

Browse source
This commit is contained in:
Martin Reinecke 2019-01-10 15:35:58 +01:00
parent 721a8e3b93
commit 14955e74f1
1 changed files with 4 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

6
libsharp/sharp_core.c
View file

@ -531,6 +531,8 @@ NOINLINE static void alm2map_spin_kernel(sxdata_v * restrict d,
d->l1p[i] = (d->cth[i]*fx10 - fx11)*d->l2p[i] - d->l1p[i];
d->l1m[i] = (d->cth[i]*fx10 + fx11)*d->l2m[i] - d->l1m[i];
d->l2p[i] = (d->cth[i]*fx20 - fx21)*d->l1p[i] - d->l2p[i];
d->l2m[i] = (d->cth[i]*fx20 + fx21)*d->l1m[i] - d->l2m[i];
d->p1pr[i] += aci2*d->l1p[i];
d->p1pi[i] -= acr2*d->l1p[i];
@ -542,8 +544,8 @@ NOINLINE static void alm2map_spin_kernel(sxdata_v * restrict d,
d->p2mr[i] += acr2*d->l1m[i];
d->p2mi[i] += aci2*d->l1m[i];
d->l2p[i] = (d->cth[i]*fx20 - fx21)*d->l1p[i] - d->l2p[i];
d->l2m[i] = (d->cth[i]*fx20 + fx21)*d->l1m[i] - d->l2m[i];
// d->l2p[i] = (d->cth[i]*fx20 - fx21)*d->l1p[i] - d->l2p[i];
// d->l2m[i] = (d->cth[i]*fx20 + fx21)*d->l1m[i] - d->l2m[i];
}
l+=2;
}