expmARPACK公司

光子晶体模型中麦克斯韦方程的Krylov子空间指数时域解。指数时间积分,即涉及矩阵指数的时间积分,是涉及麦克斯韦方程组的时域建模的一种有吸引力的工具。然而,它在实际应用中往往需要大量的数值线性代数算法知识,如Krylov子空间方法。在本文中,我们讨论指数Krylov子空间时间积分方法,并提供如何在实践中使用这些方法的简单指南。特别是针对纳米光子学的应用,我们有意保持尽可能一般的表现形式,并考虑带阻尼的全矢量麦克斯韦方程组(即,具有非零电导项)。详细讨论了Krylov移位逆方法和基于残差的停止准则等有效技术。通过数值实验验证了所讨论方法的有效性及其与网格无关的收敛性。这里描述的一些算法可以从url{http://www.math.utwente.nl/botchevma/expm/}获得Octave/Matlab代码。