From fe7e4fa52d35376da4f0326b1c02523fc22c9c09 Mon Sep 17 00:00:00 2001 From: Martin Reinecke Date: Thu, 11 Oct 2012 08:26:30 +0200 Subject: [PATCH] fix geom2 --- libsharp/sharp_bench2.c | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/libsharp/sharp_bench2.c b/libsharp/sharp_bench2.c index d637554..e86c5e2 100644 --- a/libsharp/sharp_bench2.c +++ b/libsharp/sharp_bench2.c @@ -138,21 +138,21 @@ int main(int argc, char **argv) int geom2=0; if (strcmp(argv[1],"gauss")==0) { - int nrings=lmax+1; - int ppring=geom2=atoi(argv[3]); + int nrings=geom2=lmax+1; + int ppring=atoi(argv[3]); sharp_make_gauss_geom_info (nrings, ppring, 1, ppring, &tinfo); } else if (strcmp(argv[1],"ecp")==0) { - int nrings=2*lmax+2; - int ppring=geom2=atoi(argv[3]); + int nrings=geom2=2*lmax+2; + int ppring=atoi(argv[3]); sharp_make_ecp_geom_info (nrings, ppring, 0., 1, ppring, &tinfo); } else if (strcmp(argv[1],"healpix")==0) { int nside=atoi(argv[3]); if (nside<1) nside=1; - geom2=nside; + geom2=4*nside-1; sharp_make_healpix_geom_info (nside, 1, &tinfo); } else