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.

WT1 interacts with MAD2 and regulates mitotic checkpoint function

Tumour suppressors safeguard the fidelity of the mitotic checkpoint by transcriptional regulation of genes that encode components of the mitotic checkpoint complex (MCC). Here we report a new role for the tumour suppressor and transcription factor, WT1, in the mitotic checkpoint. We show that WT1 regulates the MCC by directly interacting with the spindle assembly checkpoint protein, MAD2. WT1 colocalizes with MAD2 during mitosis and preferentially binds to the functionally active, closed-conformer, C-MAD2. Furthermore, WT1 associates with the MCC containing MAD2, BUBR1 and CDC20, resulting in prolonged inhibition of the anaphase-promoting complex/cyclosome (APC/C) and delayed degradation of its substrates SECURIN and CYCLIN B1. Strikingly, RNA interference-mediated depletion of WT1 leads to enhanced turnover of SECURIN, decreased lag time to anaphase and defects in chromosome segregation. Our findings identify WT1 as a regulator of the mitotic checkpoint and chromosomal stability.

Recruiting and retaining plasmapheresis donors: A critical belief analysis

This paper identifies critical beliefs underpinning intentions to commence and continue plasmapheresis donation. Whole blood (n = 624) and first-time plasmapheresis (n = 460) donors completed a cross-sectional survey assessing the belief-base of the theory of planned behaviour and rated their plasmapheresis donation intentions. While the idea of red blood cells being returned was a key deterrent for all donors, critical beliefs underlying commencement and continuation in the plasmapheresis donor panel differed and varied as a function of blood donation history. Findings will assist the development of targeted persuasion messages to optimise recruitment and retention of plasmapheresis donors in a non-remunerated context.

IT consumerization and its effects on IT business value, IT capabilities, and the IT function

IT consumerization is both a major opportunity and significant challenge for organizations. However, IS research has hardly discussed the implications for IT management so far. In this paper we address this topic by empirically identifying organizational themes for IT consumerization and conceptually exploring the direct and indirect effects on the business value of IT, IT capabilities, and the IT function. More specifically, based on two case studies, we identify eight organizational themes: consumer IT strategy, policy development and responsibilities, consideration of private life of employees, user involvement into IT-related processes, individualization, updated IT infrastructure, end user support, and data and system security. The contributions of this paper are: (1) the identification of organizational themes for IT consumerization; (2) the proposed effects on the business value of IT, IT capabilities and the IT function, and; (3) combining empirical insights into IT consumerization with managerial theories in the IS discipline.

Understanding the function and mechanisms of intestinal cell kinase in the growth and survival of prostate cancer cells

This project was a step forward in discovering the potential role of intestinal cell kinase in prostate cancer development. Intestinal cell kinase was shown to be upregulated in prostate cancer cells and altered expression led to changes in key cell survival proteins. This study used in vitro experiments to monitor changes in cell growth, protein and RNA expression.

Cross-cultural adaptation and validation of the Foot Function Index to Spanish

Background The purpose of this study was to adapt and validate the Foot Function Index to the Spanish (FFI-Sp) following the guidelines of the American Academy of Orthopaedic Surgeons. Methods A cross-sectional study 80 participants with some foot pathology. A statistical analysis was made, including a correlation study with other questionnaires (the Foot Health Status Questionnaire, EuroQol 5-D, Visual Analogue Pain Scale, and the Short Form SF-12 Health Survey). Data analysis included reliability, construct and criterion-related validity and factor analyses. Results The principal components analysis with varimax rotation produced 3 principal factors that explained 80% of the variance. The confirmatory factor analysis showed an acceptable fit with a comparative fit index of 0.78. The FFI-Sp demonstrated excellent internal consistency on the three subscales: pain 0.95; disability 0.96; and activity limitation 0.69, the subscale that scored lowest. The correlation between the FFI-Sp and the other questionnaires was high to moderate. Conclusions The Spanish version of the Foot Function Index (FFI-Sp) is a tool that is a valid and reliable tool with a very good internal consistency for use in the assessment of pain, disability and limitation of the function of the foot, for use both in clinic and research.

A kinematic analysis of the hand function

