From 228e3383244aa55e29b91a6f80cde112ecb54008 Mon Sep 17 00:00:00 2001 From: btracey Date: Mon, 6 Jan 2014 15:26:37 -0800 Subject: [PATCH] further simplification to Drotm --- referenceblas/level1double.go | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/referenceblas/level1double.go b/referenceblas/level1double.go index 29a20cdd..1e499ab3 100644 --- a/referenceblas/level1double.go +++ b/referenceblas/level1double.go @@ -457,15 +457,14 @@ func (Blas) Drotm(n int, x []float64, incX int, y []float64, incY int, p blas.Dr iy = (-n + 1) * incY } for i := 0; i < n; i++ { - w := x[ix] - z := y[iy] - x[ix] = w*h11 + z*h12 - y[iy] = w*h21 + z*h22 + x[ix], y[iy] = x[ix]*h11+y[iy]*h12, x[ix]*h21+y[iy]*h22 ix += incX iy += incY } /* + More complicated version of the above, though possibly faster + if flag == -2 { return }