Debugged `luDecomp' and added its complete pivoting version

Sorry, I found that my implementation of `luDecomp` doesn't generate triangular matrices for non-square inputs.
This pull request fixes that, and also implemented `luDecomp'`, the complete pivoting version of `luDecomp`.