GPU Acceleration of Partial Differential Equation Solvers

Differential equations are often directly solvable by analytical means only in their one dimensional version. Partial differential equations are generally not solvable by analytical means in two and three dimensions, with the exception of few special cases. In all other cases, numerical approximation methods need to be utilized. One of the most popular methods is the finite element method. The main areas of focus, here, are the Poisson heat equation and the plate bending equation. The purpose of this paper is to provide a quick walkthrough of the various approaches that the authors followed in pursuit of creating optimal solvers, accelerated with the use of graphical processing units, and comparing them in terms of accuracy and time efficiency with existing or self-made non-accelerated solvers.

Discussion and debate: Reflections on Pleistocene Island Occupation

Editorial Note: Within a 48-hour period during January 2014, the JMA co-editors received two papers—those of Curtis Runnels and Thomas P. Leppard printed above—that, quite fortuitously, each addressed the topic of Mediterranean island colonization by archaic hominins, albeit from radically different perspectives. Neither author was aware of the other’s paper, nor has either article subsequently been revised to take account of the other. Realizing the widespread current interest in this subject and the possibility for productive debate prompted by such variant approaches, we commissioned three sets of comments and invited Runnels and Leppard to respond. We are pleased to publish this discussion around questions of great importance for our understanding of the earliest insular prehistory of the Mediterranean, and with significant implications reaching well beyond it.

Journal of Mediterranean Archaeology 27.2 (2014) 255-278

Ecological dynamics of extinct species in empty habitat networks. 1. The role of habitat pattern and quantity, stochasticity and dispersal

We examined a remnant host plant (Primula veris L.) habitat network that was last inhabited by the rare butterfly Hamearis lucina L. in north Wales in 1943, to assess the relative contribution of several spatial parameters to its regional extinction. We first examined relationships between P. veris characteristics and H. lucina eggs in surviving H. lucina populations, and used these to predict the suitability and potential carrying capacity of the habitat network in north Wales. This resulted in an estimate of roughly 4500 eggs (ca 227 adults). We developed a discrete space, discrete time metapopulation model to evaluate the relative contribution of dispersal distance, habitat and environmental stochasticity as possible causes of extinction. We simulated the potential persistence of the butterfly in the current network as well as in three artificial (historical and present) habitat networks that differed in the quantity (current and X3) and fragmentation of the habitat (current and aggregated). We identified that reduced habitat quantity and increased isolation would have increased the probability of regional extinction, in conjunction with environmental stochasticity and H. lucina's dispersal distance. This general trend did not change in a qualitative manner when we modified the ability of dispersing females to stay in, and find suitable habitats (by changing the size of the grid cells used in the model). Contrary to most metapopulation model predictions, system persistence declined with increasing migration rate, suggesting that the mortality of migrating individuals in fragmented landscapes may pose significant risks to system-wide persistence. Based on model predictions for the present landscape we argue that a major programme of habitat restoration would be required for a re-established metapopulation to persist for > 100 years.

Ecological dynamics of extinct species in empty habitat networks. 2. The role of host plant dynamics

This paper explores the relative effects of host plant dynamics and butterfly-related parameters on butterfly persistence. It considers an empty habitat network where a rare butterfly (Cupido minimus) became extinct in 1939 in part of its historical range in north Wales, UK. Surviving populations of the butterfly in southern Britain were visited to assess use of its host plant (Anthyllis vulneraria) in order to calibrate habitat suitability and carrying capacity in the empty network in north Wales. These data were used to deduce that only a portion ( similar to 19%) of the host plant network from north Wales was likely to be highly suitable for oviposition. Nonetheless, roughly 65,460 eggs (3273 adult equivalents) could be expected to be laid in north Wales, were the empty network to be populated at the same levels as observed on comparable plants in surviving populations elsewhere. Simulated metapopulations of C. minimus in the empty network revealed that time to extinction and patch occupancy were significantly influenced by carrying capacity, butterfly mean dispersal distance and environmental stochasticity, although for most reasonable parameter values, the model system persisted. Simulation outputs differed greatly when host plant dynamics was incorporated into the modelled butterfly dynamics. Cupido minimus usually went extinct when host plant were at low densities. In these simulations host plant dynamics appeared to be the most important determinant of the butterfly's regional extirpation. Modelling the outcome of a reintroduction programme to C. minimus variation at high quality locations, revealed that 65% of systems survived at least 100 years. Given the current amount of resources of the north Wales landscape, the persistence of C. minimus under a realistic reintroduction programme has a good chance of being successful, if carried out in conjunction with a host plant management programme.

Mutagenesis of RpoE-like sigma factor genes in Bdellovibrio reveals differential control of groEL and two groES genes

BACKGROUND: Bdellovibrio bacteriovorus HD100 must regulate genes in response to a variety of environmental conditions as it enters, preys upon and leaves other bacteria, or grows axenically without prey. In addition to "housekeeping" sigma factors, its genome encodes several alternate sigma factors, including 2 Group IV-RpoE-like proteins, which may be involved in the complex regulation of its predatory lifestyle.

RESULTS: We find that one sigma factor gene, bd3314, cannot be deleted from Bdellovibrio in either predatory or prey-independent growth states, and is therefore possibly essential, likely being an alternate sigma 70. Deletion of one of two Group IV-like sigma factor genes, bd0881, affects flagellar gene regulation and results in less efficient predation, although not due to motility changes; deletion of the second, bd0743, showed that it normally represses chaperone gene expression and intriguingly we find an alternative groES gene is expressed at timepoints in the predatory cycle where intensive protein synthesis at Bdellovibrio septation, prior to prey lysis, will be occurring.

