RcppArmadillo : accelerating R with high-performance C++ linear algebra

The R statistical environment and language has demonstrated particular strengths for interactive development of statistical algorithms, as well as data modelling and visualisation. Its current implementation has an interpreter at its core which may result in a performance penalty in comparison to directly executing user algorithms in the native machine code of the host CPU. In contrast, the C++ language has no built-in visualisation capabilities, handling of linear algebra or even basic statistical algorithms; however, user programs are converted to high-performance machine code, ahead of execution. A new method avoids possible speed penalties in R by using the Rcpp extension package in conjunction with the Armadillo C++ matrix library. In addition to the inherent performance advantages of compiled code, Armadillo provides an easy-to-use template-based meta-programming framework, allowing the automatic pooling of several linear algebra operations into one, which in turn can lead to further speedups. With the aid of Rcpp and Armadillo, conversion of linear algebra centered algorithms from R to C++ becomes straightforward. The algorithms retains the overall structure as well as readability, all while maintaining a bidirectional link with the host R environment. Empirical timing comparisons of R and C++ implementations of a Kalman filtering algorithm indicate a speedup of several orders of magnitude.

Analytical Solution of a Walters’ Liquid B Flow Over a Linear Stretching Sheet in a Porous Medium

This chapter represents the analytical solution of two-dimensional linear stretching sheet problem involving a non-Newtonian liquid and suction by (a) invoking the boundary layer approximation and (b) using this result to solve the stretching sheet problem without using boundary layer approximation. The basic boundary layer equations for momentum, which are non-linear partial differential equations, are converted into non-linear ordinary differential equations by means of similarity transformation. The results reveal a new analytical procedure for solving the boundary layer equations arising in a linear stretching sheet problem involving a non-Newtonian liquid (Walters’ liquid B). The present study throws light on the analytical solution of a class of boundary layer equations arising in the stretching sheet problem.

Heteroscedastic probabilistic linear discriminant analysis for manifold learning in video-based face recognition

To recognize faces in video, face appearances have been widely modeled as piece-wise local linear models which linearly approximate the smooth yet non-linear low dimensional face appearance manifolds. The choice of representations of the local models is crucial. Most of the existing methods learn each local model individually meaning that they only anticipate variations within each class. In this work, we propose to represent local models as Gaussian distributions which are learned simultaneously using the heteroscedastic probabilistic linear discriminant analysis (PLDA). Each gallery video is therefore represented as a collection of such distributions. With the PLDA, not only the within-class variations are estimated during the training, the separability between classes is also maximized leading to an improved discrimination. The heteroscedastic PLDA itself is adapted from the standard PLDA to approximate face appearance manifolds more accurately. Instead of assuming a single global within-class covariance, the heteroscedastic PLDA learns different within-class covariances specific to each local model. In the recognition phase, a probe video is matched against gallery samples through the fusion of point-to-model distances. Experiments on the Honda and MoBo datasets have shown the merit of the proposed method which achieves better performance than the state-of-the-art technique.

Image based visual servo control for fixed wing UAVs tracking linear infrastructure in wind

This paper presents an Image Based Visual Servo control design for Fixed Wing Unmanned Aerial Vehicles tracking locally linear infrastructure in the presence of wind using a body fixed imaging sensor. Visual servoing offers improved data collection by posing the tracking task as one of controlling a feature as viewed by the inspection sensor, although is complicated by the introduction of wind as aircraft heading and course angle no longer align. In this work it is shown that the effects of wind alter the desired line angle required for continuous tracking to equal the wind correction angle as would be calculated to set a desired course. A control solution is then sort by linearizing the interaction matrix about the new feature pose such that kinematics of the feature can be augmented with the lateral dynamics of the aircraft, from which a state feedback control design is developed. Simulation results are presented comparing no compensation, integral control and the proposed controller using the wind correction angle, followed by an assessment of response to atmospheric disturbances in the form of turbulence and wind gusts

Board age and gender diversity : a test of competing linear and curvilinear predictions

We test competing linear and curvilinear predictions between board diversity and performance. The predictions were tested using archival data on 288 organizations listed on the Australian Securities Exchange. The findings provide additional evidence on the business case for board gender diversity and refine the business case for board age diversity.

A longitudinal study of grieving in family caregivers of people with dementia

