update_magnetic_fields.m

来自「The Finite Difference Time Domain Method」· M 代码 · 共 19 行

% update magnetic fields

current_time  = current_time + dt/2;
    
Hx = Chxh.*Hx+Chxey.*(Ey(1:nxp1,1:ny,2:nzp1)-Ey(1:nxp1,1:ny,1:nz)) ...
    + Chxez.*(Ez(1:nxp1,2:nyp1,1:nz)-Ez(1:nxp1,1:ny,1:nz)); 
                  
Hy = Chyh.*Hy+Chyez.*(Ez(2:nxp1,1:nyp1,1:nz)-Ez(1:nx,1:nyp1,1:nz)) ...
    + Chyex.*(Ex(1:nx,1:nyp1,2:nzp1)-Ex(1:nx,1:nyp1,1:nz)); 
               
Hz = Chzh.*Hz+Chzex.*(Ex(1:nx,2:nyp1,1:nzp1)-Ex(1:nx,1:ny,1:nzp1))  ...
    + Chzey.*(Ey(2:nxp1,1:ny,1:nzp1)-Ey(1:nx,1:ny,1:nzp1)); 

if incident_plane_wave.enabled
    Hx = Hx + Chxhic .* Hxic + Chxhip .* Hxip; 
    Hy = Hy + Chyhic .* Hyic + Chyhip .* Hyip; 
    Hz = Hz + Chzhic .* Hzic + Chzhip .* Hzip; 
end