Towards an algorithmic skeleton framework for programming the Intel R Xeon PhiTM processor

The Intel R Xeon PhiTM is the first processor based on Intel’s MIC (Many Integrated Cores) architecture. It is a co-processor specially tailored for data-parallel computations, whose basic architectural design is similar to the ones of GPUs (Graphics Processing Units), leveraging the use of many integrated low computational cores to perform parallel computations. The main novelty of the MIC architecture, relatively to GPUs, is its compatibility with the Intel x86 architecture. This enables the use of many of the tools commonly available for the parallel programming of x86-based architectures, which may lead to a smaller learning curve. However, programming the Xeon Phi still entails aspects intrinsic to accelerator-based computing, in general, and to the MIC architecture, in particular. In this thesis we advocate the use of algorithmic skeletons for programming the Xeon Phi. Algorithmic skeletons abstract the complexity inherent to parallel programming, hiding details such as resource management, parallel decomposition, inter-execution flow communication, thus removing these concerns from the programmer’s mind. In this context, the goal of the thesis is to lay the foundations for the development of a simple but powerful and efficient skeleton framework for the programming of the Xeon Phi processor. For this purpose we build upon Marrow, an existing framework for the orchestration of OpenCLTM computations in multi-GPU and CPU environments. We extend Marrow to execute both OpenCL and C++ parallel computations on the Xeon Phi. We evaluate the newly developed framework, several well-known benchmarks, like Saxpy and N-Body, will be used to compare, not only its performance to the existing framework when executing on the co-processor, but also to assess the performance on the Xeon Phi versus a multi-GPU environment.

Creation of an Age-Adjusted, Dual-Energy X-ray Absorptiometry-Derived Trabecular Bone Score Curve for the Lumbar Spine in Non-Hispanic US White Women.

The trabecular bone score (TBS, Med-Imaps, Pessac, France) is an index of bone microarchitecture texture extracted from anteroposterior dual-energy X-ray absorptiometry images of the spine. Previous studies have documented the ability of TBS of the spine to differentiate between women with and without fractures among age- and areal bone mineral density (aBMD)-matched controls, as well as to predict future fractures. In this cross-sectional analysis of data collected from 3 geographically dispersed facilities in the United States, we investigated age-related changes in the microarchitecture of lumbar vertebrae as assessed by TBS in a cohort of non-Hispanic US white American women. All subjects were 30 yr of age and older and had an L1-L4aBMDZ-score within ±2 SD of the population mean. Individuals were excluded if they had fractures, were on any osteoporosis treatment, or had any illness that would be expected to impact bone metabolism. All data were extracted from Prodigy dual-energy X-ray absorptiometry devices (GE-Lunar, Madison, WI). Cross-calibrations between the 3 participating centers were performed for TBS and aBMD. aBMD and TBS were evaluated for spine L1-L4 but also for all other possible vertebral combinations. To validate the cohort, a comparison between the aBMD normative data of our cohort and US non-Hispanic white Lunar data provided by the manufacturer was performed. A database of 619 non-Hispanic US white women, ages 30-90 yr, was created. aBMD normative data obtained from this cohort were not statistically different from the non-Hispanic US white Lunar normative data provided by the manufacturer (p = 0.30). This outcome thereby indirectly validates our cohort. TBS values at L1-L4 were weakly inversely correlated with body mass index (r = -0.17) and weight (r = -0.16) and not correlated with height. TBS values for all lumbar vertebral combinations decreased significantly with age. There was a linear decrease of 16.0% (-2.47 T-score) in TBS at L1-L4 between 45 and 90 yr of age (vs. -2.34 for aBMD). Microarchitectural loss rate increased after age 65 by 50% (-0.004 to -0.006). Similar results were obtained for other combinations of lumbar vertebra. TBS, an index of bone microarchitectural texture, decreases with advancing age in non-Hispanic US white women. Little change in TBS is observed between ages 30 and 45. Thereafter, a progressive decrease is observed with advancing age. The changes we observed in these American women are similar to that previously reported for a French population of white women (r(2) > 0.99). This reference database will facilitate the use of TBS to assess bone microarchitectural deterioration in clinical practice.

Reducing the learning curve for the treatment of morphoeic (sclerosing) basal cell carcinoma of the face.

