1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50

HDF5 file output routines for frequencydependent permittivity and permeability.
Skewed bloch periodic boundary conditions for nonorthogonal unit cells
Discrete rotational symmetry in cylindrical coordinates.
Gyromagnetic materials via general anisotropic conductivities.
Secondorder accurate subpixel smoothing scheme for
dispersive materials involving complex permittivity.
Finish implementation of active materials (multilevel atoms and
saturable absorption/gain).
Support 3d (r,phi,z) cell, where m is "Bloch wavenumber"
 gives 6fold (and nfold) symmetry as side effect
 allows computation of bent photoniccrystal waveguides
Implement transparent boundary conditions to back PML layers, so that
glancingangle waves can be absorbed.
Reentrant functions: foo(vec&) arguments should be foo(vec&, void*)
where the second argument can be used to pass state (instead of global
vars).
User proper PML for cylindrical coordinates, rather than quasiPML.
Use morestable algorithm for dispersive media.
Make sure epsilon and other material properties respect symmetry,
periodic boundaries, etc.?
Remove monitor point class (redundant with HDF5 output, DFT volumes?).
Rename fields::initialize to fields::add_to_fields or something like that,
and allow specifying a geometric_volume.
Support Pade approximates (libpadespectrum) for Fourier and modal analysis.
Check sensitivity to rounding error and if there is a better way to handle:
structure.cpp:555 (PML region boundaries)
vec.cpp:560 (interpolation weights)
anisotropic_averaging: 79 (magnitude of normal vector)

Python & Scheme Interface
Need a way to reset epsilon when restarting, and also phasing support.
Support arbitraryshaped sources.
