changelog

采用FDTD时域有限差分法计算电磁波的传播问题等。

  * initialize mpi(argc,argv) had to be modified to pass those arguments as *references*, since MPI_Init may modify the argument list (fortunately, no changes are required to user code)

    M ./src/meep/mympi.hpp -1 +1
    M ./src/mympi.cpp -1 +1

Mon Jul 25 20:30:59 EDT 2005  [e-mail]
  * typo (wrong return type)

    M ./src/bicgstab.cpp -1 +1

Mon Jul 25 20:30:49 EDT 2005  [e-mail]
  * wrong format string now that m is double

    M ./src/bands.cpp -1 +1

Mon Jul 25 17:43:23 EDT 2005  [e-mail]
  * bugfix in HDF5+MPI (called H5Fopen instead of H5Fcreate)

    M ./src/h5file.cpp -2 +2

Mon Jul 25 17:42:25 EDT 2005  [e-mail]
  * HDF5 check needs to use MPICXX in case HDF was compiled with MPI (and needs mpi.h for header file)

    M ./configure.ac +5

Mon Jul 18 18:53:51 EDT 2005  [e-mail]
  * allow non-integer values of "m" in cylindrical coordinates (you probably should only do this for cells that don't include r=0)

    M ./src/fields.cpp -7 +6
    M ./src/initialize.cpp -5 +7
    M ./src/meep.hpp -4 +6
    M ./src/step_d.cpp -3 +3
    M ./src/step_h.cpp -3 +3

Fri Jul  8 17:19:02 EDT 2005  [e-mail]
  * slightly weakened test (now that we are truly including multiple source frequencies)

    M ./tests/two_dimensional.cpp -1 +1

Fri Jul  8 16:36:36 EDT 2005  [e-mail]
  * noticed *another* apparent bug in nonlinear code for symmetry, grrr (??); disable test for now

    M ./tests/symmetry.cpp -2 +6

Fri Jul  8 16:36:04 EDT 2005  [e-mail]
  * add more multiple-source tests to symmetry.cpp

    M ./tests/symmetry.cpp +6

Thu Jul  7 20:21:19 EDT 2005  [e-mail]
  * fix by Mihai for case with multiple source times

    M ./src/sources.cpp -1 +1

Wed Jun  1 21:19:31 EDT 2005  [e-mail]
  * another potential bug exposed by my previous patch - polariton_ex test also fails when symmetry is included, so revert with comment for now

    M ./tests/symmetry.cpp -2 +6

Wed Jun  1 20:16:42 EDT 2005  [e-mail]
  * can't use "darcs changes" to generate ChangeLog during 'make dist', because then 'darcs check' will fail (doesn't run in darcs dir); instead, add new 'make darcs-dist' target that creates the ChangeLog

    M ./ChangeLog -1 +1
    M ./Makefile.am -3 +7

Wed Jun  1 20:10:55 EDT 2005  [e-mail]
  * silence warning

    M ./src/monitor.cpp +1

Wed Jun  1 20:10:31 EDT 2005  [e-mail]
  * silence C++ warning

    M ./src/bicgstab.cpp -2 +4

Wed Jun  1 19:46:28 EDT 2005  [e-mail]
  * Generate ChangeLog automatically during 'make dist'

    M ./ChangeLog -1 +1
    M ./Makefile.am +5

Wed Jun  1 14:15:23 EDT 2005  [e-mail]
  * add equal_layout method to check whether two fields/structures have same
  chunking, etcetera; testcase in symmetry.cpp
  
  also add == for symmetry and signed_direction, and != for various things
  that already have ==.
  
  also fixed bug in symmetry.cpp (some fields were being compared to
  themselves rather than to symmetrized version)
  
  NOTE: THERE IS STILL A BUG IN symmetry.cpp ... nonlinear_ex() compares
  the field to itself rather than to a symmetrized version, and fails if
  the correct thing is done, which means that there may be a bug in the
  nonlinear code for symmetry

    M ./src/fields.cpp +16
    M ./src/meep.hpp +3
    M ./src/meep/vec.hpp +7
    M ./src/structure.cpp +13
    M ./src/vec.cpp +11
    M ./tests/symmetry.cpp -3 +36

