On Aug 11, 2008, at 7:58 AM, Ryan Hao wrote: > If I only want the period in Y direction: > > (set-param! k-point (vector3 0 0 0)) > ; firstly set the X,Y,Z directions into periodic boundary, because I > saw many people do this in this mail list > (set! pml-layers (list (make pml (direction X)(thickness dpml)) > (make pml (direction Z)(thickness dpml)) )) > ; secondly change the X,Z directions into PML
The second line doesn't "change" the X and Z directions into PML. You are have the common misunderstanding that PML has something to do with the boundary condition. It doesn't. PML is not a boundary condition. You can use PML with any boundary condition, including periodic. What the second line does is to specify that the periodic X and Z directions have PML absorbing layers adjacent to those periodic boundaries. (Those directions are still periodic but the periodicity is irrelevant because waves are absorbed before they propagate back into the interior of the cell around the periodic boundaries.) > When I do this, MEEP trace an error: "infinite period for "set- > param! k-point"" Are you sure that is the *exact* error message? I searched through the source code of Meep, libctl, and Guile, and found no instances of the words "infinite period", so I don't see how that error message could have been generated. Regards, Steven G. Johnson _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