The treatment of morphoeic (or sclerosing) basal cell carcinoma (mBCC) of the face is associated with high rates of incomplete excision and recurrence. A principal risk factor for incomplete resection is the grade of surgeon. We did a prospective, randomised study of 40 consecutive patients with mBCC of the face. The extent of the tumour was assessed under standard conditions by consultant surgeons and compared with assessments by resident surgeons with the help of the Varioscope, a combination of microscope and loupe glasses with strong illumination and a maximal magnification of 7x. The data from a former retrospective study of all excisions of mBCC of the face during a five-year period at the hospital served as control. Residents with the support of the Varioscope achieved a rate of incomplete excisions similar to that of consultants under standard conditions. There was a significant reduction of the rate of incomplete resections by resident surgeons thanks to high magnification and good lighting (p=0.02). High magnification and good lighting were useful in learning how to recognise skin changes associated with mBCC of the face and achieving a low rate of incomplete excisions.

Generation and validation of a normative, age-specific reference curve for lumbar spine trabecular bone score (TBS) in French women.

Age-related changes in lumbar vertebral microarchitecture are evaluated, as assessed by trabecular bone score (TBS), in a cohort of 5,942 French women. The magnitude of TBS decline between 45 and 85 years of age is piecewise linear in the spine and averaged 14.5 %. TBS decline rate increases after 65 years by 50 %. INTRODUCTION: This study aimed to evaluate age-related changes in lumbar vertebral microarchitecture, as assessed by TBS, in a cohort of French women aged 45-85 years. METHODS: An all-comers cohort of French Caucasian women was selected from two clinical centers. Data obtained from these centers were cross-calibrated for TBS and bone mineral density (BMD). BMD and TBS were evaluated at L1-L4 and for all lumbar vertebrae combined using GE-Lunar Prodigy densitometer images. Weight, height, and body mass index (BMI) also were determined. To validate our all-comers cohort, the BMD normative data of our cohort and French Prodigy data were compared. RESULTS: A cohort of 5,942 French women aged 45 to 85 years was created. Dual-energy X-ray absorptiometry normative data obtained for BMD from this cohort were not significantly different from French prodigy normative data (p = 0.15). TBS values at L1-L4 were poorly correlated with BMI (r = -0.17) and weight (r = -0.14) and not correlated with height. TBS values obtained for all lumbar vertebra combined (L1, L2, L3, L4) decreased with age. The magnitude of TBS decline at L1-L4 between 45 and 85 years of age was piecewise linear in the spine and averaged 14.5 %, but this rate increased after 65 years by 50 %. Similar results were obtained for other region of interest in the lumbar spine. As opposed to BMD, TBS was not affected by spinal osteoarthrosis. CONCLUSION: The age-specific reference curve for TBS generated here could therefore be used to help clinicians to improve osteoporosis patient management and to monitor microarchitectural changes related to treatment or other diseases in routine clinical practice.

Study of Type 'A' Asphaltic Concrete Surface Mix Design, R-232, 1968

In some asphaltic concrete mixes asphalt absorption in field mixes is difficult to predict by the routine mix design tests presently being used. Latent or slow absorption in hot mixes is hard to compensate for in field control due to aggregate gradations being near maximum density. If critical asphalt need could be changed by increasing voids in the mineral aggregate so that more freedom could be exercised in compensating for the absorption, this may aid in design. The voids in the mineral aggregate can be related to composite gradation of total aggregate in a mixture, i.e. if a composite gradation of aggregate is finer than that of maximum density curve, the V.M.A. will be greater than that of a mix of maximum density. The typical gradation of Iowa Type 'A' mixes is finer than a gradation which is near the centerline of the specification at sieves larger than the No. 30 and coarser at the lower sieve sizes. The mixes of the typical gradation will have higher V.M.A. than those of the near centerline mixes. By studying properties of the mixes of the typical gradation and comparing them with those of the mixes of maximum density, it may aid in the modification and simplification of our present testing methods and specification requirements while still maintaining control of quality of the mix by controlling voids, stability, gradation and asphalt content.

ASTRAL-R score predicts non-recanalisation after intravenous thrombolysis in acute ischaemic stroke.

