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

changed memory layout for phase array for alm2map

Browse source
This commit is contained in:
Martin Reinecke 2013-02-08 18:25:22 +01:00
parent 4b3da2e693
commit 7c8039d34e
1 changed files with 12 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

15
libsharp/sharp.c
View file

@ -437,9 +437,18 @@ static void init_output (sharp_job *job)
static void alloc_phase (sharp_job *job, int nm, int ntheta) static void alloc_phase (sharp_job *job, int nm, int ntheta)
{ {
if ((nm&1023)==0) nm+=3; // hack to avoid critical strides if (job->type==SHARP_MAP2ALM)
job->s_m=2*job->ntrans*job->nmaps; {
job->s_th=job->s_m*nm; if ((nm&1023)==0) nm+=3; // hack to avoid critical strides
job->s_m=2*job->ntrans*job->nmaps;
job->s_th=job->s_m*nm;
}
else
{
if ((ntheta&1023)==0) ntheta+=3; // hack to avoid critical strides
job->s_th=2*job->ntrans*job->nmaps;
job->s_m=job->s_th*ntheta;
}
job->phase=RALLOC(dcmplx,2*job->ntrans*job->nmaps*nm*ntheta); job->phase=RALLOC(dcmplx,2*job->ntrans*job->nmaps*nm*ntheta);
} }