Challenges in using GPUs for the real-time reconstruction of digital hologram images

This is the pre-print version of the final published paper that is available from the link below.In-line holography has recently made the transition from silver-halide based recording media, with laser reconstruction, to recording with large-area pixel detectors and computer-based reconstruction. This form of holographic imaging is an established technique for the study of fine particulates, such as cloud or fuel droplets, marine plankton and alluvial sediments, and enables a true 3D object field to be recorded at high resolution over a considerable depth. The move to digital holography promises rapid, if not instantaneous, feedback as it avoids the need for the time-consuming chemical development of plates or film film and a dedicated replay system, but with the growing use of video-rate holographic recording, and the desire to reconstruct fully every frame, the computational challenge becomes considerable. To replay a digital hologram a 2D FFT must be calculated for every depth slice desired in the replayed image volume. A typical hologram of ~100 μm particles over a depth of a few hundred millimetres will require O(10^3) 2D FFT operations to be performed on a hologram of typically a few million pixels. In this paper we discuss the technical challenges in converting our existing reconstruction code to make efficient use of NVIDIA CUDA-based GPU cards and show how near real-time video slice reconstruction can be obtained with holograms as large as 4096 by 4096 pixels. Our performance to date for a number of different NVIDIA GPU running under both Linux and Microsoft Windows is presented. The recent availability of GPU on portable computers is discussed and a new code for interactive replay of digital holograms is presented

Mosaicking with cosmic microwave background interferometers

Measurements of cosmic microwave background (CMB) anisotropies by interferometers offer several advantages over single-dish observations. The formalism for analyzing interferometer CMB data is well developed in the flat-sky approximation, valid for small fields of view. As the area of sky is increased to obtain finer spectral resolution, this approximation needs to be relaxed. We extend the formalism for CMB interferometry, including both temperature and polarization, to mosaics of observations covering arbitrarily large areas of the sky, with each individual pointing lying within the flat-sky approximation. We present a method for computing the correlation between visibilities with arbitrary pointing centers and baselines and illustrate the effects of sky curvature on the l-space resolution that can be obtained from a mosaic.Comment: 9 pages; submitted to Ap

Systematic Errors in Cosmic Microwave Background Interferometry

Cosmic microwave background (CMB) polarization observations will require superb control of systematic errors in order to achieve their full scientific potential, particularly in the case of attempts to detect the B modes that may provide a window on inflation. Interferometry may be a promising way to achieve these goals. This paper presents a formalism for characterizing the effects of a variety of systematic errors on interferometric CMB polarization observations, with particular emphasis on estimates of the B-mode power spectrum. The most severe errors are those that couple the temperature anisotropy signal to polarization; such errors include cross-talk within detectors, misalignment of polarizers, and cross-polarization. In a B mode experiment, the next most serious category of errors are those that mix E and B modes, such as gain fluctuations, pointing errors, and beam shape errors. The paper also indicates which sources of error may cause circular polarization (e.g., from foregrounds) to contaminate the cosmologically interesting linear polarization channels, and conversely whether monitoring of the circular polarization channels may yield useful information about the errors themselves. For all the sources of error considered, estimates of the level of control that will be required for both E and B mode experiments are provided. Both experiments that interfere linear polarizations and those that interfere circular polarizations are considered. The fact that circular experiments simultaneously measure both linear polarization Stokes parameters in each baseline mitigates some sources of error.Comment: 19 pages, 9 figures, submitted to Phys. Rev.

Bayes-X: a Bayesian inference tool for the analysis of X-ray observations of galaxy clusters

We present the first public release of our Bayesian inference tool, Bayes-X, for the analysis of X-ray observations of galaxy clusters. We illustrate the use of Bayes-X by analysing a set of four simulated clusters at z=0.2-0.9 as they would be observed by a Chandra-like X-ray observatory. In both the simulations and the analysis pipeline we assume that the dark matter density follows a spherically-symmetric Navarro, Frenk and White (NFW) profile and that the gas pressure is described by a generalised NFW (GNFW) profile. We then perform four sets of analyses. By numerically exploring the joint probability distribution of the cluster parameters given simulated Chandra-like data, we show that the model and analysis technique can robustly return the simulated cluster input quantities, constrain the cluster physical parameters and reveal the degeneracies among the model parameters and cluster physical parameters. We then analyse Chandra data on the nearby cluster, A262, and derive the cluster physical profiles. To illustrate the performance of the Bayesian model selection, we also carried out analyses assuming an Einasto profile for the matter density and calculated the Bayes factor. The results of the model selection analyses for the simulated data favour the NFW model as expected. However, we find that the Einasto profile is preferred in the analysis of A262. The Bayes-X software, which is implemented in Fortran 90, is available at http://www.mrao.cam.ac.uk/facilities/software/bayesx/.Comment: 22 pages, 11 figure