A longitudinal study of grieving in family caregivers of people with dementia Recent research into dementia has identified the long term impact that the role of care giving for a relative with dementia has on family members This is largely due to the cognitive decline that characterises dementia and the losses that can be directly attributed to this. These losses include loss of memories, relationships and intimacy, and are often ambiguous so that the grief that accompanies them is commonly not recognised or acknowledged. The role and impact of pre-death or anticipatory grief has not previously been widely considered as a factor influencing health and well-being of family caregivers. Studies of grief in caregivers of a relative with dementia have concluded that grief is one of the greatest barriers to care giving and is a primary determinant of caregiver well-being. The accumulation of losses, in conjunction with experiences unique to dementia care giving, place family caregivers at risk of complicated grief. This occurs when integration of the death does not take place following bereavement and has been associated with a range of negative health outcomes. The aim of this research was to determine the influence of grief, in addition to other factors representing both positive and negative aspects of the role, on the health related quality of life of family caregivers of people with dementia, prior to and following the death of their relative with dementia. An exploratory research project underpinned by a conceptual framework of caregivers’ adaptation in the context of subjective appraisal of the strains and gains in their role was undertaken. The research comprised three studies. Study 1 was a scoping study that involved a series of semi-structured interviews with thirteen participants who were family caregivers of people with severe dementia or whose relative with dementia had died in the previous twelve months. The results of this study in conjunction with factors identified in the literature informed data collection for the further studies. Study 2 was a cross sectional survey of fifty caregivers recruited when their relative was in the moderate to severe stage of dementia. This study provided the baseline data for Study 3, a prospective cohort follow up study. Study 3 consisted of seventeen participants followed up at two time points after the death of their relative with dementia: six weeks and then six months following the death of the relative with dementia. The scoping study indicated that differences in appraisal of the care giving role and encounters with health professionals were related to levels of grief of caregivers prior to and following the death of the relative with dementia. This was supported in the baseline and follow up studies. In the baseline study, after adjusting for all variables in multivariate regression models, subjective appraisal of burden was found to make a significant contribution (p<.05) to mental health related quality of life. The two dependent variables, anticipatory grief and mental health related quality of life, were significantly (p<.01) correlated at a bivariate level. In the follow up study, linear mixed modelling and multiple regression analysis of data found that subjective appraisal of burden and resilience were significantly associated (p<.05 and p<.01, respectively) with mental health related quality of life over time. In addition, bereavement and complicated grief were significantly associated (p<.05) with mental health following the death of the relative. In this study social support and satisfaction with end of life care were found to be statistically associated (p<.05) with physical health related quality of life over time. The strong relationship between grief of caregivers and their health related quality of life over the entire care giving trajectory and period following the death of their relative highlights the urgent need for further research and interventions in this area. Overall results indicate that addressing the risk and protective factors including subjective appraisal of their care giving role, resilience, social support and satisfaction with end of life care of their relative, has the potential to both ameliorate negative health outcomes and to promote improved health for these caregivers. This research provides important information for development of targeted and appropriate interventions that aim to promote resilience and reduce the personal burden on caregivers of people with dementia.

Improved earthwork allocation planning for fuel consumption and emissions reduction in linear infrastructure construction – case study data

This document provides data for the case study presented in our recent earthwork planning papers. Some results are also provided in a graphical format using Excel.

A class B switch-mode assisted linear amplifier

A switch-mode assisted linear amplifier (SMALA) combining a linear (Class B) and a switch-mode (Class D) amplifier is presented. The usual single hysteretic controlled half-bridge current dumping stage is replaced by two parallel buck converter stages, in a parallel voltage controlled topology. These operate independently: one buck converter sources current to assist the upper Class B output device, and a complementary converter sinks current to assist the lower device. This topology lends itself to a novel control approach of a dead-band at low power levels where neither class D amplifier assists, allowing the class B amplifier to supply the load without interference, ensuring high fidelity. A 20 W implementation demonstrates 85% efficiency, with distortion below 0.08% measured across the full audio bandwidth at 15 W. The class D amplifier begins assisting at 2 W, and below this value, the distortion was below 0.03%. Complete circuitry is given, showing the simplicity of the additional class D amplifier and its corresponding control circuitry.

Transcription and expression of Plasmodium falciparum histidine-rich proteins in different stages and strains : implications for rapid diagnostic tests

Background: Although rapid diagnostic tests (RDTs) for Plasmodium falciparum infection that target histidine rich protein 2 (PfHRP2) are generally sensitive, their performance has been reported to be variable. One possible explanation for variable test performance is differences in expression level of PfHRP in different parasite isolates. Methods: Total RNA and protein were extracted from synchronised cultures of 7 P. falciparum lines over 5 time points of the life cycle, and from synchronised ring stages of 10 falciparum lines. Using quantitative real-time polymerase chain reaction, Western blot analysis and ELISA we investigated variations in the transcription and protein levels of pfhrp2, pfhrp3 and PfHRP respectively in the different parasite lines, over the parasite intraerythrocytic life cycle. Results: Transcription of pfhrp2 and pfhrp3 in different parasite lines over the parasite life cycle was observed to vary relative to the control parasite K1. In some parasite lines very low transcription of these genes was observed. The peak transcription was observed in ring-stage parasites. Pfhrp2 transcription was observed to be consistently higher than pfhrp3 transcription within parasite lines. The intraerythrocytic lifecycle stage at which the peak level of protein was present varied across strains. Total protein levels were more constant relative to total mRNA transcription, however a maximum 24 fold difference in expression at ring-stage parasites relative to the K1 strain was observed. Conclusions: The levels of transcription of pfhrp2 and pfhrp3, and protein expression of PfHRP varied between different P. falciparum strains. This variation may impact on the detection sensitivity of PfHRP2-detecting RDTs.

Non-linear position control for hover and automatic landing of unmanned aerial vehicles

