Include LQD-based solver in Solve

A type switch is included to allow for the fact that (*LQ).Solve will
not overwrite the matrix parameter, but if Solve is handed a non-*Dense
matrix, it must be converted.

Solution is confirmed by octave - the test should probably just confirm
that AX == B.