+++ /dev/null
-// Copyright 2013 The Gonum Authors. All rights reserved.
-// Use of this code is governed by a BSD-style
-// license that can be found in the LICENSE file
-
-package floats
-
-import (
- "math"
- "strconv"
- "testing"
-
- "golang.org/x/exp/rand"
-)
-
-const (
- EqTolerance = 1E-14
- Small = 10
- Medium = 1000
- Large = 100000
- Huge = 10000000
-)
-
-func AreSlicesEqual(t *testing.T, truth, comp []float64, str string) {
- if !EqualApprox(comp, truth, EqTolerance) {
- t.Errorf(str+". Expected %v, returned %v", truth, comp)
- }
-}
-
-func Panics(fun func()) (b bool) {
- defer func() {
- err := recover()
- if err != nil {
- b = true
- }
- }()
- fun()
- return
-}
-
-func TestAdd(t *testing.T) {
- a := []float64{1, 2, 3}
- b := []float64{4, 5, 6}
- c := []float64{7, 8, 9}
- truth := []float64{12, 15, 18}
- n := make([]float64, len(a))
-
- Add(n, a)
- Add(n, b)
- Add(n, c)
- AreSlicesEqual(t, truth, n, "Wrong addition of slices new receiver")
- Add(a, b)
- Add(a, c)
- AreSlicesEqual(t, truth, n, "Wrong addition of slices for no new receiver")
-
- // Test that it panics
- if !Panics(func() { Add(make([]float64, 2), make([]float64, 3)) }) {
- t.Errorf("Did not panic with length mismatch")
- }
-}
-
-func TestAddTo(t *testing.T) {
- a := []float64{1, 2, 3}
- b := []float64{4, 5, 6}
- truth := []float64{5, 7, 9}
- n1 := make([]float64, len(a))
-
- n2 := AddTo(n1, a, b)
- AreSlicesEqual(t, truth, n1, "Bad addition from mutator")
- AreSlicesEqual(t, truth, n2, "Bad addition from returned slice")
-
- // Test that it panics
- if !Panics(func() { AddTo(make([]float64, 2), make([]float64, 3), make([]float64, 3)) }) {
- t.Errorf("Did not panic with length mismatch")
- }
- if !Panics(func() { AddTo(make([]float64, 3), make([]float64, 3), make([]float64, 2)) }) {
- t.Errorf("Did not panic with length mismatch")
- }
-
-}
-
-func TestAddConst(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- c := 6.0
- truth := []float64{9, 10, 7, 13, 11}
- AddConst(c, s)
- AreSlicesEqual(t, truth, s, "Wrong addition of constant")
-}
-
-func TestAddScaled(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- alpha := 6.0
- dst := []float64{1, 2, 3, 4, 5}
- ans := []float64{19, 26, 9, 46, 35}
- AddScaled(dst, alpha, s)
- if !EqualApprox(dst, ans, EqTolerance) {
- t.Errorf("Adding scaled did not match")
- }
- short := []float64{1}
- if !Panics(func() { AddScaled(dst, alpha, short) }) {
- t.Errorf("Doesn't panic if s is smaller than dst")
- }
- if !Panics(func() { AddScaled(short, alpha, s) }) {
- t.Errorf("Doesn't panic if dst is smaller than s")
- }
-}
-
-func TestAddScaledTo(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- alpha := 6.0
- y := []float64{1, 2, 3, 4, 5}
- dst1 := make([]float64, 5)
- ans := []float64{19, 26, 9, 46, 35}
- dst2 := AddScaledTo(dst1, y, alpha, s)
- if !EqualApprox(dst1, ans, EqTolerance) {
- t.Errorf("AddScaledTo did not match for mutator")
- }
- if !EqualApprox(dst2, ans, EqTolerance) {
- t.Errorf("AddScaledTo did not match for returned slice")
- }
- AddScaledTo(dst1, y, alpha, s)
- if !EqualApprox(dst1, ans, EqTolerance) {
- t.Errorf("Reusing dst did not match")
- }
- short := []float64{1}
- if !Panics(func() { AddScaledTo(dst1, y, alpha, short) }) {
- t.Errorf("Doesn't panic if s is smaller than dst")
- }
- if !Panics(func() { AddScaledTo(short, y, alpha, s) }) {
- t.Errorf("Doesn't panic if dst is smaller than s")
- }
- if !Panics(func() { AddScaledTo(dst1, short, alpha, s) }) {
- t.Errorf("Doesn't panic if y is smaller than dst")
- }
-}
-
-func TestArgsort(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- inds := make([]int, len(s))
-
- Argsort(s, inds)
-
- sortedS := []float64{1, 3, 4, 5, 7}
- trueInds := []int{2, 0, 1, 4, 3}
-
- if !Equal(s, sortedS) {
- t.Error("elements not sorted correctly")
- }
- for i := range trueInds {
- if trueInds[i] != inds[i] {
- t.Error("inds not correct")
- }
- }
-
- inds = []int{1, 2}
- if !Panics(func() { Argsort(s, inds) }) {
- t.Error("does not panic if lengths do not match")
- }
-}
-
-func TestCount(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- f := func(v float64) bool { return v > 3.5 }
- truth := 3
- n := Count(f, s)
- if n != truth {
- t.Errorf("Wrong number of elements counted")
- }
-}
-
-func TestCumProd(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- receiver := make([]float64, len(s))
- result := CumProd(receiver, s)
- truth := []float64{3, 12, 12, 84, 420}
- AreSlicesEqual(t, truth, receiver, "Wrong cumprod mutated with new receiver")
- AreSlicesEqual(t, truth, result, "Wrong cumprod result with new receiver")
- CumProd(receiver, s)
- AreSlicesEqual(t, truth, receiver, "Wrong cumprod returned with reused receiver")
-
- // Test that it panics
- if !Panics(func() { CumProd(make([]float64, 2), make([]float64, 3)) }) {
- t.Errorf("Did not panic with length mismatch")
- }
-
- // Test empty CumProd
- emptyReceiver := make([]float64, 0)
- truth = []float64{}
- CumProd(emptyReceiver, emptyReceiver)
- AreSlicesEqual(t, truth, emptyReceiver, "Wrong cumprod returned with empty receiver")
-
-}
-
-func TestCumSum(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- receiver := make([]float64, len(s))
- result := CumSum(receiver, s)
- truth := []float64{3, 7, 8, 15, 20}
- AreSlicesEqual(t, truth, receiver, "Wrong cumsum mutated with new receiver")
- AreSlicesEqual(t, truth, result, "Wrong cumsum returned with new receiver")
- CumSum(receiver, s)
- AreSlicesEqual(t, truth, receiver, "Wrong cumsum returned with reused receiver")
-
- // Test that it panics
- if !Panics(func() { CumSum(make([]float64, 2), make([]float64, 3)) }) {
- t.Errorf("Did not panic with length mismatch")
- }
-
- // Test empty CumSum
- emptyReceiver := make([]float64, 0)
- truth = []float64{}
- CumSum(emptyReceiver, emptyReceiver)
- AreSlicesEqual(t, truth, emptyReceiver, "Wrong cumsum returned with empty receiver")
-
-}
-
-func TestDistance(t *testing.T) {
- norms := []float64{1, 2, 4, math.Inf(1)}
- slices := []struct {
- s []float64
- t []float64
- }{
- {
- nil,
- nil,
- },
- {
- []float64{8, 9, 10, -12},
- []float64{8, 9, 10, -12},
- },
- {
- []float64{1, 2, 3, -4, -5, 8},
- []float64{-9.2, -6.8, 9, -3, -2, 1},
- },
- }
-
- for j, test := range slices {
- tmp := make([]float64, len(test.s))
- for i, L := range norms {
- dist := Distance(test.s, test.t, L)
- copy(tmp, test.s)
- Sub(tmp, test.t)
- norm := Norm(tmp, L)
- if dist != norm { // Use equality because they should be identical
- t.Errorf("Distance does not match norm for case %v, %v. Expected %v, Found %v.", i, j, norm, dist)
- }
- }
- }
-
- if !Panics(func() { Distance([]float64{}, norms, 1) }) {
- t.Errorf("Did not panic with unequal lengths")
- }
-
-}
-
-func TestDiv(t *testing.T) {
- s1 := []float64{5, 12, 27}
- s2 := []float64{1, 2, 3}
- ans := []float64{5, 6, 9}
- Div(s1, s2)
- if !EqualApprox(s1, ans, EqTolerance) {
- t.Errorf("Mul doesn't give correct answer")
- }
- s1short := []float64{1}
- if !Panics(func() { Div(s1short, s2) }) {
- t.Errorf("Did not panic with unequal lengths")
- }
- s2short := []float64{1}
- if !Panics(func() { Div(s1, s2short) }) {
- t.Errorf("Did not panic with unequal lengths")
- }
-}
-
-func TestDivTo(t *testing.T) {
- s1 := []float64{5, 12, 27}
- s1orig := []float64{5, 12, 27}
- s2 := []float64{1, 2, 3}
- s2orig := []float64{1, 2, 3}
- dst1 := make([]float64, 3)
- ans := []float64{5, 6, 9}
- dst2 := DivTo(dst1, s1, s2)
- if !EqualApprox(dst1, ans, EqTolerance) {
- t.Errorf("DivTo doesn't give correct answer in mutated slice")
- }
- if !EqualApprox(dst2, ans, EqTolerance) {
- t.Errorf("DivTo doesn't give correct answer in returned slice")
- }
- if !EqualApprox(s1, s1orig, EqTolerance) {
- t.Errorf("S1 changes during multo")
- }
- if !EqualApprox(s2, s2orig, EqTolerance) {
- t.Errorf("s2 changes during multo")
- }
- DivTo(dst1, s1, s2)
- if !EqualApprox(dst1, ans, EqTolerance) {
- t.Errorf("DivTo doesn't give correct answer reusing dst")
- }
- dstShort := []float64{1}
- if !Panics(func() { DivTo(dstShort, s1, s2) }) {
- t.Errorf("Did not panic with s1 wrong length")
- }
- s1short := []float64{1}
- if !Panics(func() { DivTo(dst1, s1short, s2) }) {
- t.Errorf("Did not panic with s1 wrong length")
- }
- s2short := []float64{1}
- if !Panics(func() { DivTo(dst1, s1, s2short) }) {
- t.Errorf("Did not panic with s2 wrong length")
- }
-}
-
-func TestDot(t *testing.T) {
- s1 := []float64{1, 2, 3, 4}
- s2 := []float64{-3, 4, 5, -6}
- truth := -4.0
- ans := Dot(s1, s2)
- if ans != truth {
- t.Errorf("Dot product computed incorrectly")
- }
-
- // Test that it panics
- if !Panics(func() { Dot(make([]float64, 2), make([]float64, 3)) }) {
- t.Errorf("Did not panic with length mismatch")
- }
-}
-
-func TestEquals(t *testing.T) {
- s1 := []float64{1, 2, 3, 4}
- s2 := []float64{1, 2, 3, 4}
- if !Equal(s1, s2) {
- t.Errorf("Equal slices returned as unequal")
- }
- s2 = []float64{1, 2, 3, 4 + 1e-14}
- if Equal(s1, s2) {
- t.Errorf("Unequal slices returned as equal")
- }
- if Equal(s1, []float64{}) {
- t.Errorf("Unequal slice lengths returned as equal")
- }
-}
-
-func TestEqualApprox(t *testing.T) {
- s1 := []float64{1, 2, 3, 4}
- s2 := []float64{1, 2, 3, 4 + 1e-10}
- if EqualApprox(s1, s2, 1e-13) {
- t.Errorf("Unequal slices returned as equal for absolute")
- }
- if !EqualApprox(s1, s2, 1e-5) {
- t.Errorf("Equal slices returned as unequal for absolute")
- }
- s1 = []float64{1, 2, 3, 1000}
- s2 = []float64{1, 2, 3, 1000 * (1 + 1e-7)}
- if EqualApprox(s1, s2, 1e-8) {
- t.Errorf("Unequal slices returned as equal for relative")
- }
- if !EqualApprox(s1, s2, 1e-5) {
- t.Errorf("Equal slices returned as unequal for relative")
- }
- if EqualApprox(s1, []float64{}, 1e-5) {
- t.Errorf("Unequal slice lengths returned as equal")
- }
-}
-
-func TestEqualFunc(t *testing.T) {
- s1 := []float64{1, 2, 3, 4}
- s2 := []float64{1, 2, 3, 4}
- eq := func(x, y float64) bool { return x == y }
- if !EqualFunc(s1, s2, eq) {
- t.Errorf("Equal slices returned as unequal")
- }
- s2 = []float64{1, 2, 3, 4 + 1e-14}
- if EqualFunc(s1, s2, eq) {
- t.Errorf("Unequal slices returned as equal")
- }
- if EqualFunc(s1, []float64{}, eq) {
- t.Errorf("Unequal slice lengths returned as equal")
- }
-}
-
-func TestEqualsRelative(t *testing.T) {
- var equalityTests = []struct {
- a, b float64
- tol float64
- equal bool
- }{
- {1000000, 1000001, 0, true},
- {1000001, 1000000, 0, true},
- {10000, 10001, 0, false},
- {10001, 10000, 0, false},
- {-1000000, -1000001, 0, true},
- {-1000001, -1000000, 0, true},
- {-10000, -10001, 0, false},
- {-10001, -10000, 0, false},
- {1.0000001, 1.0000002, 0, true},
- {1.0000002, 1.0000001, 0, true},
- {1.0002, 1.0001, 0, false},
- {1.0001, 1.0002, 0, false},
- {-1.000001, -1.000002, 0, true},
- {-1.000002, -1.000001, 0, true},
- {-1.0001, -1.0002, 0, false},
- {-1.0002, -1.0001, 0, false},
- {0.000000001000001, 0.000000001000002, 0, true},
- {0.000000001000002, 0.000000001000001, 0, true},
- {0.000000000001002, 0.000000000001001, 0, false},
- {0.000000000001001, 0.000000000001002, 0, false},
- {-0.000000001000001, -0.000000001000002, 0, true},
- {-0.000000001000002, -0.000000001000001, 0, true},
- {-0.000000000001002, -0.000000000001001, 0, false},
- {-0.000000000001001, -0.000000000001002, 0, false},
- {0, 0, 0, true},
- {0, -0, 0, true},
- {-0, -0, 0, true},
- {0.00000001, 0, 0, false},
- {0, 0.00000001, 0, false},
- {-0.00000001, 0, 0, false},
- {0, -0.00000001, 0, false},
- {0, 1e-310, 0.01, true},
- {1e-310, 0, 0.01, true},
- {1e-310, 0, 0.000001, false},
- {0, 1e-310, 0.000001, false},
- {0, -1e-310, 0.1, true},
- {-1e-310, 0, 0.1, true},
- {-1e-310, 0, 0.00000001, false},
- {0, -1e-310, 0.00000001, false},
- {math.Inf(1), math.Inf(1), 0, true},
- {math.Inf(-1), math.Inf(-1), 0, true},
- {math.Inf(-1), math.Inf(1), 0, false},
- {math.Inf(1), math.MaxFloat64, 0, false},
- {math.Inf(-1), -math.MaxFloat64, 0, false},
- {math.NaN(), math.NaN(), 0, false},
- {math.NaN(), 0, 0, false},
- {-0, math.NaN(), 0, false},
- {math.NaN(), -0, 0, false},
- {0, math.NaN(), 0, false},
- {math.NaN(), math.Inf(1), 0, false},
- {math.Inf(1), math.NaN(), 0, false},
- {math.NaN(), math.Inf(-1), 0, false},
- {math.Inf(-1), math.NaN(), 0, false},
- {math.NaN(), math.MaxFloat64, 0, false},
- {math.MaxFloat64, math.NaN(), 0, false},
- {math.NaN(), -math.MaxFloat64, 0, false},
- {-math.MaxFloat64, math.NaN(), 0, false},
- {math.NaN(), math.SmallestNonzeroFloat64, 0, false},
- {math.SmallestNonzeroFloat64, math.NaN(), 0, false},
- {math.NaN(), -math.SmallestNonzeroFloat64, 0, false},
- {-math.SmallestNonzeroFloat64, math.NaN(), 0, false},
- {1.000000001, -1.0, 0, false},
- {-1.0, 1.000000001, 0, false},
- {-1.000000001, 1.0, 0, false},
- {1.0, -1.000000001, 0, false},
- {10 * math.SmallestNonzeroFloat64, 10 * -math.SmallestNonzeroFloat64, 0, true},
- {1e11 * math.SmallestNonzeroFloat64, 1e11 * -math.SmallestNonzeroFloat64, 0, false},
- {math.SmallestNonzeroFloat64, -math.SmallestNonzeroFloat64, 0, true},
- {-math.SmallestNonzeroFloat64, math.SmallestNonzeroFloat64, 0, true},
- {math.SmallestNonzeroFloat64, 0, 0, true},
- {0, math.SmallestNonzeroFloat64, 0, true},
- {-math.SmallestNonzeroFloat64, 0, 0, true},
- {0, -math.SmallestNonzeroFloat64, 0, true},
- {0.000000001, -math.SmallestNonzeroFloat64, 0, false},
- {0.000000001, math.SmallestNonzeroFloat64, 0, false},
- {math.SmallestNonzeroFloat64, 0.000000001, 0, false},
- {-math.SmallestNonzeroFloat64, 0.000000001, 0, false},
- }
- for _, ts := range equalityTests {
- if ts.tol == 0 {
- ts.tol = 1e-5
- }
- if equal := EqualWithinRel(ts.a, ts.b, ts.tol); equal != ts.equal {
- t.Errorf("Relative equality of %g and %g with tolerance %g returned: %v. Expected: %v",
- ts.a, ts.b, ts.tol, equal, ts.equal)
- }
- }
-}
-
-func nextAfterN(x, y float64, n int) float64 {
- for i := 0; i < n; i++ {
- x = math.Nextafter(x, y)
- }
- return x
-}
-
-func TestEqualsULP(t *testing.T) {
- if f := 67329.242; !EqualWithinULP(f, nextAfterN(f, math.Inf(1), 10), 10) {
- t.Errorf("Equal values returned as unequal")
- }
- if f := 67329.242; EqualWithinULP(f, nextAfterN(f, math.Inf(1), 5), 1) {
- t.Errorf("Unequal values returned as equal")
- }
- if f := 67329.242; EqualWithinULP(nextAfterN(f, math.Inf(1), 5), f, 1) {
- t.Errorf("Unequal values returned as equal")
- }
- if f := nextAfterN(0, math.Inf(1), 2); !EqualWithinULP(f, nextAfterN(f, math.Inf(-1), 5), 10) {
- t.Errorf("Equal values returned as unequal")
- }
- if !EqualWithinULP(67329.242, 67329.242, 10) {
- t.Errorf("Equal float64s not returned as equal")
- }
- if EqualWithinULP(1, math.NaN(), 10) {
- t.Errorf("NaN returned as equal")
- }
-
-}
-
-func TestEqualLengths(t *testing.T) {
- s1 := []float64{1, 2, 3, 4}
- s2 := []float64{1, 2, 3, 4}
- s3 := []float64{1, 2, 3}
- if !EqualLengths(s1, s2) {
- t.Errorf("Equal lengths returned as unequal")
- }
- if EqualLengths(s1, s3) {
- t.Errorf("Unequal lengths returned as equal")
- }
- if !EqualLengths(s1) {
- t.Errorf("Single slice returned as unequal")
- }
- if !EqualLengths() {
- t.Errorf("No slices returned as unequal")
- }
-}
-
-func eqIntSlice(one, two []int) string {
- if len(one) != len(two) {
- return "Length mismatch"
- }
- for i, val := range one {
- if val != two[i] {
- return "Index " + strconv.Itoa(i) + " mismatch"
- }
- }
- return ""
-}
-
-func TestFind(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- f := func(v float64) bool { return v > 3.5 }
- allTrueInds := []int{1, 3, 4}
-
- // Test finding first two elements
- inds, err := Find(nil, f, s, 2)
- if err != nil {
- t.Errorf("Find first two: Improper error return")
- }
- trueInds := allTrueInds[:2]
- str := eqIntSlice(inds, trueInds)
- if str != "" {
- t.Errorf("Find first two: " + str)
- }
-
- // Test finding no elements with non nil slice
- inds = []int{1, 2, 3, 4, 5, 6}
- inds, err = Find(inds, f, s, 0)
- if err != nil {
- t.Errorf("Find no elements: Improper error return")
- }
- str = eqIntSlice(inds, []int{})
- if str != "" {
- t.Errorf("Find no non-nil: " + str)
- }
-
- // Test finding first two elements with non nil slice
- inds = []int{1, 2, 3, 4, 5, 6}
- inds, err = Find(inds, f, s, 2)
- if err != nil {
- t.Errorf("Find first two non-nil: Improper error return")
- }
- str = eqIntSlice(inds, trueInds)
- if str != "" {
- t.Errorf("Find first two non-nil: " + str)
- }
-
- // Test finding too many elements
- inds, err = Find(inds, f, s, 4)
- if err == nil {
- t.Errorf("Request too many: No error returned")
- }
- str = eqIntSlice(inds, allTrueInds)
- if str != "" {
- t.Errorf("Request too many: Does not match all of the inds: " + str)
- }
-
- // Test finding all elements
- inds, err = Find(nil, f, s, -1)
- if err != nil {
- t.Errorf("Find all: Improper error returned")
- }
- str = eqIntSlice(inds, allTrueInds)
- if str != "" {
- t.Errorf("Find all: Does not match all of the inds: " + str)
- }
-}
-
-func TestHasNaN(t *testing.T) {
- for i, test := range []struct {
- s []float64
- ans bool
- }{
- {},
- {
- s: []float64{1, 2, 3, 4},
- },
- {
- s: []float64{1, math.NaN(), 3, 4},
- ans: true,
- },
- {
- s: []float64{1, 2, 3, math.NaN()},
- ans: true,
- },
- } {
- b := HasNaN(test.s)
- if b != test.ans {
- t.Errorf("HasNaN mismatch case %d. Expected %v, Found %v", i, test.ans, b)
- }
- }
-}
-
-func TestLogSpan(t *testing.T) {
- receiver1 := make([]float64, 6)
- truth := []float64{0.001, 0.01, 0.1, 1, 10, 100}
- receiver2 := LogSpan(receiver1, 0.001, 100)
- tst := make([]float64, 6)
- for i := range truth {
- tst[i] = receiver1[i] / truth[i]
- }
- comp := make([]float64, 6)
- for i := range comp {
- comp[i] = 1
- }
- AreSlicesEqual(t, comp, tst, "Improper logspace from mutator")
-
- for i := range truth {
- tst[i] = receiver2[i] / truth[i]
- }
- AreSlicesEqual(t, comp, tst, "Improper logspace from returned slice")
-
- if !Panics(func() { LogSpan(nil, 1, 5) }) {
- t.Errorf("Span accepts nil argument")
- }
- if !Panics(func() { LogSpan(make([]float64, 1), 1, 5) }) {
- t.Errorf("Span accepts argument of len = 1")
- }
-}
-
-func TestLogSumExp(t *testing.T) {
- s := []float64{1, 2, 3, 4, 5}
- val := LogSumExp(s)
- // http://www.wolframalpha.com/input/?i=log%28exp%281%29+%2B+exp%282%29+%2B+exp%283%29+%2B+exp%284%29+%2B+exp%285%29%29
- truth := 5.4519143959375933331957225109748087179338972737576824
- if math.Abs(val-truth) > EqTolerance {
- t.Errorf("Wrong logsumexp for many values")
- }
- s = []float64{1, 2}
- // http://www.wolframalpha.com/input/?i=log%28exp%281%29+%2B+exp%282%29%29
- truth = 2.3132616875182228340489954949678556419152800856703483
- val = LogSumExp(s)
- if math.Abs(val-truth) > EqTolerance {
- t.Errorf("Wrong logsumexp for two values. %v expected, %v found", truth, val)
- }
- // This case would normally underflow
- s = []float64{-1001, -1002, -1003, -1004, -1005}
- // http://www.wolframalpha.com/input/?i=log%28exp%28-1001%29%2Bexp%28-1002%29%2Bexp%28-1003%29%2Bexp%28-1004%29%2Bexp%28-1005%29%29
- truth = -1000.54808560406240666680427748902519128206610272624
- val = LogSumExp(s)
- if math.Abs(val-truth) > EqTolerance {
- t.Errorf("Doesn't match for underflow case. %v expected, %v found", truth, val)
- }
- // positive infinite case
- s = []float64{1, 2, 3, 4, 5, math.Inf(1)}
- val = LogSumExp(s)
- truth = math.Inf(1)
- if val != truth {
- t.Errorf("Doesn't match for pos Infinity case. %v expected, %v found", truth, val)
- }
- // negative infinite case
- s = []float64{1, 2, 3, 4, 5, math.Inf(-1)}
- val = LogSumExp(s)
- truth = 5.4519143959375933331957225109748087179338972737576824 // same as first case
- if math.Abs(val-truth) > EqTolerance {
- t.Errorf("Wrong logsumexp for values with negative infinity")
- }
-
-}
-
-func TestMaxAndIdx(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- ind := MaxIdx(s)
- val := Max(s)
- if val != 7 {
- t.Errorf("Wrong value returned")
- }
- if ind != 3 {
- t.Errorf("Wrong index returned")
- }
-}
-
-func TestMinAndIdx(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- ind := MinIdx(s)
- val := Min(s)
- if val != 1 {
- t.Errorf("Wrong value returned")
- }
- if ind != 2 {
- t.Errorf("Wrong index returned")
- }
-}
-
-func TestMul(t *testing.T) {
- s1 := []float64{1, 2, 3}
- s2 := []float64{1, 2, 3}
- ans := []float64{1, 4, 9}
- Mul(s1, s2)
- if !EqualApprox(s1, ans, EqTolerance) {
- t.Errorf("Mul doesn't give correct answer")
- }
- s1short := []float64{1}
- if !Panics(func() { Mul(s1short, s2) }) {
- t.Errorf("Did not panic with unequal lengths")
- }
- s2short := []float64{1}
- if !Panics(func() { Mul(s1, s2short) }) {
- t.Errorf("Did not panic with unequal lengths")
- }
-}
-
-func TestMulTo(t *testing.T) {
- s1 := []float64{1, 2, 3}
- s1orig := []float64{1, 2, 3}
- s2 := []float64{1, 2, 3}
- s2orig := []float64{1, 2, 3}
- dst1 := make([]float64, 3)
- ans := []float64{1, 4, 9}
- dst2 := MulTo(dst1, s1, s2)
- if !EqualApprox(dst1, ans, EqTolerance) {
- t.Errorf("MulTo doesn't give correct answer in mutated slice")
- }
- if !EqualApprox(dst2, ans, EqTolerance) {
- t.Errorf("MulTo doesn't give correct answer in returned slice")
- }
- if !EqualApprox(s1, s1orig, EqTolerance) {
- t.Errorf("S1 changes during multo")
- }
- if !EqualApprox(s2, s2orig, EqTolerance) {
- t.Errorf("s2 changes during multo")
- }
- MulTo(dst1, s1, s2)
- if !EqualApprox(dst1, ans, EqTolerance) {
- t.Errorf("MulTo doesn't give correct answer reusing dst")
- }
- dstShort := []float64{1}
- if !Panics(func() { MulTo(dstShort, s1, s2) }) {
- t.Errorf("Did not panic with s1 wrong length")
- }
- s1short := []float64{1}
- if !Panics(func() { MulTo(dst1, s1short, s2) }) {
- t.Errorf("Did not panic with s1 wrong length")
- }
- s2short := []float64{1}
- if !Panics(func() { MulTo(dst1, s1, s2short) }) {
- t.Errorf("Did not panic with s2 wrong length")
- }
-}
-
-func TestNaNWith(t *testing.T) {
- tests := []struct {
- payload uint64
- bits uint64
- }{
- {0, math.Float64bits(0 / func() float64 { return 0 }())}, // Hide the division by zero from the compiler.
- {1, math.Float64bits(math.NaN())},
- {1954, 0x7ff80000000007a2}, // R NA.
- }
-
- for _, test := range tests {
- nan := NaNWith(test.payload)
- if !math.IsNaN(nan) {
- t.Errorf("expected NaN value, got:%f", nan)
- }
-
- bits := math.Float64bits(nan)
-
- // Strip sign bit.
- const sign = 1 << 63
- bits &^= sign
- test.bits &^= sign
-
- if bits != test.bits {
- t.Errorf("expected NaN bit representation: got:%x want:%x", bits, test.bits)
- }
- }
-}
-
-func TestNaNPayload(t *testing.T) {
- tests := []struct {
- f float64
- payload uint64
- ok bool
- }{
- {0 / func() float64 { return 0 }(), 0, true}, // Hide the division by zero from the compiler.
-
- // The following two line are written explicitly to defend against potential changes to math.Copysign.
- {math.Float64frombits(math.Float64bits(math.NaN()) | (1 << 63)), 1, true}, // math.Copysign(math.NaN(), -1)
- {math.Float64frombits(math.Float64bits(math.NaN()) &^ (1 << 63)), 1, true}, // math.Copysign(math.NaN(), 1)
-
- {NaNWith(1954), 1954, true}, // R NA.
-
- {math.Copysign(0, -1), 0, false},
- {0, 0, false},
- {math.Inf(-1), 0, false},
- {math.Inf(1), 0, false},
-
- {math.Float64frombits(0x7ff0000000000001), 0, false}, // Signalling NaN.
- }
-
- for _, test := range tests {
- payload, ok := NaNPayload(test.f)
- if payload != test.payload {
- t.Errorf("expected NaN payload: got:%x want:%x", payload, test.payload)
- }
- if ok != test.ok {
- t.Errorf("expected NaN status: got:%t want:%t", ok, test.ok)
- }
- }
-}
-
-func TestNearest(t *testing.T) {
- s := []float64{6.2, 3, 5, 6.2, 8}
- ind := Nearest(s, 2.0)
- if ind != 1 {
- t.Errorf("Wrong index returned when value is less than all of elements")
- }
- ind = Nearest(s, 9.0)
- if ind != 4 {
- t.Errorf("Wrong index returned when value is greater than all of elements")
- }
- ind = Nearest(s, 3.1)
- if ind != 1 {
- t.Errorf("Wrong index returned when value is greater than closest element")
- }
- ind = Nearest(s, 3.1)
- if ind != 1 {
- t.Errorf("Wrong index returned when value is greater than closest element")
- }
- ind = Nearest(s, 2.9)
- if ind != 1 {
- t.Errorf("Wrong index returned when value is less than closest element")
- }
- ind = Nearest(s, 3)
- if ind != 1 {
- t.Errorf("Wrong index returned when value is equal to element")
- }
- ind = Nearest(s, 6.2)
- if ind != 0 {
- t.Errorf("Wrong index returned when value is equal to several elements")
- }
- ind = Nearest(s, 4)
- if ind != 1 {
- t.Errorf("Wrong index returned when value is exactly between two closest elements")
- }
-}
-
-func TestNearestWithinSpan(t *testing.T) {
- if !Panics(func() { NearestWithinSpan(10, 8, 2, 4.5) }) {
- t.Errorf("Did not panic when upper bound is lower than greater bound")
- }
- for i, test := range []struct {
- length int
- lower float64
- upper float64
- value float64
- idx int
- }{
- {
- length: 13,
- lower: 7,
- upper: 8.2,
- value: 6,
- idx: -1,
- },
- {
- length: 13,
- lower: 7,
- upper: 8.2,
- value: 10,
- idx: -1,
- },
- {
- length: 13,
- lower: 7,
- upper: 8.2,
- value: 7.19,
- idx: 2,
- },
- {
- length: 13,
- lower: 7,
- upper: 8.2,
- value: 7.21,
- idx: 2,
- },
- {
- length: 13,
- lower: 7,
- upper: 8.2,
- value: 7.2,
- idx: 2,
- },
- {
- length: 13,
- lower: 7,
- upper: 8.2,
- value: 7.151,
- idx: 2,
- },
- {
- length: 13,
- lower: 7,
- upper: 8.2,
- value: 7.249,
- idx: 2,
- },
- } {
- if idx := NearestWithinSpan(test.length, test.lower, test.upper, test.value); test.idx != idx {
- t.Errorf("Case %v mismatch: Want: %v, Got: %v", i, test.idx, idx)
- }
- }
-}
-
-func TestNorm(t *testing.T) {
- s := []float64{-1, -3.4, 5, -6}
- val := Norm(s, math.Inf(1))
- truth := 6.0
- if math.Abs(val-truth) > EqTolerance {
- t.Errorf("Doesn't match for inf norm. %v expected, %v found", truth, val)
- }
- // http://www.wolframalpha.com/input/?i=%28%28-1%29%5E2+%2B++%28-3.4%29%5E2+%2B+5%5E2%2B++6%5E2%29%5E%281%2F2%29
- val = Norm(s, 2)
- truth = 8.5767126569566267590651614132751986658027271236078592
- if math.Abs(val-truth) > EqTolerance {
- t.Errorf("Doesn't match for inf norm. %v expected, %v found", truth, val)
- }
- // http://www.wolframalpha.com/input/?i=%28%28%7C-1%7C%29%5E3+%2B++%28%7C-3.4%7C%29%5E3+%2B+%7C5%7C%5E3%2B++%7C6%7C%5E3%29%5E%281%2F3%29
- val = Norm(s, 3)
- truth = 7.2514321388020228478109121239004816430071237369356233
- if math.Abs(val-truth) > EqTolerance {
- t.Errorf("Doesn't match for inf norm. %v expected, %v found", truth, val)
- }
-
- //http://www.wolframalpha.com/input/?i=%7C-1%7C+%2B+%7C-3.4%7C+%2B+%7C5%7C%2B++%7C6%7C
- val = Norm(s, 1)
- truth = 15.4
- if math.Abs(val-truth) > EqTolerance {
- t.Errorf("Doesn't match for inf norm. %v expected, %v found", truth, val)
- }
-}
-
-func TestProd(t *testing.T) {
- s := []float64{}
- val := Prod(s)
- if val != 1 {
- t.Errorf("Val not returned as default when slice length is zero")
- }
- s = []float64{3, 4, 1, 7, 5}
- val = Prod(s)
- if val != 420 {
- t.Errorf("Wrong prod returned. Expected %v returned %v", 420, val)
- }
-}
-
-func TestReverse(t *testing.T) {
- for _, s := range [][]float64{
- {0},
- {1, 0},
- {2, 1, 0},
- {3, 2, 1, 0},
- {9, 8, 7, 6, 5, 4, 3, 2, 1, 0},
- } {
- Reverse(s)
- for i, v := range s {
- if v != float64(i) {
- t.Errorf("unexpected values for element %d: got:%v want:%v", i, v, i)
- }
- }
- }
-}
-
-func TestRound(t *testing.T) {
- for _, test := range []struct {
- x float64
- prec int
- want float64
- }{
- {x: 0, prec: 1, want: 0},
- {x: math.Inf(1), prec: 1, want: math.Inf(1)},
- {x: math.NaN(), prec: 1, want: math.NaN()},
- {x: func() float64 { var f float64; return -f }(), prec: 1, want: 0},
- {x: math.MaxFloat64 / 2, prec: 1, want: math.MaxFloat64 / 2},
- {x: 1 << 64, prec: 1, want: 1 << 64},
- {x: 454.4445, prec: 3, want: 454.445},
- {x: 454.44445, prec: 4, want: 454.4445},
- {x: 0.42499, prec: 4, want: 0.425},
- {x: 0.42599, prec: 4, want: 0.426},
- {x: 0.424999999999993, prec: 2, want: 0.42},
- {x: 0.425, prec: 2, want: 0.43},
- {x: 0.425000000000001, prec: 2, want: 0.43},
- {x: 123.4244999999999, prec: 3, want: 123.424},
- {x: 123.4245, prec: 3, want: 123.425},
- {x: 123.4245000000001, prec: 3, want: 123.425},
-
- {x: 454.45, prec: 0, want: 454},
- {x: 454.45, prec: 1, want: 454.5},
- {x: 454.45, prec: 2, want: 454.45},
- {x: 454.45, prec: 3, want: 454.45},
- {x: 454.445, prec: 0, want: 454},
- {x: 454.445, prec: 1, want: 454.4},
- {x: 454.445, prec: 2, want: 454.45},
- {x: 454.445, prec: 3, want: 454.445},
- {x: 454.445, prec: 4, want: 454.445},
- {x: 454.55, prec: 0, want: 455},
- {x: 454.55, prec: 1, want: 454.6},
- {x: 454.55, prec: 2, want: 454.55},
- {x: 454.55, prec: 3, want: 454.55},
- {x: 454.455, prec: 0, want: 454},
- {x: 454.455, prec: 1, want: 454.5},
- {x: 454.455, prec: 2, want: 454.46},
- {x: 454.455, prec: 3, want: 454.455},
- {x: 454.455, prec: 4, want: 454.455},
-
- // Negative precision.
- {x: 454.45, prec: -1, want: 450},
- {x: 454.45, prec: -2, want: 500},
- {x: 500, prec: -3, want: 1000},
- {x: 500, prec: -4, want: 0},
- {x: 1500, prec: -3, want: 2000},
- {x: 1500, prec: -4, want: 0},
- } {
- for _, sign := range []float64{1, -1} {
- got := Round(sign*test.x, test.prec)
- want := sign * test.want
- if want == 0 {
- want = 0
- }
- if (got != want || math.Signbit(got) != math.Signbit(want)) && !(math.IsNaN(got) && math.IsNaN(want)) {
- t.Errorf("unexpected result for Round(%g, %d): got: %g, want: %g", sign*test.x, test.prec, got, want)
- }
- }
- }
-}
-
-func TestRoundEven(t *testing.T) {
- for _, test := range []struct {
- x float64
- prec int
- want float64
- }{
- {x: 0, prec: 1, want: 0},
- {x: math.Inf(1), prec: 1, want: math.Inf(1)},
- {x: math.NaN(), prec: 1, want: math.NaN()},
- {x: func() float64 { var f float64; return -f }(), prec: 1, want: 0},
- {x: math.MaxFloat64 / 2, prec: 1, want: math.MaxFloat64 / 2},
- {x: 1 << 64, prec: 1, want: 1 << 64},
- {x: 454.4445, prec: 3, want: 454.444},
- {x: 454.44445, prec: 4, want: 454.4444},
- {x: 0.42499, prec: 4, want: 0.425},
- {x: 0.42599, prec: 4, want: 0.426},
- {x: 0.424999999999993, prec: 2, want: 0.42},
- {x: 0.425, prec: 2, want: 0.42},
- {x: 0.425000000000001, prec: 2, want: 0.43},
- {x: 123.4244999999999, prec: 3, want: 123.424},
- {x: 123.4245, prec: 3, want: 123.424},
- {x: 123.4245000000001, prec: 3, want: 123.425},
-
- {x: 454.45, prec: 0, want: 454},
- {x: 454.45, prec: 1, want: 454.4},
- {x: 454.45, prec: 2, want: 454.45},
- {x: 454.45, prec: 3, want: 454.45},
- {x: 454.445, prec: 0, want: 454},
- {x: 454.445, prec: 1, want: 454.4},
- {x: 454.445, prec: 2, want: 454.44},
- {x: 454.445, prec: 3, want: 454.445},
- {x: 454.445, prec: 4, want: 454.445},
- {x: 454.55, prec: 0, want: 455},
- {x: 454.55, prec: 1, want: 454.6},
- {x: 454.55, prec: 2, want: 454.55},
- {x: 454.55, prec: 3, want: 454.55},
- {x: 454.455, prec: 0, want: 454},
- {x: 454.455, prec: 1, want: 454.5},
- {x: 454.455, prec: 2, want: 454.46},
- {x: 454.455, prec: 3, want: 454.455},
- {x: 454.455, prec: 4, want: 454.455},
-
- // Negative precision.
- {x: 454.45, prec: -1, want: 450},
- {x: 454.45, prec: -2, want: 500},
- {x: 500, prec: -3, want: 0},
- {x: 500, prec: -4, want: 0},
- {x: 1500, prec: -3, want: 2000},
- {x: 1500, prec: -4, want: 0},
- } {
- for _, sign := range []float64{1, -1} {
- got := RoundEven(sign*test.x, test.prec)
- want := sign * test.want
- if want == 0 {
- want = 0
- }
- if (got != want || math.Signbit(got) != math.Signbit(want)) && !(math.IsNaN(got) && math.IsNaN(want)) {
- t.Errorf("unexpected result for RoundEven(%g, %d): got: %g, want: %g", sign*test.x, test.prec, got, want)
- }
- }
- }
-}
-
-func TestSame(t *testing.T) {
- s1 := []float64{1, 2, 3, 4}
- s2 := []float64{1, 2, 3, 4}
- if !Same(s1, s2) {
- t.Errorf("Equal slices returned as unequal")
- }
- s2 = []float64{1, 2, 3, 4 + 1e-14}
- if Same(s1, s2) {
- t.Errorf("Unequal slices returned as equal")
- }
- if Same(s1, []float64{}) {
- t.Errorf("Unequal slice lengths returned as equal")
- }
- s1 = []float64{1, 2, math.NaN(), 4}
- s2 = []float64{1, 2, math.NaN(), 4}
- if !Same(s1, s2) {
- t.Errorf("Slices with matching NaN values returned as unequal")
- }
- s1 = []float64{1, 2, math.NaN(), 4}
- s2 = []float64{1, math.NaN(), 3, 4}
- if !Same(s1, s2) {
- t.Errorf("Slices with unmatching NaN values returned as equal")
- }
-}
-
-func TestScale(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- c := 5.0
- truth := []float64{15, 20, 5, 35, 25}
- Scale(c, s)
- AreSlicesEqual(t, truth, s, "Bad scaling")
-}
-
-func TestSpan(t *testing.T) {
- receiver1 := make([]float64, 5)
- truth := []float64{1, 2, 3, 4, 5}
- receiver2 := Span(receiver1, 1, 5)
- AreSlicesEqual(t, truth, receiver1, "Improper linspace from mutator")
- AreSlicesEqual(t, truth, receiver2, "Improper linspace from returned slice")
- receiver1 = make([]float64, 6)
- truth = []float64{0, 0.2, 0.4, 0.6, 0.8, 1.0}
- Span(receiver1, 0, 1)
- AreSlicesEqual(t, truth, receiver1, "Improper linspace")
- if !Panics(func() { Span(nil, 1, 5) }) {
- t.Errorf("Span accepts nil argument")
- }
- if !Panics(func() { Span(make([]float64, 1), 1, 5) }) {
- t.Errorf("Span accepts argument of len = 1")
- }
-}
-
-func TestSub(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- v := []float64{1, 2, 3, 4, 5}
- truth := []float64{2, 2, -2, 3, 0}
- Sub(s, v)
- AreSlicesEqual(t, truth, s, "Bad subtract")
- // Test that it panics
- if !Panics(func() { Sub(make([]float64, 2), make([]float64, 3)) }) {
- t.Errorf("Did not panic with length mismatch")
- }
-}
-
-func TestSubTo(t *testing.T) {
- s := []float64{3, 4, 1, 7, 5}
- v := []float64{1, 2, 3, 4, 5}
- truth := []float64{2, 2, -2, 3, 0}
- dst1 := make([]float64, len(s))
- dst2 := SubTo(dst1, s, v)
- AreSlicesEqual(t, truth, dst1, "Bad subtract from mutator")
- AreSlicesEqual(t, truth, dst2, "Bad subtract from returned slice")
- // Test that all mismatch combinations panic
- if !Panics(func() { SubTo(make([]float64, 2), make([]float64, 3), make([]float64, 3)) }) {
- t.Errorf("Did not panic with dst different length")
- }
- if !Panics(func() { SubTo(make([]float64, 3), make([]float64, 2), make([]float64, 3)) }) {
- t.Errorf("Did not panic with subtractor different length")
- }
- if !Panics(func() { SubTo(make([]float64, 3), make([]float64, 3), make([]float64, 2)) }) {
- t.Errorf("Did not panic with subtractee different length")
- }
-}
-
-func TestSum(t *testing.T) {
- s := []float64{}
- val := Sum(s)
- if val != 0 {
- t.Errorf("Val not returned as default when slice length is zero")
- }
- s = []float64{3, 4, 1, 7, 5}
- val = Sum(s)
- if val != 20 {
- t.Errorf("Wrong sum returned")
- }
-}
-
-func TestWithin(t *testing.T) {
- for i, test := range []struct {
- s []float64
- v float64
- idx int
- panics bool
- }{
- {
- s: []float64{1, 2, 5, 9},
- v: 1,
- idx: 0,
- },
- {
- s: []float64{1, 2, 5, 9},
- v: 9,
- idx: -1,
- },
- {
- s: []float64{1, 2, 5, 9},
- v: 1.5,
- idx: 0,
- },
- {
- s: []float64{1, 2, 5, 9},
- v: 2,
- idx: 1,
- },
- {
- s: []float64{1, 2, 5, 9},
- v: 2.5,
- idx: 1,
- },
- {
- s: []float64{1, 2, 5, 9},
- v: -3,
- idx: -1,
- },
- {
- s: []float64{1, 2, 5, 9},
- v: 15,
- idx: -1,
- },
- {
- s: []float64{1, 2, 5, 9},
- v: math.NaN(),
- idx: -1,
- },
- {
- s: []float64{5, 2, 6},
- panics: true,
- },
- {
- panics: true,
- },
- {
- s: []float64{1},
- panics: true,
- },
- } {
- var idx int
- panics := Panics(func() { idx = Within(test.s, test.v) })
- if panics {
- if !test.panics {
- t.Errorf("Case %v: bad panic", i)
- }
- continue
- }
- if test.panics {
- if !panics {
- t.Errorf("Case %v: did not panic when it should", i)
- }
- continue
- }
- if idx != test.idx {
- t.Errorf("Case %v: Idx mismatch. Want: %v, got: %v", i, test.idx, idx)
- }
- }
-
-}
-
-func randomSlice(l int) []float64 {
- s := make([]float64, l)
- for i := range s {
- s[i] = rand.Float64()
- }
- return s
-}
-
-func benchmarkMin(b *testing.B, size int) {
- s := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- Min(s)
- }
-}
-func BenchmarkMinSmall(b *testing.B) { benchmarkMin(b, Small) }
-func BenchmarkMinMed(b *testing.B) { benchmarkMin(b, Medium) }
-func BenchmarkMinLarge(b *testing.B) { benchmarkMin(b, Large) }
-func BenchmarkMinHuge(b *testing.B) { benchmarkMin(b, Huge) }
-
-func benchmarkAdd(b *testing.B, size int) {
- s1 := randomSlice(size)
- s2 := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- Add(s1, s2)
- }
-}
-func BenchmarkAddSmall(b *testing.B) { benchmarkAdd(b, Small) }
-func BenchmarkAddMed(b *testing.B) { benchmarkAdd(b, Medium) }
-func BenchmarkAddLarge(b *testing.B) { benchmarkAdd(b, Large) }
-func BenchmarkAddHuge(b *testing.B) { benchmarkAdd(b, Huge) }
-
-func benchmarkAddTo(b *testing.B, size int) {
- s1 := randomSlice(size)
- s2 := randomSlice(size)
- dst := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- AddTo(dst, s1, s2)
- }
-}
-func BenchmarkAddToSmall(b *testing.B) { benchmarkAddTo(b, Small) }
-func BenchmarkAddToMed(b *testing.B) { benchmarkAddTo(b, Medium) }
-func BenchmarkAddToLarge(b *testing.B) { benchmarkAddTo(b, Large) }
-func BenchmarkAddToHuge(b *testing.B) { benchmarkAddTo(b, Huge) }
-
-func benchmarkCumProd(b *testing.B, size int) {
- s := randomSlice(size)
- dst := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- CumProd(dst, s)
- }
-}
-func BenchmarkCumProdSmall(b *testing.B) { benchmarkCumProd(b, Small) }
-func BenchmarkCumProdMed(b *testing.B) { benchmarkCumProd(b, Medium) }
-func BenchmarkCumProdLarge(b *testing.B) { benchmarkCumProd(b, Large) }
-func BenchmarkCumProdHuge(b *testing.B) { benchmarkCumProd(b, Huge) }
-
-func benchmarkCumSum(b *testing.B, size int) {
- s := randomSlice(size)
- dst := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- CumSum(dst, s)
- }
-}
-func BenchmarkCumSumSmall(b *testing.B) { benchmarkCumSum(b, Small) }
-func BenchmarkCumSumMed(b *testing.B) { benchmarkCumSum(b, Medium) }
-func BenchmarkCumSumLarge(b *testing.B) { benchmarkCumSum(b, Large) }
-func BenchmarkCumSumHuge(b *testing.B) { benchmarkCumSum(b, Huge) }
-
-func benchmarkDiv(b *testing.B, size int) {
- s := randomSlice(size)
- dst := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- Div(dst, s)
- }
-}
-func BenchmarkDivSmall(b *testing.B) { benchmarkDiv(b, Small) }
-func BenchmarkDivMed(b *testing.B) { benchmarkDiv(b, Medium) }
-func BenchmarkDivLarge(b *testing.B) { benchmarkDiv(b, Large) }
-func BenchmarkDivHuge(b *testing.B) { benchmarkDiv(b, Huge) }
-
-func benchmarkDivTo(b *testing.B, size int) {
- s1 := randomSlice(size)
- s2 := randomSlice(size)
- dst := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- DivTo(dst, s1, s2)
- }
-}
-func BenchmarkDivToSmall(b *testing.B) { benchmarkDivTo(b, Small) }
-func BenchmarkDivToMed(b *testing.B) { benchmarkDivTo(b, Medium) }
-func BenchmarkDivToLarge(b *testing.B) { benchmarkDivTo(b, Large) }
-func BenchmarkDivToHuge(b *testing.B) { benchmarkDivTo(b, Huge) }
-
-func benchmarkSub(b *testing.B, size int) {
- s1 := randomSlice(size)
- s2 := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- Sub(s1, s2)
- }
-}
-func BenchmarkSubSmall(b *testing.B) { benchmarkSub(b, Small) }
-func BenchmarkSubMed(b *testing.B) { benchmarkSub(b, Medium) }
-func BenchmarkSubLarge(b *testing.B) { benchmarkSub(b, Large) }
-func BenchmarkSubHuge(b *testing.B) { benchmarkSub(b, Huge) }
-
-func benchmarkSubTo(b *testing.B, size int) {
- s1 := randomSlice(size)
- s2 := randomSlice(size)
- dst := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- SubTo(dst, s1, s2)
- }
-}
-func BenchmarkSubToSmall(b *testing.B) { benchmarkSubTo(b, Small) }
-func BenchmarkSubToMed(b *testing.B) { benchmarkSubTo(b, Medium) }
-func BenchmarkSubToLarge(b *testing.B) { benchmarkSubTo(b, Large) }
-func BenchmarkSubToHuge(b *testing.B) { benchmarkSubTo(b, Huge) }
-
-func benchmarkLogSumExp(b *testing.B, size int) {
- s := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- LogSumExp(s)
- }
-}
-func BenchmarkLogSumExpSmall(b *testing.B) { benchmarkLogSumExp(b, Small) }
-func BenchmarkLogSumExpMed(b *testing.B) { benchmarkLogSumExp(b, Medium) }
-func BenchmarkLogSumExpLarge(b *testing.B) { benchmarkLogSumExp(b, Large) }
-func BenchmarkLogSumExpHuge(b *testing.B) { benchmarkLogSumExp(b, Huge) }
-
-func benchmarkDot(b *testing.B, size int) {
- s1 := randomSlice(size)
- s2 := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- Dot(s1, s2)
- }
-}
-func BenchmarkDotSmall(b *testing.B) { benchmarkDot(b, Small) }
-func BenchmarkDotMed(b *testing.B) { benchmarkDot(b, Medium) }
-func BenchmarkDotLarge(b *testing.B) { benchmarkDot(b, Large) }
-func BenchmarkDotHuge(b *testing.B) { benchmarkDot(b, Huge) }
-
-func benchmarkAddScaledTo(b *testing.B, size int) {
- dst := randomSlice(size)
- y := randomSlice(size)
- s := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i++ {
- AddScaledTo(dst, y, 2.3, s)
- }
-}
-func BenchmarkAddScaledToSmall(b *testing.B) { benchmarkAddScaledTo(b, Small) }
-func BenchmarkAddScaledToMedium(b *testing.B) { benchmarkAddScaledTo(b, Medium) }
-func BenchmarkAddScaledToLarge(b *testing.B) { benchmarkAddScaledTo(b, Large) }
-func BenchmarkAddScaledToHuge(b *testing.B) { benchmarkAddScaledTo(b, Huge) }
-
-func benchmarkScale(b *testing.B, size int) {
- dst := randomSlice(size)
- b.ResetTimer()
- for i := 0; i < b.N; i += 2 {
- Scale(2.0, dst)
- Scale(0.5, dst)
- }
-}
-func BenchmarkScaleSmall(b *testing.B) { benchmarkScale(b, Small) }
-func BenchmarkScaleMedium(b *testing.B) { benchmarkScale(b, Medium) }
-func BenchmarkScaleLarge(b *testing.B) { benchmarkScale(b, Large) }
-func BenchmarkScaleHuge(b *testing.B) { benchmarkScale(b, Huge) }