Background The evaluation of the hand function is an essential element within the clinical practice. The usual assessments are focus on the ability to perform activities of daily life. The inclusion of instruments to measure kinematic variables provides a new approach to the assessment. Inertial sensors adapted to the hand could be used as a complementary instrument to the traditional assessment. Material: clinimetric assessment (Upper Limb Functional Index, Quick Dash), antrophometric variables (eight and weight), dynamometry (palm preasure) was taken. Functional analysis was made with Acceleglove system for the right hand and computer system. The glove has six acceleration sensor, one on each finger and another one on the reverse palm. Method Analytic, transversal approach. Ten healthy subject made six task on evaluation table (tripod pinch, lateral pinch and tip pinch, extension grip, spherical grip and power grip). Each task was made and measure three times, the second one was analyze for the results section. A Matlab script was created for the analysis of each movement and detection phase based on module vector. Results The module acceleration vector offers useful information of the hand function. The data analysis obtained during the performance of functional gestures allows to identify five different phases within the movement, three static phase and tow dynamic, each module vector was allied to one task. Conclusion Module vector variables could be used for the analysis of the different task made by the hand. Inertial sensor could be use as a complement for the traditional assessment system.

Brain activity during the encoding, retention, and retrieval of stimulus representations

Studies of delayed nonmatching-to-sample (DNMS) performance following lesions of the monkey cortex have revealed a critical circuit of brain regions involved in forming memories and retaining and retrieving stimulus representations. Using event-related functional magnetic resonance imaging (fMRI), we measured brain activity in 10 healthy human participants during performance of a trial-unique visual DNMS task using novel barcode stimuli. The event-related design enabled the identification of activity during the different phases of the task (encoding, retention, and retrieval). Several brain regions identified by monkey studies as being important for successful DNMS performance showed selective activity during the different phases, including the mediodorsal thalamic nucleus (encoding), ventrolateral prefrontal cortex (retention), and perirhinal cortex (retrieval). Regions showing sustained activity within trials included the ventromedial and dorsal prefrontal cortices and occipital cortex. The present study shows the utility of investigating performance on tasks derived from animal models to assist in the identification of brain regions involved in human recognition memory.

Probabilistic multi-tensor estimation using the tensor distribution function

Diffusion weighted magnetic resonance (MR) imaging is a powerful tool that can be employed to study white matter microstructure by examining the 3D displacement profile of water molecules in brain tissue. By applying diffusion-sensitized gradients along a minimum of 6 directions, second-order tensors can be computed to model dominant diffusion processes. However, conventional DTI is not sufficient to resolve crossing fiber tracts. Recently, a number of high-angular resolution schemes with greater than 6 gradient directions have been employed to address this issue. In this paper, we introduce the Tensor Distribution Function (TDF), a probability function defined on the space of symmetric positive definite matrices. Here, fiber crossing is modeled as an ensemble of Gaussian diffusion processes with weights specified by the TDF. Once this optimal TDF is determined, the diffusion orientation distribution function (ODF) can easily be computed by analytic integration of the resulting displacement probability function.

A novel measure of fractional anisotropy based on the tensor distribution function

Fractional anisotropy (FA), a very widely used measure of fiber integrity based on diffusion tensor imaging (DTI), is a problematic concept as it is influenced by several quantities including the number of dominant fiber directions within each voxel, each fiber's anisotropy, and partial volume effects from neighboring gray matter. With High-angular resolution diffusion imaging (HARDI) and the tensor distribution function (TDF), one can reconstruct multiple underlying fibers per voxel and their individual anisotropy measures by representing the diffusion profile as a probabilistic mixture of tensors. We found that FA, when compared with TDF-derived anisotropy measures, correlates poorly with individual fiber anisotropy, and may sub-optimally detect disease processes that affect myelination. By contrast, mean diffusivity (MD) as defined in standard DTI appears to be more accurate. Overall, we argue that novel measures derived from the TDF approach may yield more sensitive and accurate information than DTI-derived measures.

Analyzing multi-fiber reconstruction in high angular resolution diffusion imaging using the tensor distribution function

