skip to content

Splitting methods for oscillatory non-autonomous linear systems

Tuesday 3rd July 2007 - 11:30 to 12:00
INI Seminar Room 1

Several relevant physical phenomena are modelled (totally or partially) by linear PDEs that, once spatially discretized, give rise to systems of coupled harmonic oscillators. To diagonalise these systems is usually prohibitively expensive. Then, to use splitting methods (involving matrix-vector products and possibly FFTs) is a valid alternative. A theoretical analysis about the stability and accuracy of splitting methods on the harmonic oscillator allows us to build new methods which outperform the existing methods from the literature.

Non-autonomous problems are also of great interest but, the most efficient methods for the autonomous case are not valid in this setting. From the Magnus series expansion (as a formal solution to the non-autonomous problem) we show how to adapt these methods by treating the "time" separately from the coordinates. This technique allows us to build new methods whose performance is tested on the Schrödinger equation with time-dependent potentials.

Related Links

Presentation Material: 
University of Cambridge Research Councils UK
    Clay Mathematics Institute London Mathematical Society NM Rothschild and Sons