Jan's OpenCL FDTD electromagnetic simulation in 3 dimensions
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
Go to file
jan cbbc4e294a
gitignore changes
7 years ago
fdtd Initial commit 7 years ago
.gitignore gitignore changes 7 years ago
LICENSE.txt Add AGPL license 7 years ago
README.md Instructions for runnning 7 years ago
fdtd.py Initial commit 7 years ago
pcgen.py Initial commit 7 years ago
requirements.txt Initial commit 7 years ago



opencl-fdtd is a python package for running 3D time-domain electromagnetic simulations on parallel compute hardware (mainly GPUs).

Performance highly depends on what hardware you have available:

  • A 395x345x73 cell simulation (~10 million points, 8-cell absorbing boundaries) runs at around 42 iterations/sec. on my Nvidia GTX 580.
  • On my laptop (Nvidia 940M) the same simulation achieves ~8 iterations/sec.
  • An L3 photonic crystal cavity ringdown simulation (1550nm source, 40nm discretization, 8000 steps) takes about 5 minutes on my laptop.

Capabilities are currently pretty minimal:

  • Absorbing boundaries (CPML)
  • Conducting boundaries (PMC)
  • Anisotropic media (eps_xx, eps_yy, eps_zz)
  • Direct access to fields (eg., you can trivially add a soft or hard current source with just sim.E[1] += sin(f0 * t), or save any portion of a field to a file)



  • python 3 (written and tested with 3.5)
  • numpy
  • pyopencl
  • h5py (for file output)
  • gridlock
  • masque

You can install the requirements and their dependencies easily with

pip intall -r requirements.txt


The root directory contains fdtd.py, which sets up and runs a sample simulation (cavity ringdown).

python3 fdtd.py