6 "gonum.org/v1/gonum/blas"
7 "gonum.org/v1/gonum/floats"
10 type Dtrmmer interface {
11 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)
14 func DtrmmTest(t *testing.T, blasser Dtrmmer) {
15 for i, test := range []struct {
798 aFlat := flatten(test.a)
799 bFlat := flatten(test.b)
800 ansFlat := flatten(test.ans)
801 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]))
802 if !floats.EqualApprox(ansFlat, bFlat, 1e-14) {
803 t.Errorf("Case %v. Want %v, got %v.", i, ansFlat, bFlat)