Weighted decomposition in high-performance lattice-Boltzmann simulations: Are some lattice sites more equal than others?

Obtaining a good load balance is a significant challenge in scaling up lattice-Boltzmann simulations of realistic sparse problems to the exascale. Here we analyze the effect of weighted decomposition on the performance of the HemeLB lattice-Boltzmann simulation environment, when applied to sparse domains. Prior to domain decomposition, we assign wall and in/outlet sites with increased weights which reflect their increased computational cost. We combine our weighted decomposition with a second optimization, which is to sort the lattice sites according to a space filling curve. We tested these strategies on a sparse bifurcation and very sparse aneurysm geometry, and find that using weights reduces calculation load imbalance by up to 85 %, although the overall communication overhead is higher than some of our runs.This work has received funding from the CRESTA and MAPPER projects within the EC-FP7 (ICT-2011.9.13) under Grant Agreements nos. 287703 and 261507, and from EPSRC Grants EP/I017909/1 (www.2020science.net) and EP/I034602/1

Chaste: an open source C++ library for computational physiology and biology

Chaste - Cancer, Heart And Soft Tissue Environment - is an open source C++ library for the computational simulation of mathematical models developed for physiology and biology. Code development has been driven by two initial applications: cardiac electrophysiology and cancer development. A large number of cardiac electrophysiology studies have been enabled and performed, including high performance computational investigations of defibrillation on realistic human cardiac geometries. New models for the initiation and growth of tumours have been developed. In particular, cell-based simulations have provided novel insight into the role of stem cells in the colorectal crypt. Chaste is constantly evolving and is now being applied to a far wider range of problems. The code provides modules for handling common scientific computing components, such as meshes and solvers for ordinary and partial differential equations (ODEs/PDEs). Re-use of these components avoids the need for researchers to "re-invent the wheel" with each new project, accelerating the rate of progress in new applications. Chaste is developed using industrially-derived techniques, in particular test-driven development, to ensure code quality, re-use and reliability. In this article we provide examples that illustrate the types of problems Chaste can be used to solve, which can be run on a desktop computer. We highlight some scientific studies that have used or are using Chaste, and the insights they have provided. The source code, both for specific releases and the development version, is available to download under an open source Berkeley Software Distribution (BSD) licence at http://www.cs.ox.ac.uk/chaste, together with details of a mailing list and links to documentation and tutorials

Analyzing and Modeling the Performance of the HemeLB Lattice-Boltzmann Simulation Environment

We investigate the performance of the HemeLB lattice-Boltzmann simulator for cerebrovascular blood flow, aimed at providing timely and clinically relevant assistance to neurosurgeons. HemeLB is optimised for sparse geometries, supports interactive use, and scales well to 32,768 cores for problems with ∼81 million lattice sites. We obtain a maximum performance of 29.5 billion site updates per second, with only an 11% slowdown for highly sparse problems (5% fluid fraction). We present steering and visualisation performance measurements and provide a model which allows users to predict the performance, thereby determining how to run simulations with maximum accuracy within time constraints

Coalesced communication: a design pattern for complex parallel scientific software

We present a new design pattern for high-performance parallel scientific software, named coalesced communication. This pattern allows for a structured way to improve the communication performance through coalescence of multiple communication needs using two communication management components. We apply the design pattern to several simulations of a lattice-Boltzmann blood flow solver with streaming visualisation which engenders a reduction in the communication overhead of approximately 40%

Flexible composition and execution of high performance, high fidelity multiscale biomedical simulations

Multiscale simulations are essential in the biomedical domain to accurately model human physiology. We present a modular approach for designing, constructing and executing multiscale simulations on a wide range of resources, from laptops to petascale supercomputers, including combinations of these. Our work features two multiscale applications, in-stent restenosis and cerebrovascular bloodflow, which combine multiple existing single-scale applications to create a multiscale simulation. These applications can be efficiently coupled, deployed and executed on computers up to the largest (peta) scale, incurring a coupling overhead of 1–10% of the total execution time

Ten Simple Rules for Effective Computational Research

<p>Ten Simple Rules for Effective Computational Research</p

Mice lacking endoglin in macrophages show an impaired immune response

