@karlnapf I should change this to
P*c*P^-1=L*U, i.e. c=LP*UP, where LP=P^-1*L and UP=U*P