Wed Jun  1 12:46:06 EDT 2005  [e-mail]
  * pass chunk index to chunkloop function, so that in the future we can easily write functions that e.g. compare two different fields objects

    M ./src/dft.cpp -1 +1
    M ./src/energy_and_flux.cpp -1 +1
    M ./src/h5fields.cpp -2 +2
    M ./src/integrate.cpp -1 +1
    M ./src/loop_in_chunks.cpp -1 +1
    M ./src/meep.hpp -1 +1
    M ./src/sources.cpp -1 +1

Wed Jun  1 00:50:53 EDT 2005  [e-mail]
  * current unit fix: field update should add current*dt, not current, to fields

    M ./src/step.cpp -2 +2

Wed Jun  1 00:33:49 EDT 2005  [e-mail]
  * minor update to bicgstab comments

    M ./src/bicgstab.cpp -2 +2

Wed Jun  1 00:27:39 EDT 2005  [e-mail]
  * whoops, forgot to return a value from solve_cw

    M ./src/cw_fields.cpp -1 +1

Tue May 31 20:20:39 EDT 2005  [e-mail]
  * need string.h in bicgstab.cpp for memset()

    M ./src/bicgstab.cpp +1

Tue May 31 20:05:00 EDT 2005  [e-mail]
  * to get current(t) from dipole(t), use forward rather than backward difference (makes slightly more sense with solve_cw since many currents turn on at t = 0, and backwards difference would then give undesired spike at t=0)

    M ./src/meep.hpp -1 +1
    M ./src/meep_internals.hpp -2

Tue May 31 19:55:04 EDT 2005  [e-mail]
  * get sources directly in solve_cw; add new step_d_source to sidestep current hacks

    M ./src/cw_fields.cpp -10 +5
    M ./src/meep.hpp -1 +3
    M ./src/step.cpp -5 +25

Tue May 31 19:24:53 EDT 2005  [e-mail]
  * use solve_cw in physical.cpp instead of timestepping

    M ./tests/physical.cpp -53 +21

Tue May 31 19:24:16 EDT 2005  [e-mail]
  * put -J on right hand side of solve_cw

    M ./src/cw_fields.cpp -3 +3

Tue May 31 18:36:51 EDT 2005  [e-mail]
  * yikes, silly bug in fields::solve_cw

    M ./src/cw_fields.cpp -1 +1

Tue May 31 18:17:47 EDT 2005  [e-mail]
  * allow CW sources to have complex frequencies (e.g. exponentially decaying)

    M ./src/meep.hpp -7 +9
    M ./src/sources.cpp -11 +3

Tue May 31 18:07:17 EDT 2005  [e-mail]
  * increase maxiters default to 10000 in solve_cw

    M ./src/meep.hpp -1 +1

Tue May 31 18:05:13 EDT 2005  [e-mail]
  * add solve_cw() variant that gets frequency from sources (and has reasonable defaults for the tolerance and maximum iterations)

    M ./src/cw_fields.cpp +15
    M ./src/meep.hpp +4

Tue May 31 17:52:24 EDT 2005  [e-mail]
  * don't output so often in bicgstab

    M ./src/bicgstab.cpp -1 +7

Tue May 31 17:49:00 EDT 2005  [e-mail]
  * slight tweak in bicgstab

    M ./src/bicgstab.cpp +5

Tue May 31 17:26:11 EDT 2005  [e-mail]
  * in fields::solve_cw, use the BiCGSTAB(L) algorithm instead of BiCGSTAB.  With L >= 2, this avoids a known convergence failure in BiCGSTAB for operators (like ours) where the eigenvalues are nearly purely imaginary.  Now, solve_cw converges to arbitrary tolerance (limited by numerical precision, of course) even in > 1d!  Hooray, it works!

    M ./configure.ac -1 +1
    M ./src/bicgstab.cpp -50 +167
    M ./src/bicgstab.hpp -6 +24
    M ./src/cw_fields.cpp -6 +18
    M ./src/meep.hpp -1 +1

Tue May 31 01:12:22 EDT 2005  [e-mail]
  * whoops, forgot to add cw_fields.cpp

    M ./src/bicgstab.cpp -1 +2
    M ./src/bicgstab.hpp -1 +2
    A ./src/cw_fields.cpp

Mon May 30 22:35:05 EDT 2005  [e-mail]
  * remove (commented-out) stabilized dnrm2, which doesn't parallelize anyway

    M ./src/bicgstab.cpp -22 +1

