On Aug 5, 2008, at 12:32 PM, Zhu, Lei wrote: > Therefore, c (the speed of light) is cancelled for my code, whose > value will not affect the result. > If my NTF code is consistent for the meep, I do need to divide the > parameters (omiga and rk) of my code by 2*pi. > However, the phase change for sin(wt) or exp(ikr) is 1 per period. > Could you tell me how you deal with this problem? I tried to find it > in the source but I just touch the c++ for a short time.
No, we multiply the frequency and wavevector by 2*pi before computing phases. The 2*pi factors are "there", you just don't include them when you specify them in the interface. So, for example, if you specify a frequency 0.3 in the Meep interface, then the time-dependence of the source is exp(-2*pi*i*0.3*t), where time t is also in Meep units. i.e. what you are specifying is frequency, not angular frequency. And when you specify the k-point for Bloch-periodic boundary conditions, you are also specifying it without the 2*pi factor, so if you specify k = (0.4,0,0) in the interface, then the phase factor between adjacent unit cells with period L in the x direction is exp(2*pi*i*0.4*L). Steven _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
