GPU accelerated algorithms for computing matrix function vector products with applications to exponential integrators and fractional diffusion

The efficient computation of matrix function vector products has become an important area of research in recent times, driven in particular by two important applications: the numerical solution of fractional partial differential equations and the integration of large systems of ordinary differential equations. In this work we consider a problem that combines these two applications, in the form of a numerical solution algorithm for fractional reaction diffusion equations that after spatial discretisation, is advanced in time using the exponential Euler method. We focus on the efficient implementation of the algorithm on Graphics Processing Units (GPU), as we wish to make use of the increased computational power available with this hardware. We compute the matrix function vector products using the contour integration method in [N. Hale, N. Higham, and L. Trefethen. Computing Aα, log(A), and related matrix functions by contour integrals. SIAM J. Numer. Anal., 46(5):2505–2523, 2008]. Multiple levels of preconditioning are applied to reduce the GPU memory footprint and to further accelerate convergence. We also derive an error bound for the convergence of the contour integral method that allows us to pre-determine the appropriate number of quadrature points. Results are presented that demonstrate the effectiveness of the method for large two-dimensional problems, showing a speedup of more than an order of magnitude compared to a CPU-only implementation.

Prediction of maximal surface electromyographically based voluntary contractions of erector spinae muscles from sonographic measurements during isometric contractions

Objectives Currently, there are no studies combining electromyography (EMG) and sonography to estimate the absolute and relative strength values of erector spinae (ES) muscles in healthy individuals. The purpose of this study was to establish whether the maximum voluntary contraction (MVC) of the ES during isometric contractions could be predicted from the changes in surface EMG as well as in fiber pennation and thickness as measured by sonography. Methods Thirty healthy adults performed 3 isometric extensions at 45° from the vertical to calculate the MVC force. Contractions at 33% and 100% of the MVC force were then used during sonographic and EMG recordings. These measurements were used to observe the architecture and function of the muscles during contraction. Statistical analysis was performed using bivariate regression and regression equations. Results The slope for each regression equation was statistically significant (P < .001) with R2 values of 0.837 and 0.986 for the right and left ES, respectively. The standard error estimate between the sonographic measurements and the regression-estimated pennation angles for the right and left ES were 0.10 and 0.02, respectively. Conclusions Erector spinae muscle activation can be predicted from the changes in fiber pennation during isometric contractions at 33% and 100% of the MVC force. These findings could be essential for developing a regression equation that could estimate the level of muscle activation from changes in the muscle architecture.

The moderating role of task characteristics in determining responses to a stressful work simulation

Two experimental studies were conducted to examine whether the stress-buffering effects of behavioral control on work task responses varied as a function of procedural information. Study 1 manipulated low and high levels of task demands, behavioral control, and procedural information for 128 introductory psychology students completing an in-basket activity. ANOVA procedures revealed a significant three-way interaction among these variables in the prediction of subjective task performance and task satisfaction. It was found that procedural information buffered the negative effects of task demands on ratings of performance and satisfaction only under conditions of low behavioral control. This pattern of results suggests that procedural information may have a compensatory effect when the work environment is characterized by a combination of high task demands and low behavioral control. Study 2 (N=256) utilized simple and complex versions of the in-basket activity to examine the extent to which the interactive relationship among task demands, behavioral control, and procedural information varied as a function of task complexity. There was further support for the stress-buffering role of procedural information on work task responses under conditions of low behavioral control. This effect was, however, only present when the in-basket activity was characterized by high task complexity, suggesting that the interactive relationship among these variables may depend on the type of tasks performed at work.

Suitability of diodes for point dose measurements in IMRT/VMAT beams

This study investigated a potential source of inaccuracy for diode measurements in modulated beams; the effect of diode housing asymmetry on measurement results. The possible effects of diode housing asymmetry on the measurement of steep dose gradients were evaluated by measuring 5x5 cm2 beam profiles, with three cylindrical diodes and two commonly used ionization chambers, with each dosimeter positioned in a 3D scanning water tank with its stem perpendicular to the beam axis (horizontal) and parallel to the direction of scanning. The resulting profiles were used to compare the penumbrae measured with the diode stem pointing into (equivalent to a “stem-first” setup) and out of the field (equivalent to a “stem-last” setup) in order to evaluate the effects of dosimeter alignment and thereby identify the effects of dosimeter asymmetry. The stem-first and stem-last orientations resulted in differences of up to 0.2 mm in the measured 20-80% penumbra widths and differences of up to 0.4 mm in the off axis position of the 90% isodose. These differences, which are smaller than previously reported for older model dosimeters, were apparent in the profile results for both diodes and small volume ionization chambers. As an extension to this study, the practical use of all five dosimeters was exemplified by measuring point doses in IMRT test beams. These measurements showed good agreement (within 2%) between the diodes and the small volume ionization chamber, with all of these dosimeters being able to identify a region 3% under-dosage which was not identified by a larger volume (6 mm diameter) ionization chamber. The results of this work should help to remove some of the barriers to the use of diodes for modulated radiotherapy dosimetry in the future.

