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.

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.

The tensor distribution function

Diffusion weighted magnetic resonance 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 six directions, second-order tensors (represented by three-by-three positive definite matrices) can be computed to model dominant diffusion processes. However, conventional DTI is not sufficient to resolve more complicated white matter configurations, e.g., crossing fiber tracts. Recently, a number of high-angular resolution schemes with more than six gradient directions have been employed to address this issue. In this article, we introduce the tensor distribution function (TDF), a probability function defined on the space of symmetric positive definite matrices. Using the calculus of variations, we solve the TDF that optimally describes the observed data. 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 orientation distribution function (ODF) can easily be computed by analytic integration of the resulting displacement probability function. Moreover, a tensor orientation distribution function (TOD) may also be derived from the TDF, allowing for the estimation of principal fiber directions and their corresponding eigenvalues.

Early mathematical learning: Number processing skills and executive function at 5 and 8 years of age

This research investigated differences and associations in performance in number processing and executive function for children attending primary school in a large Australian metropolitan city. In a cross-sectional study, performance of 25 children in the first full-time year of school, (Prep; mean age = 5.5 years) and 21 children in Year 3 (mean age = 8.5 years) completed three number processing tasks and three executive function tasks. Year 3 children consistently outperformed the Prep year children on measures of accuracy and reaction time, on the tasks of number comparison, calculation, shifting, and inhibition but not on number line estimation. The components of executive function (shifting, inhibition, and working memory) showed different patterns of correlation to performance on number processing tasks across the early years of school. Findings could be used to enhance teachers’ understanding about the role of the cognitive processes employed by children in numeracy learning, and so inform teachers’ classroom practices.

The MS risk allele of CD40 is associated with reduced cell-membrane bound expression in antigen presenting cells: Implications for gene function

Human genetic and animal studies have implicated the costimulatory molecule CD40 in the development of multiple sclerosis (MS). We investigated the cell specific gene and protein expression variation controlled by the CD40 genetic variant(s) associated with MS, i.e. the T-allele at rs1883832. Previously we had shown that the risk allele is expressed at a lower level in whole blood, especially in people with MS. Here, we have defined the immune cell subsets responsible for genotype and disease effects on CD40 expression at the mRNA and protein level. In cell subsets in which CD40 is most highly expressed, B lymphocytes and dendritic cells, the MS-associated risk variant is associated with reduced CD40 cell-surface protein expression. In monocytes and dendritic cells, the risk allele additionally reduces the ratio of expression of full-length versus truncated CD40 mRNA, the latter encoding secreted CD40. We additionally show that MS patients, regardless of genotype, express significantly lower levels of CD40 cell-surface protein compared to unaffected controls in B lymphocytes. Thus, both genotype-dependent and independent down-regulation of cell-surface CD40 is a feature of MS. Lower expression of a co-stimulator of T cell activation, CD40, is therefore associated with increased MS risk despite the same CD40 variant being associated with reduced risk of other inflammatory autoimmune diseases. Our results highlight the complexity and likely individuality of autoimmune pathogenesis, and could be consistent with antiviral and/or immunoregulatory functions of CD40 playing an important role in protection from MS. © 2015 Field et al.

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.

Crystallization kinetics of poly(l-lactic acid)

The morphology and crystal growth of poly(l-lactic acid), PLLA have been studied from the melt as a function of undercooling and molecular weight using hot stage microscopy. Attention has been given to the application of growth rate equation on the growth rate data of PLLA and thus various nucleation parameters have been calculated. The criteria of Regime I and Regime II types of crystallization has been applied for the evaluation of substrate lengths.

Age-related changes in hepatic function: An update on implications for drug therapy

The accumulation of deficits with increasing age results in a decline in the functional capacity of multiple organs and systems. These changes can have a significant influence on the pharmacokinetics and pharmacodynamics of prescribed drugs. Although alterations in body composition and worsening renal clearance are important considerations, for most drugs the liver has the greatest effect on metabolism. Age-related change in hepatic function thereby causes much of the variability in older people’s responses to medication. In this review, we propose that a decline in the ability of the liver to inactivate toxins may contribute to a proinflammatory state in which frailty can develop. Since inflammation also downregulates drug metabolism, medication prescribed to frail older people in accordance with disease-specific guidelines may undergo reduced systemic clearance, leading to adverse drug reactions, further functional decline and increasing polypharmacy, exacerbating rather than ameliorating frailty status. We also describe how increasing chronological age and frailty status impact liver size, blood flow and protein binding and enzymes of drug metabolism. This is used to contextualise our discussion of appropriate prescribing practices. For example, while the general axiom of ‘start low, go slow’ should underpin the initiation of medication (titrating to a defined therapeutic goal), it is important to consider whether drug clearance is flow or capacity-limited. By summarising the effect of age-related changes in hepatic function on medications commonly used in older people, we aim to provide a guide that will have high clinical utility for practising geriatricians.

