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

explanatory comment

Browse source
This commit is contained in:
Martin Reinecke 2013-01-17 11:45:19 +01:00
parent 2353da0042
commit 3027924ce4
1 changed files with 2 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

2
libsharp/sharp_mpi.c
View file

@ -226,6 +226,8 @@ static void sharp_execute_job_mpi (sharp_job *job, MPI_Comm comm)
for (int isub=0; isub<nsub; ++isub) for (int isub=0; isub<nsub; ++isub)
{ {
sharp_job ljob=*job; sharp_job ljob=*job;
// When creating a_lm, every sub-job produces a complete set of
// coefficients; they need to be added up.
if ((isub>0)&&(job->type==SHARP_MAP2ALM)) ljob.flags|=SHARP_ADD; if ((isub>0)&&(job->type==SHARP_MAP2ALM)) ljob.flags|=SHARP_ADD;
sharp_geom_info lginfo; sharp_geom_info lginfo;
lginfo.pair=RALLOC(sharp_ringpair,(job->ginfo->npairs/nsub)+1); lginfo.pair=RALLOC(sharp_ringpair,(job->ginfo->npairs/nsub)+1);