2 回答
TA贡献1770条经验 获得超3个赞
通过 cgo 调用例如 cblas 很容易:
package main
// #include <cblas.h>
// #cgo LDFLAGS: -L/usr/lib64/atlas -lcblas
import "C"
import "fmt"
type matrix struct {
rows int
cols int
elems []float32
}
func (a matrix) cblasmul(b matrix) (c matrix) {
c = matrix{a.rows, b.cols, make([]float32, a.rows*b.cols)}
C.cblas_sgemm(
C.CblasRowMajor, C.CblasNoTrans, C.CblasNoTrans,
C.int(a.rows), C.int(b.cols), C.int(a.cols),
1.0,
(*C.float)(&a.elems[0]), C.int(a.cols),
(*C.float)(&b.elems[0]), C.int(b.cols),
0.0,
(*C.float)(&c.elems[0]), C.int(c.cols))
return c
}
func main() {
a := matrix{100, 100, make([]float32, 100*100)}
b := matrix{100, 100, make([]float32, 100*100)}
// ...
c := a.cblasmul(b)
fmt.Println(c)
}
TA贡献1777条经验 获得超3个赞
- 2 回答
- 0 关注
- 204 浏览
添加回答
举报