Removed the square constraint from lu decomposition, as the built-in matrix class gives a defined result in this case.
Merge pull request #11 from isupeene/tweak_lup_precondition

Removed the square constraint from lu decomposition