share/lapack/rtest_dgemm.mac

   1 (if not ?fboundp('dgemm) then load("lapack"),
   2   'done);
   3 'done$
   4
   5 (float_approx_equal_tolerance: 1e-14, 'done);
   6 'done$
   7
   8 (A : matrix([1,2,3],[4,5,6],[7,8,9]),
   9   B : matrix([-1,-2,-3],[-4,-5,-6],[-7,-8,-9]),
  10   C : matrix([3,2,1],[6,5,4],[9,8,7]),
  11   'done);
  12 'done$
  13
  14 dgemm(A, B);
  15 matrix([-30.0, -36.0, -42.0], [-66.0, -81.0, -96.0], [-102.0, -126.0, -150.0]);
  16
  17 dgemm(A, B, transpose_a = true);
  18 matrix([-66.0, -78.0, -90.0], [-78.0, -93.0, -108.0], [-90.0, -108.0, -126.0]);
  19
  20 dgemm(A, B, transpose_b = true);
  21 matrix([-14.0,-32.0,-50.0],[-32.0,-77.0,-122.0],[-50.0,-122.0,-194.0]);
  22
  23 dgemm(A, B, c=C, beta=2);
  24 matrix([-24.0, -32.0, -40.0], [-54.0, -71.0, -88.0], [-84.0, -110.0, -136.0]);
  25
  26 dgemm(A, B, c=C, beta=2, alpha=-1);
  27 matrix([36.0,40.0,44.0],[78.0,91.0,104.0],[120.0,142.0,164.0]);
  28
  29 dgemm(A, B, c=C, beta=2, alpha=-1, transpose_a = true);
  30 matrix([72.0,82.0,92.0],[90.0,103.0,116.0],[108.0,124.0,140.0]);