could it be due to the frequency(for the gaussion source) going into negative
values? (if frequency is minus it wont be a plane wave anymore)
cheers!
Gursoy Akguc
________________________________
From: Alex <[email protected]>
To: [email protected]
Sent: Wednesday, March 20, 2013 12:39 PM
Subject: [Meep-discuss] Transmission spectrum of free space
Dear meep users and developers.
I try to observe transmission spectrum of 1D structure. To check
correctness of fluxes output I use free space (no geometry set). I
guess spectrum of incoming flux and outcoming flux should be equal.
But meep prints different values for a same frequency in incoming and
outcoming fluxes.
Where is my mistake?
Alex Friman
MEEP code:
(set! geometry-lattice (make lattice (size 1 1 1)))
(set! resolution 32)
(define-param fcen 0.5)
(define-param df 1.5)
(set! sources (list (make source (src (make gaussian-src (frequency
fcen) (fwidth df))) (component Hz) (center 0.5 0))))
(define-param nfreq 400)
(define-param fmin (- fcen (/ df 2)))
(define-param fmax (+ fcen (/ df 2)))
(define outcome
(add-flux (* 0.5 (+ fmin fmax)) (- fmax fmin) nfreq
(make flux-region (center 0.9 0.5 0.5) (size 0 1 1) )))
(define income
(add-flux (* 0.5 (+ fmin fmax)) (- fmax fmin) nfreq
(make flux-region (center 0 0.5 0.5) (size 0 1 1) )))
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ey)
(center 0.0 0.5 0.5)
(size 0 1 1))))
(run-sources+ 20)
(display-fluxes income outcome)
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
