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

fix geom2

Browse source
This commit is contained in:
Martin Reinecke 2012-10-11 08:26:30 +02:00
parent 7852214981
commit fe7e4fa52d
1 changed files with 5 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

10
libsharp/sharp_bench2.c
View file

@ -138,21 +138,21 @@ int main(int argc, char **argv)
int geom2=0; int geom2=0;
if (strcmp(argv[1],"gauss")==0) if (strcmp(argv[1],"gauss")==0)
{ {
int nrings=lmax+1; int nrings=geom2=lmax+1;
int ppring=geom2=atoi(argv[3]); int ppring=atoi(argv[3]);
sharp_make_gauss_geom_info (nrings, ppring, 1, ppring, &tinfo); sharp_make_gauss_geom_info (nrings, ppring, 1, ppring, &tinfo);
} }
else if (strcmp(argv[1],"ecp")==0) else if (strcmp(argv[1],"ecp")==0)
{ {
int nrings=2*lmax+2; int nrings=geom2=2*lmax+2;
int ppring=geom2=atoi(argv[3]); int ppring=atoi(argv[3]);
sharp_make_ecp_geom_info (nrings, ppring, 0., 1, ppring, &tinfo); sharp_make_ecp_geom_info (nrings, ppring, 0., 1, ppring, &tinfo);
} }
else if (strcmp(argv[1],"healpix")==0) else if (strcmp(argv[1],"healpix")==0)
{ {
int nside=atoi(argv[3]); int nside=atoi(argv[3]);
if (nside<1) nside=1; if (nside<1) nside=1;
geom2=nside; geom2=4*nside-1;
sharp_make_healpix_geom_info (nside, 1, &tinfo); sharp_make_healpix_geom_info (nside, 1, &tinfo);
} }
else else