def Mat.solve

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

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

245# X = A^-1 B; Fail with the zero pivot as the code when A is singular
246def Mat.solve(+n: Nat, +nrhs: Nat, a: Mat<n, n>, b: Mat<n, nrhs>) -> IO(Result<&1, &1, U32 & String, Mat<n, nrhs>>):
247  MatA{x} = a
248  MatA{y} = b
249  IO.bind(Result<&1, &1, U32 & String, Array<F32>>, Result<&1, &1, U32 & String, Mat<n, nrhs>>, L.Lapack.sgesv(U32.from_nat(n), U32.from_nat(nrhs), x, y), r => Mat.wrap(n, nrhs, r))