Parameters read by syn_alm_cxx: nlmax (integer): maximum order of l nmmax (integer): maximum order of m (must not be larger than nlmax, default=nlmax) infile (string): input file containing the CMB power spectrum outfile (string): output file name for the calculated a_lm rand_seed (integer): random-number seed fwhm_arcmin (real): FWHM (in arcmin) of a Gaussian beam, which is used to smooth the resulting sky a_lm (default=0) polarisation (bool): if false, only the intensity a_lm are generated, if true, T, G and C a_lm are generated if (polarisation) full_ps (bool, default=false): if false, only TT, GG, CC and TG components of the spectrum are used if true, also TC and GC components are used endif double_precision (bool, default=false): if false, the a_lm are created in single precision, otherwise in double precision.
Parameters read by alm2map_cxx: nlmax (integer): maximum order of l nmmax (integer): maximum order of m (must not be larger than nlmax, default=nlmax) infile (string): input file containing the a_lm outfile (string): output file name for the Healpix map(s) nside (integer): nside parameter for the output map(s) polarisation (bool): if false, only the intensity map is generated, if true, maps for I, Q and U are generated fwhm_arcmin (double, default=0): FWHM in arc minutes of a Gaussian beam, which is used to smooth the a_lm pixel_window (bool, default=false): if true, the a_lm are multiplied by the appropriate pixel window function if (pixel_window) healpix_data (string): directory containing the Healpix data files endif double_precision (bool, default=false): if false, a_lm and maps are read/written in single precision, otherwise in double precision.
Anafast_cxx reads a file containing the map(s) and produces a file containing the temperature power spectrum and, if requested, also the polarisation power spectra
,
and
. The
coefficients computed during the execution also can be written to a file if requested.
Anafast_cxx executes an approximate, discrete point-set quadrature on a sphere sampled at the HEALPix pixel centers. Spherical harmonic transforms are computed using recurrence relations for Legendre polynomials on co-latitude () and Fast Fourier transforms on longitude (
).
Anafast permits two execution options which allow a significant improvement of accuracy of the approximate quadrature performed by this facility:
Parameters read by anafast_cxx: nlmax (integer): maximum order of l nmmax (integer): maximum order of m (must not be larger than nlmax, default=nlmax) infile (string): input file containing the Healpix map outfile (string, default=""): output file name for power spectrum; if empty, no spectrum is written outfile_alms (string, default=""): output file name for the a_lm; if empty, no a_lm are written polarisation (bool): if false, only the intensity a_lm are generated, if true, a_lm for T, G and C component are generated weighted (bool, default=false): if true, weighted quadrature is used if (weighted) healpix_data (string): directory containing the Healpix data files endif iter_order (integer, default=0) number of iterations for the analysis (0: standard analysis) double_precision (bool, default=false): if false, maps and a_lm are read/written in single precision, otherwise in double precision. if (polarisation==true && outfile!="") full_powerspectrum (bool, default=false): if true, write a 6-column power spectrum; if false, write a 4-column power spectrum. endif remove_monopole (bool, default=false): if true, subtract the average pixel value from the temperature map, and add the corresponding value to almT(0,0) after the map2alm transform.
Usage: map2tga <parameter file> or: map2tga <input file> <output file> [-sig <int>] [-pal <int>] [-xsz <int>] [-bar] [-log] [-asinh] [-lon <float>] [-lat <float>] [-mul <float>] [-add <float>] [-min <float>] [-max <float>] [-res <float>] [-title <string>] [-flippal] [-gnomonic] [-interpol] [-equalize] [-viewer <viewer>] Parameters read by map2tga: infile (string): input file containing the Healpix map outfile (string): output TGA file sig (integer, default=1): column number of the requested Healpix map pal (integer, default=4): number of the color palette flippal (bool, default=false): whether the palette should be flipped xsz (integer, default=1024): number of image pixels in x direction bar (logical, default=false): whether a color bar should be added to the image log (logical, default=false): whether the logarithm of the map values should be displayed equalize (logical, default=false): whether histogram equalisation should be performed asinh (logical, default=false): whether the hyperbolic arcsine of the map values should be displayed lon (double, default=0): the longitude (in degrees) of the image center lat (double, default=0): the latitude (in degrees) of the image center mul (double, default=1): scale factor applied to the data add (double, default=0): offset added to the data (before multiplication) min (double, optional): if specified, this value is used as minimum of the color scale max (double, optional): if specified, this value is used as maximum of the color scale res (double, default=1): only for gnomonic projection: the size (in arcmin) of an image pixel title (string, optional): if specified, this value is used as the image title viewer (string, optional): if specified, this executable is used to show the resulting image pro (string, default="mollw"): if this is equal to "gno", gnomonic projection is used, else Mollweide interpol (bool, default=false): false: no interpolation true : bilinear interpolation
Usage: rotalm_cxx <infile> <outfile> <itransform> <pol> Transform 1: Equatorial (2000) -> Galactic (2000) 2: Galactic (2000) -> Equatorial (2000) 3: Equatorial (2000) -> Ecliptic (2000) 4: Ecliptic (2000) -> Equatorial (2000) 5: Ecliptic (2000) -> Galactic (2000) 6: Galactic (2000) -> Ecliptic (2000) 7: Equatorial (1950) -> Galactic (1950) 8: Galactic (1950) -> Equatorial (1950) 9: Equatorial (1950) -> Ecliptic (1950) 10: Ecliptic (1950) -> Equatorial (1950) 11: Ecliptic (1950) -> Galactic (1950) 12: Galactic (1950) -> Ecliptic (1950)
Parameters read by mult_alm: infile (string): input file containing the a_lm outfile (string): output file name for the calculated a_lm fwhm_arcmin_in (real, default=0): FWHM (in arcmin) of a Gaussian beam, which will be _removed_ from the input a_lm fwhm_arcmin_out (real, default=0): FWHM (in arcmin) of a Gaussian beam, which will be used to smoothe the output a_lm nside_pixwin_in (integer, default=0): If this is different from 0, the corresponding pixel window function will be _removed_ from the input a_lm nside_pixwin_out (integer, default=0): If this is different from 0, the corresponding pixel window function will be applied to the output a_lm if ((nside_pixwin_in>0)||(nside_pixwin_out>0)) healpix_data (string): directory containing the Healpix data files endif cl_in (string, default=""): if supplied, the power spectrum from this file will be _removed_ from the input a_lm NOTE: currently only supported for unpolarised a_lm cl_out (string, default=""): if supplied, the power pectrum from this file will be applied to the output a_lm NOTE: currently only supported for unpolarised a_lm polarisation (bool): if false, only the intensity a_lm are generated, if true, T, G and C a_lm are generated double_precision (bool, default=false): if false, the a_lm are read/written in single precision, otherwise in double precision.
Parameters read by smoothing_cxx: fwhm_arcmin (double): FWHM (in arcmin) of the Gaussian beam used for smoothing. If the value is negative, the smoothing of a Gaussian beam with an FWHM of -fwhm_arcmin will be _removed_ from the input map. nlmax (integer): maximum order of l infile (string): input file containing the Healpix map outfile (string, default=""): output file for the smoothed Healpix map polarisation (bool): if false, only an intensity map is smoothed if true, an IQU map is smoothed weighted (bool, default=false): if true, weighted quadrature is used if (weighted) healpix_data (string): directory containing the Healpix data files endif iter_order (integer, default=0) number of iterations for the analysis (0: standard analysis) double_precision (bool, default=false): if false, a_lm are read/written in single precision, otherwise in double precision.
Usage: calc_powspec <almfile1> [<almfile2>] <powspec_file>
Usage: median_filter_cxx <input map> <output map> <radius in arcmin>