jan
/
opencl_fdtd
2
2
Fork 1
Code Issues Pull Requests Releases Wiki Activity
OpenCL FDTD electromagnetic simulation in 3 dimensions
You cannot select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
56 Commits
6 Branches
0 Tags
60 MiB
Python 73.8%
C 26.2%
 
 
ongoing
master
new_poynting
image
array_pml
lookup_tables
Branches Tags
${ item.name }
Create tag ${ searchTerm }
Create branch ${ searchTerm }
from 'master'
${ noResults }
Go to file
Jan Petykiewicz d5fd78d493 Revert "Add restrict keyword to pointers (not sharing the same memory for multiple fields)" 
This reverts commit 60b70bb332.

It appears to have minimal performance impact, and fails to compile on
nvidia cards.
 		4 years ago
opencl_fdtd Revert "Add restrict keyword to pointers (not sharing the same memory for multiple fields)" 4 years ago
.gitignore expand gitignores 4 years ago
LICENSE.md Use markdown for license 7 years ago
README.md move code to new location 5 years ago
fdtd.py Update example with bloch fields 5 years ago
pcgen.py Fix triangular lattice code in pgcgen 7 years ago
requirements.txt move code to new location 5 years ago
setup.py Use readme as long_description 5 years ago

README.md

opencl_fdtd

opencl_fdtd is a python application 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 91 iterations/sec. on my AMD RX480.
  • On an Nvidia GTX 580, it runs at 66 iterations/sec
  • On my laptop (Nvidia 940M) the same simulation achieves ~12 iterations/sec.
  • An L3 photonic crystal cavity ringdown simulation (1550nm source, 40nm discretization, 8000 steps) takes about 3 minutes on my laptop.

Capabilities are currently pretty minimal:

  • Absorbing boundaries (CPML)
  • Perfect electrical conductors (PECs; to use set epsilon to inf)
  • Anisotropic media (eps_xx, eps_yy, eps_zz, mu_xx, ...)
  • Direct access to fields (eg., you can trivially add a soft or hard current source with just sim.E[ind] += sin(f0 * t), or save any portion of a field to a file)

Installation

Requirements:

  • python 3 (written and tested with 3.5)
  • numpy
  • pyopencl
  • jinja2
  • fdfd_tools

Optional (used for examples):

To get the code, just clone this repository:

git clone https://mpxd.net/code/jan/opencl_fdtd.git

You can install the requirements and their dependencies easily with

pip install -r requirements.txt

Running

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

python3 fdtd.py