Hi Arthur, Your advice for implementing stop-when-fields-decayed has been extremely helpful. Thank you so much for that.
I am trying to implement "(run-k-points T k-points)". Thus I need to first implement "(run-k-point T k)" and I am stuck at how to implement these two lines. (change-k-point! k) (restart-fields) >From what I understand, the vector "k-point" defines the direction in which >the wave propagates through the waveguide. So I guess in C++ I should >re-initialize the fields using the new k-point in each iteration of the loop >in "(run-k-points T k-points)". I am not sure how to do it. In fact, I >searched through fieds.cpp, the method "initialize_field" is never called. >When I write my own C++ program, I also never had to explicitly initialize the >field. My question is: 1) how to change k-point parameter? 2) how to initialize the field with the effect of the change in k-point? Thanks so much! -- 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
