add overload for plu_decompose which outputs an array represents the permutation matrix.
add some test utility functions
add test for new plu_decompose which outputs permutation array.