Mango (Mangifera indica L.) peel extract fractions from different cultivars differentially affect lipid accumulation in 3T3-L1 adipocyte cells

Plant phytochemicals are increasingly recognised as sources of bioactive molecules which may have potential benefit in many health conditions. In mangoes, peel extracts from different cultivars exhibit varying effects on adipogenesis in the 3T3-L1 adipocyte cell line. In this study, the effects of preparative HPLC fractions of methanol peel extracts from Irwin, Nam Doc Mai and Kensington Pride mangoes were evaluated. Fraction 1 contained the most hydrophilic components while subsequent fractions contained increasingly more hydrophobic components. High content imaging was used to assess mango peel fraction effects on lipid accumulation, nuclei count and nuclear area in differentiating 3T3-L1 cells. For all three mango cultivars, the more hydrophilic peel fractions 1-3 inhibited lipid accumulation with greater potency than the more hydrophobic peel fractions 4. For all three cultivars, the more lipophilic fraction 4 had concentrations that enhanced lipid accumulation greater than fractions 1-3 as assessed by lipid droplet integrated intensity. The potency of this fraction 4 varied significantly between cultivars. Using mass spectrometry, five long chain free fatty acids were detected in fraction 4; these were not present in any other peel extract fractions. Total levels varied between cultivars, with Irwin fraction 4 containing the highest levels of these free fatty acids. Lipophilic components appear to be responsible for the lipid accumulation promoting effects of some mango extracts and are the likely cause of the diverse effects of peel extracts from different mango cultivars on lipid accumulation.

Study of the structure and function of the sheath of tropically adapted bulls and risk factors for preputial eversion and preputial prolapse

Preputial prolapse is an obvious condition affecting bulls from many breeds. Unfortunately, the losses in production and welfare concerns associated with preputial prolapse can remain undetected for long periods of time in the extensive beef areas of northern Australia where the bulls are not inspected regularly. Thus, there is a critical need to identify the structural factors predisposing to preputial prolapse in young bulls so that they can be culled early. Despite there being no firm scientific evidence of an association between preputial eversion and preputial prolapse, it seems logical that the increased exposure of the sensitive prepuce as a consequence of preputial eversion may increase the risk of bulls developing preputial pathology, in particular preputial prolapse. This may be particularly relevant in Bos indicus bulls as they have a more pendulous sheath and thus eversion of the prepuce may be associated with a greater risk of injury to the prepuce compared to that in Bos taurus bulls. Further, studies of preputial eversion in Bos taurus bulls have concluded that there is an association between polledness and increased prevalence and severity (length of everted prepuce and duration of eversion) of preputial eversion due primarily to the absence or poor development of the caudal preputial muscles. No similar definitive work in Bos indicus bulls has been conducted and thus anatomical studies reported in this thesis were conducted to determine if a similar association occurred in Bos indicus bulls. A survey of a sample of large beef breeding herds in northern Australia found that preputial prolapse is a significant problem in Bos indicus and Bos indicus derived bulls and affected both young and older bulls. The importance of preputial prolapse confirmed the value of further research into the causes of this problem. A series of anatomical studies confirmed that preputial eversion in Bos indicus derived bulls was not more prevalent in polled bulls than horned bulls and was not associated with deficiency of the caudal preputial muscles as was established in Bos taurus bulls. An anatomical study of Bos indicus derived bulls with preputial prolapse found that preputial prolapse occurred in horned bulls of varying ages and these bulls did not have any evidence of deficiency in the caudal preputial muscles. However, preputial prolapse was observed in young polled bulls that had poorly developed or absent caudal preputial muscles. It was concluded that deficiency of the caudal preputial muscles in polled Bos indicus derived bulls may predispose to preputial prolapse at an early age, but no predisposing anatomical factors were found for horned Bos indicus derived bulls. In these studies, preputial eversion and preputial prolapse were found in horned Bos indicus derived bulls that did not have any preputial muscle deficiency and it was noted that preputial eversion was not related to the length of the prepuce. Further studies confirmed that preputial eversion was linearly and consistently associated with position of the glans penis within the sheath in Bos indicus derived bulls, and movement of the glans penis towards the preputial orifice consistently resulted in preputial eversion in these bulls. A method to objectively measure the relationship between movement of the glans penis within the sheath and preputial eversion was developed. Studies in humans have linked function of some abdominal muscles to function of the pelvic organs. This relationship was investigated in Bos indicus derived bulls to determine whether the function of specific abdominal muscles affected position of the penis in the sheath. Using the method developed to objectively measure the relationship between penis movement and preputial eversion, the abdominal muscles that potentially were associated with movement of the glans penis or preputial eversion were examined but no significant relationships were observed. In the anatomical study of Bos indicus derived bulls not affected with preputial prolapse a more pendulous sheath was associated with increased prevalence of preputial eversion. This relationship was confirmed for horned and polled bulls in the penis movement studies. Bos indicus derived bulls with more pendulous sheaths evert their prepuces more than bulls with less pendulous sheaths thus increasing the risk of damage to the prepuce either from the environment, other bulls, or from them inadvertently stepping on the everted prepuce when they get to their feet. Culling Bos indicus derived bulls with more pendulous sheaths should reduce the incidence of preputial eversion and possibly preputial prolapse. The anatomical study of Bos indicus derived bulls that did not have preputial prolapse demonstrates that there are herds of bulls where the polled bulls do not have any evidence of deficiency of the caudal preputial iv muscles. There is a need to develop a practical and cost effective test to identify polled Bos indicus bulls that have a deficiency in their caudal preputial muscles.