Mon May 30 22:34:10 EDT 2005  [e-mail]
  * whoops

    M ./src/bicgstab.cpp -1 +1

Mon May 30 22:23:05 EDT 2005  [e-mail]
  * during phasing, don't forget to update the E boundaries as well as E

    M ./src/step.cpp -4 +7

Mon May 30 22:20:06 EDT 2005  [e-mail]
  * add new fields::solve_cw method that, instead of timestepping, treats the whole FDTD system as a big linear equation (specified implicitly by the fields::step) and solves it by the stabilized bi-conjugate gradient method to get the CW response at a particular frequency.  It seems to work reasonably well, except for some reason I can't fathom you can't set the tolerance very low (less then 1e-2) or it never converges in 2+ dimensions

    M ./configure.ac -1 +1
    M ./src/Makefile.am -7 +7
    A ./src/bicgstab.cpp
    A ./src/bicgstab.hpp
    M ./src/fields.cpp +2
    M ./src/meep.hpp -5 +13
    M ./src/sources.cpp -8 +12
    M ./src/step.cpp -1 +21
    M ./src/update_e_from_d.cpp -1 +5

Thu May 19 15:21:05 EDT 2005  [e-mail]
  * add FIXME comment to unnecessary allocation

    M ./src/fields.cpp +1

Tue May 17 16:57:48 EDT 2005  [e-mail]
  * fixed slight bug ... volume::nowned() really returned the minimum nowned over all components, which is slightly less than the actual nowned for some components in cylindrical coords; changed the name to nowned_min() to reflect this, added a new nowned(component) that gives the exact nowned, and added a regression test

    M ./src/meep/vec.hpp -1 +2
    M ./src/vec.cpp -4 +11
    M ./tests/integrate.cpp +3

Tue May 17 16:19:33 EDT 2005  [e-mail]
  * rename all C++ headers to .hpp, for consistency with .cpp and for clarity

     ./src/meep.h -> ./src/meep.hpp
     ./src/meep/mympi.h -> ./src/meep/mympi.hpp
     ./src/meep/vec.h -> ./src/meep/vec.hpp
     ./src/meep_internals.h -> ./src/meep_internals.hpp
     ./src/ran.h -> ./src/ran.hpp
     ./src/threevec.h -> ./src/threevec.hpp
    M ./doc/complicated.tex -1 +1
    M ./doc/energy_cons.tex -1 +1
    M ./doc/energy_cons_1d.tex -1 +1
    M ./doc/epsilon_polariton_1d.tex -1 +1
    M ./doc/example.cpp -1 +1
    M ./doc/lossgain_epsilon.tex -1 +1
    M ./doc/nonlinear.tex -1 +1
    M ./doc/omniguide.tex -1 +1
    M ./doc/polaritonbands.tex -1 +1
    M ./doc/simple.tex -1 +1
    M ./doc/simplebands.tex -1 +1
    M ./src/Makefile.am -12 +13
    M ./src/anisotropic_averaging.cpp -2 +2
    M ./src/bands.cpp -2 +2
    M ./src/boundaries.cpp -2 +2
    M ./src/control_c.cpp -1 +1
    M ./src/dft.cpp -2 +2
    M ./src/energy_and_flux.cpp -2 +2
    M ./src/fields.cpp -2 +2
    M ./src/grace.cpp -2 +2
    M ./src/h5fields.cpp -1 +1
    M ./src/h5file.cpp -1 +1
    M ./src/initialize.cpp -2 +2
    M ./src/integrate.cpp -2 +2
    M ./src/loop_in_chunks.cpp -2 +2
    M ./src/meep.hpp -3 +2
    M ./src/meep_internals.hpp -2 +1
    M ./src/monitor.cpp -2 +2
    M ./src/mympi.cpp -1 +1
    M ./src/output_directory.cpp -1 +1
    M ./src/polarization.cpp -2 +2
    M ./src/ran.cpp -1 +1
    M ./src/ran.hpp -1
    M ./src/slices.cpp -2 +2
    M ./src/sources.cpp -2 +2
    M ./src/step.cpp -2 +2
    M ./src/step_d.cpp -4 +4
    M ./src/step_h.cpp -4 +4
    M ./src/structure.cpp -2 +2
    M ./src/threevec.hpp -1
    M ./src/time.cpp -1 +1
    M ./src/update_e_from_d.cpp -2 +2
    M ./src/update_e_from_d_prepare.cpp -4 +4
    M ./src/update_e_from_d_sources.cpp -3 +3
    M ./src/update_e_from_d_update.cpp -3 +3
    M ./src/update_from_e.cpp -4 +4
    M ./src/vec.cpp -1 +1
    M ./tests/2D_convergence.cpp -1 +1
    M ./tests/bench.cpp -1 +1
    M ./tests/bragg_transmission.cpp -1 +1
    M ./tests/convergence_cyl_waveguide.cpp -1 +1
    M ./tests/cylindrical.cpp -1 +1
    M ./tests/flux.cpp -1 +1
    M ./tests/h5test.cpp -2 +2
    M ./tests/integrate.cpp -1 +1
    M ./tests/known_results.cpp -1 +1
    M ./tests/one_dimensional.cpp -1 +1
    M ./tests/physical.cpp -1 +1
    M ./tests/symmetry.cpp -1 +1
    M ./tests/three_d.cpp -1 +1
    M ./tests/two_dimensional.cpp -1 +1

