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.
13 "gonum.org/v1/gonum/blas"
16 func newHermitianFrom(a HermitianCols) Hermitian {
20 Data: make([]complex128, a.N*a.N),
27 func (m Hermitian) n() int { return m.N }
28 func (m Hermitian) at(i, j int) complex128 {
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 Hermitian) uplo() blas.Uplo { return m.Uplo }
39 func newHermitianColsFrom(a Hermitian) HermitianCols {
43 Data: make([]complex128, a.N*a.N),
50 func (m HermitianCols) n() int { return m.N }
51 func (m HermitianCols) at(i, j int) complex128 {
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 HermitianCols) uplo() blas.Uplo { return m.Uplo }
62 type hermitian interface {
64 at(i, j int) complex128
68 func sameHermitian(a, b hermitian) 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 hermitianTests = []Hermitian{
88 {N: 3, Stride: 3, Data: []complex128{
93 {N: 3, Stride: 5, Data: []complex128{
100 func TestConvertHermitian(t *testing.T) {
101 for _, test := range hermitianTests {
102 for _, uplo := range []blas.Uplo{blas.Upper, blas.Lower} {
104 colmajor := newHermitianColsFrom(test)
105 if !sameHermitian(colmajor, test) {
106 t.Errorf("unexpected result for row major to col major conversion:\n\tgot: %#v\n\tfrom:%#v",
109 rowmajor := newHermitianFrom(colmajor)
110 if !sameHermitian(rowmajor, test) {
111 t.Errorf("unexpected result for col major to row major conversion:\n\tgot: %#v\n\twant:%#v",
117 func newHermitianBandFrom(a HermitianBandCols) HermitianBand {
122 Data: make([]complex128, a.N*(a.K+1)),
125 for i := range t.Data {
126 t.Data[i] = math.NaN()
132 func (m HermitianBand) n() (n int) { return m.N }
133 func (m HermitianBand) at(i, j int) complex128 {
135 Rows: m.N, Cols: m.N,
141 panic("cblas128: bad BLAS uplo")
155 func (m HermitianBand) bandwidth() (k int) { return m.K }
156 func (m HermitianBand) uplo() blas.Uplo { return m.Uplo }
158 func newHermitianBandColsFrom(a HermitianBand) HermitianBandCols {
159 t := HermitianBandCols{
163 Data: make([]complex128, a.N*(a.K+1)),
166 for i := range t.Data {
167 t.Data[i] = math.NaN()
173 func (m HermitianBandCols) n() (n int) { return m.N }
174 func (m HermitianBandCols) at(i, j int) complex128 {
176 Rows: m.N, Cols: m.N,
182 panic("cblas128: bad BLAS uplo")
196 func (m HermitianBandCols) bandwidth() (k int) { return m.K }
197 func (m HermitianBandCols) uplo() blas.Uplo { return m.Uplo }
199 type hermitianBand interface {
201 at(i, j int) complex128
206 func sameHermitianBand(a, b hermitianBand) 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 hermitianBandTests = []HermitianBand{
231 {N: 3, K: 0, Stride: 1, Uplo: blas.Upper, Data: []complex128{
236 {N: 3, K: 0, Stride: 1, Uplo: blas.Lower, Data: []complex128{
241 {N: 3, K: 1, Stride: 2, Uplo: blas.Upper, Data: []complex128{
246 {N: 3, K: 1, Stride: 2, Uplo: blas.Lower, Data: []complex128{
251 {N: 3, K: 2, Stride: 3, Uplo: blas.Upper, Data: []complex128{
256 {N: 3, K: 2, Stride: 3, Uplo: blas.Lower, Data: []complex128{
262 {N: 3, K: 0, Stride: 5, Uplo: blas.Upper, Data: []complex128{
267 {N: 3, K: 0, Stride: 5, Uplo: blas.Lower, Data: []complex128{
272 {N: 3, K: 1, Stride: 5, Uplo: blas.Upper, Data: []complex128{
277 {N: 3, K: 1, Stride: 5, Uplo: blas.Lower, Data: []complex128{
282 {N: 3, K: 2, Stride: 5, Uplo: blas.Upper, Data: []complex128{
287 {N: 3, K: 2, Stride: 5, Uplo: blas.Lower, Data: []complex128{
294 func TestConvertHermBand(t *testing.T) {
295 for _, test := range hermitianBandTests {
296 colmajor := newHermitianBandColsFrom(test)
297 if !sameHermitianBand(colmajor, test) {
298 t.Errorf("unexpected result for row major to col major conversion:\n\tgot: %#v\n\tfrom:%#v",
301 rowmajor := newHermitianBandFrom(colmajor)
302 if !sameHermitianBand(rowmajor, test) {
303 t.Errorf("unexpected result for col major to row major conversion:\n\tgot: %#v\n\twant:%#v",