package/gonum
1
0
Fork 0
mirror of https://github.com/gonum/gonum.git synced 2025-10-05 23:26:52 +08:00
Code Issues Actions Packages Projects Releases Wiki Activity

lapack/testlapack: remove unused SymBand helpers

Browse Source
This commit is contained in:
Vladimir Chalupecky
2019-06-16 18:54:29 +02:00
committed by Vladimír Chalupecký
parent ce6986a678
commit d7f6bfa9af
1 changed files with 0 additions and 46 deletions
Download Patch File Download Diff File Expand all files Collapse all files

46
lapack/testlapack/general.go
View File

@@ -1076,52 +1076,6 @@ func distSymBand(uplo blas.Uplo, n, kd int, a []float64, lda int, b []float64, l
return dist return dist
} }
// symToSymBand takes the data in a Symmetric matrix and returns a
// SymmetricBanded matrix.
func symToSymBand(ul blas.Uplo, a []float64, n, lda, kb, ldab int) []float64 {
if ul == blas.Upper {
band := make([]float64, (n-1)*ldab+kb+1)
for i := 0; i < n; i++ {
for j := i; j < min(i+kb+1, n); j++ {
band[i*ldab+j-i] = a[i*lda+j]
}
}
return band
}
band := make([]float64, (n-1)*ldab+kb+1)
for i := 0; i < n; i++ {
for j := max(0, i-kb); j <= i; j++ {
band[i*ldab+j-i+kb] = a[i*lda+j]
}
}
return band
}
// symBandToSym takes a banded symmetric matrix and returns the same data as
// a Symmetric matrix.
func symBandToSym(ul blas.Uplo, band []float64, n, kb, ldab int) blas64.Symmetric {
sym := make([]float64, n*n)
if ul == blas.Upper {
for i := 0; i < n; i++ {
for j := 0; j < min(kb+1+i, n)-i; j++ {
sym[i*n+i+j] = band[i*ldab+j]
}
}
} else {
for i := 0; i < n; i++ {
for j := kb - min(i, kb); j < kb+1; j++ {
sym[i*n+i-kb+j] = band[i*ldab+j]
}
}
}
return blas64.Symmetric{
N: n,
Stride: n,
Data: sym,
Uplo: ul,
}
}
// eye returns an identity matrix of given order and stride. // eye returns an identity matrix of given order and stride.
func eye(n, stride int) blas64.General { func eye(n, stride int) blas64.General {
ans := nanGeneral(n, n, stride) ans := nanGeneral(n, n, stride)