// void dpotrs_(char *uplo, int *n, int *nrhs, double *A,
// int *lda, double *B, int *ldb, int *info);
func dpotrs(uplo string, N, Nrhs int, A []float64, lda int, B []float64, ldb int) int {
var info int = 0
cuplo := C.CString(uplo)
defer C.free(unsafe.Pointer(cuplo))
C.dpotrs_(cuplo,
(*C.int)(unsafe.Pointer(&N)),
(*C.int)(unsafe.Pointer(&Nrhs)),
(*C.double)(unsafe.Pointer(&A[0])),
(*C.int)(unsafe.Pointer(&lda)),
(*C.double)(unsafe.Pointer(&B[0])),
(*C.int)(unsafe.Pointer(&ldb)),
(*C.int)(unsafe.Pointer(&info)))
return info
}
// DPOTRS: (Double-precision) POsitive-definite TRiangular factor Solve
//
// http://www.netlib.org/lapack/double/dpotrs.f
func dpotrs(uplo Triangle, n, nrhs int, a []float64, lda int, b []float64, ldb int) error {
var (
uplo_ = uploChar(uplo)
n_ = C.integer(n)
nrhs_ = C.integer(nrhs)
a_ = ptrFloat64(a)
lda_ = C.integer(lda)
b_ = ptrFloat64(b)
ldb_ = C.integer(ldb)
)
var info_ C.integer
C.dpotrs_(&uplo_, &n_, &nrhs_, a_, &lda_, b_, &ldb_, &info_)
info := int(info_)
switch {
case info < 0:
return errInvalidArg(-info)
case info == 0:
return nil
default:
panic(errUnknown(info))
}
}