Mon May 16 18:52:39 EDT 2005  [e-mail]
  * don't re-allocate d_minus_p array at each timestep

    M ./src/fields.cpp +5
    M ./src/meep.h -3 +6
    M ./src/update_e_from_d.cpp -11 +6
    M ./src/update_e_from_d_prepare.cpp -2 +1
    M ./src/update_e_from_d_sources.cpp -2 +1
    M ./src/update_e_from_d_update.cpp -2 +1

Mon May 16 18:03:20 EDT 2005  [e-mail]
  * add possibility to change PML strength, e.g. pml(X) * 2.0 will increase PML conductivity by 2 (squaring attenuation for same thickness, in theory, albeit increasing numerical reflection)

    M ./src/meep.h -10 +17
    M ./src/structure.cpp -6 +11

Mon May 16 15:45:01 EDT 2005  [e-mail]
  * comment regarding accuracy of modal volume calcs

    M ./src/energy_and_flux.cpp -2 +14

Mon May 16 14:58:27 EDT 2005  [e-mail]
  * Darcs is now right-handed (there was a bug in 1d/2d/3d that caused the curls to have the wrong signs, making the system "left-handed" and giving fluxes the wrong signs, although it didn't affect most other calculations)

    M ./src/fields.cpp -1 +1
    M ./tests/flux.cpp -6 +6

Mon May 16 13:36:37 EDT 2005  [e-mail]
  * work around datatype change in HDF5 1.6.4 (grrr)

    M ./src/h5file.cpp -2 +13

Mon May 16 13:29:19 EDT 2005  [e-mail]
  * eliminate unnecessary allocation of PML fields

    M ./src/energy_and_flux.cpp -8 +14
    M ./src/fields.cpp -9 +13

Wed Apr 27 15:33:47 EDT 2005  [e-mail]
  * whoops, fix compilation failure for --without-hdf5

    M ./src/h5file.cpp +2

Mon Apr 25 17:25:20 EDT 2005  [e-mail]
  * fixed incorrect types in hdf5 calls

    M ./src/h5file.cpp -2 +2

Mon Apr 25 16:33:34 EDT 2005  [e-mail]
  * print a little more information from "aack, odd number of grid points"

    M ./src/structure.cpp -1 +2

Mon Apr 25 16:31:35 EDT 2005  [e-mail]
  * when phasing epsilon, use update_e_from_d to keep E in synch, rather than "manual" division and multiplication by inveps diagonal"

    M ./src/step.cpp -31 +1

Mon Apr 25 16:30:55 EDT 2005  [e-mail]
  * order spherical quadrature points to maximize spatial separation (to make interface detection more efficient, in future)

    M ./src/sphere-quad.cpp -7 +49

Mon May  2 10:44:24 EDT 2005  [e-mail]
  * fix bug when running in cylindrical coordinates with real fields.

    M ./src/step_d.cpp -1 +3
    M ./src/step_h.cpp -1 +3

Sat Feb  5 13:44:55 EST 2005  [e-mail]
  * output s/step along with time-step status

    M ./src/meep.h +1
    M ./src/step.cpp -2 +8
    M ./src/time.cpp -1 +1