*DEPRECATED* Tools for optical simulations
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.
Go to file
Jan Petykiewicz 03fc9e6d70 deprecate 4 years ago
examples add bloch example 7 years ago
fdfd_tools fix operator test 6 years ago
.gitignore initial development 8 years ago
LICENSE.md add license 8 years ago
README.md deprecate 4 years ago
setup.py bump version number 7 years ago




The functionality in this module is now provided by meanas.

fdfd_tools is a python package containing utilities for creating and analyzing 2D and 3D finite-difference frequency-domain (FDFD) electromagnetic simulations.


  • Library of sparse matrices for representing the electromagnetic wave equation in 3D, as well as auxiliary matrices for conversion between fields
  • Waveguide mode solver and waveguide mode operators
  • Stretched-coordinate PML boundaries (SCPML)
  • Functional versions of most operators
  • Anisotropic media (eps_xx, eps_yy, eps_zz, mu_xx, ...)
  • Arbitrary distributions of perfect electric and magnetic conductors (PEC / PMC)

This package does not provide a fast matrix solver, though by default fdfd_tools.solvers.generic(...) will call scipy.sparse.linalg.qmr(...) to perform a solve. For 2D problems this should be fine; likewise, the waveguide mode solver uses scipy's eigenvalue solver, with reasonable results.

For solving large (or 3D) problems, I recommend a GPU-based iterative solver, such as opencl_fdfd or those included in MAGMA). Your solver will need the ability to solve complex symmetric (non-Hermitian) linear systems, ideally with double precision.



  • python 3 (written and tested with 3.5)
  • numpy
  • scipy

Install with pip, via git:

pip install git+https://mpxd.net/gogs/jan/fdfd_tools.git@release


See examples/test.py for some simple examples; you may need additional packages such as gridlock to run the examples.