1 // Copyright ©2016 The Gonum Authors. All rights reserved.
2 // Use of this source code is governed by a BSD-style
3 // license that can be found in the LICENSE file.
10 "golang.org/x/exp/rand"
12 "gonum.org/v1/gonum/blas/blas64"
13 "gonum.org/v1/gonum/lapack"
16 var resultGeneral blas64.General
18 func DgeevBenchmark(b *testing.B, impl Dgeever) {
19 rnd := rand.New(rand.NewSource(1))
20 benchmarks := []struct {
24 {"AntisymRandom3", NewAntisymRandom(3, rnd).Matrix()},
25 {"AntisymRandom4", NewAntisymRandom(4, rnd).Matrix()},
26 {"AntisymRandom5", NewAntisymRandom(5, rnd).Matrix()},
27 {"AntisymRandom10", NewAntisymRandom(10, rnd).Matrix()},
28 {"AntisymRandom50", NewAntisymRandom(50, rnd).Matrix()},
29 {"AntisymRandom100", NewAntisymRandom(100, rnd).Matrix()},
30 {"AntisymRandom200", NewAntisymRandom(200, rnd).Matrix()},
31 {"AntisymRandom500", NewAntisymRandom(500, rnd).Matrix()},
32 {"Circulant3", Circulant(3).Matrix()},
33 {"Circulant4", Circulant(4).Matrix()},
34 {"Circulant5", Circulant(5).Matrix()},
35 {"Circulant10", Circulant(10).Matrix()},
36 {"Circulant50", Circulant(50).Matrix()},
37 {"Circulant100", Circulant(100).Matrix()},
38 {"Circulant200", Circulant(200).Matrix()},
39 {"Circulant500", Circulant(500).Matrix()},
41 for _, bm := range benchmarks {
46 wr := make([]float64, n)
47 wi := make([]float64, n)
48 work := make([]float64, 1)
49 impl.Dgeev(lapack.ComputeLeftEV, lapack.ComputeRightEV, n, nil, n, nil, nil, nil, n, nil, n, work, -1)
50 work = make([]float64, int(work[0]))
51 b.Run(bm.name, func(b *testing.B) {
52 for i := 0; i < b.N; i++ {
56 impl.Dgeev(lapack.ComputeLeftEV, lapack.ComputeRightEV, n, a.Data, a.Stride, wr, wi,
57 vl.Data, vl.Stride, vr.Data, vr.Stride, work, len(work))