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

add test for Fejer2 quadrature

Browse source
This commit is contained in:
Martin Reinecke 2012-12-03 11:11:18 +01:00
parent 0d8d82677d
commit 7cf0c8c282
1 changed files with 14 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

16
libsharp/sharp_test.c
View file

@ -193,7 +193,8 @@ int main(int argc, char **argv)
MPI_Init(NULL,NULL); MPI_Init(NULL,NULL);
#endif #endif
sharp_module_startup("sharp_test",argc,7, sharp_module_startup("sharp_test",argc,7,
"<healpix|ecp|gauss|hw> <lmax> <nside|nphi> <niter> <spin> <ntrans>",1); "<healpix|ecp|gauss|hw|fejer2> <lmax> <nside|nphi> <niter> <spin> <ntrans>",
1);
int lmax=atoi(argv[2]); int lmax=atoi(argv[2]);
int niter=atoi(argv[4]); int niter=atoi(argv[4]);
@ -217,7 +218,7 @@ int main(int argc, char **argv)
} }
else if (strcmp(argv[1],"ecp")==0) else if (strcmp(argv[1],"ecp")==0)
{ {
int nrings=2*lmax+2; int nrings=2*lmax+1;
int ppring=atoi(argv[3]); int ppring=atoi(argv[3]);
ptrdiff_t npix=(ptrdiff_t)nrings*ppring; ptrdiff_t npix=(ptrdiff_t)nrings*ppring;
printf("\nTesting ECP grid (%d rings, %d pixels/ring, %ld pixels)\n", printf("\nTesting ECP grid (%d rings, %d pixels/ring, %ld pixels)\n",
@ -248,6 +249,17 @@ int main(int argc, char **argv)
check_accuracy(tinfo,lmax,lmax,npix,spin,ntrans,niter); check_accuracy(tinfo,lmax,lmax,npix,spin,ntrans,niter);
sharp_destroy_geom_info(tinfo); sharp_destroy_geom_info(tinfo);
} }
else if (strcmp(argv[1],"fejer2")==0)
{
int nrings=2*lmax+1;
int ppring=atoi(argv[3]);
ptrdiff_t npix=(ptrdiff_t)nrings*ppring;
printf("\nTesting Fejer2 grid (%d rings, %d pixels/ring, %ld pixels)\n",
nrings,ppring,(long)npix);
sharp_make_fejer2_geom_info (nrings, ppring, 0., 1, ppring, &tinfo);
check_accuracy(tinfo,lmax,lmax,npix,spin,ntrans,niter);
sharp_destroy_geom_info(tinfo);
}
else else
UTIL_FAIL("unknown grid geometry"); UTIL_FAIL("unknown grid geometry");