Computational Modelling Group


The GNU Scientific Library (GSL) is a numerical library for C and C++ programmers. It is free software under the GNU General Public License.

The library provides a wide range of mathematical routines such as random number generators, special functions and least-squares fitting. There are over 1000 functions in total with an extensive test suite.

Unlike the licenses of proprietary numerical libraries the license of GSL does not restrict scientific cooperation. It allows you to share your programs freely with others.

For more information, follow this link:

For queries about this topic, contact Sam Dolan.


B-meson coupling with relativistic heavy quarks

Jonathan Flynn (Investigator), Ben Samways, Dirk Broemmel, Patrick Fritzsch

We non-perturbatively compute the coupling between B* and B pi meson states relying on relativistic heavy quarks and domain wall light fermions. The coupling is of importance for an effective description of hadronic heavy meson decays.

Hadronic structure on the computer

Jonathan Flynn (Investigator), Dirk Broemmel, Thomas Rae, Ben Samways

In experiments at the Large Hadron Collider (LHC) at CERN, Geneva, the interactions that occur between the colliding particles (protons in this case) can be factorised into a simple scattering between two constituent particles, called quarks, followed by a hadronisation process, which describes the dynamics of forming the bound proton states. Quarks are particles within the proton that bind to form composite particles (hadrons) such as a proton. The scattering process can be computed relatively easily, but hadronisation is intrinsically non-perturbative and hard to calculate. Lattice QCD (computer simulation of QCD on a discrete space-time lattice) provides our only known first-principles and systematically-improvable method to address problems like hadronisation. This project uses Iridis to extract parton distribution amplitudes which are experimentally inaccessible, but needed to describe the quark structure of hadrons.

Kaon to two pion decays in lattice QCD

Jonathan Flynn (Investigator), Elaine Goode, Dirk Broemmel

We calculate kaon decay amplitudes on the lattice so we may compare the Standard Model to experiment.

Non-Perturbative Renormalisation on the Lattice

Jonathan Flynn (Investigator), Dirk Broemmel, Thomas Rae

In this project we compute renormalisation factors for various physical observables in a non-perturbative lattice framework. Renormalisation hereby arises due to a fundamental scale dependence of the physical processes.

Self-Force and Black Hole Inspirals

Sam Dolan (Investigator)

We use IRIDIS to compute the self-force acting on a solar-mass black hole orbiting a supermassive black hole.

Structured low-rank approximation

Ivan Markovsky

Today's state-of-the-art methods for data processing are model based. We propose a fundamentally new approach that does not depend on an explicit model representation and can be used for model-free data processing. From a theoretical point of view, the prime advantage of the newly proposed paradigm is conceptual unification of existing methods. From a practical point of view, the proposed paradigm opens new possibilities for development of computational methods for data processing.


Jonathan Flynn
Professor, Physics & Astronomy (FPAS)
Ivan Markovsky
Lecturer, Electronics and Computer Science (FPAS)
Dirk Broemmel
Research Fellow, Physics & Astronomy (FPAS)
Sam Dolan
Research Fellow, Mathematics (FSHS)
Elaine Goode
Postgraduate Research Student, Physics & Astronomy (FPAS)
Thomas Rae
Postgraduate Research Student, Physics & Astronomy (FPAS)
Ben Samways
Postgraduate Research Student, Physics & Astronomy (FPAS)
Jess Jones
Technical Staff, iSolutions
Petrina Butler
Administrative Staff, Research and Innovation Services
Andrew Penner
Alumnus, Mathematics (FSHS)