Mass-matrix orthogonality of normal modes (source code)

= Mass-matrix orthogonality of normal modes

For symmetric stiffness matrix $K$ and positive-definite mass matrix $M$, eigenvectors $a_i,a_j$ of
$$
Ka=\omega^2Ma
$$
with distinct squared frequencies satisfy
$$
a_i^TMa_j=0.
$$
Thus normal modes are orthogonal in the kinetic-energy inner product.