CONCLUSIONS: We have taken the first step in understanding how alternate sigma factors regulate different processes in the predatory lifecycle of Bdellovibrio and discovered that alternate chaperones regulated by one of them are expressed at different stages of the lifecycle.

Differential evolution combined with clonal selection for dynamic economic dispatch

Dynamic economic load dispatch (DELD) is one of the most important steps in power system operation. Various optimisation algorithms for solving the problem have been developed; however, due to the non-convex characteristics and large dimensionality of the problem, it is necessary to explore new methods to further improve the dispatch results and minimise the costs. This article proposes a hybrid differential evolution (DE) algorithm, namely clonal selection-based differential evolution (CSDE), to solve the problem. CSDE is an artificial intelligence technique that can be applied to complex optimisation problems which are for example nonlinear, large scale, non-convex and discontinuous. This hybrid algorithm combines the clonal selection algorithm (CSA) as the local search technique to update the best individual in the population, which enhances the diversity of the solutions and prevents premature convergence in DE. Furthermore, we investigate four mutation operations which are used in CSA as the hyper-mutation operations. Finally, an efficient solution repair method is designed for DELD to satisfy the complicated equality and inequality constraints of the power system to guarantee the feasibility of the solutions. Two benchmark power systems are used to evaluate the performance of the proposed method. The experimental results show that the proposed CSDE/best/1 approach significantly outperforms nine other variants of CSDE and DE, as well as most other published methods, in terms of the quality of the solution and the convergence characteristics.

Differential ecological impacts of an invader and a native under environmental change are revealed by comparative functional responses.

Predicting the ecological impacts of damaging invasive species under relevant environmental contexts is a major challenge, for which comparative functional responses (the relationship between resource availability and consumer uptake rate) have great potential. Here, the functional responses of Gammarus pulex, an ecologically damaging invader in freshwaters in Ireland and other islands, were compared with those of a native trophic equivalent Gammarus duebeni celticus. Experiments were conducted at two dissolved oxygen concentrations (80 and 50 % saturation), representative of anthropogenic water quality changes, using two larval prey, blackfly (Simuliidae spp.) and mayfly (Baetis rhodani). Overall, G. pulex had higher Type II functional responses and hence predatory impacts than G. d. celticus and the functional responses of both predators were reduced by lowered oxygen concentration. However, this reduction was of lower magnitude for the invader as compared to the native. Further, the invader functional response at low oxygen was comparable to that of the native at high oxygen. Attack rates of the two predators were similar, with low oxygen reducing these attack rates, but this effect occurred more strongly for blackfly than mayfly prey. Handling times were significantly lower for the invader compared with the native, and significantly higher at low oxygen, however, the effect of lowered oxygen on handling times was minimal for the invader and pronounced for the native. Maximum feeding rates were significantly greater for the invader compared with the native, and significantly reduced at low oxygen, with this effect again lesser for the invader as compared to the native. The greater functional responses of the invader corroborate with its impacts on recipient macroinvertebrate communities when it replaces the native. Further, our experiments predict that the impact of the invader will be less affected than the native under altered oxygen regimes driven by anthropogenic influences.

Is differential phytochelatin production related to decreased arsenate influx in arsenate tolerant Holcus lanatus?

Phytochelatins (PCs) are required for arsenic (As) detoxification in nontolerant plants. In addition, a role for PCs in arsenate tolerance has recently been proven, with tolerant plants able to accumulate significantly higher concentrations of As-PC complexes at equivalent levels of stress than nontolerant plants. The relationship between arsenate influx and PC production in tolerant and non-tolerant Holcus lanatus plants was determined in this study, along with an investigation of the effect of inhibition of PC synthesis by buthionine sulfoximine (BSO) on arsenate tolerance. A strong correlation between PC production and arsenate influx was demonstrated in arsenate-tolerant plants. In addition, inhibition of PC synthesis by BSO in tolerant plants increased arsenate sensitivity to that of the nontolerant clone. This dramatic reduction in tolerance proves that PC production is an essential component of the arsenate tolerance mechanism in H. lanatus. This study proposes that while there is a single major gene for arsenate tolerance, hypostatic modifiers are also in operation, affecting the expression of the tolerance character. © New Phytologist (2002).

Copper- and arsenate-induced oxidative stress in Holcus lanatus L. clones with differential sensitivity

The biochemical responses of Holcus lanatus L. to copper and arsenate exposure were investigated in arsenate-tolerant and -non-tolerant plants from uncontaminated and arsenic/copper-contaminated sites. Increases in lipid peroxidation, superoxide dismutase (SOD) activity and phytochelatin (PC) production were correlated with increasing copper and arsenate exposure. In addition, significant differences in biochemical responses were observed between arsenate-tolerant and -non-tolerant plants. Copper and arsenate exposure led to the production of reactive oxygen species, resulting in significant lipid peroxidation in non-tolerant plants. However, SOD activity was suppressed upon metal exposure, possibly due to interference with metallo-enzymes. It was concluded that in non-tolerant plants, rapid arsenate influx resulted in PC production, glutathione depletion and lipid peroxidation. This process would also occur in tolerant plants, but by decreasing the rate of influx, they were able to maintain their constitutive functions, detoxify the metals though PC production and quench reactive oxygen species by SOD activity.