A microbenchmark measuring the performance of various programming languages and parallelisation strategies by solving a matrix-free linear system using an iterative Jacobi method. Essentially solves the diffusion equation. Zig and C are bloody fast.

#C++ #Zig #C #Julia #Python

A pseudospectral solver framework designed to simulate 2-dimensional advection-diffusion problems. It is written entirely in Python and leverages CuPy for GPU parallelisation. It has been used to simulate thermohaline staircase formation, Rayleigh-Bénard convection, and the resistive tearing instability.

#Python #Fluid-Dynamics #CuPy #GPU #Simulation #Spectral-Method

Custom renderer which uses radiosity to implement global illumination. Uses OpenGL. Developed at ICHEC, Dublin as part of the PRACE Summer of HPC.

#C++ #OpenGL

Starting code for a parallel programming challenge I ran in 2022. Includes simulations of a general N-body problem, and a fully implicit fluid solver simulating Rayleight-Benard convection. The event taught students how to program GPUs and multi-core CPUs in Python.

#Python #Simulation #Fluid-Dynamics #Finite-Difference

Tool to calculate the integration of various quantities along magnetic field lines using Mayavi. Consumes SDF files.

#Python #Research

Implementation of Stam's Real-Time Fluid Dynamics for Games in OpenCL. Designed to run on FPGAs but only CPU and GPU were tested.

#C++ #Opencl #Fpga #Fluid-Dynamics #Simulation #Finite-Difference

Implementation of a finite-volume Lagrangian-Remap scheme in 1D with shock capturing. Used as a toy model for the numerical methods section of my PhD thesis.

#C++ #Research #Fluid-Dynamics #Simulation

My version of the finite-volume magnetohydrodynamics code Lare3d including various implementations of anisotropic viscosity.

#Finite-Volume #Fortran #Research #Simulation #Fluid-Dynamics

Earlier version of Melvin.py written in C++ and parallelised with OpenMP and CUDA.

#C++ #Fluid-Dynamics #CUDA #OpenMP #GPU #Simulation #Spectral-Method

Implementation of the remove-restore algorithm for combining bathymetry and topography datasets. Developed with Dr Alessandro Felder at UCL.

#Python #Research

Simple bytebeat synth written in C++ with multiple filters and waveforms. Bytebeat synths write 8-bit values to stdout which can be consumed by, for example, `aplay`.

#Music #C++

Simulate the interaction between molecules GFP-CRIB and Rac in 2 or 3 dimensions. Core numerical solver is a custom finite-element method. Developed for the Beatson Institute for Cancer Research. Technically works, but not particularly readable.

#Research #Python #Simulation #Finite-Element

Trade Pokemon between Nintendo DS and a PC

#Python

Web scraper for the Journal of Open Source Software. Used in the associated blog post to analyse how many published projects use continuous testing.

#Python

Tool for manipulating SDF (self-describing format) files to make analysis easier.

#Fortran

An experiment in producing art with only 4 lines.

#Art #Javascript

Simulation of a double pendulum system in javascript with various options for creating interesting visualisations.

#Art #Javascript #Simulation

Visualisation of the Lorenz attractor with various options for making pretty pictures.

#Javascript #P5.js #Simulation #Art

Generate fractals from regular expressions in javascript.

#Javascript #Fractal #Art

Generate Julia fractals with WebGL

#Javascript #Webgl #Fractal #Art

Incomplete implementation of Lorena Barba's 12 steps to Navier-Stokes using WebGL.

#Javascript #Fluid-Dynamics #Webgl

The Latex files and Python scripts used to render my PhD thesis. Includes all analysis code.

#Research #Python #Latex

Facts about my friend Ross that may or may not be true.

#Javascript

Implementation of a numerical test of the Collatz Conjecture in C++, parallelised with OpenMP, std::threads and Boost Compute.

#C++ #OpenMP

Analysis code for my paper "The effect of anisotropic viscosity on the nonlinear kink instability"

#Python #Analysis #Research

Analysis code for my paper "Flute and kink instabilities in a dynamically twisted flux tube with anisotropic plasma viscosity".

#Research #Analysis #Python

Analysis code for my paper "Kelvin-Helmholtz instability and collapse of a twisted magnetic nullpoint with anisotropic viscosity"

#Research #Python #Analysis