Accelerating frequency-domain diffuse optical tomographic image reconstruction using graphics processing units

Diffuse optical tomographic image reconstruction uses advanced numerical models that are computationally costly to be implemented in the real time. The graphics processing units (GPUs) offer desktop massive parallelization that can accelerate these computations. An open-source GPU-accelerated linear algebra library package is used to compute the most intensive matrix-matrix calculations and matrix decompositions that are used in solving the system of linear equations. These open-source functions were integrated into the existing frequency-domain diffuse optical image reconstruction algorithms to evaluate the acceleration capability of the GPUs (NVIDIA Tesla C 1060) with increasing reconstruction problem sizes. These studies indicate that single precision computations are sufficient for diffuse optical tomographic image reconstruction. The acceleration per iteration can be up to 40, using GPUs compared to traditional CPUs in case of three-dimensional reconstruction, where the reconstruction problem is more underdetermined, making the GPUs more attractive in the clinical settings. The current limitation of these GPUs in the available onboard memory (4 GB) that restricts the reconstruction of a large set of optical parameters, more than 13, 377. (C) 2010 Society of Photo-Optical Instrumentation Engineers. DOI: 10.1117/1.3506216]

Massive black-hole binary inspirals: results from the LISA parameter estimation taskforce

The LISA Parameter Estimation Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models.

Compressible boundary-layer flow at a three-dimensional stagnation point with massive blowing

The effect of massive blowing rates on the steady laminar compressible boundary-layer flow with variable gas properties at a 3-dim. stagnation point (which includes both nodal and saddle points of attachment) has been studied. The equations governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique for nodal points of attachment but employing a parametric differentiation technique instead of quasilinearization for saddle points of attachment. It is found that the effect of massive blowing rates is to move the viscous layer away from the surface. The effect of the variation of the density- viscosity product across the boundary layer is found to be negligible for massive blowing rates but significant for moderate blowing rates. The velocity profiles in the transverse direction for saddle points of attachment in the presence of massive blowing show both the reverse flow as well as velocity overshoot.

Semi-similar solution of an unsteady compressible three-dimensional stagnation point boundary layer flow with massive blowing

A semi-similar solution of an unsteady laminar compressible three-dimensional stagnation point boundary layer flow with massive blowing has been obtained when the free stream velocity varies arbitrarily with time. The resulting partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme with a quasi-linearization technique in the nodal point region and an implicit finite-difference scheme with a parametric differentiation technique in the saddle point region. The results have been obtained for two particular unsteady free stream velocity distributions: (i) an accelerating stream and (ii) a fluctuating stream. Results show that the skin-friction and heat-transfer parameters respond significantly to the time dependent arbitrary free stream velocity. Velocity and enthalpy profiles approach their free stream values faster as time increases. There is a reverse flow in the y-wise velocity profile, and overshoot in the x-wise velocity and enthalpy profiles in the saddle point region, which increase as injection and wall temperature increase. Location of the dividing streamline increases as injection increases, but as the wall temperature and time increase, it decreases.

Magnetic Field Due to the Self-Gravity-Induced Electric Polarization of a Rotating Massive Body

he notion of the gravity-induced electric field has been applied to an entire self-gravitating massive body. The resulting electric polarization of the otherwise neutral body, when taken in conjunction with the latter's rotation, is shown to generate an axial-magnetic field of the right type and order of magnitude for certain astrophysical objects. In the present treatment the electric polarization is calculated in the ion-continuum Thomas-Fermi approximation while the electrodynamics of the continuous medium is treated in the nonrelativistic approximation.

A parametric differentiation version with finite-difference scheme applicable to a class of problems in boundary layer flow with massive blowing

A numerical procedure, based on the parametric differentiation and implicit finite difference scheme, has been developed for a class of problems in the boundary-layer theory for saddle-point regions. Here, the results are presented for the case of a three-dimensional stagnation-point flow with massive blowing. The method compares very well with other methods for particular cases (zero or small mass blowing). Results emphasize that the present numerical procedure is well suited for the solution of saddle-point flows with massive blowing, which could not be solved by other methods.

Parallelization of an unstructured data based cell centre finite volume code, HIFUN-3D

This work describes the parallelization of High Resolution flow solver on unstructured meshes, HIFUN-3D, an unstructured data based finite volume solver for 3-D Euler equations. For mesh partitioning, we use METIS, a software based on multilevel graph partitioning. The unstructured graph used for partitioning is associated with weights both on its vertices and edges. The data residing on every processor is split into four layers. Such a novel procedure of handling data helps in maintaining the effectiveness of the serial code. The communication of data across the processors is achieved by explicit message passing using the standard blocking mode feature of Message Passing Interface (MPI). The parallel code is tested on PACE++128 available in CFD Center

Nonsingular cosmological models: the massive scalar field case

The nonminimal coupling of a massive self-interacting scalar field with a gravitational field is studied. Spontaneous symmetry breaking occurs in the open universe even when the sign on the mass term is positive. In contrast to grand unified theories, symmetry breakdown is more important for the early universe and it is restored only in the limit of an infinite expansion. Symmetry breakdown is shown to occur in flat and closed universes when the mass term carries a wrong sign. The model has a naturally defined effective gravitational coupling coefficient which is rendered time-dependent due to the novel symmetry breakdown. It changes sign below a critical value of the cosmic scale factor indicating the onset of a repulsive field. The presence of the mass term severely alters the behaviour of ordinary matter and radiation in the early universe. The total energy density becomes negative in a certain domain. These features make possible a nonsingular cosm

