I wrote the data simulation in R, but I couldn't find out the function for "Complete pivoting LU decomposition (which is equivalent to FullPivLU in C++)"
If we need to manually calculate this, this link would be helpful.
Complete pivoting is for making more stable the decomposition when having pivots. Since we are filling the matrixes with norm(0,1) then there shouldn't be any pivots and hence is not necessary.
I wrote the data simulation in R, but I couldn't find out the function for "Complete pivoting LU decomposition (which is equivalent to FullPivLU in C++)"
If we need to manually calculate this, this link would be helpful.