High-angular resolution diffusion imaging (HARDI) can reconstruct fiber pathways in the brain with extraordinary detail, identifying anatomical features and connections not seen with conventional MRI. HARDI overcomes several limitations of standard diffusion tensor imaging, which fails to model diffusion correctly in regions where fibers cross or mix. As HARDI can accurately resolve sharp signal peaks in angular space where fibers cross, we studied how many gradients are required in practice to compute accurate orientation density functions, to better understand the tradeoff between longer scanning times and more angular precision. We computed orientation density functions analytically from tensor distribution functions (TDFs) which model the HARDI signal at each point as a unit-mass probability density on the 6D manifold of symmetric positive definite tensors. In simulated two-fiber systems with varying Rician noise, we assessed how many diffusionsensitized gradients were sufficient to (1) accurately resolve the diffusion profile, and (2) measure the exponential isotropy (EI), a TDF-derived measure of fiber integrity that exploits the full multidirectional HARDI signal. At lower SNR, the reconstruction accuracy, measured using the Kullback-Leibler divergence, rapidly increased with additional gradients, and EI estimation accuracy plateaued at around 70 gradients.

Numerical computation of an Evans function for travelling waves

We demonstrate a geometrically inspired technique for computing Evans functions for the linearised operators about travelling waves. Using the examples of the F-KPP equation and a Keller–Segel model of bacterial chemotaxis, we produce an Evans function which is computable through several orders of magnitude in the spectral parameter and show how such a function can naturally be extended into the continuous spectrum. In both examples, we use this function to numerically verify the absence of eigenvalues in a large region of the right half of the spectral plane. We also include a new proof of spectral stability in the appropriate weighted space of travelling waves of speed c≥sqrt(2δ) in the F-KPP equation.

Hydraulic conductivity increases in a sodic clay soil in response to gypsum applications impacts of bulk density and cation exchange

Amelioration of sodic soils is commonly achieved by applying gypsum, which increases soil hydraulic conductivity by altering soil chemistry. The magnitude of hydraulic conductivity increases expected in response to gypsum applications depends on soil properties including clay content, clay mineralogy, and bulk density. The soil analyzed in this study was a kaolinite rich sodic clay soil from an irrigated area of the Lower Burdekin coastal floodplain in tropical North Queensland, Australia. The impact of gypsum amelioration was investigated by continuously leaching soil columns with a saturated gypsum solution, until the hydraulic conductivity and leachate chemistry stabilized. Extended leaching enabled the full impacts of electrolyte effects and cation exchange to be determined. For the columns packed to 1.4 g/cm3, exchangeable sodium concentrations were reduced from 5.0 ± 0.5 mEq/100 g to 0.41 ± 0.06 mEq/100 g, exchangeable magnesium concentrations were reduced from 13.9 ± 0.3 mEq/100 g to 4.3 ± 2.12 mEq/100 g, and hydraulic conductivity increased to 0.15 ± 0.04 cm/d. For the columns packed to 1.3 g/cm3, exchangeable sodium concentrations were reduced from 5.0 ± 0.5 mEq/100 g to 0.51 ± 0.03 mEq/100 g, exchangeable magnesium concentrations were reduced from 13.9 ± 0.3 mEq/100 g to 0.55 ± 0.36 mEq/100 g, and hydraulic conductivity increased to 0.96 ± 0.53 cm/d. The results of this study highlight that both sodium and magnesium need to be taken into account when determining the suitability of water quality for irrigation of sodic soils and that soil bulk density plays a major role in controlling the extent of reclamation that can be achieved using gypsum applications.

Effect of nitrification inhibitors (DMPP and 3MP+ TZ) on soil nitrous oxide emissions from a sub-tropical vegetable production system in Queensland, Australia