Classifying LISA gravitational wave burst signals using Bayesian evidence

We consider the problem of characterisation of burst sources detected with the Laser Interferometer Space Antenna (LISA) using the multi-modal nested sampling algorithm, MultiNest. We use MultiNest as a tool to search for modelled bursts from cosmic string cusps, and compute the Bayesian evidence associated with the cosmic string model. As an alternative burst model, we consider sine-Gaussian burst signals, and show how the evidence ratio can be used to choose between these two alternatives. We present results from an application of MultiNest to the last round of the Mock LISA Data Challenge, in which we were able to successfully detect and characterise all three of the cosmic string burst sources present in the release data set. We also present results of independent trials and show that MultiNest can detect cosmic string signals with signal-to-noise ratio (SNR) as low as ~7 and sine-Gaussian signals with SNR as low as ~8. In both cases, we show that the threshold at which the sources become detectable coincides with the SNR at which the evidence ratio begins to favour the correct model over the alternative.Comment: 21 pages, 11 figures, accepted by CQG; v2 has minor changes for consistency with accepted versio

Sensitivity of a Bolometric Interferometer to the CMB power spectrum

Context. The search for B-mode polarization fluctuations in the Cosmic Microwave Background is one of the main challenges of modern cosmology. The expected level of the B-mode signal is very low and therefore requires the development of highly sensitive instruments with low systematic errors. An appealing possibility is bolometric interferometry. Aims. We compare in this article the sensitivity on the CMB angular power spectrum achieved with direct imaging, heterodyne and bolometric interferometry. Methods. Using a simple power spectrum estimator, we calculate its variance leading to the counterpart for bolometric interferometry of the well known Knox formula for direct imaging. Results. We find that bolometric interferometry is less sensitive than direct imaging. However, as expected, it is finally more sensitive than heterodyne interferometry due to the low noise of the bolometers. It therefore appears as an alternative to direct imagers with different and possibly lower systematic errors, mainly due to the absence of an optical setup in front of the horns.Comment: 5 pages, 3 figures. This last version matches the published version (Astronomy and Astrophysics 491 3 (2008) 923-927). Sensitivity of Heterodyne Interferometers modified by a factor of tw

The effect of point sources on satellite observations of the cosmic microwave background

We study the effect of extragalactic point sources on satellite observations of the cosmic microwave background (CMB). In order to separate the contributions due to different foreground components, a maximum-entropy method is applied to simulated observations by the Planck Surveyor satellite. In addition to point sources, the simulations include emission from the CMB and the kinetic and thermal Sunyaev-Zel'dovich (SZ) effects from galaxy clusters, as well as Galactic dust, free-free and synchrotron emission. We find that the main input components are faithfully recovered and, in particular, that the quality of the CMB reconstruction is only slightly reduced by the presence of point sources. In addition, we find that it is possible to recover accurate point source catalogues at each of the Planck Surveyor observing frequencies.Comment: 12 pages, 9 figures, submitted to MNRA

Filtering techniques for the detection of Sunyaev-Zel'dovich clusters in multifrequency CMB maps

The problem of detecting Sunyaev-Zel'dovich (SZ) clusters in multifrequency CMB observations is investigated using a number of filtering techniques. A multifilter approach is introduced, which optimizes the detection of SZ clusters on microwave maps. An alternative method is also investigated, in which maps at different frequencies are combined in an optimal manner so that existing filtering techniques can be applied to the single combined map. The SZ profiles are approximated by the circularly-symmetric template $\tau (x) = [1 +(x/r_c)^2]^{-\lambda}$, with $\lambda \simeq \tfrac{1}{2}$ and $x\equiv |\vec{x}|$, where the core radius $r_c$ and the overall amplitude of the effect are not fixed a priori, but are determined from the data. The background emission is modelled by a homogeneous and isotropic random field, characterized by a cross-power spectrum $P_{\nu_1 \nu_2}(q)$ with $q\equiv |\vec{q}|$. The filtering methods are illustrated by application to simulated Planck observations of a $12.8^\circ \times 12.8^\circ$ patch of sky in 10 frequency channels. Our simulations suggest that the Planck instrument should detect $\approx 10000$ SZ clusters in 2/3 of the sky. Moreover, we find the catalogue to be complete for fluxes $S > 170$ mJy at 300 GHz.Comment: 12 pages, 7 figures; Corrected figures. Submitted to MNRA