Biomedical Photonic Imaging (BMPI)

Running even larger wave propagation simulations made possible with latest developments to WaveSim published in article 

Wednesday 24 December 2025

Swapnil has done it to get his supergood WaveSim article published!! That is titled:
Domain decomposition of the modified Born series approach for large-scale wave propagation simulations.
We are all so proud!

Abstract

WaveSim is a free and open-source solver developed in-house at BMPI for simulating wave propagation in large and complex structures. We have been working on improving and accelerating our fast and accurate wave propagation solver, and recently, in December 2025, we published a significant addition to WaveSim’s capabilities. 

Using a strategy called domain decomposition to split the simulations over multiple processors or graphics processing units (GPUs), WaveSim can now simulate light propagation at much larger sizes. This new method maintains the accuracy, memory efficiency, and the guarantee of converging with every iteration of the original method without domain decomposition. 

To demonstrate the new improvements, we ran a simulation with an unprecedented 2.1 billion voxels (320×320×320 wavelengths, sampled at 4 points per wavelength) in just 45 minutes on two A40 GPUs. This was over 4500 times the size of a recent simulation using the finite difference time domain method, which is very popular in light propagation simulations. You can try WaveSim for yourself, it is already available on GitHub. 

The below image shows a small 2D slice of 100×100 wavelengths from the large 3D problem. Here the light source is a plane wave propagating from the left side of the x axis through a collection of small spheres. The image is showing 160000 voxels, just 0.007% of the total 2.1 billion that WaveSim simulated! 

You can find the link for the super article here!