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

cleanup

Browse source
This commit is contained in:
Martin Reinecke 2018-12-13 10:42:07 +01:00
parent 88a78b2fcb
commit 5057843dac
1 changed files with 9 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

9
libsharp/sharp_core.c
View file

@ -447,7 +447,16 @@ NOINLINE static void inner_loop_m2a(sharp_job *job, const int *ispair,
++ith; ++ith;
} }
if (nth>0) if (nth>0)
{
int i2=((nth+VLEN-1)/VLEN)*VLEN;
for (int i=nth; i<i2; ++i)
{
cth.s[i]=cth.s[nth-1];
sth.s[i]=sth.s[nth-1];
p1.s.r[i]=p1.s.i[i]=p2.s.r[i]=p2.s.i[i]=0.;
}
calc_map2alm(cth.b,sth.b,gen,job,&p1.b,&p2.b, nth); calc_map2alm(cth.b,sth.b,gen,job,&p1.b,&p2.b, nth);
}
} }
} }
else else