// 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))) }