Support 3d (r,phi,z) cell, where m is "Bloch wavenumber"      -- gives 6-fold (and n-fold) symmetry as side effect      -- allows computation of bended photonic-crystal waveguidesImplement transparent boundary conditions to back PML layers, so thatglancing-angle waves can be absorbed.Re-entrant functions: foo(vec&) arguments should be foo(vec&, void*)where the second argument can be used to pass state (instead of globalvars).Fix bugs: divergence in saturated polarization (tests/symmetry.cpp),          and problems with nonlinearity + symmetry (idem)User-specified anisotropic dielectric and mu functions.Unify mu and epsilon support in the code, to support dispersive mu etc.Use anisotropic PML rather than split-fields PML (former requires lessmemory).  User proper PML for dispersive media and cylindrical.Use more-stable algorithm for dispersive media.Field dump/load (i.e. restarting)?Make sure epsilon and other material properties respect symmetry,periodic boundaries, etc.?Remove monitor point class (redundant with HDF5 output, DFT volumes?).Integrated support for calling MPB mode solver for current distributions.rename fields::initialize to fields::add_to_fields or something like that,and allow specifying a geometric_volume.Dump field Fourier transforms in a readable, documented format.Support Pade approximates (libpadespectrum) for Fourier and modal analysis.-------------------LIBCTL FRONT-END:Epsilon from file support.Need a way to reset epsilon when restarting, and also phasing support.Support arbitrary-shaped sources.