Added a function to reduce the dimensionality of a linear transformation. This allows us to reduce the number of internal modes being used.

Also got the matrix matrix multiplication working as needed for now.

Looks like we are completely set for the internal wave model?