Although Meep computes Fourier transforms of field components in space
internally, it doesn't currently give you an easy way to output this
information, I'm afraid. (It is mainly used right now for computing flux
spectra.)
The Fourier transform of a region in space is computed when you all
fields::add_dft to create a "dft_chunk" data structure that is updated as
you time step. At any point in time thereafter, the dft_chunk::dft member
points to an array of the Fourier amplitudes that have been currently
accumulated for that region of space. This is complicated by the fact
that the dft_chunk in general is a linked list, corresponding to the DFTs
for each "chunk" of the fields (fields arrays are broken into chunks to
deal with different processors, PML regions, etcetera).
In short, the information is all there, but it requires dealing with the
C++ internals at a fairly low level right now. If the region of interest
is not too big, it is probably easier to just output the fields in that
region as a function of time and Fourier transform them yourself
afterwards.
On Thu, 18 May 2006, Frank Rodjius wrote:
It seems meep can give field profile in space as Ez(t,x,y) for a source.
Would you please let me know if I can get the Ez(w, x,y) from meep? How
to do that? I mean I would like to see the field profile in space for a
frequency w. Thank you very much
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
