For weird reasons I'd like to solve Ux = y where U is the Cholesky square root of a given symmetric positive definite matrix A. Is there a way to abuse the PETSc solver interface into doing this?

all the best

--cjc
--