Two recursive and numerically stable matrix algorithms for modeling layered diffraction gratings, the S-matrix algorithm and the R-matrix algorithm, are systematically presented in a form that is independent of the underlying grating models, geometries, and mountings. Many implementation variants of the algorithms are also presented. Their physical interpretations are given, and their numerical stabilities and efficiencies are discussed in detail. The single most important criterion for achieving unconditional numerical stability with both algorithms is to avoid the exponentially growing functions in every step of the matrix recursion. From the viewpoint of numerical efficiency, the S-matrix algorithm is generally preferred to the R-matrix algorithm, but exceptional cases are noted.
© 1996 Optical Society of America
Original Manuscript: July 20, 1995
Revised Manuscript: November 13, 1995
Manuscript Accepted: December 4, 1995
Published: May 1, 1996
Lifeng Li, "Formulation and comparison of two recursive matrix algorithms for modeling layered diffraction gratings," J. Opt. Soc. Am. A 13, 1024-1035 (1996)