From 14955e74f1010a89e817cfcf6dd26bfc2df9fac7 Mon Sep 17 00:00:00 2001 From: Martin Reinecke Date: Thu, 10 Jan 2019 15:35:58 +0100 Subject: [PATCH] temporary commit for future reference --- libsharp/sharp_core.c | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/libsharp/sharp_core.c b/libsharp/sharp_core.c index 4a6967e..f8bc16d 100644 --- a/libsharp/sharp_core.c +++ b/libsharp/sharp_core.c @@ -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; }