Return permutation matrix in LU decomposition.

Current LU decomposition implementation is useless without permutation matrix P because  A = L*U is incorrect, correct formula: A = P*L*U: [LU decomposition](http://en.wikipedia.org/wiki/Lu_decomposition)