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

Fast directional continuous spherical wavelet transform algorithms

We describe the construction of a spherical wavelet analysis through the inverse stereographic projection of the Euclidean planar wavelet framework, introduced originally by Antoine and Vandergheynst and developed further by Wiaux et al. Fast algorithms for performing the directional continuous wavelet analysis on the unit sphere are presented. The fast directional algorithm, based on the fast spherical convolution algorithm developed by Wandelt and Gorski, provides a saving of O(sqrt(Npix)) over a direct quadrature implementation for Npix pixels on the sphere, and allows one to perform a directional spherical wavelet analysis of a 10^6 pixel map on a personal computer.Comment: 10 pages, 3 figures, replaced to match version accepted by IEEE Trans. Sig. Pro

Markov chain Monte Carlo analysis of Bianchi VII_h models

We have extended the analysis of Jaffe et al. to a complete Markov chain Monte Carlo (MCMC) study of the Bianchi type ${\rm VII_h}$ models including a dark energy density, using 1-year and 3-year Wilkinson Microwave Anisotropy Probe (WMAP) cosmic microwave background (CMB) data. Since we perform the analysis in a Bayesian framework our entire inference is contained in the multidimensional posterior distribution from which we can extract marginalised parameter constraints and the comparative Bayesian evidence. Treating the left-handed Bianchi CMB anisotropy as a template centred upon the `cold-spot' in the southern hemisphere, the parameter estimates derived for the total energy density, `tightness' and vorticity from 3-year data are found to be: $\Omega_{tot} = 0.43\pm 0.04$, $h = 0.32^{+0.02}_{-0.13}$, $\omega = 9.7^{+1.6}_{-1.5}\times 10^{-10}$ with orientation $\gamma = {337^{\circ}}^{+17}_{-23}$). This template is preferred by a factor of roughly unity in log-evidence over a concordance cosmology alone. A Bianchi type template is supported by the data only if its position on the sky is heavily restricted. The low total energy density of the preferred template, implies a geometry that is incompatible with cosmologies inferred from recent CMB observations. Jaffe et al. found that extending the Bianchi model to include a term in $\Omega_{\Lambda}$ creates a degeneracy in the $\Omega_m - \Omega_{\Lambda}$ plane. We explore this region fully by MCMC and find that the degenerate likelihood contours do not intersect areas of parameter space that 1 or 3 year WMAP data would prefer at any significance above $2\sigma$. Thus we can confirm that a physical Bianchi ${\rm VII_h}$ model is not responsible for this signature.Comment: 8 pages, 10 figures, significant update to include more accurate results and conclusions to match version accepted by MNRA

Probing dark energy with steerable wavelets through correlation of WMAP and NVSS local morphological measures

Using local morphological measures on the sphere defined through a steerable wavelet analysis, we examine the three-year WMAP and the NVSS data for correlation induced by the integrated Sachs-Wolfe (ISW) effect. The steerable wavelet constructed from the second derivative of a Gaussian allows one to define three local morphological measures, namely the signed-intensity, orientation and elongation of local features. Detections of correlation between the WMAP and NVSS data are made with each of these morphological measures. The most significant detection is obtained in the correlation of the signed-intensity of local features at a significance of 99.9%. By inspecting signed-intensity sky maps, it is possible for the first time to see the correlation between the WMAP and NVSS data by eye. Foreground contamination and instrumental systematics in the WMAP data are ruled out as the source of all significant detections of correlation. Our results provide new insight on the ISW effect by probing the morphological nature of the correlation induced between the cosmic microwave background and large scale structure of the Universe. Given the current constraints on the flatness of the Universe, our detection of the ISW effect again provides direct and independent evidence for dark energy. Moreover, this new morphological analysis may be used in future to help us to better understand the nature of dark energy.Comment: 12 pages, 10 figures, replaced to match version accepted by MNRA

Ephemeral properties and the illusion of microscopic particles

Founding our analysis on the Geneva-Brussels approach to quantum mechanics, we use conventional macroscopic objects as guiding examples to clarify the content of two important results of the beginning of twentieth century: Einstein-Podolsky-Rosen's reality criterion and Heisenberg's uncertainty principle. We then use them in combination to show that our widespread belief in the existence of microscopic particles is only the result of a cognitive illusion, as microscopic particles are not particles, but are instead the ephemeral spatial and local manifestations of non-spatial and non-local entities

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

Detection of the ISW effect and corresponding dark energy constraints made with directional spherical wavelets

Using a directional spherical wavelet analysis we detect the integrated Sachs-Wolfe (ISW) effect, indicated by a positive correlation between the first-year Wilkinson Microwave Anisotropy Probe (WMAP) and NRAO VLA Sky Survey (NVSS) data. Detections are made using both a directional extension of the spherical Mexican hat wavelet and the spherical butterfly wavelet. We examine the possibility of foreground contamination and systematics in the WMAP data and conclude that these factors are not responsible for the signal that we detect. The wavelet analysis inherently enables us to localise on the sky those regions that contribute most strongly to the correlation. On removing these localised regions the correlation that we detect is reduced in significance, as expected, but it is not eliminated, suggesting that these regions are not the sole source of correlation between the data. This finding is consistent with predictions made using the ISW effect, where one would expect weak correlations over the entire sky. In a flat universe the detection of the ISW effect provides direct and independent evidence for dark energy. We use our detection to constrain dark energy parameters by deriving a theoretical prediction for the directional wavelet covariance statistic for a given cosmological model. Comparing these predictions with the data we place constraints on the equation-of-state parameter $w$ and the vacuum energy density $\Omega_\Lambda$. We also consider the case of a pure cosmological constant, i.e. $w=-1$. For this case we rule out a zero cosmological constant at greater than the 99.9% significance level. All parameter estimates that we obtain are consistent with the standand cosmological concordance model values.Comment: 16 pages, 13 figures; replaced to match version accepted by MNRA