LU decomposition
function lu input Real A[:, size(A, 1)]; constant Integer n = size(A, 1); output Real L[size(A, 1), size(A, 1)]; output Real U[size(A, 1), size(A, 1)]; output Real P[size(A, 1), size(A, 1)]; end lu;