1<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN"> 2 3<!--Converted with LaTeX2HTML 2018.3 (Released July 19, 2018) --> 4<HTML lang="EN"> 5<HEAD> 6<TITLE>smoothing</TITLE> 7<META NAME="description" CONTENT="smoothing"> 8<META NAME="keywords" CONTENT="facilities"> 9<META NAME="resource-type" CONTENT="document"> 10<META NAME="distribution" CONTENT="global"> 11 12<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=utf-8"> 13<META NAME="viewport" CONTENT="width=device-width, initial-scale=1.0"> 14<META NAME="Generator" CONTENT="LaTeX2HTML v2018.3"> 15 <link rel='apple-touch-icon' sizes='180x180' href='images/favicons/apple-touch-icon.png?v=2017'> 16 <link rel='icon' type='image/png' sizes='32x32' href='images/favicons/favicon-32x32.png?v=2017'> 17 <link rel='icon' type='image/png' sizes='16x16' href='images/favicons/favicon-16x16.png?v=2017'> 18 <link rel='manifest' href='images/favicons/manifest.json?v=2017'> 19 <link rel='mask-icon' href='images/favicons/safari-pinned-tab.svg?v=2017' color='#5bbad5'> 20 <link rel='shortcut icon' href='images/favicons/favicon.ico?v=2017'> 21 <meta name='apple-mobile-web-app-title' content='HEALPix'> 22 <meta name='application-name' content='HEALPix'> 23 <meta name='msapplication-config' content='images/favicons/browserconfig.xml?v=2017'> 24 <meta name='theme-color' content='#ffffff'> 25 26<LINK REL="STYLESHEET" HREF="facilities.css"> 27 28<LINK REL="next" HREF="fac_synfast.htm"> 29<LINK REL="previous" HREF="fac_sky_ng_sim.htm"> 30<LINK REL="next" HREF="fac_synfast.htm"> 31</HEAD> 32 33<body text="#000000" bgcolor="#FFFFFA"> 34 35<DIV CLASS="navigation"><!--Navigation Panel--> 36<A 37 HREF="fac_sky_ng_sim.htm"> 38<IMG WIDTH="63" HEIGHT="24" ALT="previous" SRC="prev.png"></A> 39<A 40 HREF="fac_HEALPix_F90_facilities.htm"> 41<IMG WIDTH="26" HEIGHT="24" ALT="up" SRC="up.png"></A> 42<A 43 HREF="fac_synfast.htm"> 44<IMG WIDTH="37" HEIGHT="24" ALT="next" SRC="next.png"></A> 45<A ID="tex2html131" 46 HREF="fac_TABLE_CONTENTS.htm"> 47<IMG WIDTH="65" HEIGHT="24" ALT="contents" SRC="contents.png"></A> 48<BR> 49<B> Previous:</B> <A 50 HREF="fac_sky_ng_sim.htm">sky_ng_sim</A> 51 52<B>Up:</B> <A 53 HREF="fac_HEALPix_F90_facilities.htm">HEALPix/F90 facilities</A> 54 55<B> Next:</B> <A 56 HREF="fac_synfast.htm">synfast</A> 57<B> Top:</B> <a href="main.htm">Main Page</a></DIV> 58<!--End of Navigation Panel--> 59 60<H1><A ID="SECTION1300"></A> 61<A ID="fac:smoothing"></A> 62<BR> 63 64</H1> 65 66<P> 67<b><font size=+6><FONT COLOR="#CC0000">smoothing</FONT></font></b><hr> 68<H3>This program can be used to convolve a map with a gaussian beam. 69The input map can be given in RING or NESTED scheme and the smoothed map 70is written 71to a FITS file in the RING scheme. 72<BR> 73NOTE: This automated facility is susceptible to problems with non-commutativity 74of discrete spherical harmonics transforms, described in the Recommendations 75for Users of the <A HREF="fac_anafast.htm#fac:anafast">anafast</A> facility. 76If very high accuracy of the results is 77required in the spectral regime of <!-- MATH 78 $\ell > 2\cdot nsmax$ 79 --> 80<SPAN CLASS="MATH"><IMG STYLE="height: 1.69ex; vertical-align: -0.16ex; " SRC="fac_img83.png" 81 ALT="$\ell > 2\cdot nsmax$"></SPAN>, it is recommended 82to choose an iterative computation of the <SPAN CLASS="MATH"><IMG STYLE="height: 1.46ex; vertical-align: -0.45ex; " SRC="fac_img14.png" 83 ALT="$a_{\ell m}$"></SPAN> coefficients. 84</H3> 85Location in HEALPix directory tree: <a href="https://sourceforge.net/p/healpix/code/1005/tree/trunk/src/f90/smoothing/smoothing.f90"><b>src/f90/smoothing/smoothing.f90</b></a> 86 87<P> 88<hr><h1>FORMAT </h1><blockquote><h3>% 89smoothing [options] [parameter_file] 90</h3></blockquote> 91 92<P> 93<hr><H1>COMMAND LINE OPTIONS</H1> 94 <DL COMPACT><DT> 95<B><TT>-d</TT></B> 96<DD><DT> 97<B><TT>--double</TT></B> 98<DD>double precision mode (see 99 <A HREF="fac_Using_HEALPix_Fortran_90_fa.htm#fac:subsec:ioprec">Notes on double/single precision modes</A>) 100 <DT> 101<B><TT>-s</TT></B> 102<DD><DT> 103<B><TT>--single</TT></B> 104<DD>single precision mode (default) 105 </DL> 106 107<P> 108<hr> 109<H1>QUALIFIERS</H1> 110 111 <DL COMPACT><DT> 112<B>simul_type = </B> 113<DD><A ID="fac:smoothing:simul_type"></A>Defines which map(s) to analyse, 1=temperature only, 2=temperature AND polarisation. 114(default= 1) 115 <DT> 116<B>infile = </B> 117<DD><A ID="fac:smoothing:infile"></A>Defines the filename for the FITS file containing the map to be smoothed. 118 (default= 'map.fits') 119 <DT> 120<B>nlmax = </B> 121<DD><A ID="fac:smoothing:nlmax"></A>Defines the <!-- MATH 122 $\ell_{\mathrm{max}}$ 123 --> 124<SPAN CLASS="MATH"><IMG STYLE="height: 2.04ex; vertical-align: -0.45ex; " SRC="fac_img25.png" 125 ALT="$\ell_{\mathrm{max}}$"></SPAN> value for the application. 126(default= 64) 127 <DT> 128<B>iter_order = </B> 129<DD><A ID="fac:smoothing:iter_order"></A>Defines the maximum order of quadrature 130 iteration to be used. (default=0, no iteration). 131For details, see the <A HREF="./sub_map2alm_iterative.htm#sub:map2alm_iterative"><TT>map2alm_iterative</TT></A> routine 132described in the <A ID="tex2html23" 133 HREF="subroutines.htm">”Fortran Subroutines”</A> 134document. 135 <DT> 136<B>fwhm_arcmin = </B> 137<DD><A ID="fac:smoothing:fwhm_arcmin"></A>Defines the FWHM in arcminutes of the gaussian 138beam for the convolution. (default=10) 139 <DT> 140<B>beam_file = </B> 141<DD>Defines the FITS file describing the 142 Legendre window 143 function of the circular beam to be used for the 144 simulation (see <A HREF="fac_Using_HEALPix_Fortran_90_fa.htm#fac:subsec:beamfiles">”Beam window function files” in introduction</A>). 145 If set to an existing file name, it will override the 146 <TT>fhwm_arcmin</TT> given above. default=`' 147<DT> 148<B>outfile = </B> 149<DD><A ID="fac:smoothing:outfile"></A>Defines the filename for the file that will contain 150the smoothed map. (default='map_smoothed.fits') 151 <DT> 152<B>plmfile = </B> 153<DD>Defines the name for an input file 154 containing precomputed Legendre polynomials <SPAN CLASS="MATH"><IMG STYLE="height: 2.04ex; vertical-align: -0.45ex; " SRC="fac_img51.png" 155 ALT="$P_{\ell m}$"></SPAN>. 156(default= no entry -- <FONT COLOR="#CC0000">smoothing</FONT>executes the recursive evaluation 157of <SPAN CLASS="MATH"><IMG STYLE="height: 2.04ex; vertical-align: -0.45ex; " SRC="fac_img51.png" 158 ALT="$P_{\ell m}$"></SPAN>s) 159<DT> 160<B>w8file = </B> 161<DD><A ID="fac:smoothing:w8file"></A>Defines name for an input file containing ring 162 weights in the improved quadrature mode (default= no entry -- 163the name is assumed to be 'weight_ring_n0xxxx.fits' where xxxx is nsmax) 164<DT> 165<B>w8filedir = </B> 166<DD><A ID="fac:smoothing:w8filedir"></A>Gives the directory where the weight files are 167to be found (default= no entry -- smoothing searches in the default 168directories, see <A HREF="fac_Using_HEALPix_Fortran_90_fa.htm#fac:subsec:defdir">”Default file names and directories” in introduction</A>) 169<DT> 170<B>won = </B> 171<DD><A ID="fac:smoothing:won"></A>Set this to 1 if ring-based quadrature weight files are to be used, 172 or to 2 to use pixel-based weight files instead; 173otherwise set it to 0. (default= 0) 174 </DL> 175 176<P> 177<hr> 178<H1>DESCRIPTION</H1> 179<blockquote> 180 181A FITS file containing a <b>HEALPix</b> map in RING or NESTED scheme is read in. 182When dealing with polarized signal maps, the <FONT COLOR="#CC0000">smoothing</FONT> behavior will depend on the value of the <SPAN CLASS="texttt">POLCCONV</SPAN> FITS keyword 183(see <A HREF="./intro_HEALPix_conventions.htm#intro:polcconv">note on POLCCONV</A> in <A ID="tex2html24" 184 HREF="intro.htm">The <b>HEALPix</b> Primer</A>). 185The map is analysed and smoothed in fourier space with a gaussian beam 186of a given FHWM. 187A new map is then synthesized using the smoothed <SPAN CLASS="MATH"><IMG STYLE="height: 1.46ex; vertical-align: -0.45ex; " SRC="fac_img14.png" 188 ALT="$a_{\ell m}$"></SPAN> coefficients. 189For a more accurate application, an iteration of arbitrary order can be applied. 190The output map is stored in <EM>the same scheme</EM> as the input map. 191 192</blockquote> 193 194<P> 195<hr><H1>DATASETS</H1> 196<h3>The following datasets are involved in the <b><FONT COLOR="#CC0000">smoothing</FONT></b> 197 processing.</h3> 198 199<TABLE CELLPADDING=3 BORDER="1"> 200<TR><TH ALIGN="LEFT" VALIGN="TOP" WIDTH=150><SPAN CLASS="textbf">Dataset</SPAN></TH> 201<TH ALIGN="LEFT" VALIGN="TOP" WIDTH=175><SPAN CLASS="textbf">Description</SPAN></TH> 202</TR> 203<TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=150> </TD> 204<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=175> </TD> 205</TR> 206<TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=150> 207data/weight_ring_n0xxxx.fits</TD> 208<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=175>Files containing ring weights 209 for the smoothing improved quadrature mode.</TD> 210</TR> 211<TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=150> </TD> 212<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=175> </TD> 213</TR> 214<TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=150></TD> 215<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=175> </TD> 216</TR> 217</TABLE> 218 219<P> 220<hr><H1>SUPPORT </H1><H3>This section lists those routines and facilities (including those <i>external</i> to the <b>HEALPix</b> distribution) which can assist in the utilisation of <b><FONT COLOR="#CC0000">smoothing</FONT></b>.</H3> 221 <DL COMPACT><DT> 222<B><A HREF="./sub_generate_beam.htm#sub:generate_beam">generate_beam</A></B> 223<DD>This <b>HEALPix</b> Fortran 224subroutine generates or reads the <SPAN CLASS="MATH"><IMG STYLE="height: 2.33ex; vertical-align: -0.68ex; " SRC="fac_img21.png" 225 ALT="$B(\ell)$"></SPAN> window function used in <FONT COLOR="#CC0000">smoothing</FONT> 226 <DT> 227<B><A HREF="fac_map2gif.htm#fac:map2gif">map2gif</A></B> 228<DD>This <b>HEALPix</b> Fortran facility can be used to visualise the 229 input and output maps of <FONT COLOR="#CC0000">smoothing</FONT>. 230 <DT> 231<B><A HREF="./idl_mollview.htm#idl:mollview">mollview</A></B> 232<DD>This <b>HEALPix</b> IDL facility can be used to visualise the 233 input and output maps of <FONT COLOR="#CC0000">smoothing</FONT>. 234 <DT> 235<B><A HREF="fac_synfast.htm#fac:synfast">synfast</A></B> 236<DD>This <b>HEALPix</b> facility can generate a map and also do the smoothing. 237 <DT> 238<B><A HREF="fac_anafast.htm#fac:anafast">anafast</A></B> 239<DD>This <b>HEALPix</b> facility can analyse a smoothed map. 240 </DL> 241 242<P> 243<hr><H1>EXAMPLE # 1:</H1> 244<tt><TABLE CELLPADDING=3> 245<TR><TD ALIGN="LEFT">smoothing</TD> 246<TD ALIGN="LEFT"> </TD> 247</TR> 248</TABLE></tt> 249<blockquote> 250Smoothing runs in interactive mode, self-explanatory. 251</blockquote> 252 253<P> 254 255<P> 256<hr><H1>EXAMPLE # 2:</H1> 257<tt><TABLE CELLPADDING=3> 258<TR><TD ALIGN="LEFT">smoothing filename</TD> 259<TD ALIGN="LEFT"> </TD> 260</TR> 261</TABLE></tt> 262<blockquote>When `filename' is present, smoothing enters the non-interactive mode and parses 263its inputs from the file `filename'. This has the following 264structure: the first entry is a qualifier which announces to the parser 265which input immediately follows. If this input is omitted in the 266input file, the parser assumes the default value. 267If the equality sign is omitted, then the parser ignores the entry. 268In this way comments may also be included in the file. 269In this example, the file contains the following 270qualifiers: 271<BR> 272<tt><A HREF="#fac:smoothing:simul_type">simul_type</A>= 1</tt><br> 273<tt><A HREF="#fac:smoothing:nlmax">nlmax</A>= 64</tt><br> 274<tt><A HREF="#fac:smoothing:infile">infile</A>= map.fits</tt><br> 275<tt><A HREF="#fac:smoothing:outfile">outfile</A>= map_smoothed.fits</tt><br> 276<tt><A HREF="#fac:smoothing:fwhm_arcmin">fwhm_arcmin</A>= 10.</tt><br> 277<tt><A HREF="#fac:smoothing:iter_order">iter_order</A>= 1</tt><br> 278smoothes the <b>HEALPix</b> temperature map contained in `map.fits' with 279a 10 arcmin FWHM beam. The resulting map is saved 280in `map_smoothed.fits'. The map analysis/synthesis was carried 281out using fourier coeffecients up to an <SPAN CLASS="MATH"><IMG STYLE="height: 1.69ex; vertical-align: -0.10ex; " SRC="fac_img5.png" 282 ALT="$\ell$"></SPAN> value of 64. A first 283order iteration of the quadrature was performed.</blockquote> 284 285<P> 286<hr><H1>RELEASE NOTES</H1><blockquote> 287 <DL COMPACT><DT><DD><IMG WIDTH="14" HEIGHT="14" SRC="blueball.png" ALT="*"> 288 Initial release (<b>HEALPix</b> 0.90) 289 <DT><DD><IMG WIDTH="14" HEIGHT="14" SRC="blueball.png" ALT="*"> 290 Extension to polarization and arbitrary <EM>circular</EM> beams (<b>HEALPix</b> 1.20) 291 <DT><DD><IMG WIDTH="14" HEIGHT="14" SRC="blueball.png" ALT="*"> 292 Support for pixel-based quadrature weights when <A HREF="fac_anafast.htm#fac:anafast:won"><SPAN CLASS="texttt">won=2</SPAN></A> (<b>HEALPix</b> 3.40) 293 </DL> 294</blockquote> 295 296<P> 297 298<hr><H1>MESSAGES</H1><h3>This section describes error messages generated by <b><FONT COLOR="#CC0000">smoothing</FONT></b> 299</h3> 300 301<TABLE CELLPADDING=3 BORDER="1"> 302<TR><TH ALIGN="LEFT" VALIGN="TOP" WIDTH=125><SPAN CLASS="textbf">Message</SPAN></TH> 303<TH ALIGN="LEFT" VALIGN="TOP" WIDTH=50><SPAN CLASS="textbf">Severity</SPAN></TH> 304<TH ALIGN="LEFT" VALIGN="TOP" WIDTH=175><SPAN CLASS="textbf">Text</SPAN></TH> 305</TR> 306<TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=125> </TD> 307<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=50> </TD> 308<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=175> </TD> 309</TR> 310<TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=125> 311can not allocate memory for array xxx</TD> 312<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=50>Fatal</TD> 313<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=175>You do not have 314 sufficient system resources to run this 315 facility at the map resolution you required. 316 Try a lower map resolution.</TD> 317</TR> 318<TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=125> </TD> 319<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=50> </TD> 320<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=175> </TD> 321</TR> 322<TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=125></TD> 323<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=50> </TD> 324<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=175> </TD> 325</TR> 326</TABLE> 327 328<P> 329 330<P> 331 332 333<P> 334 335<P> 336 337<DIV CLASS="navigation"><HR> 338<!--Navigation Panel--> 339<A 340 HREF="fac_sky_ng_sim.htm"> 341<IMG WIDTH="63" HEIGHT="24" ALT="previous" SRC="prev.png"></A> 342<A 343 HREF="fac_HEALPix_F90_facilities.htm"> 344<IMG WIDTH="26" HEIGHT="24" ALT="up" SRC="up.png"></A> 345<A 346 HREF="fac_synfast.htm"> 347<IMG WIDTH="37" HEIGHT="24" ALT="next" SRC="next.png"></A> 348<A ID="tex2html131" 349 HREF="fac_TABLE_CONTENTS.htm"> 350<IMG WIDTH="65" HEIGHT="24" ALT="contents" SRC="contents.png"></A> 351<BR> 352<B> Previous:</B> <A 353 HREF="fac_sky_ng_sim.htm">sky_ng_sim</A> 354 355<B>Up:</B> <A 356 HREF="fac_HEALPix_F90_facilities.htm">HEALPix/F90 facilities</A> 357 358<B> Next:</B> <A 359 HREF="fac_synfast.htm">synfast</A> 360<B> Top:</B> <a href="main.htm">Main Page</a></DIV> 361<!--End of Navigation Panel--> 362<ADDRESS> 363Version 3.50, 2018-12-10 364</ADDRESS> 365</BODY> 366</HTML> 367