Re: [Meep-discuss] meep-discuss Digest, Vol 37, Issue 3

2009-03-05 Thread Ian Buss

Hi Nathan,

I had a similar problem when simulating real metals - it turns out  
that Meep will not accept negative permittivities without blowing up.   
You have to sneak up on the negative permittivity by using the  
polarizability feature of materials (see http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial/Material_dispersion) 
.  If you set the infinite permittivity to 1 you can create the  
negative part using the Drude model (see Skinner et al., Appl. Optics  
Vol. 45, No. 9, p 1943, 2006).  Another member of the mailing list  
helped me out with this.


Regards,

Ian


On 5 Mar 2009, at 17:00, meep-discuss-requ...@ab-initio.mit.edu wrote:


Send meep-discuss mailing list submissions to
meep-discuss@ab-initio.mit.edu

To subscribe or unsubscribe via the World Wide Web, visit
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
or, via email, send a message with subject or body 'help' to
meep-discuss-requ...@ab-initio.mit.edu

You can reach the person managing the list at
meep-discuss-ow...@ab-initio.mit.edu

When replying, please edit your Subject line so it is more specific
than Re: Contents of meep-discuss digest...


Today's Topics:

  1. meep installation problem (guile related) (Hua Bao)
  2. Complex field in meep (adrian)
  3.  meep initial field (Michal Martin?)
  4. regarding single pulse source (prabhat behere)
  5. Unstable simulation with evanescent fields in active   medium
 (Nathan)


--

Message: 1
Date: Wed, 04 Mar 2009 22:20:42 -0500
From: Hua Bao h...@purdue.edu
Subject: [Meep-discuss] meep installation problem (guile related)
To: meep-discuss@ab-initio.mit.edu
Message-ID: 49af450a.9070...@purdue.edu
Content-Type: text/plain; charset=ISO-8859-1; format=flowed

Dear All,

I am a new user of meep.
I was trying to build up a meep executable but it shows error at the
configure step.

checking for guile-config... no
checking for readline in -lreadline... no
checking for dlopen in -ldl... yes
checking for scm_lt_dlinit in -lguile-ltdl... no
checking for dlopen in -lltdl... no
checking for gh_eval_str in -lguile... no
configure: error: linking to guile failed

The linux cluster (amd64) has guile installed. There is no problem to
run guile --version. Also, there is a libguile.so under /usr/lib64.

/usr/lib64/libguile-srfi-srfi-1-v-3.so
/usr/lib64/libguile-srfi-srfi-1-v-3.so.3
/usr/lib64/libguile-srfi-srfi-1-v-3.so.3.0.0
/usr/lib64/libguile-srfi-srfi-13-14-v-3.so
/usr/lib64/libguile-srfi-srfi-13-14-v-3.so.3
/usr/lib64/libguile-srfi-srfi-13-14-v-3.so.3.0.0
/usr/lib64/libguile-srfi-srfi-4-v-3.so
/usr/lib64/libguile-srfi-srfi-4-v-3.so.3
/usr/lib64/libguile-srfi-srfi-4-v-3.so.3.0.0
/usr/lib64/libguile-srfi-srfi-60-v-2.so
/usr/lib64/libguile-srfi-srfi-60-v-2.so.2
/usr/lib64/libguile-srfi-srfi-60-v-2.so.2.0.0
/usr/lib64/libguile.so.17
/usr/lib64/libguile.so.17.0.0
/usr/lib64/libguilereadline-v-17.so
/usr/lib64/libguilereadline-v-17.so.17
/usr/lib64/libguilereadline-v-17.so.17.0.0
**

I tried

setenv LDFLAGS -L/usr/lib64

It is still not working.

Also, it can not find lapack and blas on the cluster, either. And I am
very sure they are appropriately installed.

Can anyone give me some suggestions on this problem?

Thanks,
Hua




--

Message: 2
Date: Thu, 5 Mar 2009 13:12:48 +0800 (CST)
From: adrian hiaquar...@yahoo.com.tw
Subject: [Meep-discuss] Complex field in meep
To: meep-discuss@ab-initio.mit.edu
Message-ID: 825480.96270...@web73705.mail.tp2.yahoo.com
Content-Type: text/plain; charset=utf-8

? Hi Meep users and Steven:
? ? I know this question has been discussed but I can't find the  
answer to previous poster.

? ? My question is:

??? 1. what is the definition of complex field in FDTD. Is it just  
E0*exp(j(wt - kr)) and the real field and imaginary field correspond  
to real part and imaginary part, respectively?
? If so, if we set force-complex-fields? false, then the real  
field we normally got is E0 times cos(wt - kr)?
? ? ? I try to reach Complex field in meep? but it requires  
password. so I can't find the definition of complex fields in meep.


 The purpose for getting complex field is to get the far field  
patterns from the near field patterns calculated by meep


Thanks in advanced

Adrian
??



  
__

?Yahoo!??2.0? http://tw.mg0.mail.yahoo.com/dc/landing
-- next part --
An HTML attachment was scrubbed...
URL: http://ab-initio.mit.edu/pipermail/meep-discuss/attachments/20090305/9896d4d6/attachment.html 



--

Message: 3
Date: Thu, 5 Mar 2009 08:33:59 +0100
From: Michal Martin? cmelak9...@centrum.cz
Subject: [Meep-discuss]  meep initial field
To: meep-discuss@ab-initio.mit.edu
Message-ID: 01c99d64$bf8dee10$6f01a...@michalp4

Re: [Meep-discuss] meep-discuss Digest, Vol 37, Issue 3

2009-03-05 Thread Steven G. Johnson


On Mar 5, 2009, at 12:12 PM, Ian Buss wrote:
I had a similar problem when simulating real metals - it turns out  
that Meep will not accept negative permittivities without blowing up.


This is not a property of Meep specifically.  FDTD is unstable for non- 
dispersive negative epsilon (which also violates the Kramers-Kronig  
relations).


___
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss


Re: [Meep-discuss] meep-discuss Digest, Vol 37, Issue 3

2009-03-05 Thread Nathan

Dr. Johnson,

Is this true when just the complex part of epsilon is negative?  I was 
aware of the instability for non-dispersive negative epsilon, but 
thought it was only for the real part. I know the FD-TD method can model 
non-dispersive gain situations, but is there a limit on the value of the 
gain?  I have found some work that suggests the time step must be 
reduced when modeling complex epsilon to ensure stability, but that 
didn't seem to fix my issue.  I tried cutting the Courant number in half 
and the results were the same in my simulation. 


Regards,
Nathan

Steven G. Johnson wrote:

On Mar 5, 2009, at 12:12 PM, Ian Buss wrote:
  
I had a similar problem when simulating real metals - it turns out  
that Meep will not accept negative permittivities without blowing up.



This is not a property of Meep specifically.  FDTD is unstable for non- 
dispersive negative epsilon (which also violates the Kramers-Kronig  
relations).


___
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
  


___
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss