Previous: plmgen Up: Healpix/F90 facilities Next: sky_ng_sim Top: Main Page

smoothing

This program can be used to convolve a map with a gaussian beam. The input map can be given in RING or NESTED scheme and the smoothed map is written to a FITS file in the RING scheme.

NOTE: This automated facility is susceptible to problems with non-commutativity of discrete spherical harmonics transforms, described in the Recommendations for Users of the anafast facility. If very high accuracy of the results is required in the spectral regime of , it is recommended to choose an iterative computation of the a_lm coefficients.

FORMAT

% smoothing [options] [parameter_file]

COMMAND LINE OPTIONS

-d

--double

double precision mode (see Notes on double/single precision modes on page )

-s

--single

single precision mode (default)

QUALIFIERS

simul_type =

Defines which map(s) to analyse, 1=temperature only, 2=temperature AND polarisation. (default= 1)

infile =

Defines the filename for the FITS file containing the map to be smoothed. (default= 'map.fits')

nlmax =

Defines the l_max value for the application. (default= 64)

iter_order =

Defines the maximum order of quadrature iteration to be used. (default=0, no iteration)

fwhm_arcmin =

Defines the FWHM in arcminutes of the gaussian beam for the convolution. (default=10)

beam_file =

Defines the FITS file describing the Legendre window function of the circular beam to be used for the simulation. If set to an existing file name, it will override the fhwm_arcmin given above. default=`'

outfile =

Defines the filename for the file that will contain the smoothed map. (default='map_smoothed.fits')

plmfile =

Defines the name for an input file containing precomputed Legendre polynomials P_lm. (default= no entry -- anafast executes the recursive evaluation of P_lms)

w8file =

Defines name for an input file containing ring weights in the improved quadrature mode (default= no entry -- the name is assumed to be 'weight_ring_n0xxxx.fits' where xxxx is nsmax)

w8filedir =

Gives the directory where the weight files are to be found (default= no entry -- smoothing searches in the default directories, see introduction)

won =

Set this to 1 if weight files are to be used, otherwise set it to 0 (or 2). (default= 0)

DESCRIPTION

A FITS file containing a HEALPix map in RING or NESTED scheme is read in. The map is analysed and smoothed in fourier space with a gaussian beam of a given FHWM. A new map is then synthesized using the smoothed a_lm coeffecients. For a more accurate application, an iteration of arbitrary order can be applied. The output map is stored in the same scheme as the input map.

DATASETS

The following datasets are involved in the smoothing processing.

SUPPORT

This section lists those routines and facilities (including those external to the Healpix distribution) which can assist in the utilisation of smoothing.

synfast

This HEALPix facility can generate a map and also do the smoothing.

anafast

This HEALPix facility can analyse a smoothed map.

EXAMPLE # 1:

Smoothing runs in interactive mode, self-explanatory.

EXAMPLE # 2:

When `filename' is present, smoothing enters the non-interactive mode and parses its inputs from the file `filename'. This has the following structure: the first entry is a qualifier which announces to the parser which input immediately follows. If this input is omitted in the input file, the parser assumes the default value. If the equality sign is omitted, then the parser ignores the entry. In this way comments may also be included in the file. In this example, the file contains the following qualifiers:
simul_type = 1
nlmax = 64
infile = map.fits
outfile = map_smoothed.fits
fwhm_arcmin = 10.
iter_order = 1

RELEASE NOTES

Initial release (HEALPix 0.90)

Extension to polarization and arbitrary circular beams (HEALPix 1.20)

MESSAGES

This section describes error messages generated by smoothing

---

Previous: plmgen Up: Healpix/F90 facilities Next: sky_ng_sim Top: Main Page