24 p.-9 fig.-1 tab. Ojeda Fernández, Luisa et al.Endoglin is an auxiliary receptor for members of the TGF-β superfamily and plays an important role in the homeostasis of the vessel wall. Mutations in endoglin gene (ENG) or in the closely related TGF-β receptor type I ACVRL1/ALK1 are responsible for a rare dominant vascular dysplasia, the Hereditary Hemorrhagic Telangiectasia (HHT), or Rendu-OslerWeber syndrome. Endoglin is also expressed in human macrophages, but its role in macrophage function remains unknown. In this work, we show that endoglin expression is triggered during the monocyte-macrophage differentiation process, both in vitro and during the in vivo differentiation of blood monocytes recruited to foci of inflammation in wild-type C57BL/6 mice. To analyze the role of endoglin in macrophages in vivo, an endoglin myeloid lineage specific knock-out mouse line (Engfl/flLysMCre) was generated. These mice show a predisposition to develop spontaneous infections by opportunistic bacteria. Engfl/flLysMCre mice also display increased survival following LPS-induced peritonitis, suggesting a delayed immune response. Phagocytic activity is impaired in peritoneal macrophages, altering one of the main functions of macrophages which contributes to the initiation of the immune response. We also observed altered expression of TGF-β1 target genes in endoglin deficient peritoneal macrophages. Overall, the altered immune activity of endoglin deficient macrophages could help to explain the higher rate of infectious diseases seen in HHT1 patients.This work was funded by: Ministerio de Economía y Competitividad of Spain (SAF2011-23475 to LMB; SAF2013-43421-R and SAF2010- 19222 to CB.Peer reviewe

Genotype-phenotype relationships of truncating mutations, p.E297G and p.D482G in bile salt export pump deficiency

Background & Aims: Bile salt export pump (BSEP) deficiency frequently necessitates liver transplantation in childhood. Homozygous p.D482G or p.E297G mutations are associated with relatively mild phenotypes, responsive to surgical interruption of the enterohepatic circulation (siEHC), in contrast to patients with two predicted protein truncating mutations (PPTM). The phenotype of patients with a compound heterozygous genotype of one p.D482G or p.E297G mutation and one PPTM has remained unclear. We aimed to assess their genotype-phenotype relationship. Methods: From the NAPPED database, we selected patients with homozygous p.D482G or p.E297G mutations (BSEP1/1; n=31), with one p.D482G or p.E297G, and one PPTM (BSEP1/3; n=30), and with two PPTMs (BSEP3/3; n=77). We compared presentation, native liver survival (NLS), and effect of siEHC on NLS. Results: The groups had a similar median age at presentation (0.7-1.3 years). Overall NLS at age 10 years was 21% in BSEP1/3 vs. 75% in BSEP1/1 and 23% in BSEP3/3 (P<0.001). Without siEHC in their follow-up, NLS of BSEP1/3 was similar to BSEP3/3 patients, but considerably lower than BSEP1/1 patients (at age 10 years: 38%, 30%, and 71%, resp; P=0.003). After siEHC, BSEP1/3 and BSEP3/3 patients had similarly low NLS, while this was much higher in BSEP1/1 patients (10 years after siEHC, 27%, 14%, and 92%, resp.; P<0.001). Conclusions: BSEP deficiency patients with one p.E297G or p.D482G mutation and one PPTM have a similarly severe disease course and low responsiveness to siEHC as patients with two PPTMs. This identifies a considerable subgroup of patients who are unlikely to benefit from interruption of the enterohepatic circulation by either surgical or ileal bile acid transporter inhibitor treatment

Dynamic endothelial cell rearrangements drive developmental vessel regression

Patterning of functional blood vessel networks is achieved by pruning of superfluous connections. The cellular and molecular principles of vessel regression are poorly understood. Here we show that regression is mediated by dynamic and polarized migration of endothelial cells, representing anastomosis in reverse. Establishing and analyzing the first axial polarity map of all endothelial cells in a remodeling vascular network, we propose that balanced movement of cells maintains the primitive plexus under low shear conditions in a metastable dynamic state. We predict that flow-induced polarized migration of endothelial cells breaks symmetry and leads to stabilization of high flow/shear segments and regression of adjacent low flow/shear segments.status: publishe