+++ /dev/null
-package testblas
-
-import (
- "testing"
-
- "gonum.org/v1/gonum/blas"
- "gonum.org/v1/gonum/floats"
-)
-
-type Dtrmmer interface {
- Dtrmm(s blas.Side, ul blas.Uplo, tA blas.Transpose, d blas.Diag, m, n int, alpha float64, a []float64, lda int, b []float64, ldb int)
-}
-
-func DtrmmTest(t *testing.T, blasser Dtrmmer) {
- for i, test := range []struct {
- s blas.Side
- ul blas.Uplo
- tA blas.Transpose
- d blas.Diag
- m int
- n int
- alpha float64
- a [][]float64
- b [][]float64
- ans [][]float64
- }{
- {
- s: blas.Left,
- ul: blas.Upper,
- tA: blas.NoTrans,
- d: blas.NonUnit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3, 4},
- {0, 5, 6, 7},
- {0, 0, 8, 9},
- {0, 0, 0, 10},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {320, 340, 360},
- {588, 624, 660},
- {598, 632, 666},
- {380, 400, 420},
- },
- },
- {
- s: blas.Left,
- ul: blas.Upper,
- tA: blas.NoTrans,
- d: blas.NonUnit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2},
- {0, 5},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {72, 78, 84},
- {130, 140, 150},
- },
- },
- {
- s: blas.Left,
- ul: blas.Upper,
- tA: blas.NoTrans,
- d: blas.Unit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3, 4},
- {0, 5, 6, 7},
- {0, 0, 8, 9},
- {0, 0, 0, 10},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {320, 340, 360},
- {484, 512, 540},
- {374, 394, 414},
- {38, 40, 42},
- },
- },
- {
- s: blas.Left,
- ul: blas.Upper,
- tA: blas.NoTrans,
- d: blas.Unit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2},
- {0, 5},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {72, 78, 84},
- {26, 28, 30},
- },
- },
- {
- s: blas.Left,
- ul: blas.Lower,
- tA: blas.NoTrans,
- d: blas.NonUnit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0, 0},
- {2, 5, 0, 0},
- {3, 6, 8, 0},
- {4, 7, 9, 10},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {20, 22, 24},
- {170, 184, 198},
- {472, 506, 540},
- {930, 990, 1050},
- },
- },
- {
- s: blas.Left,
- ul: blas.Lower,
- tA: blas.NoTrans,
- d: blas.NonUnit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0},
- {2, 5},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {20, 22, 24},
- {170, 184, 198},
- },
- },
- {
- s: blas.Left,
- ul: blas.Lower,
- tA: blas.NoTrans,
- d: blas.Unit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0, 0},
- {2, 5, 0, 0},
- {3, 6, 8, 0},
- {4, 7, 9, 10},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {20, 22, 24},
- {66, 72, 78},
- {248, 268, 288},
- {588, 630, 672},
- },
- },
- {
- s: blas.Left,
- ul: blas.Lower,
- tA: blas.NoTrans,
- d: blas.Unit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0},
- {2, 5},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {20, 22, 24},
- {66, 72, 78},
- },
- },
- {
- s: blas.Left,
- ul: blas.Upper,
- tA: blas.Trans,
- d: blas.NonUnit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3, 4},
- {0, 5, 6, 7},
- {0, 0, 8, 9},
- {0, 0, 0, 10},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {20, 22, 24},
- {170, 184, 198},
- {472, 506, 540},
- {930, 990, 1050},
- },
- },
- {
- s: blas.Left,
- ul: blas.Upper,
- tA: blas.Trans,
- d: blas.NonUnit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2},
- {0, 5},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {20, 22, 24},
- {170, 184, 198},
- },
- },
- {
- s: blas.Left,
- ul: blas.Upper,
- tA: blas.Trans,
- d: blas.Unit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3, 4},
- {0, 5, 6, 7},
- {0, 0, 8, 9},
- {0, 0, 0, 10},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {20, 22, 24},
- {66, 72, 78},
- {248, 268, 288},
- {588, 630, 672},
- },
- },
- {
- s: blas.Left,
- ul: blas.Upper,
- tA: blas.Trans,
- d: blas.Unit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2},
- {0, 5},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {20, 22, 24},
- {66, 72, 78},
- },
- },
- {
- s: blas.Left,
- ul: blas.Lower,
- tA: blas.Trans,
- d: blas.NonUnit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0, 0},
- {2, 5, 0, 0},
- {3, 6, 8, 0},
- {4, 7, 9, 10},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {320, 340, 360},
- {588, 624, 660},
- {598, 632, 666},
- {380, 400, 420},
- },
- },
- {
- s: blas.Left,
- ul: blas.Lower,
- tA: blas.Trans,
- d: blas.NonUnit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0},
- {2, 5},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {72, 78, 84},
- {130, 140, 150},
- },
- },
- {
- s: blas.Left,
- ul: blas.Lower,
- tA: blas.Trans,
- d: blas.Unit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0, 0},
- {2, 5, 0, 0},
- {3, 6, 8, 0},
- {4, 7, 9, 10},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {320, 340, 360},
- {484, 512, 540},
- {374, 394, 414},
- {38, 40, 42},
- },
- },
- {
- s: blas.Left,
- ul: blas.Lower,
- tA: blas.Trans,
- d: blas.Unit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0},
- {2, 5},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {72, 78, 84},
- {26, 28, 30},
- },
- },
- {
- s: blas.Right,
- ul: blas.Upper,
- tA: blas.NoTrans,
- d: blas.NonUnit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3},
- {0, 4, 5},
- {0, 0, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {20, 128, 314},
- {26, 164, 398},
- {32, 200, 482},
- {38, 236, 566},
- },
- },
- {
- s: blas.Right,
- ul: blas.Upper,
- tA: blas.NoTrans,
- d: blas.NonUnit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3},
- {0, 4, 5},
- {0, 0, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {20, 128, 314},
- {26, 164, 398},
- },
- },
- {
- s: blas.Right,
- ul: blas.Upper,
- tA: blas.NoTrans,
- d: blas.Unit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3},
- {0, 4, 5},
- {0, 0, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {20, 62, 194},
- {26, 80, 248},
- {32, 98, 302},
- {38, 116, 356},
- },
- },
- {
- s: blas.Right,
- ul: blas.Upper,
- tA: blas.NoTrans,
- d: blas.Unit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3},
- {0, 4, 5},
- {0, 0, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {20, 62, 194},
- {26, 80, 248},
- },
- },
- {
- s: blas.Right,
- ul: blas.Lower,
- tA: blas.NoTrans,
- d: blas.NonUnit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0},
- {2, 4, 0},
- {3, 5, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {136, 208, 144},
- {172, 262, 180},
- {208, 316, 216},
- {244, 370, 252},
- },
- },
- {
- s: blas.Right,
- ul: blas.Lower,
- tA: blas.NoTrans,
- d: blas.NonUnit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0},
- {2, 4, 0},
- {3, 5, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {136, 208, 144},
- {172, 262, 180},
- },
- },
- {
- s: blas.Right,
- ul: blas.Lower,
- tA: blas.NoTrans,
- d: blas.Unit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0},
- {2, 4, 0},
- {3, 5, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {136, 142, 24},
- {172, 178, 30},
- {208, 214, 36},
- {244, 250, 42},
- },
- },
- {
- s: blas.Right,
- ul: blas.Lower,
- tA: blas.NoTrans,
- d: blas.Unit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0},
- {2, 4, 0},
- {3, 5, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {136, 142, 24},
- {172, 178, 30},
- },
- },
-
- {
- s: blas.Right,
- ul: blas.Upper,
- tA: blas.Trans,
- d: blas.NonUnit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3},
- {0, 4, 5},
- {0, 0, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {136, 208, 144},
- {172, 262, 180},
- {208, 316, 216},
- {244, 370, 252},
- },
- },
- {
- s: blas.Right,
- ul: blas.Upper,
- tA: blas.Trans,
- d: blas.NonUnit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3},
- {0, 4, 5},
- {0, 0, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {136, 208, 144},
- {172, 262, 180},
- },
- },
- {
- s: blas.Right,
- ul: blas.Upper,
- tA: blas.Trans,
- d: blas.Unit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3},
- {0, 4, 5},
- {0, 0, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {136, 142, 24},
- {172, 178, 30},
- {208, 214, 36},
- {244, 250, 42},
- },
- },
- {
- s: blas.Right,
- ul: blas.Upper,
- tA: blas.Trans,
- d: blas.Unit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 2, 3},
- {0, 4, 5},
- {0, 0, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {136, 142, 24},
- {172, 178, 30},
- },
- },
-
- {
- s: blas.Right,
- ul: blas.Lower,
- tA: blas.Trans,
- d: blas.NonUnit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0},
- {2, 4, 0},
- {3, 5, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {20, 128, 314},
- {26, 164, 398},
- {32, 200, 482},
- {38, 236, 566},
- },
- },
- {
- s: blas.Right,
- ul: blas.Lower,
- tA: blas.Trans,
- d: blas.NonUnit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0},
- {2, 4, 0},
- {3, 5, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {20, 128, 314},
- {26, 164, 398},
- },
- },
- {
- s: blas.Right,
- ul: blas.Lower,
- tA: blas.Trans,
- d: blas.Unit,
- m: 4,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0},
- {2, 4, 0},
- {3, 5, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- {16, 17, 18},
- {19, 20, 21},
- },
- ans: [][]float64{
- {20, 62, 194},
- {26, 80, 248},
- {32, 98, 302},
- {38, 116, 356},
- },
- },
- {
- s: blas.Right,
- ul: blas.Lower,
- tA: blas.Trans,
- d: blas.Unit,
- m: 2,
- n: 3,
- alpha: 2,
- a: [][]float64{
- {1, 0, 0},
- {2, 4, 0},
- {3, 5, 6},
- },
- b: [][]float64{
- {10, 11, 12},
- {13, 14, 15},
- },
- ans: [][]float64{
- {20, 62, 194},
- {26, 80, 248},
- },
- },
- } {
- aFlat := flatten(test.a)
- bFlat := flatten(test.b)
- ansFlat := flatten(test.ans)
- blasser.Dtrmm(test.s, test.ul, test.tA, test.d, test.m, test.n, test.alpha, aFlat, len(test.a[0]), bFlat, len(test.b[0]))
- if !floats.EqualApprox(ansFlat, bFlat, 1e-14) {
- t.Errorf("Case %v. Want %v, got %v.", i, ansFlat, bFlat)
- }
- }
-}