diff --git a/mat/band.go b/mat/band.go
index b8994b4e..8b26f02b 100644
--- a/mat/band.go
+++ b/mat/band.go
@@ -322,13 +322,13 @@ func (b *BandDense) MulVecTo(dst *VecDense, trans bool, x Vector) {
 			blas64.Gbmv(t, 1, b.mat, xVec.mat, 0, dst.mat)
 		} else {
 			xCopy := getWorkspaceVec(n, false)
-			xCopy.CloneVec(xVec)
+			xCopy.CloneFromVec(xVec)
 			blas64.Gbmv(t, 1, b.mat, xCopy.mat, 0, dst.mat)
 			putWorkspaceVec(xCopy)
 		}
 	} else {
 		xCopy := getWorkspaceVec(n, false)
-		xCopy.CloneVec(x)
+		xCopy.CloneFromVec(x)
 		blas64.Gbmv(t, 1, b.mat, xCopy.mat, 0, dst.mat)
 		putWorkspaceVec(xCopy)
 	}
diff --git a/mat/symband.go b/mat/symband.go
index ceefb091..be3f2c2f 100644
--- a/mat/symband.go
+++ b/mat/symband.go
@@ -267,13 +267,13 @@ func (s *SymBandDense) MulVecTo(dst *VecDense, _ bool, x Vector) {
 			blas64.Sbmv(1, s.mat, xVec.mat, 0, dst.mat)
 		} else {
 			xCopy := getWorkspaceVec(n, false)
-			xCopy.CloneVec(xVec)
+			xCopy.CloneFromVec(xVec)
 			blas64.Sbmv(1, s.mat, xCopy.mat, 0, dst.mat)
 			putWorkspaceVec(xCopy)
 		}
 	} else {
 		xCopy := getWorkspaceVec(n, false)
-		xCopy.CloneVec(x)
+		xCopy.CloneFromVec(x)
 		blas64.Sbmv(1, s.mat, xCopy.mat, 0, dst.mat)
 		putWorkspaceVec(xCopy)
 	}
diff --git a/mat/vector.go b/mat/vector.go
index c3f8b1b6..d3cbc721 100644
--- a/mat/vector.go
+++ b/mat/vector.go
@@ -194,9 +194,9 @@ func (v *VecDense) Zero() {
 	}
 }
 
-// CloneVec makes a copy of a into the receiver, overwriting the previous value
+// CloneFromVec makes a copy of a into the receiver, overwriting the previous value
 // of the receiver.
-func (v *VecDense) CloneVec(a Vector) {
+func (v *VecDense) CloneFromVec(a Vector) {
 	if v == a {
 		return
 	}
@@ -218,7 +218,7 @@ func (v *VecDense) CloneVec(a Vector) {
 // VecDenseCopyOf returns a newly allocated copy of the elements of a.
 func VecDenseCopyOf(a Vector) *VecDense {
 	v := &VecDense{}
-	v.CloneVec(a)
+	v.CloneFromVec(a)
 	return v
 }
 
diff --git a/optimize/bfgs.go b/optimize/bfgs.go
index 4ab6d204..b44ef81e 100644
--- a/optimize/bfgs.go
+++ b/optimize/bfgs.go
@@ -91,8 +91,8 @@ func (b *BFGS) InitDirection(loc *Location, dir []float64) (stepSize float64) {
 
 	x := mat.NewVecDense(dim, loc.X)
 	grad := mat.NewVecDense(dim, loc.Gradient)
-	b.x.CloneVec(x)
-	b.grad.CloneVec(grad)
+	b.x.CloneFromVec(x)
+	b.grad.CloneFromVec(grad)
 
 	b.y.Reset()
 	b.s.Reset()