@Test
public void testGemmDenseDenseDense() {

    // C := alpha * A * B + beta * C
    double alpha = 1.0;
    DenseMatrix a = new DenseMatrix(new double[][] {{10.0, 11.0}, {0.0, 1.0}});
    DenseMatrix b = new DenseMatrix(new double[][] {{1.0, 0.3}, {0.0, 1.0}});
    double beta = 0.0;
    DenseMatrix c = new DenseMatrix(new double[][] {{1.0, 2.0}, {2.0, 3.0}});
    DenseMatrix exp = (DenseMatrix)a.times(b);//.times(alpha).plus(c.times(beta));
    Blas.gemm(alpha, a, b, beta, c);
    Assert.assertEquals(exp, c);
}