Electroanatomic mapping systems collect increasingly large quantities of spatially-distributed electrical data which may be potentially further scrutinized post-operatively to expose mechanistic properties which sustain and perpetuate atrial fibrillation. We describe a modular software platform, developed to post-process and rapidly analyse data exported from electroanatomic mapping systems using a range of existing and novel algorithms. Imaging […]
An Automated Algorithm for Determining Conduction Velocity, Wavefront Direction and Focal Source Location using a Multipolar Catheter
Determining locations of focal arrhythmia sources and quantifying myocardial conduction velocity (CV) are two major challenges in clinical catheter ablation cases. CV, wavefront direction and focal source location can be estimated from multipolar catheter data, but currently available methods are time-consuming, limited to specific electrode configurations, and can be inaccurate. We developed automated algorithms to […]
From h to p efficiently: Optimal implementation strategies for explicit time-dependent problems using the spectral/hp element method
We investigate the relative performance of a second-order Adams-Bashforth scheme and second- and fourth-order Runge-Kutta schemes when time-stepping a 2D linear advection problem discretised using a spectral/hp element technique for a range of different mesh sizes and polynomial orders. Numerical experiments explore the effects of short (2 wavelengths) and long (32 wavelengths) time integration for […]
Simplifying the Development, Use and Sustainability of HPC Software
Developing software to undertake complex, compute-intensive scientific processes requires a challenging combination of both specialist domain knowledge and software development skills to convert this knowledge into efficient code. As computational platforms become increasingly heterogeneous and newer types of platform such as Infrastructure-as-a-Service (IaaS) cloud computing become more widely accepted for HPC computations, scientists require more […]
High-order spectral/hp element discretisation for reaction-diffusion problems on surfaces: Application to cardiac electrophysiology
We present a numerical discretisation of an embedded two-dimensional manifold using high-order continuous Galerkin spectral/hp elements, which provide exponential convergence of the solution with increasing polynomial order, while retaining geometric flexibility in the representation of the domain. Our work is motivated by applications in cardiac electrophysiology where sharp gradients in the solution benefit from the […]
Creating histograms in gnuplot
I needed to produce a histogram in gnuplot of multiple sets of raw data (i.e. not sumarised). The data in each dataset required binning and displaying on the same plot.
Nekkloud: A software environment for high-order finite element analysis on clusters and clouds
As the capabilities and diversity of computational platforms continue to grow, scientific software is becoming ever more complex in order to target resources effectively. In the libhpc project we are developing a suite of tools and services to simplify job description and execution on heterogeneous infrastructures.
Severe gap junctional uncoupling causes electrogram fractionation in a novel 1D numerical model of discrete myocytes
Determining the resolution of data required for tracking rotors in a downsampled computational model of atrial fibrillation
Finite element assembly strategies on multi-core and many-core architectures
We demonstrate that radically differing implementations of finite element methods (FEMs) are needed on multi-core (CPU) and many-core (GPU) architectures, if their respective performance potential is to be realised. Our numerical investigations using a finite element advection–diffusion solver show that increased performance on each architecture can only be achieved by committing to specific and diverse […]