The C++ code that steps H and E sources is different although they should be treated in same way (symmetry between electric and magnetic current density) ?
for example, if the source was a certain function f(t), for E sources: E(i,j,k) -= f (t) (update_e_from_d.cpp) for H sources: H(i,j,k) -=d/dt[f(t)] * dt (step.cpp) Why do we have this difference ? Thanks Samer _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
