Return permutation matrix in LU decomposition.
Merge pull request #136 from nbeloglazov/lu-fix

Return permutation matrix in LU decomposition.