Three parallel programming paradigms: Comparisons on an archetypal PDE computation

Three paradigms for distributed-memory parallel computation that free the application programmer from the details of message passing are compared for an archetypal structured scientific computation -- a nonlinear, structured-grid partial differential equation boundary value problem -- using the same algorithm on the same hardware. All of the paradigms -- parallel languages represented by the Portland Group's HPF, (semi-)automated serial-to-parallel source-to-source translation represented by CAP-Tools from the University of Greenwich, and parallel libraries represented by Argonne's PETSc -- are found to be easy to use for this problem class, and all are reasonably effective in exploiting concurrency after a short learning curve. The level of involvement required by the application programmer under any paradigm includes specification of the data partitioning, corresponding to a geometrically simple decomposition of the domain of the PDE. Programming in SPMD style for the PETSc library requires writing only the routines that discretize the PDE and its Jacobian, managing subdomain-to-processor mappings (affine global-to-local index mappings), and interfacing to library solver routines. Programming for HPF requires a complete sequential implementation of the same algorithm as a starting point, introduction of concurrency through subdomain blocking (a task similar to the index mapping), and modest experimentation with rewriting loops to elucidate to the compiler the latent concurrency. Programming with CAPTools involves feeding the same sequential implementation to the CAPTools interactive parallelization system, and guiding the source-to-source code transformation by responding to various queries about quantities knowable only at runtime. Results representative of "the state of the practice" for a scaled sequence of structured grid problems are given on three of the most important contemporary high-performance platforms: the IBM SP, the SGI Origin 2000, and the CRAYY T3E.

A new parameterization of Johnson's SB distribution with application to fitting forest tree diameter data

The SB distributional model of Johnson's 1949 paper was introduced by a transformation to normality, that is, z ~ N(0, 1), consisting of a linear scaling to the range (0, 1), a logit transformation, and an affine transformation, z = γ + δu. The model, in its original parameterization, has often been used in forest diameter distribution modelling. In this paper, we define the SB distribution in terms of the inverse transformation from normality, including an initial linear scaling transformation, u = γ′ + δ′z (δ′ = 1/δ and γ′ = �γ/δ). The SB model in terms of the new parameterization is derived, and maximum likelihood estimation schema are presented for both model parameterizations. The statistical properties of the two alternative parameterizations are compared empirically on 20 data sets of diameter distributions of Changbai larch (Larix olgensis Henry). The new parameterization is shown to be statistically better than Johnson's original parameterization for the data sets considered here.

Kinematics and EMG analysis of expert pole vaulter's lower limb during take off phase (P218)

Introduction: The critical phase, in jumping events in track and field, appears to be between touchdown and take-off. Since obvious similarities exist between the take off phase in both long jump and pole vault, numerous 3D kinematics and electromyographic studies have only looked at long jump. Currently there are few detailed kinematics electromyographic data on the pole vault take-off phase. The aim of this study was therefore to characterise kinematics and electromyographic variables during the take-off phase to provide a better understanding of this phase in pole vaulting and its role in performance outcome. Material and methods: Six pole-vaulters took part in the study. Kinematics data were captured with retro reflective markers fixed on the body. Hip, knee and ankle angle were calculated. Differential bipolar surface electrodes were placed on the following muscles of the take-off leg: tibialis anterior, lateral gastrocnemius, vastus lateralis, rectus femoris, bicep femoris and gluteus maximus. EMG activity was synchronously acquired with the kinematic data. EMG data were rectified and smoothed using a second order low pass Butterworth Bidirectional filter (resulting in a 4th order filter) with a cut-off frequency of 14 Hz. Results: Evolution of hip, knee and ankle angle show no significant differences during the last step before touchdown, the take-off phase and the beginning of fly phase. Meanwhile, strong differences in EMG signal are noted inter and intra pole vaulter. However for a same subject the EMG activities seem to converge to some phase locked point. Discussion: All pole vaulters have approximately the same visible coordination This coordination reflects a different muscular control among pole vaulters but also for a considered pole vaulter. These phase locked point could be considered as invariant of motor control i.e. a prerequisite for a normal sequence of the movement and performance realization.