// For symmetric matrix and vector X solve
// Y = alpha * A * X + beta * Y
func dsymv(uplo string, N int, alpha float64,
A []float64, lda int, X []float64, incX int, beta float64,
Y []float64, incY int) {
cuplo := C.CString(uplo)
defer C.free(unsafe.Pointer(cuplo))
C.dsymv_(cuplo, (*C.int)(unsafe.Pointer(&N)),
(*C.double)(unsafe.Pointer(&alpha)),
(*C.double)(unsafe.Pointer(&A[0])),
(*C.int)(unsafe.Pointer(&lda)),
(*C.double)(unsafe.Pointer(&X[0])),
(*C.int)(unsafe.Pointer(&incX)),
(*C.double)(unsafe.Pointer(&beta)),
(*C.double)(unsafe.Pointer(&Y[0])),
(*C.int)(unsafe.Pointer(&incY)))
}
// For symmetric matrix and vector X solve
// Y = alpha * A * X + beta * Y
func dsymv(uplo string, N int, alpha float64,
A []float64, lda int, X []float64, incX int, beta float64,
Y []float64, incY int) {
cuplo := C.CString(uplo)
defer C.free(unsafe.Pointer(cuplo))
// ?? TODO: protect against index out of bounds panics.
C.dsymv_(cuplo, (*C.int)(unsafe.Pointer(&N)),
(*C.double)(unsafe.Pointer(&alpha)),
(*C.double)(unsafe.Pointer(&A[0])),
(*C.int)(unsafe.Pointer(&lda)),
(*C.double)(unsafe.Pointer(&X[0])),
(*C.int)(unsafe.Pointer(&incX)),
(*C.double)(unsafe.Pointer(&beta)),
(*C.double)(unsafe.Pointer(&Y[0])),
(*C.int)(unsafe.Pointer(&incY)))
}