Statistical analysis of network data and evolution on GPUs: High-performance statistical computing

Network analysis typically involves as set of repetitive tasks that are particularly amenable to poor-man's parallelization. This is therefore an ideal application are for GPU architectures, which help to alleviate the tedium inherent to statistically sound analysis of network data. Here we will illustrate the use of GPUs in a range of applications, which include percolation processes on networks, the evolution of protein-protein interaction networks, and the fusion of different types of biomedical and disease data in the context of molecular interaction networks. We will pay particular attention to the numerical performance of different routines that are frequently invoked in network analysis problems. We conclude with a review over recent developments in the generation of random numbers that address the specific requirements posed by GPUs and high-performance computing needs

General-Relativistic MHD for the Numerical Construction of Dynamical Spacetimes

We assemble the equations of general relativistic magnetohydrodynamics (MHD) in 3+1 form. These consist of the complete coupled set of Maxwell equations for the electromagnetic field, Einstein's equations for the gravitational field, and the equations of relativistic MHD for a perfectly conducting ideal gas. The adopted form of the equations is suitable for evolving numerically a relativistic MHD fluid in a dynamical spacetime characterized by a strong gravitational field.Comment: 8 pages; scheduled for March 10 issue of Ap

Corrections and Comments on the Multipole Moments of Axisymmetric Electrovacuum Spacetimes

Following the method of Hoenselaers and Perj\'{e}s we present a new corrected and dimensionally consistent set of multipole gravitational and electromagnetic moments for stationary axisymmetric spacetimes. Furthermore, we use our results to compute the multipole moments, both gravitational and electromagnetic, of a Kerr-Newman black hole.Comment: This is a revised and corrected versio

Gravitational Wavetrains in the Quasi-Equilibrium Approximation: A Model Problem in Scalar Gravitation

A quasi-equilibrium (QE) computational scheme was recently developed in general relativity to calculate the complete gravitational wavetrain emitted during the inspiral phase of compact binaries. The QE method exploits the fact that the the gravitational radiation inspiral timescale is much longer than the orbital period everywhere outside the ISCO. Here we demonstrate the validity and advantages of the QE scheme by solving a model problem in relativistic scalar gravitation theory. By adopting scalar gravitation, we are able to numerically track without approximation the damping of a simple, quasi-periodic radiating system (an oscillating spherical matter shell) to final equilibrium, and then use the exact numerical results to calibrate the QE approximation method. In particular, we calculate the emitted gravitational wavetrain three different ways: by integrating the exact coupled dynamical field and matter equations, by using the scalar-wave monopole approximation formula (corresponding to the quadrupole formula in general relativity), and by adopting the QE scheme. We find that the monopole formula works well for weak field cases, but fails when the fields become even moderately strong. By contrast, the QE scheme remains quite reliable for moderately strong fields, and begins to breakdown only for ultra-strong fields. The QE scheme thus provides a promising technique to construct the complete wavetrain from binary inspiral outside the ISCO, where the gravitational fields are strong, but where the computational resources required to follow the system for more than a few orbits by direct numerical integration of the exact equations are prohibitive.Comment: 15 pages, 14 figure

Implementing an apparent-horizon finder in three dimensions

