mat: don't clobber optimisation for TriDense matrix input
MODIFY
mat/vector.go
mat/vector.go
@@ -629,6 +629,7 @@ func (v *VecDense) MulVec(a Matrix, b Vector) {
return
}
case *TriDense:
+ if fast {
v.CopyVec(b)
aU.checkOverlap(v.asGeneral())
ta := blas.NoTrans
@@ -636,6 +637,8 @@ func (v *VecDense) MulVec(a Matrix, b Vector) {
ta = blas.Trans
}
blas64.Trmv(ta, aU.mat, v.mat)
+ return
+ }
case *Dense:
if fast {
aU.checkOverlap(v.asGeneral())

mat: don't clobber optimisation for TriDense matrix input