Hi,
Suppose I have a test.ctl like that:
(define-param dt 0.5)
(define-param sx 10)
(define-param sy 2)
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! resolution 10)
(run-until 60 (to-appended "ex" (at-every dt output-efield-x)))
I should get 60*.5=120 data points and that what I get:
meep test.ctl ; h5ls test-ex.h5
-----------
Initializing structure...
Working in 2D dimensions.
time for set_epsilon = 0.010647 s
-----------
creating output file "./test-ex.h5"...
run 0 finished at t = 60.0 (1200 timesteps)
Elapsed run time = 0.612535 s
ex Dataset {100, 20, 120/Inf}
But with dt=0.1 for example, I don't get 60*.1=600 data points:
meep dt=0.1 test.ctl ; h5ls test-ex.h5
command-line param: dt=0.1
-----------
Initializing structure...
Working in 2D dimensions.
time for set_epsilon = 0.00977708 s
-----------
creating output file "./test-ex.h5"...
run 0 finished at t = 60.0 (1200 timesteps)
Elapsed run time = 1.91347 s
ex Dataset {100, 20, 507/Inf}
Why 507 only ? I don't understand.
Cheers.
--
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
