Hi guys,

Is there any way I can implement this scheme syntax in C++?

(run-sources+ (stop-when-fields-decayed
               50 Ey
               (vector3 (- (* 0.5 sx) dpml 0.5) 0)
               1e-3)
              (during-sources
               (in-volume (volume (center 0 0) (size sx 0))
                (to-appended "hz-slice" (at-every 0.4 output-hfield-z)))))

It seems like usually the while loop is as follows:
---------------------------
while ( f.time() < f.last_source_time()+400)
    {
        f.step();
        f.output_hdf5(Ey, vol.surroundings());
    }
---------------------------
How do I increment the time in steps of 50 time units till Ey has decayed to 
1/1000 on the plane indicated by  "(vector3 (- (* 0.5 sx) dpml 0.5) 0)"? 
How do I output hfield on that slice indicated by "(volume (center 0 0) (size 
sx 0))" every 0.4 time unit? It seems like f.output_hdf5() can only output the 
field data in the entire waveguide.

Thanks a lot in advance.
--
Bin Huang (Bryan)
MIT Graduate Student '10
Computation for Design and Optimization
(+65)98947649
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Reply via email to