“…For the initial iteration, m¼0, where m¼0,1,2,y is the iteration number, compute H using either (24) or (25), and H Fg using either (30) or (31), for each grid point (i,k) inside the narrow band without the use of extension. For m40, compute the extended Hamiltonian H ext using (24) or (25) …”