From cae2fd60d5582ba4e55b148e865ade006b44dd16 Mon Sep 17 00:00:00 2001 From: Vladimir Chalupecky Date: Mon, 18 Apr 2016 11:14:48 +0900 Subject: [PATCH] zero index remaining occurences --- lapack.go | 4 ++-- native/dlatrd.go | 6 +++--- 2 files changed, 5 insertions(+), 5 deletions(-) diff --git a/lapack.go b/lapack.go index a4134559..c099ed2e 100644 --- a/lapack.go +++ b/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. diff --git a/native/dlatrd.go b/native/dlatrd.go index 653aee37..8f929388 100644 --- a/native/dlatrd.go +++ b/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