A Bayesian latent class safety performance function for identifying motor vehicle crash blackspots

The current state of the practice in Blackspot Identification (BSI) utilizes safety performance functions based on total crash counts to identify transport system sites with potentially high crash risk. This paper postulates that total crash count variation over a transport network is a result of multiple distinct crash generating processes including geometric characteristics of the road, spatial features of the surrounding environment, and driver behaviour factors. However, these multiple sources are ignored in current modelling methodologies in both trying to explain or predict crash frequencies across sites. Instead, current practice employs models that imply that a single underlying crash generating process exists. The model mis-specification may lead to correlating crashes with the incorrect sources of contributing factors (e.g. concluding a crash is predominately caused by a geometric feature when it is a behavioural issue), which may ultimately lead to inefficient use of public funds and misidentification of true blackspots. This study aims to propose a latent class model consistent with a multiple crash process theory, and to investigate the influence this model has on correctly identifying crash blackspots. We first present the theoretical and corresponding methodological approach in which a Bayesian Latent Class (BLC) model is estimated assuming that crashes arise from two distinct risk generating processes including engineering and unobserved spatial factors. The Bayesian model is used to incorporate prior information about the contribution of each underlying process to the total crash count. The methodology is applied to the state-controlled roads in Queensland, Australia and the results are compared to an Empirical Bayesian Negative Binomial (EB-NB) model. A comparison of goodness of fit measures illustrates significantly improved performance of the proposed model compared to the NB model. The detection of blackspots was also improved when compared to the EB-NB model. In addition, modelling crashes as the result of two fundamentally separate underlying processes reveals more detailed information about unobserved crash causes.

Identifying the Function of Sorghum's Drought Tolerance Stay-Green QTL

The goal of this research is to understand the function of allelic variation of genes underpinning the stay-green drought adaptation trait in sorghum in order to enhance yield in water-limited environments. Stay-green, a delayed leaf senescence phenotype in sorghum, is primarily an emergent consequence of the improved balance between the supply and demand of water. Positional and functional fine-mapping of candidate genes associated with stay-green in sorghum is the focus of an international research partnership between Australian (UQ/DAFFQ) and US (Texas A&M University) scientists. Stay-green was initially mapped to four chromosomal regions (Stg1, Stg2, Stg3, and Stg4) by a number of research groups in the US and Australia. Physiological dissection of near-isolines containing single introgressions of Stg QTL (Stg1-4) indicate that these QTL reduce water demand before flowering by constricting the size of the canopy, thereby increasing water availability during grain filling and, ultimately, grain yield. Stg and root angle QTL are also co-located and, together with crop water use data, suggest the role of roots in the stay-green phenomenon. Candidate genes have been identified in Stg1-4, including genes from the PIN family of auxin efflux carriers in Stg1 and Stg2, with 10 of 11 PIN genes in sorghum co-locating with Stg QTL. Modified gene expression in some of these PIN candidates in the stay-green compared with the senescent types has been found in preliminary RNA expression profiling studies. Further proof-of-function studies are underway, including comparative genomics, SNP analysis to assess diversity at candidate genes, reverse genetics and transformation.

Dirichlet problem at infinity for p-harmonic functions on negatively curved spaces

The object of this dissertation is to study globally defined bounded p-harmonic functions on Cartan-Hadamard manifolds and Gromov hyperbolic metric measure spaces. Such functions are constructed by solving the so called Dirichlet problem at infinity. This problem is to find a p-harmonic function on the space that extends continuously to the boundary at inifinity and obtains given boundary values there. The dissertation consists of an overview and three published research articles. In the first article the Dirichlet problem at infinity is considered for more general A-harmonic functions on Cartan-Hadamard manifolds. In the special case of two dimensions the Dirichlet problem at infinity is solved by only assuming that the sectional curvature has a certain upper bound. A sharpness result is proved for this upper bound. In the second article the Dirichlet problem at infinity is solved for p-harmonic functions on Cartan-Hadamard manifolds under the assumption that the sectional curvature is bounded outside a compact set from above and from below by functions that depend on the distance to a fixed point. The curvature bounds allow examples of quadratic decay and examples of exponential growth. In the final article a generalization of the Dirichlet problem at infinity for p-harmonic functions is considered on Gromov hyperbolic metric measure spaces. Existence and uniqueness results are proved and Cartan-Hadamard manifolds are considered as an application.