* Overwrites A with [L\U] and returns P, st PLU=A. L is considered to have 1s in the diagonal. */
()
source not stored for this graph (policy: none)