From 022eeec0ef663019fd325c3ed595b1010af1b42c Mon Sep 17 00:00:00 2001 From: Vladimir Chalupecky Date: Wed, 18 Jan 2017 14:05:47 +0100 Subject: [PATCH] native: pass correct length ipiv from Dgetrf to Dlaswp --- native/dgetrf.go | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/native/dgetrf.go b/native/dgetrf.go index dd23ff10..d3200b51 100644 --- a/native/dgetrf.go +++ b/native/dgetrf.go @@ -51,9 +51,9 @@ func (impl Implementation) Dgetrf(m, n int, a []float64, lda int, ipiv []int) (o for i := j; i <= min(m-1, j+jb-1); i++ { ipiv[i] = j + ipiv[i] } - impl.Dlaswp(j, a, lda, j, j+jb-1, ipiv, 1) + impl.Dlaswp(j, a, lda, j, j+jb-1, ipiv[:j+jb], 1) if j+jb < n { - impl.Dlaswp(n-j-jb, a[j+jb:], lda, j, j+jb-1, ipiv, 1) + impl.Dlaswp(n-j-jb, a[j+jb:], lda, j, j+jb-1, ipiv[:j+jb], 1) bi.Dtrsm(blas.Left, blas.Lower, blas.NoTrans, blas.Unit, jb, n-j-jb, 1, a[j*lda+j:], lda,