zero index remaining occurences

This commit is contained in:
Vladimir Chalupecky
2016-04-18 11:14:48 +09:00
parent 7c9410be91
commit cae2fd60d5
2 changed files with 5 additions and 5 deletions

View File

@@ -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.

View File

@@ -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