zero index remaining occurences

Author: vladimir-ch

lapack.go

@@ -48,8 +48,8 @@ type Float64 interface {
 type Direct byte
 
 const (
-	Forward  Direct = 'F' // Reflectors are right-multiplied, H(0) * H(1) * ... * H(k).
-	Backward Direct = 'B' // Reflectors are left-multiplied, H(k) * ... * H(1) * H(0).
+	Forward  Direct = 'F' // Reflectors are right-multiplied, H(0) * H(1) * ... * H(k-1).
+	Backward Direct = 'B' // Reflectors are left-multiplied, H(k-1) * ... * H(1) * H(0).
 )
 
 // Sort is the sorting order.

native/dlatrd.go

@@ -53,12 +53,12 @@ import (
 // H has the form
 //  I - tau * v * v^T
 // If uplo == blas.Upper,
-//  Q = H(n) * H(n-1) * ... * H(n-nb+1)
+//  Q = H(n-1) * H(n-2) * ... * H(n-nb)
 // where v[:i-1] is stored in A[:i-1,i], v[i-1] = 1, and v[i:n] = 0.
 //
 // If uplo == blas.Lower,
-//  Q = H(1) * H(2) * ... H(nb)
-// where v[1:i+1] = 0, v[i+1] = 1, and v[i+2:n] is stored in A[i+2:n,i].
+//  Q = H(0) * H(1) * ... * H(nb-1)
+// where v[:i+1] = 0, v[i+1] = 1, and v[i+2:n] is stored in A[i+2:n,i].
 //
 // The vectors v form the n×nb matrix V which is used with W to apply a
 // symmetric rank-2 update to the unreduced part of A