Intravenous thrombolysis (IVT) as treatment in acute ischaemic strokes may be insufficient to achieve recanalisation in certain patients. Predicting probability of non-recanalisation after IVT may have the potential to influence patient selection to more aggressive management strategies. We aimed at deriving and internally validating a predictive score for post-thrombolytic non-recanalisation, using clinical and radiological variables. In thrombolysis registries from four Swiss academic stroke centres (Lausanne, Bern, Basel and Geneva), patients were selected with large arterial occlusion on acute imaging and with repeated arterial assessment at 24 hours. Based on a logistic regression analysis, an integer-based score for each covariate of the fitted multivariate model was generated. Performance of integer-based predictive model was assessed by bootstrapping available data and cross validation (delete-d method). In 599 thrombolysed strokes, five variables were identified as independent predictors of absence of recanalisation: Acute glucose > 7 mmol/l (A), significant extracranial vessel STenosis (ST), decreased Range of visual fields (R), large Arterial occlusion (A) and decreased Level of consciousness (L). All variables were weighted 1, except for (L) which obtained 2 points based on β-coefficients on the logistic scale. ASTRAL-R scores 0, 3 and 6 corresponded to non-recanalisation probabilities of 18, 44 and 74 % respectively. Predictive ability showed AUC of 0.66 (95 %CI, 0.61-0.70) when using bootstrap and 0.66 (0.63-0.68) when using delete-d cross validation. In conclusion, the 5-item ASTRAL-R score moderately predicts non-recanalisation at 24 hours in thrombolysed ischaemic strokes. If its performance can be confirmed by external validation and its clinical usefulness can be proven, the score may influence patient selection for more aggressive revascularisation strategies in routine clinical practice.

Dose-response curve to soil applied herbicides and susceptibility evaluation of different amaranthus species using model identity

Greenhouse studies were conducted in 2008-2009 with the objective of adjusting dose-response curves of the main soil-applied herbicides currently used in cotton for the control of Amaranthus viridis, A. hybridus, A. spinosus, A. lividus, as well as comparing susceptibility among different species, using the identity test models. Thirty six individual experiments were simultaneously carried out in greenhouse, in a sandy clay loam soil (21% clay, 2.36% OM) combining increasing doses of the herbicides alachlor, clomazone, diuron, oxyfluorfen, pendimethalin, prometryn, S-metolachlor, and trifluralin applied to each species. Dose-response curves were adjusted for visual weed control at 28 days after herbicide application and doses required for 80% (C80) and 95% (C95) control were calculated. All herbicides, except clomazone and trifluralin, provided efficient control of most Amaranthus species, but substantial differences in susceptibility to herbicides were found. In general, A. lividus was the least sensitive species, whereas A. spinosus demonstrated the highest sensitivity to herbicides. Alachlor, diuron, oxyfluorfen, pendimethalin, S-metolachlor, and prometryn are efficient alternatives to control Amaranthus spp. in a range of doses that are currently lower than those recommended to cotton.

Stability of relative oxygen pulse curve during repeated maximal cardiopulmonary testing in professional soccer players

During cardiopulmonary exercise testing (CPET), stroke volume can be indirectly assessed by O2 pulse profile. However, for a valid interpretation, the stability of this variable over time should be known. The objective was to analyze the stability of the O2 pulse curve relative to body mass in elite athletes. VO2, heart rate (HR), and relative O2 pulse were compared at every 10% of the running time in two maximal CPETs, from 2005 to 2010, of 49 soccer players. Maximal values of VO2 (63.4 ± 0.9 vs 63.5 ± 0.9 mL O2•kg-1•min-1), HR (190 ± 1 vs188 ± 1 bpm) and relative O2 pulse (32.9 ± 0.6 vs 32.6 ± 0.6 mL O2•beat-1•kg-1) were similar for the two CPETs (P > 0.05), while the final treadmill velocity increased from 18.5 ± 0.9 to 18.9 ± 1.0 km/h (P < 0.01). Relative O2 pulse increased linearly and similarly in both evaluations (r² = 0.64 and 0.63) up to 90% of the running time. Between 90 and 100% of the running time, the values were less stable, with up to 50% of the players showing a tendency to a plateau in the relative O2 pulse. In young healthy men in good to excellent aerobic condition, the morphology of the relative O2 pulse curve is consistent up to close to the peak effort for a CPET repeated within a 1-year period. No increase in relative O2pulse at peak effort could represent a physiologic stroke volume limitation in these athletes.

