Hi all,
I looked at meep.hpp and meep_wrap.cxx files, but couldn't figure out
how to call do-harminv directly. can someone let me know how I can use
do-harminv in a *.ctl file to perform the same function as (given, let's
say, Hz at (0,0,0))
(run-k-point T k)
, except that Q_thresh is different from 50.
Thank you,
Changsoon
Date: Tue, 14 Nov 2006 00:13:02 -0500 (EST)
From: "Steven G. Johnson" <[EMAIL PROTECTED]>
Subject: Re: [Meep-discuss] -Q option for harminv
To: [email protected]
Message-ID: <[EMAIL PROTECTED]>
Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed
On Mon, 13 Nov 2006, Changsoon Kim wrote:
I am calculating a band structure for a photonic crystal slab, and the
results do not show some of high frequency leaky guided modes. I was
wondering whether there is a way to set -Q option for harminv from a *.ctl
file.
It is possible, by calling do-harminv directly and passing a different
Q_thresh parameter (the default is 50) (see do_harminv in meep.hpp). Or
you could edit meep.hpp to change the default Q_thresh parameter for
do_harminv.
However, I would be suspicious of Q < 50 modes in general; not only is the
lifetime too short to extract them accurately, but even the leaky-mode
approximation of the modes as perfectly exponentially decaying begins to
break down.
Steven
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
