## expmARPACK

swMATH ID: | 13396 |

Software Authors: | Botchev, Mikhail A. |

Description: | Krylov subspace exponential time domain solution of Maxwell’s equations in photonic crystal modeling. The exponential time integration, i.e., time integration which involves the matrix exponential, is an attractive tool for time domain modeling involving Maxwell’s equations. However, its application in practice often requires a substantial knowledge of numerical linear algebra algorithms, such as Krylov subspace methods. In this note we discuss exponential Krylov subspace time integration methods and provide a simple guide on how to use these methods in practice. While specifically aiming at nanophotonics applications, we intentionally keep the presentation as general as possible and consider full vector Maxwell’s equations with damping (i.e., with nonzero conductivity terms). Efficient techniques such as the Krylov shift-and-invert method and residual-based stopping criteria are discussed in detail. Numerical experiments are presented to demonstrate the efficiency of the discussed methods and their mesh independent convergence. Some of the algorithms described here are available as Octave/Matlab codes from url{http://www.math.utwente.nl/ botchevma/expm/}. |

Homepage: | http://wwwhome.math.utwente.nl/~botchevma/expm/ |

Keywords: | finite difference time domain (FDTD) method; matrix exponential; Maxwell’s equations; Krylov subspace methods; photonic crystals |

Related Software: | Expokit; mftoolbox; phipm; Algorithm 919; Expint; MATLAB expm; PARAEXP; Matlab; UMFPACK; KIOPS; IFISS; MPAS-Ocean; RKC; RODAS; AGMG; Octave |

Cited in: | 9 Publications |

all
top 5

### Cited by 13 Authors

all
top 5

### Cited in 7 Serials

### Cited in 5 Fields

9 | Numerical analysis (65-XX) |

3 | Partial differential equations (35-XX) |

2 | Linear and multilinear algebra; matrix theory (15-XX) |

1 | Statistical mechanics, structure of matter (82-XX) |

1 | Geophysics (86-XX) |