The association between ambient temperature and children's lung function in Baotou, China

The objective of this study is to examine the association between ambient temperature and children’s lung function in Baotou, China. We recruited 315 children (8–12 years) from Baotou, China in the spring of 2004, 2005, and 2006. They performed three successive forced expiratory measurements three times daily (morning, noon, and evening) for about 5 weeks. The highest peak expiratory flow (PEF) was recorded for each session. Daily data on ambient temperature, relative humidity, and air pollution were monitored during the same period. Mixed models with a distributed lag structure were used to examine the effects of temperature on lung function while adjusting for individual characteristics and environmental factors. Low temperatures were significantly associated with decreases in PEF. The effects lasted for lag 0–2 days. For all participants, the cumulative effect estimates (lag 0–2 days) were −1.44 (−1.93, −0.94) L/min, −1.39 (−1.92, −0.86) L/min, −1.40 (−1.97, −0.82) L/min, and −1.28 (−1.69, −0.88) L/min for morning, noon, evening, and daily mean PEF, respectively, associated with 1 °C decrease in daily mean temperature. Generally, the effects of temperature were slightly stronger in boys than in girls for noon, evening, and daily mean PEF, while the effects were stronger in girls for morning PEF. PM2.5 had joint effects with temperature on children’s PEF. Higher PM2.5 increased the impacts of low temperature. Low ambient temperatures are associated with lower lung function in children in Baotou, China. Preventive health policies will be required for protecting children from the cold weather.

Frequency response function based damage identification using principal component analysis and pattern recognition technique

Pattern recognition is a promising approach for the identification of structural damage using measured dynamic data. Much of the research on pattern recognition has employed artificial neural networks (ANNs) and genetic algorithms as systematic ways of matching pattern features. The selection of a damage-sensitive and noise-insensitive pattern feature is important for all structural damage identification methods. Accordingly, a neural networks-based damage detection method using frequency response function (FRF) data is presented in this paper. This method can effectively consider uncertainties of measured data from which training patterns are generated. The proposed method reduces the dimension of the initial FRF data and transforms it into new damage indices and employs an ANN method for the actual damage localization and quantification using recognized damage patterns from the algorithm. In civil engineering applications, the measurement of dynamic response under field conditions always contains noise components from environmental factors. In order to evaluate the performance of the proposed strategy with noise polluted data, noise contaminated measurements are also introduced to the proposed algorithm. ANNs with optimal architecture give minimum training and testing errors and provide precise damage detection results. In order to maximize damage detection results, the optimal architecture of ANN is identified by defining the number of hidden layers and the number of neurons per hidden layer by a trial and error method. In real testing, the number of measurement points and the measurement locations to obtain the structure response are critical for damage detection. Therefore, optimal sensor placement to improve damage identification is also investigated herein. A finite element model of a two storey framed structure is used to train the neural network. It shows accurate performance and gives low error with simulated and noise-contaminated data for single and multiple damage cases. As a result, the proposed method can be used for structural health monitoring and damage detection, particularly for cases where the measurement data is very large. Furthermore, it is suggested that an optimal ANN architecture can detect damage occurrence with good accuracy and can provide damage quantification with reasonable accuracy under varying levels of damage.

Effects of the sensorimotor training volume on sensorimotor function in patients following lower limb arthroplasty

Background Sensorimotor function is degraded in patients after lower limb arthroplasty. Sensorimotor training is thought to improve sensorimotor skills, however, the optimal training stimulus with regard to volume, frequency, duration, and intensity is still unknown. The aim of this study, therefore, was to firstly quantify the progression of sensorimotor function after total hip (THA) or knee (TKA) arthroplasty and, as second step, to evaluate effects of different sensorimotor training volumes. Methods 58 in-patients during their rehabilitation after THA or TKA participated in this prospective cohort study. Sensorimotor function was assessed using a test battery including measures of stabilization capacity, static balance, proprioception, and gait, along with a self-reported pain and function. All participants were randomly assigned to one of three intervention groups performing sensorimotor training two, four, or six times per week. Outcome measures were taken at three instances, at baseline (pre), after 1.5 weeks (mid) and at the conclusion of the 3 week program (post). Results All measurements showed significant improvements over time, with the exception of proprioception and static balance during quiet bipedal stance which showed no significant main effects for time or intervention. There was no significant effect of sensorimotor training volume on any of the outcome measures. Conclusion We were able to quantify improvements in measures of dynamic, but not static, sensorimotor function during the initial three weeks of rehabilitation following TKA/THA. Although sensorimotor improvements were independent of the training volume applied in the current study, long-term effects of sensorimotor training volume need to be investigated to optimize training stimulus recommendations.

Estimating cell diffusivity and cell proliferation rate by interpreting IncuCyte ZOOM™ assay data using the Fisher-Kolmogorov model

