https://kitchingroup.cheme.cmu.edu/blog/2013/04/01/Computing-determinants-from-matrix-decompositions/
1. A = PLU,
det(L) = 1 det(U) = Product of diagonal elements det(P) = (-1)^p