This study presents a disturbance attenuation controller for horizontal position stabilisation for hover and automatic landings of a rotary-wing unmanned aerial vehicle (RUAV) operating close to the landing deck in rough seas. Based on a helicopter model representing aerodynamics during the landing phase, a non-linear state feedback H∞ controller is designed to achieve rapid horizontal position tracking in a gusty environment. Practical constraints including flapping dynamics, servo dynamics and time lag effect are considered. A high-fidelity closed-loop simulation using parameters of the Vario XLC gas-turbine helicopter verifies performance of the proposed horizontal position controller. The proposed controller not only increases the disturbance attenuation capability of the RUAV, but also enables rapid position response when gusts occur. Comparative studies show that the H∞ controller exhibits performance improvement and can be applied to ship/RUAV landing systems.

Hysteretic multilevel NPC converter control for a Switch Mode Assisted Linear Amplifier

A Switch-Mode Assisted Linear Amplifier (SMALA) combines the high quality of a linear amplifier required for audio applications with the high efficiency of a switch-mode amplifier. The careful choice of current sense point and switch placement allows a simple non-isolated hysteresis current controller for the switch-mode section. This paper explains the extension of the hysteresis current controller for the control of a three level Neutral Point Clamped (NPC) converter, with simulations as proof of concept. The NPC topology allows the use of lower voltage switches and lower switching frequencies to implement high power audio amplifiers using the SMALA topology.

Board age and gender diversity : a test of competing linear and curvilinear predictions

The inconsistent findings of past board diversity research demand a test of competing linear and curvilinear diversity–performance predictions. This research focuses on board age and gender diversity, and presents a positive linear prediction based on resource dependence theory, a negative linear prediction based on social identity theory, and an inverted U-shaped curvilinear prediction based on the integration of resource dependence theory with social identity theory. The predictions were tested using archival data on 288 large organizations listed on the Australian Securities Exchange, with a 1-year time lag between diversity (age and gender) and performance (employee productivity and return on assets). The results indicate a positive linear relationship between gender diversity and employee productivity, a negative linear relationship between age diversity and return on assets, and an inverted U-shaped curvilinear relationship between age diversity and return on assets. The findings provide additional evidence on the business case for board gender diversity and refine the business case for board age diversity.

A simplified approach for calculating the moments of action for linear reaction-diffusion equations

The mean action time is the mean of a probability density function that can be interpreted as a critical time, which is a finite estimate of the time taken for the transient solution of a reaction-diffusion equation to effectively reach steady state. For high-variance distributions, the mean action time under-approximates the critical time since it neglects to account for the spread about the mean. We can improve our estimate of the critical time by calculating the higher moments of the probability density function, called the moments of action, which provide additional information regarding the spread about the mean. Existing methods for calculating the nth moment of action require the solution of n nonhomogeneous boundary value problems which can be difficult and tedious to solve exactly. Here we present a simplified approach using Laplace transforms which allows us to calculate the nth moment of action without solving this family of boundary value problems and also without solving for the transient solution of the underlying reaction-diffusion problem. We demonstrate the generality of our method by calculating exact expressions for the moments of action for three problems from the biophysics literature. While the first problem we consider can be solved using existing methods, the second problem, which is readily solved using our approach, is intractable using previous techniques. The third problem illustrates how the Laplace transform approach can be used to study coupled linear reaction-diffusion equations.

Block models for improved earthwork allocation planning in linear infrastructure construction

Earthwork planning has been considered in this article and a generic block partitioning and modelling approach has been devised to provide strategic plans of various levels of detail. Conceptually this approach is more accurate and comprehensive than others, for instance those that are section based. In response to environmental concerns the metric for decision making was fuel consumption and emissions. Haulage distance and gradient are also included as they are important components of these metrics. Advantageously the fuel consumption metric is generic and captures the physical difficulties of travelling over inclines of different gradients, that is consistent across all hauling vehicles. For validation, the proposed models and techniques have been applied to a real world road project. The numerical investigations have demonstrated that the models can be solved with relatively little CPU time. The proposed block models also result in solutions of superior quality, i.e. they have reduced fuel consumption and cost. Furthermore the plans differ considerably from those based solely upon a distance based metric thus demonstrating a need for industry to reflect upon their current practices.

Higher-order non-linear analysis of steel structures. Part I : elastic second-order formulation

This paper presents a higher-order beam-column formulation that can capture the geometrically non-linear behaviour of steel framed structures which contain a multiplicity of slender members. Despite advances in computational frame software, analyses of large frames can still be problematic from a numerical standpoint and so the intent of the paper is to fulfil a need for versatile, reliable and efficient non-linear analysis of general steel framed structures with very many members. Following a comprehensive review of numerical frame analysis techniques, a fourth-order element is derived and implemented in an updated Lagrangian formulation, and it is able to predict flexural buckling, snap-through buckling and large displacement post-buckling behaviour of typical structures whose responses have been reported by independent researchers. The solutions are shown to be efficacious in terms of a balance of accuracy and computational expediency. The higher-order element forms a basis for augmenting the geometrically non-linear approach with material non-linearity through the refined plastic hinge methodology described in the companion paper.