Background: Standard methods for quantifying IncuCyte ZOOM™ assays involve measurements that quantify how rapidly the initially-vacant area becomes re-colonised with cells as a function of time. Unfortunately, these measurements give no insight into the details of the cellular-level mechanisms acting to close the initially-vacant area. We provide an alternative method enabling us to quantify the role of cell motility and cell proliferation separately. To achieve this we calibrate standard data available from IncuCyte ZOOM™ images to the solution of the Fisher-Kolmogorov model. Results: The Fisher-Kolmogorov model is a reaction-diffusion equation that has been used to describe collective cell spreading driven by cell migration, characterised by a cell diffusivity, D, and carrying capacity limited proliferation with proliferation rate, λ, and carrying capacity density, K. By analysing temporal changes in cell density in several subregions located well-behind the initial position of the leading edge we estimate λ and K. Given these estimates, we then apply automatic leading edge detection algorithms to the images produced by the IncuCyte ZOOM™ assay and match this data with a numerical solution of the Fisher-Kolmogorov equation to provide an estimate of D. We demonstrate this method by applying it to interpret a suite of IncuCyte ZOOM™ assays using PC-3 prostate cancer cells and obtain estimates of D, λ and K. Comparing estimates of D, λ and K for a control assay with estimates of D, λ and K for assays where epidermal growth factor (EGF) is applied in varying concentrations confirms that EGF enhances the rate of scratch closure and that this stimulation is driven by an increase in D and λ, whereas K is relatively unaffected by EGF. Conclusions: Our approach for estimating D, λ and K from an IncuCyte ZOOM™ assay provides more detail about cellular-level behaviour than standard methods for analysing these assays. In particular, our approach can be used to quantify the balance of cell migration and cell proliferation and, as we demonstrate, allow us to quantify how the addition of growth factors affects these processes individually.

Complex permittivity measurements of RF plasma polymerized polyterpenol organic thin films employing split post dielectric resonator

In the fields of organic electronics and biotechnology, applications for organic polymer thin films fabricated using low-temperature non-equilibrium plasma techniques are gaining significant attention because of the physical and chemical stability of thin films and the low cost of production. Polymer thin films were fabricated from non-synthetic terpinen-4-ol using radiofrequency polymerization (13.56 MHz) on low loss dielectric substrates and their permittivity properties were ascertained to determine potential applications for these organic films. Real and imaginary parts of permittivity as a function of frequency were measured using the variable angle spectroscopic ellipsometer. The real part of permittivity (k) was found to be between 2.34 and 2.65 in the wavelength region of 400–1100 nm, indicating a potential low-k material. These permittivity values were confirmed at microwave frequencies. Dielectric properties of polyterpenol films were measured by means of split post dielectric resonators (SPDRs) operating at frequencies of 10 GHz and 20 GHz. Permittivity increased for samples deposited at higher RF energy – from 2.65 (25 W) to 2.83 (75 W) measured by a 20-GHz SPDR and from 2.32 (25 W) to 2.53 (100 W) obtained using a 10-GHz SPDR. The error in permittivity measurement was predominantly attributed to the uncertainty in film thickness measurement.

Intergenerational Perceptions and Conflicts in Multi-Age and Multigenerational Work Environments

One underappreciated consequence of the aging population phenomenon is that we are now experiencing what is arguably the most age-diverse workforce in modern history (Hanks & Icenogle, 2001; Newton, 2006; Toossi, 2004). As our workforce continues to age, shifts in the age demographic composition (i.e., the age diversity) of organizations and their subunits will become more apparent (Roth, Wegge, & Schmidt, 2007). Several factors have influenced and will continue to drive this trend. For example, in Western countries, younger people entering the workforce are more educated than ever before (Hussar & Bailey, 2013; Ryan & Siebens, 2012; Stoops, 2003) and could feasibly rise to positions of power in organizations more quickly than others have in the past (e.g., promotion rates vary as a function of age) (Rosenbaum, 1979; see also Clemens, 2012 conceptualization of the "fast track effect"). Furthermore, older workers are increasingly delaying retirement beyond the normative retirement age (Baltes & Rudolph, 2012; Burtless, 2012; Flynn, 2010), and already retired individuals are seeking re-employment in bridge employment roles in higher numbers than before (e.g., Adams & Rau, 2004; Kim & Feldman, 2000; Weckerle & Shultz, 1999).

The effect of internal audit function quality and internal audit contribution to external audit on audit fees

Corporate governance mandates and listing rules identify internal audit functions (IAF) as a central internal control mechanism. External audits are expected to assess the quality of IAF before placing reliance on its work. We provide evidence on the effect of IAF quality and IAF contribution to external audit on audit fees. Using data from a matched survey of both external and internal audits, we extend prior research which is based mainly on internal audits' assessment and conducted predominantly in highly developed markets. We find a positive relationship between IAF quality and audit fees as well as a reduction in audit fees as a result of external auditors' reliance on IAF. The interaction between IAF quality and IAF contribution to external audit suggests that high quality IAF induces greater external auditor reliance on internal auditors' work and thus result in lower external audit fees.