Locating apparent horizons is not only important for a complete understanding of numerically generated spacetimes, but it may also be a crucial component of the technique for evolving black-hole spacetimes accurately. A scheme proposed by Libson et al., based on expanding the location of the apparent horizon in terms of symmetric trace-free tensors, seems very promising for use with three-dimensional numerical data sets. In this paper, we generalize this scheme and perform a number of code tests to fully calibrate its behavior in black-hole spacetimes similar to those we expect to encounter in solving the binary black-hole coalescence problem. An important aspect of the generalization is that we can compute the symmetric trace-free tensor expansion to any order. This enables us to determine how far we must carry the expansion to achieve results of a desired accuracy. To accomplish this generalization, we describe a new and very convenient set of recurrence relations which apply to symmetric trace-free tensors.Comment: 14 pages (RevTeX 3.0 with 3 figures

Radiative Falloff in Neutron Star Spacetimes

We systematically study late-time tails of scalar waves propagating in neutron star spacetimes. We consider uniform density neutron stars, for which the background spacetime is analytic and the compaction of the star can be varied continously between the Newtonian limit 2M/R << 1 and the relativistic Buchdahl limit 2M/R = 8/9. We study the reflection of a finite wave packet off neutron stars of different compactions 2M/R and find that a Newtonian, an intermediate, and a highly relativistic regime can be clearly distinguished. In the highly relativistic regime, the reflected signal is dominated by quasi-periodic peaks, which originate from the wave packet bouncing back and forth between the center of the star and the maximum of the background curvature potential at R ~ 3 M. Between these peaks, the field decays according to a power-law. In the Buchdahl limit 2M/R -> 8/9 the light travel time between the center and the maximum or the curvature potential grows without bound, so that the first peak arrives only at infinitely late time. The modes of neutron stars can therefore no longer be excited in the ultra-relativistic limit, and it is in this sense that the late-time radiative decay from neutron stars looses all its features and gives rise to power-law tails reminiscent of Schwarzschild black holes.Comment: 10 pages, 7 figures, to appear in PR

Computing the Complete Gravitational Wavetrain from Relativistic Binary Inspiral

We present a new method for generating the nonlinear gravitational wavetrain from the late inspiral (pre-coalescence) phase of a binary neutron star system by means of a numerical evolution calculation in full general relativity. In a prototype calculation, we produce 214 wave cycles from corotating polytropes, representing the final part of the inspiral phase prior to reaching the ISCO. Our method is based on the inequality that the orbital decay timescale due to gravitational radiation is much longer than an orbital period and the approximation that gravitational radiation has little effect on the structure of the stars. We employ quasi-equilibrium sequences of binaries in circular orbit for the matter source in our field evolution code. We compute the gravity-wave energy flux, and, from this, the inspiral rate, at a discrete set of binary separations. From these data, we construct the gravitational waveform as a continuous wavetrain. Finally, we discuss the limitations of our current calculation, planned improvements, and potential applications of our method to other inspiral scenarios.Comment: 4 pages, 4 figure

Evolving Einstein's Field Equations with Matter: The ``Hydro without Hydro'' Test

We include matter sources in Einstein's field equations and show that our recently proposed 3+1 evolution scheme can stably evolve strong-field solutions. We insert in our code known matter solutions, namely the Oppenheimer-Volkoff solution for a static star and the Oppenheimer-Snyder solution for homogeneous dust sphere collapse to a black hole, and evolve the gravitational field equations. We find that we can evolve stably static, strong-field stars for arbitrarily long times and can follow dust sphere collapse accurately well past black hole formation. These tests are useful diagnostics for fully self-consistent, stable hydrodynamical simulations in 3+1 general relativity. Moreover, they suggest a successive approximation scheme for determining gravitational waveforms from strong-field sources dominated by longitudinal fields, like binary neutron stars: approximate quasi-equilibrium models can serve as sources for the transverse field equations, which can be evolved without having to re-solve the hydrodynamical equations (``hydro without hydro'').Comment: 4 postscript figures. Submitted to Phys. Rev. D15 as a Brief Repor

Insights from 20 years of temperature parallel measurements in Mauritius around the turn of the 20th century

There is considerable import in creating more complete, better understood holdings of early meteorological data. Such data permit an improved understanding of climate variability and long-term changes. Early records are particularly incomplete in the tropics, with implications for estimates of global and regional temperature. There is also a relatively low level of scientific understanding of how these early measurements were made and, as a result, of their homogeneity and comparability to more modern techniques and measurements. Herein we describe and analyse a newly rescued set of long-term, up to six-way parallel measurements undertaken over 1884–1903 in Mauritius, an island situated in the southern Indian Ocean. Data include (i) measurements from a well-ventilated room, (ii) a shaded thermograph, (iii) instruments housed in a manner broadly equivalent to a modern Stevenson screen, (iv) a set of measurements by a hygrometer mounted in a Stevenson screen, and for a much shorter period (v) two additional Stevenson screen configurations. All measurements were undertaken within an ∼ 80 m radius of each other. To our knowledge this is the first such multidecadal multi-instrument assessment of meteorological instrument transition impacts ever undertaken, providing potentially unique insights. The intercomparison also considers the impact of different ways of deriving daily and monthly averages. The long-term comparison is sufficient to robustly characterize systematic offsets between all the instruments and seasonally varying impacts. Differences between all techniques range from tenths of a degree Celsius to more than 1 ∘C and are considerably larger for maximum and minimum temperatures than for means or averages. Systematic differences of several tenths of a degree Celsius also exist for the different ways of deriving average and mean temperatures. All differences, except two average temperature series pairs, are significant at the 0.01 level using a paired t test. Given that all thermometers were regularly calibrated against a primary Kew standard thermometer maintained by the observatory, this analysis highlights significant impacts of instrument exposure, housing, siting, and measurement practices in early meteorological records. These results reaffirm the importance of thoroughly assessing the homogeneity of early meteorological records