Program 11 1D FDTD Vector update with Sparse matrix update

1D unitless FDTD with vector and sparse matrix update

Youtube Links :-
http://www.youtube.com/watch?v=cOFA-V7DhKg
http://www.youtube.com/watch?v=EWajv11FlDs
http://www.youtube.com/watch?v=2Ui09lMi2Zw
http://www.youtube.com/watch?v=1XdMmI8NVTw

Written for Course :-
Computational Electromagnetics, Fall 2011
Department of Electrical Engineering
Indian Institute of Technology Madras (IITM)
Chennai - 600036, India

Authors :-
Sathya Swaroop Ganta, B.Tech., M.Tech. Electrical Engg.
Kayatri, M.S. Engineering Design
Pankaj, M.S. Electrical Engg.
Sumantra Chaudhuri, M.S. Electrical Engg.
Projesh Basu, M.S. Electrical Engg.
Nikhil Kumar CS, M.S. Electrical Engg.

Instructor :-
Dr. Ananth Krishnan
Assistant Professor
Department of Electrical Engineering
Indian Institute of Technology Madras
http://www.ee.iitm.ac.in/~ananthk

Description :-
The 1D TEM wave is x-directed z-polarized TEM wave containing the y-directed magnetic field Hy and z-directed electric field Ez. The time update in Yee Algorithm is done using Leapfrog time-stepping. Here, the H fields are updated every half time-step and E fileds are updated every full time-step. This is shown by two alternating vector updates spanning entire spatial grid inside a main for-loop for time update spanning the entire time-grid. Yes, there are No For Loops except for time steps. The vector updates span only a part of spatial grid where the wave, starting from source, has reached at that particular time instant (exploiting sparse vectors) avoiding field updates at all points in the grid which is unnecessary at that time instant.