1 // Code generated by "go generate gonum.org/v1/gonum/blas”; DO NOT EDIT.
3 // Copyright ©2015 The Gonum Authors. All rights reserved.
4 // Use of this source code is governed by a BSD-style
5 // license that can be found in the LICENSE file.
10 math "gonum.org/v1/gonum/internal/math32"
13 "gonum.org/v1/gonum/blas"
16 func newSymmetricFrom(a SymmetricCols) Symmetric {
20 Data: make([]float32, a.N*a.N),
27 func (m Symmetric) n() int { return m.N }
28 func (m Symmetric) at(i, j int) float32 {
29 if m.Uplo == blas.Lower && i < j && j < m.N {
32 if m.Uplo == blas.Upper && i > j {
35 return m.Data[i*m.Stride+j]
37 func (m Symmetric) uplo() blas.Uplo { return m.Uplo }
39 func newSymmetricColsFrom(a Symmetric) SymmetricCols {
43 Data: make([]float32, a.N*a.N),
50 func (m SymmetricCols) n() int { return m.N }
51 func (m SymmetricCols) at(i, j int) float32 {
52 if m.Uplo == blas.Lower && i < j {
55 if m.Uplo == blas.Upper && i > j && i < m.N {
58 return m.Data[i+j*m.Stride]
60 func (m SymmetricCols) uplo() blas.Uplo { return m.Uplo }
62 type symmetric interface {
68 func sameSymmetric(a, b symmetric) bool {
74 if a.uplo() != b.uplo() {
77 for i := 0; i < an; i++ {
78 for j := 0; j < an; j++ {
79 if a.at(i, j) != b.at(i, j) || math.IsNaN(a.at(i, j)) != math.IsNaN(b.at(i, j)) {
87 var symmetricTests = []Symmetric{
88 {N: 3, Stride: 3, Data: []float32{
93 {N: 3, Stride: 5, Data: []float32{
100 func TestConvertSymmetric(t *testing.T) {
101 for _, test := range symmetricTests {
102 for _, uplo := range []blas.Uplo{blas.Upper, blas.Lower} {
104 colmajor := newSymmetricColsFrom(test)
105 if !sameSymmetric(colmajor, test) {
106 t.Errorf("unexpected result for row major to col major conversion:\n\tgot: %#v\n\tfrom:%#v",
109 rowmajor := newSymmetricFrom(colmajor)
110 if !sameSymmetric(rowmajor, test) {
111 t.Errorf("unexpected result for col major to row major conversion:\n\tgot: %#v\n\twant:%#v",
117 func newSymmetricBandFrom(a SymmetricBandCols) SymmetricBand {
122 Data: make([]float32, a.N*(a.K+1)),
125 for i := range t.Data {
126 t.Data[i] = math.NaN()
132 func (m SymmetricBand) n() (n int) { return m.N }
133 func (m SymmetricBand) at(i, j int) float32 {
135 Rows: m.N, Cols: m.N,
141 panic("blas32: bad BLAS uplo")
155 func (m SymmetricBand) bandwidth() (k int) { return m.K }
156 func (m SymmetricBand) uplo() blas.Uplo { return m.Uplo }
158 func newSymmetricBandColsFrom(a SymmetricBand) SymmetricBandCols {
159 t := SymmetricBandCols{
163 Data: make([]float32, a.N*(a.K+1)),
166 for i := range t.Data {
167 t.Data[i] = math.NaN()
173 func (m SymmetricBandCols) n() (n int) { return m.N }
174 func (m SymmetricBandCols) at(i, j int) float32 {
176 Rows: m.N, Cols: m.N,
182 panic("blas32: bad BLAS uplo")
196 func (m SymmetricBandCols) bandwidth() (k int) { return m.K }
197 func (m SymmetricBandCols) uplo() blas.Uplo { return m.Uplo }
199 type symmetricBand interface {
206 func sameSymmetricBand(a, b symmetricBand) bool {
212 if a.uplo() != b.uplo() {
220 for i := 0; i < an; i++ {
221 for j := 0; j < an; j++ {
222 if a.at(i, j) != b.at(i, j) || math.IsNaN(a.at(i, j)) != math.IsNaN(b.at(i, j)) {
230 var symmetricBandTests = []SymmetricBand{
231 {N: 3, K: 0, Stride: 1, Uplo: blas.Upper, Data: []float32{
236 {N: 3, K: 0, Stride: 1, Uplo: blas.Lower, Data: []float32{
241 {N: 3, K: 1, Stride: 2, Uplo: blas.Upper, Data: []float32{
246 {N: 3, K: 1, Stride: 2, Uplo: blas.Lower, Data: []float32{
251 {N: 3, K: 2, Stride: 3, Uplo: blas.Upper, Data: []float32{
256 {N: 3, K: 2, Stride: 3, Uplo: blas.Lower, Data: []float32{
262 {N: 3, K: 0, Stride: 5, Uplo: blas.Upper, Data: []float32{
267 {N: 3, K: 0, Stride: 5, Uplo: blas.Lower, Data: []float32{
272 {N: 3, K: 1, Stride: 5, Uplo: blas.Upper, Data: []float32{
277 {N: 3, K: 1, Stride: 5, Uplo: blas.Lower, Data: []float32{
282 {N: 3, K: 2, Stride: 5, Uplo: blas.Upper, Data: []float32{
287 {N: 3, K: 2, Stride: 5, Uplo: blas.Lower, Data: []float32{
294 func TestConvertSymBand(t *testing.T) {
295 for _, test := range symmetricBandTests {
296 colmajor := newSymmetricBandColsFrom(test)
297 if !sameSymmetricBand(colmajor, test) {
298 t.Errorf("unexpected result for row major to col major conversion:\n\tgot: %#v\n\tfrom:%#v",
301 rowmajor := newSymmetricBandFrom(colmajor)
302 if !sameSymmetricBand(rowmajor, test) {
303 t.Errorf("unexpected result for col major to row major conversion:\n\tgot: %#v\n\twant:%#v",