lapack/gonum: don't loop unnecessarily in Dlanv2

gonum · Oct 4, 2020 · 4380c62 · 4380c62
1 parent de0df9a
commit 4380c62
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/lapack/gonum/dlanv2.go b/lapack/gonum/dlanv2.go
@@ -69,6 +69,7 @@ func (impl Implementation) Dlanv2(a, b, c, d float64) (aa, bb, cc, dd float64, r
 			safmn2 := math.Pow(dlamchB, math.Log(dlamchS/dlamchE)/math.Log(dlamchB)/2)
 			safmx2 := 1 / safmn2
 			sigma := b + c
+		loop:
 			for iter := 0; iter < 20; iter++ {
 				scale = math.Max(math.Abs(temp), math.Abs(sigma))
 				switch {
@@ -78,6 +79,8 @@ func (impl Implementation) Dlanv2(a, b, c, d float64) (aa, bb, cc, dd float64, r
 				case scale <= safmn2:
 					sigma *= safmx2
 					temp *= safmx2
+				default:
+					break loop
 				}
 			}
 			p = temp / 2