The effects of R(+)-lipoic acid supplementation on regulation of human skeletal muscle pyruvate dehydrogenase

This thesis investigated whole body glucose disposal and the adaptive changes in skeletal muscle carbohydrate metabolism following 28 d of supplementation with 1000 mg R(+)-lipoic acid in young sedentary males (age, 22.1 ± 0.67 yr, body mass, 78.7 ± 10.3 kg, n=9). In certain individuals, lipoic acid decreased the 180-min area under the glucose concentration and insulin concentration curve during an oral glucose tolerance test (OGTT) (n=4). In the same individuals, lipoic acid supplementation decreased pyruvate dehydrogenase kinase activity (PDK) (0.09 ± 0.024 min"^ vs. 0.137 ± 0.023 min'\ n=4). The fasting levels of the activated form of pyruvate dehydrogenase (PDHa) were decreased following lipoic acid (0.42 ± 0.13 mmol-min'kg'^ vs. 0.82 ± 0.32 mmolrnin'^kg"\ n=4), yet increased to a greater extent during the OGTT (1.21 ± 0.34 mmol-min'kg"' vs. 0.81 ±0.13 mmolmin"'kg'\ n=4) following hpoic acid supplementation. No changes were demonstrated in the remaining subjects (n=5). It was concluded that improved glucose clearance during an OGTT following lipoic acid supplementation is assisted by increased muscle glucose oxidation through increased PDHa activation and decreased PDK activity in certain individuals.

Studies on scattering and thermodynamics of black holes in f(R) theory and Einstein's theory

One of the interesting consequences of Einstein's General Theory of Relativity is the black hole solutions. Until the observation made by Hawking in 1970s, it was believed that black holes are perfectly black. The General Theory of Relativity says that black holes are objects which absorb both matter and radiation crossing the event horizon. The event horizon is a surface through which even light is not able to escape. It acts as a one sided membrane that allows the passage of particles only in one direction i.e. towards the center of black holes. All the particles that are absorbed by black hole increases the mass of the black hole and thus the size of event horizon also increases. Hawking showed in 1970s that when applying quantum mechanical laws to black holes they are not perfectly black but they can emit radiation. Thus the black hole can have temperature known as Hawking temperature. In the thesis we have studied some aspects of black holes in f(R) theory of gravity and Einstein's General Theory of Relativity. The scattering of scalar field in this background space time studied in the first chapter shows that the extended black hole will scatter scalar waves and have a scattering cross section and applying tunneling mechanism we have obtained the Hawking temperature of this black hole. In the following chapter we have investigated the quasinormal properties of the extended black hole. We have studied the electromagnetic and scalar perturbations in this space-time and find that the black hole frequencies are complex and show exponential damping indicating the black hole is stable against the perturbations. In the present study we show that not only the black holes exist in modified gravities but also they have similar properties of black hole space times in General Theory of Relativity. 2 + 1 black holes or three dimensional black holes are simplified examples of more complicated four dimensional black holes. Thus these models of black holes are known as toy models of black holes in four dimensional black holes in General theory of Relativity. We have studied some properties of these types of black holes in Einstein model (General Theory of Relativity). A three dimensional black hole known as MSW is taken for our study. The thermodynamics and spectroscopy of MSW black hole are studied and obtained the area spectrum which is equispaced and different thermo dynamical properties are studied. The Dirac perturbation of this three dimensional black hole is studied and the resulting quasinormal spectrum of this three dimensional black hole is obtained. The different quasinormal frequencies are tabulated in tables and these values show an exponential damping of oscillations indicating the black hole is stable against the mass less Dirac perturbation. In General Theory of Relativity almost all solutions contain singularities. The cosmological solution and different black hole solutions of Einstein's field equation contain singularities. The regular black hole solutions are those which are solutions of Einstein's equation and have no singularity at the origin. These solutions possess event horizon but have no central singularity. Such a solution was first put forward by Bardeen. Hayward proposed a similar regular black hole solution. We have studied the thermodynamics and spectroscopy of Hay-ward regular black holes. We have also obtained the different thermodynamic properties and the area spectrum. The area spectrum is a function of the horizon radius. The entropy-heat capacity curve has a discontinuity at some value of entropy showing a phase transition.

The crime kuznets curve

We document the existence of a Crime Kuznets Curve in US states since the 1970s. As income levels have risen, crime has followed an inverted U-shaped pattern, first increasing and then dropping. The Crime Kuznets Curve is not explained by income inequality. In fact, we show that during the sample period inequality has risen monotonically with income, ruling out the traditional Kuznets Curve. Our finding is robust to adding a large set of controls that are used in the literature to explain the incidence of crime, as well as to controlling for state and year fixed effects. The Curve is also revealed in nonparametric specifications. The Crime Kuznets Curve exists for property crime and for some categories of violent crime.

Fitting models for the joint action of two drugs using SAS(R)

Combinations of drugs are increasingly being used for a wide variety of diseases and conditions. A pre-clinical study may allow the investigation of the response at a large number of dose combinations. In determining the response to a drug combination, interest may lie in seeking evidence of synergism, in which the joint action is greater than the actions of the individual drugs, or of antagonism, in which it is less. Two well-known response surface models representing no interaction are Loewe additivity and Bliss independence, and Loewe or Bliss synergism or antagonism is defined relative to these. We illustrate an approach to fitting these models for the case in which the marginal single drug dose-response relationships are represented by four-parameter logistic curves with common upper and lower limits, and where the response variable is normally distributed with a common variance about the dose-response curve. When the dose-response curves are not parallel, the relative potency of the two drugs varies according to the magnitude of the desired effect and the models for Loewe additivity and synergism/antagonism cannot be explicitly expressed. We present an iterative approach to fitting these models without the assumption of parallel dose-response curves. A goodness-of-fit test based on residuals is also described. Implementation using the SAS NLIN procedure is illustrated using data from a pre-clinical study. Copyright © 2007 John Wiley & Sons, Ltd.

Switching curve control of functional electrical stimulation assisted rowing exercise in paraplegia

An indoor rowing machine has been modified for functional electrical stimulation (FES) assisted rowing exercise in paraplegia. To perform the rowing manoeuvre successfully, however, the voluntarily controlled upper body movements must be co-ordinated with the movements of the electrically stimulated paralysed legs. To achieve such co-ordination, an automatic FES controller was developed that employs two levels of hierarchy. At the upper level, a finite state controller identifies the state or phase of the rowing cycle and activates the appropriate lower-level controller, in which electrical stimulation to the paralysed leg muscles is applied with reference to switching curves representing the desired seat velocity as a function of the seat position. In a pilot study, the hierarchical control of FES rowing was shown to be intuitive, reliable and easy to use. Compared with open-loop control of stimulation, all three variants of the closed-loop switching curve controllers used less muscle stimulation per rowing cycle (73% of the open-loop control on average). Further, the closed-loop controller that used switching curves derived from normal rowing kinematics used the lowest muscle stimulation (65% of the open-loop control) and was the most convenient to use for the client.

The Carnegie curve

The Earth’s fair weather atmospheric electric field shows, in clean air, an average daily variation which follows universal time, globally independent of the measurement position. This single diurnal cycle variation (maximum around 19UT and minimum around 03UT) is widely known as the Carnegie curve, after the geophysical survey vessel of the Carnegie Institution of Washington on which the original measurement campaigns demonstrating the universal time variation were undertaken. The Carnegie curve’s enduring importance is in providing a reference variation against which atmospheric electricity measurements are still compared; it is believed to originate from regular daily variations in atmospheric electrification associated with the different global disturbed weather regions. Details of the instrumentation, measurement principles and data obtained on the Carnegie’s seventh and final cruise are reviewed here, also deriving new harmonic coefficients allowing calculation of the Carnegie curve for different seasons. The additional harmonic analysis now identifies changes in the phasing of the maximum and minimum in the Carnegie curve, which shows a systematic seasonal variation, linked to the solstices and equinoxes, respectively.

Computing the Cassels–Tate pairing on the 3-Selmer group of an elliptic curve

We extend the method of Cassels for computing the Cassels-Tate pairing on the 2-Selmer group of an elliptic curve, to the case of 3-Selmer groups. This requires significant modifications to both the local and global parts of the calculation. Our method is practical in sufficiently small examples, and can be used to improve the upper bound for the rank of an elliptic curve obtained by 3-descent.