The use of nitrification inhibitors, in combination with ammonium based fertilisers, has been promoted recently as an effective method to reduce nitrous oxide (N2O) emissions from fertilised agricultural fields, whilst increasing yield and nitrogen use efficiency. Vegetable cropping systems are often characterised by high inputs of nitrogen fertiliser and consequently elevated emissions of nitrous oxide (N2O) can be expected. However, to date only limited data is available on the use of nitrification inhibitors in sub-tropical vegetable systems. A field experiment investigated the effect of the nitrification inhibitors (DMPP & 3MP+TZ) on N2O emissions and yield from a typical vegetable production system in sub-tropical Australia. Soil N2O fluxes were monitored continuously over an entire year with a fully automated system. Measurements were taken from three subplots for each treatment within a randomized complete blocks design. There was a significant inhibition effect of DMPP and 3MP+TZ on N2O emissions and soil mineral N content directly following the application of the fertiliser over the vegetable cropping phase. However this mitigation was offset by elevated N2O emissions from the inhibitor treatments over the post-harvest fallow period. Cumulative annual N2O emissions amounted to 1.22 kg-N/ha, 1.16 kg-N/ha, 1.50 kg-N/ha and 0.86 kg-N/ha in the conventional fertiliser (CONV), the DMPP treatment, the 3MP+TZ treatment and the zero fertiliser (0N) respectively. Corresponding fertiliser induced emission factors (EFs) were low with only 0.09 - 0.20% of the total applied fertiliser lost as N2O. There was no significant effect of the nitrification inhibitors on yield compared to the CONV treatment for the three vegetable crops (green beans, broccoli, lettuce) grown over the experimental period. This study highlights that N2O emissions from such vegetable cropping system are primarily controlled by post-harvest emissions following the incorporation of vegetable crop residues into the soil. It also shows that the use of nitrification inhibitors can lead to elevated N2O emissions by storing N in the soil profile that is available to soil microbes during the decomposition of the vegetable residues over the post-harvest phase. Hence the use of nitrification inhibitors in vegetable systems has to be treated carefully and fertiliser rates need to be adjusted to avoid excess soil nitrogen during the postharvest phase.

A novel 15N tracer approach for the quantification of N2 and N2O emissions from soil incubations in a completely automated laboratory set up

The microbial mediated production of nitrous oxide (N2O) and its reduction to dinitrogen (N2) via denitrification represents a loss of nitrogen (N) from fertilised agro-ecosystems to the atmosphere. Although denitrification has received great interest by biogeochemists in the last decades, the magnitude of N2lossesand related N2:N2O ratios from soils still are largely unknown due to methodical constraints. We present a novel 15N tracer approach, based on a previous developed tracer method to study denitrification in pure bacterial cultures which was modified for the use on soil incubations in a completely automated laboratory set up. The method uses a background air in the incubation vessels that is replaced with a helium-oxygen gas mixture with a 50-fold reduced N2 background (2 % v/v). This method allows for a direct and sensitive quantification of the N2 and N2O emissions from the soil with isotope-ratio mass spectrometry after 15N labelling of denitrification N substrates and minimises the sensitivity to the intrusion of atmospheric N2 at the same time. The incubation set up was used to determine the influence of different soil moisture levels on N2 and N2O emissions from a sub-tropical pasture soil in Queensland/Australia. The soil was labelled with an equivalent of 50 μg-N per gram dry soil by broadcast application of KNO3solution (4 at.% 15N) and incubated for 3 days at 80% and 100% water filled pore space (WFPS), respectively. The headspace of the incubation vessel was sampled automatically over 12hrs each day and 3 samples (0, 6, and 12 hrs after incubation start) of headspace gas analysed for N2 and N2O with an isotope-ratio mass spectrometer (DELTA V Plus, Thermo Fisher Scientific, Bremen, Germany(. In addition, the soil was analysed for 15N NO3- and NH4+ using the 15N diffusion method, which enabled us to obtain a complete N balance. The method proved to be highly sensitive for N2 and N2O emissions detecting N2O emissions ranging from 20 to 627 μN kg-1soil-1hr-1and N2 emissions ranging from 4.2 to 43 μN kg-1soil-1hr-1for the different treatments. The main end-product of denitrification was N2O for both water contents with N2 accounting for 9% and 13% of the total denitrification losses at 80% and 100%WFPS, respectively. Between 95-100% of the added 15N fertiliser could be recovered. Gross nitrification over the 3 days amounted to 8.6 μN g-1 soil-1 and 4.7 μN g-1 soil-1, denitrification to 4.1 μN g-1 soil-1 and 11.8 μN g-1 soil-1at 80% and 100%WFPS, respectively. The results confirm that the tested method allows for a direct and highly sensitive detection of N2 and N2O fluxes from soils and hence offers a sensitive tool to study denitrification and N turnover in terrestrial agro-ecosystems.