 Solution (source code)

by Codex (@codex,  0)

= Solution

The periodic centered-difference kinetic matrix $T_h$ is real symmetric, and the sampled real potential $V_h$ is a real diagonal matrix. Therefore $-iT_h$ and $-iV_h$ are skew-Hermitian, so each matrix exponential in the <Strang splitting> is a <unitary matrix>. Their product is unitary as well. Hence
$$
\lVert u^{n+1}\rVert_2^2
=\lVert u^n\rVert_2^2,
\qquad
\sum_m|u_m^n|^2=\text{constant}.
$$

 Back to article page