Hypersonic stagnation‐point boundary layers with massive blowing in the presence of a magnetic field

The effect of massive blowing rates on the steady laminar hypersonic boundary‐layer flow of an electrically conducting fluid in the stagnation region of an axisymmetric body with an applied magnetic field has been studied. The governing equations have been solved numerically by combining the implicit finite‐difference scheme with the quasi‐linearization technique. It is observed that the effect of massive blowing rates is to remove the viscous layer away from the boundary, whereas the effect of the magnetic field is just the opposite. It is also found that the velocity overshoot increases with blowing rates and also with magnetic field. The effect of the variation of the density‐viscosity product across the boundary layer is strong only when the blowing rate is small, but for the massive blowing rate the effect is negligible.

Infrared divergences and confinement in massive (2+1)-dimensional QED

In (2+1)-dimensional quantum electrodynamics with massless photons and massive matter fields, it is shown that the mass renormalization of the latter is infrared divergent at one loop. This result remains unchanged at two loops. A simple argument based on a similar divergence of the Coulomb potential leads us to conjecture that charged states are not observable in this model. This argument holds in 1+1 dimensions also.

An automatic parallelization framework for multicomputers

Several researchers have looked into various issues related to automatic parallelization of sequential programs for multicomputers. But there is a need for a coherent framework which encompasses all these issues. In this paper we present a such a framework which takes best advantage of the multicomputer architecture. We resort to tiling transformation for iteration space partitioning and propose a scheme of automatic data partitioning and dynamic data distribution. We have tried a simple implementation of our scheme on a transputer based multicomputer [1] and the results are encouraging.

New Massive Gravity and AdS(4) Counterterms

We show that the recently proposed Dirac-Born-Infeld extension of new massive gravity emerges naturally as a counterterm in four-dimensional anti-de Sitter space (AdS(4)). The resulting on-shell Euclidean action is independent of the cutoff at zero temperature. We also find that the same choice of counterterm gives the usual area law for the AdS(4) Schwarzschild black hole entropy in a cutoff-independent manner. The parameter values of the resulting counterterm action correspond to a c = 0 theory in the context of the duality between AdS(3) gravity and two-dimensional conformal field theory. We rewrite this theory in terms of the gauge field that is used to recast 3D gravity as a Chern-Simons theory.

Unsteady laminar compressible swirling flow with massive blowing

THE study of swirling boundary layers is of considerable importance in many rotodynamic machines such as rockets, jet engines, swirl generators, swirl atomizers, arc heaters, etc. For example, the introduction of swirl in a flow acceleration device such as a nozzle in a rocket engine promises efficient mass flow control. In nuclear rockets, swirl is used to retain the uranium atoms in the rocket chamber. With these applications in mind, Back1 and Muthanna and Nath2 have obtained the similarity solutions for a low-speed three-dimensional steady laminar compressible boundary layer with swirl inside an axisymmetric surface of variable cross section. The aim of the present analysis is to study the effect of massive blowing rates on the unsteady laminar swirling compressible boundary-layer flow of an axisymmetric body of arbitrary cross section when the freestream velocity and blowing rate vary with time. The type of swirl considered here is that of a free vortex superimposed on the longitudinal flow of a compressible fluid with variable properties. The analysis is applicable to external flow over a body as well as internal flow along a surface. For the case of external flow, strong blowing can have significant use in cooling the surface of hypervelocity vehicles, particularly when ablation occurs under large aerodynamic or radiative heating, but there may not be such an important application of strong blowing in the case of internal flow. The governing partial differential equations have been solved numerically using an implicit finite difference scheme with a quasilinearization technique.3 High temperature gas effects, such as radiation, dissociation, and ionization, etc., are not investigated. The nomenclature is usually that of Ref. 4 and is listed in the full paper.

Online Unsupervised Pattern Discovery in Speech Using Parallelization

Segmental dynamic time warping (DTW) has been demonstrated to be a useful technique for finding acoustic similarity scores between segments of two speech utterances. Due to its high computational requirements, it had to be computed in an offline manner, limiting the applications of the technique. In this paper, we present results of parallelization of this task by distributing the workload in either a static or dynamic way on an 8-processor cluster and discuss the trade-offs among different distribution schemes. We show that online unsupervised pattern discovery using segmental DTW is plausible with as low as 8 processors. This brings the task within reach of today's general purpose multi-core servers. We also show results on a 32-processor system, and discuss factors affecting scalability of our methods.

Integrated parallelization of computations and visualization for large-scale applications

Critical applications like cyclone tracking and earthquake modeling require simultaneous high-performance simulations and online visualization for timely analysis. Faster simulations and simultaneous visualization enable scientists provide real-time guidance to decision makers. In this work, we have developed an integrated user-driven and automated steering framework that simultaneously performs numerical simulations and efficient online remote visualization of critical weather applications in resource-constrained environments. It considers application dynamics like the criticality of the application and resource dynamics like the storage space, network bandwidth and available number of processors to adapt various application and resource parameters like simulation resolution, simulation rate and the frequency of visualization. We formulate the problem of finding an optimal set of simulation parameters as a linear programming problem. This leads to 30% higher simulation rate and 25-50% lesser storage consumption than a naive greedy approach. The framework also provides the user control over various application parameters like region of interest and simulation resolution. We have also devised an adaptive algorithm to reduce the lag between the simulation and visualization times. Using experiments with different network bandwidths, we find that our adaptive algorithm is able to reduce lag as well as visualize the most representative frames.