Parallelized Inference for Gravitational-Wave Astronomy

Bayesian inference is the workhorse of gravitational-wave astronomy, for example, determining the mass and spins of merging black holes, revealing the neutron star equation of state, and unveiling the population properties of compact binaries. The science enabled by these inferences comes with a computational cost that can limit the questions we are able to answer. This cost is expected to grow. As detectors improve, the detection rate will go up, allowing less time to analyze each event. Improvement in low-frequency sensitivity will yield longer signals, increasing the number of computations per event. The growing number of entries in the transient catalog will drive up the cost of population studies. While Bayesian inference calculations are not entirely parallelizable, key components are embarrassingly parallel: calculating the gravitational waveform and evaluating the likelihood function. Graphical processor units (GPUs) are adept at such parallel calculations. We report on progress porting gravitational-wave inference calculations to GPUs. Using a single code - which takes advantage of GPU architecture if it is available - we compare computation times using modern GPUs (NVIDIA P100) and CPUs (Intel Gold 6140). We demonstrate speed-ups of $\sim 50 \times$ for compact binary coalescence gravitational waveform generation and likelihood evaluation and more than $100\times$ for population inference within the lifetime of current detectors. Further improvement is likely with continued development. Our python-based code is publicly available and can be used without familiarity with the parallel computing platform, CUDA.Comment: 5 pages, 4 figures, submitted to PRD, code can be found at https://github.com/ColmTalbot/gwpopulation https://github.com/ColmTalbot/GPUCBC https://github.com/ADACS-Australia/ADACS-SS18A-RSmith Add demonstration of improvement in BNS spi

Elliptical Distance Transforms And Object Splitting

The classical morphological method to separate fused objects in binary images is to use the watershed transform on the complement of the distance transform of the binary image. This method assumes roughly disk-like objects and cannot separate objects when they are fused together beyond a certain point. In this paper we revisit the issue by assuming that fused objects are unions of ellipses rather than mere disks. The problem is recast in terms of finding the constituent primary grains given a boolean model of ellipses. To this end, we modify the well-known pseudo-Euclidean distance transform algorithm to generate arbitrary elliptical distance transforms to reduce the dimension of the problem and we present a goodness-of-fit measure that allows us to select ellipses. The results of the methods are given on both synthetic sample boolean models and real data

Open Medical Gesture: An Open-Source Experiment in Naturalistic Physical Interactions for Mixed and Virtual Reality Simulations

Mixed Reality (MR) and Virtual Reality (VR) simulations are hampered by requirements for hand controllers or attempts to perseverate in use of two-dimensional computer interface paradigms from the 1980s. From our efforts to produce more naturalistic interactions for combat medic training for the military, USC has developed an open-source toolkit that enables direct hand controlled responsive interactions that is sensor independent and can function with depth sensing cameras, webcams or sensory gloves. Natural approaches we have examined include the ability to manipulate virtual smart objects in a similar manner to how they are used in the real world. From this research and review of current literature, we have discerned several best approaches for hand-based human computer interactions which provide intuitive, responsive, useful, and low frustration experiences for VR users.Comment: AHFE 202

Aerosol chemical composition in Asian continental outflow during the TRACE-P campaign: Comparison with PEM-West B

Aerosol associated soluble ions and the radionuclide tracers 7Be and 210Pb were quantified in 414 filter samples collected in spring 2001 from the DC-8 during the Transport and Chemical Evolution over the Pacific (TRACE-P) campaign. Binning the data into near Asia (flights from Hong Kong and Japan) and remote Pacific (all other flights) revealed large enhancements of NO3−, SO4=, C2O4=, NH4+, K+, Mg2+, and Ca2+ near Asia. The boundary layer and lower troposphere were most strongly influenced by continental outflow, and the largest enhancements were seen in Ca2+ (a dust tracer) and NO3− (reflecting uptake of HNO3 onto the dust). Comparing the TRACE P near Asia bin with earlier results from the same region during PEM-West B (in 1994) shows at least twofold enhancements during TRACE P in most of the ions listed above. Calcium and NO3− were most enhanced in this comparison as well (more than sevenfold higher in the boundary layer and threefold higher in the lower troposphere). Independent estimation of Asian emissions of gaseous precursors of the aerosol-associated ions suggest only small changes between the two missions, and precipitation fields do not suggest any significant difference in the efficiency of the primary sink, precipitation scavenging. It thus appears that with the possible exception of dust, the enhancements of aerosol-associated species during TRACE P cannot be explained by stronger sources or weaker sinks. We argue that the enhancements largely reflect the fact that TRACE P focused on characterizing Asian outflow, and thus the DC-8 was more frequently flown into regions that were influenced by well-organized flow off the continent

Large-scale distributions of tropospheric nitric, formic, and acetic acids over the western Pacific basin during wintertime

We report here measurements of the acidic gases nitric (HNO3), formic (HCOOH), and acetic (CH3COOH) over the western Pacific basin during the February-March 1994 Pacific Exploratory Mission-West (PEM-West B). These data were obtained aboard the NASA DC-8 research aircraft as it flew missions in the altitude range of 0.3–12.5 km over equatorial regions near Guam and then further westward encompassing the entire Pacific Rim arc. Aged marine air over the equatorial Pacific generally exhibited mixing ratios of acidic gases \u3c100 parts per trillion by volume (pptv). Near the Asian continent, discrete plumes encountered below 6 km altitude contained up to 8 parts per billion by volume (ppbv) HNO3 and 10 ppbv HCOOH and CH3COOH. Overall there was a general correlation between mixing ratios of acidic gases with those of CO, C2H2, and C2Cl4, indicative of emissions from combustion and industrial sources. The latitudinal distributions of HNO3 and CO showed that the largest mixing ratios were centered around 15°N, while HCOOH, CH3COOH, and C2Cl4 peaked at 25°N. The mixing ratios of HCOOH and CH3COOH were highly correlated (r2 = 0.87) below 6 km altitude, with a slope (0.89) characteristic of the nongrowing season at midlatitudes in the northern hemisphere. Above 6 km altitude, HCOOH and CH3COOH were marginally correlated (r2 = 0.50), and plumes well defined by CO, C2H2, and C2Cl4 were depleted in acidic gases, most likely due to scavenging during vertical transport of air masses through convective cloud systems over the Asian continent. In stratospheric air masses, HNO3 mixing ratios were several parts per billion by volume (ppbv), yielding relationships with O3 and N2O consistent with those previously reported for NOy

A mathematical model of a single main rotor helicopter for piloted simulation

A mathematical model, suitable for piloted simulation of the flying qualities of helicopters, is a nonlinear, total force and moment model of a single main rotor helicopter. The model has ten degrees of freedom: six rigid body, three rotor flapping, and the rotor rotational degrees of freedom. The rotor model assumes rigid blades with rotor forces and moments radially integrated and summed about the azimuth. The fuselage aerodynamic model uses a detailed representation over a nominal angle of attack and sideslip range of + or - 15 deg., as well as a simplified curve fit at large angles of attack or sideslip. Stabilizing surface aerodynamics are modeled with a lift curve slope between stall limits and a general curve fit for large angles of attack. A generalized stability and control augmentation system is described. Additional computer subroutines provide options for a simplified engine/governor model, atmospheric turbulence, and a linearized six degree of freedom dynamic model for stability and control analysis

Photochemistry in the arctic free troposphere: NOx budget and the role of odd nitrogen reservoir recycling

The budget of nitrogen oxides (NOx) in the arctic free troposphere is calculated with a constrained photochemical box model using aircraft observations from the Tropospheric O3 Production about the Spring Equinox (TOPSE) campaign between February and May. Peroxyacetic nitric anhydride (PAN) was observed to be the dominant odd nitrogen species (NOy) in the arctic free troposphere and showed a pronounced seasonal increase in mixing ratio. When constrained to observed acetaldehyde (CH3CHO) mixing ratios, the box model calculates unrealistically large net NOx losses due to PAN formation (62pptv/day for May, 1-3km). Thus, given our current understanding of atmospheric chemistry, these results cast doubt on the robustness of the CH3CHO observations during TOPSE. When CH3CHO was calculated to steady state in the box model, the net NOx loss to PAN was of comparable magnitude to the net NOx loss to HNO3 (NO2 reaction with OH) for spring conditions. During the winter, net NOx loss due to N2O5 hydrolysis dominates other NOx loss processes and is near saturation with respect to further increases in aerosol surface area concentration. NOx loss due to N2O5 hydrolysis is sensitive to latitude and month due to changes in diurnal photolysis (sharp day-night transitions in winter to continuous sun in spring for the arctic). Near NOx sources, HNO4 is a net sink for NOx; however, for more aged air masses HNO4 is a net source for NOx, largely countering the NOx loss to PAN, N2O5 and HNO3. Overall, HNO4 chemistry impacts the timing of NOx decay and O3 production; however, the cumulative impact on O3 and NOx mixing ratios after a 20-day trajectory is minimal. © 2003 Elsevier Science Ltd. All rights reserved