On Thu, 21 Sep 2006, matt wrote:
If you're using run-k-points on a unit cell, you want to be able to tell how
many time steps is long enough until the characteristic modes have settled.
For this purpose, it would be nice to watch an animation of the DFT evolution
for a given field point.
The DFT may not be the best thing for this. The reason is that Meep does
not use the DFT in order to compute the modes; instead, it uses harminv's
filter-diagonalization implementation, which has better resolution than
the DFT.
If you want to look at the convergence of harminv's frequencies, the best
thing to do is to look at the harminv output itself.
To do this, I can imagine adding the output of (display-fluxes) to an h5
file every certain number of time steps, just because h5 files are
convenient for multidimensional data. The question then is how you get
matlab to turn the h5 data into an animation of the spectrum.
There are plenty of animation methods in Matlab; google "matlab animation"
to find them. If you have multidimensional data, however, I usually find
it easier to convert them to png files with h5topng and then to use
convert (from ImageMagick) to combine the png files into an mpeg or
whatever.
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
