[Meep-discuss] bug in interpolate procedure

2009-01-30 Thread Antoine Monmayrant
Hi everyone, I think I found something fishy with the procedure interpolate. It seems there is a rounding bug somewhere: meep (interpolate 9 (list 0.5 0)) (0.5 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0) which is what I expected. However, if I interpolate the other way round I get this: meep

Re: [Meep-discuss] bug in interpolate procedure

2009-01-30 Thread matt
This is caused by floating point arithmetic. Have a look at this thread: http://www.mail-archive.com/meep-discuss@ab-initio.mit.edu/msg00978.html Best, Matt On Fri, 30 Jan 2009, Antoine Monmayrant wrote: Hi everyone, I think I found something fishy with the procedure interpolate.

[Meep-discuss] question: timestep in Meep

2009-01-30 Thread Ajayram S K S M
Hello! I am trying to simulate em fields in the vicinity of nanoparticles. parameters: a = 1 micron How do we find the size of the timestep in seconds? How do we know how many wavelengths will pass at a particular point in two timesteps? Here the last part of the program: (set!

Re: [Meep-discuss] question: timestep in Meep

2009-01-30 Thread Paweł Smoliński
How do we find the size of the timestep in seconds? How do we know how many wavelengths will pass at a particular point in two timesteps? Hi, timestep (dt) is set when you are define your structure and it depends on two parameters: structure resolution r and Courant parameter (dt = r/Courant,