Abstract
We recently found that the coordinate transformation method (the C method) equipped with well-established recursive algorithms for solving the system of linear equations is numerically instable when it is applied to thinly coated gratings. The origin of this new kind of numerical instability is not the exponential dependence of the field in the coated layers but the ill condition of the eigenvector matrix of the C method when the truncation number is high. Two simple and effective methods to circumvent the new instability are recommended. We also found that the popular recursive matrix algorithms have different (poor) immunities to the new instability, and they all perform inferiorly to the full matrix (nonrecursive) algorithm.
© 2020 Optical Society of America
Full Article | PDF ArticleMore Like This
Xihong Xu and Lifeng Li
J. Opt. Soc. Am. A 34(6) 881-891 (2017)
R. H. Morf
J. Opt. Soc. Am. A 12(5) 1043-1056 (1995)
Lifeng Li
J. Opt. Soc. Am. A 11(11) 2816-2828 (1994)