如何在Python中实现部分旋转的LU分解?

I want to implement my own LU decomposition P,L,U = my_lu(A), so that given a matrix A, computes the LU decomposition with partial pivoting. But I only know how to do it without pivoting. Can anyone help to do the partial pivoting?I want to implement my own LU decomposition P,L