def Mat.cholesky_solve

bend-blas-lapack@0.0.0.1 · mat.bend:262

import Base
import bend-blas-lapack@0.0.0.1/lapack.bend as L

262# X = A^-1 B given L from Mat.cholesky
263def Mat.cholesky_solve(+n: Nat, +nrhs: Nat, l: Mat<n, n>, b: Mat<n, nrhs>) -> IO(Result<&1, &1, U32 & String, Mat<n, nrhs>>):
264  MatA{x} = l
265  MatA{y} = b
266  IO.bind(Result<&1, &1, U32 & String, Array<F32>>, Result<&1, &1, U32 & String, Mat<n, nrhs>>, L.Lapack.spotrs(U32.from_nat(n), U32.from_nat(nrhs), x, y), r => Mat.wrap(n, nrhs, r))