blas/testblas: Added 7x7 dgemv test case.
MODIFY
blas/testblas/dgemv.go
blas/testblas/dgemv.go
@@ -567,6 +567,52 @@ var DgemvCases = []DgemvCase{
},
},
},
+ {
+ Name: "M_eq_N_Lg_IncNot1_Trans",
+ tA: blas.Trans,
+ m: 7,
+ n: 7,
+ A: [][]float64{
+ {4.1, 6.2, 8.1, 2.5, 3.3, 7.4, 9.3},
+ {9.6, 3.5, 9.1, 1.2, 5.4, 4.8, 8.7},
+ {10, 7, 3, 2, 4, 1, 12},
+ {9.6, 3.5, 9.1, 1.2, 5.4, 4.8, 8.7},
+ {4.1, 6.2, 8.1, 2.5, 3.3, 7.4, 9.3},
+ {10, 7, 3, 2, 4, 1, 12},
+ {9.6, 3.5, 9.1, 1.2, 5.4, 4.8, 8.7},
+ },
+ incX: 2,
+ incY: 3,
+ x: []float64{1, 105, 2, 150, 3, 200, 4, 300, 5, 600, 6, 700, 7},
+ y: []float64{7, 200, 600, 8, -400, -500, 9, 600, 700, 10, 500, 400, 11, 200, 300, 12, 100, 200, 13, 300, 400, 14},
+
+ Subcases: []DgemvSubcase{
+ {
+ alpha: 8,
+ beta: -6,
+ ans: []float64{1873.1999999999998, 200, 600, 1117.6, -400, -500, 1497.1999999999998, 600, 700, 328.8, 500, 400, 942, 200, 300, 854.4000000000001, 100, 200, 2137.2, 300, 400, 14},
+ },
+ {
+ mulXNeg1: true,
+ alpha: 8,
+ beta: -6,
+ ans: []float64{1690.8, 200, 600, 1148, -400, -500, 1562.8, 600, 700, 357.6, 500, 400, 897.2, 200, 300, 998.4, 100, 200, 2103.6000000000004, 300, 400, 14},
+ },
+ {
+ mulYNeg1: true,
+ alpha: 8,
+ beta: -6,
+ ans: []float64{2173.2, 200, 600, 878.4000000000001, -400, -500, 954, 600, 700, 328.8, 500, 400, 1485.1999999999998, 200, 300, 1093.6, 100, 200, 1837.1999999999998, 300, 400, 14},
+ },
+ {
+ mulXNeg1: true,
+ mulYNeg1: true,
+ alpha: 8,
+ beta: -6,
+ ans: []float64{2139.6, 200, 600, 1022.4, -400, -500, 909.2, 600, 700, 357.6, 500, 400, 1550.8, 200, 300, 1124, 100, 200, 1654.8, 300, 400, 14},
+ },
+ },
+ },
// TODO: A can be longer than mxn. Add cases where it is longer
// TODO: x and y can also be longer. Add tests for these

blas/testblas: Added 7x7 dgemv test case.