C 代码 构成 3 级基本线性代数子程序 （BLAS） 对于矩阵-矩阵运算， 使用双精度实数算术.rar

19 January 2020 08:48:26 AM BLAS3_D_TEST C version Test the BLAS3_D library. DGEMM_TEST DGEMM carries out matrix multiplications for double precision real matrices. 1: C = alpha * A * B + beta * C; 2: C = alpha * A' * B + beta * C; 3: C = alpha * A * B' + beta * C; 4: C = alpha * A' * B' + beta * C; We carry out all four calculations, but in each case, we choose our input matrices so that we get the same result. C = alpha * A * B + beta * C: Col: 0 1 2 3 4 Row 0: 1585 1660 1735 1810 1885 1: 2875 3010 3145 3280 3415 2: 4165 4360 4555 4750 4945 3: 5455 5710 5965 6220 6475 C = alpha * A' * B + beta * C: Col: 0 1 2 3 4 Row 0: 1585 1660 1735 1810 1885 1: 2875 3010 3145 3280 3415 2: 4165 4360 4555 4750 4945 3: 5455 5710 5965 6220 6475 C = alpha * A * B' + beta * C: Col: 0 1 2 3 4 Row 0: 1585 1660 1735 1810 1885 1: 2875 3010 3145 3280 3415 2: 4165 4360 4555 4750 4945 3: 5455 5710 5965 6220 6475 C = alpha * A' * B' + beta * C: Col: 0 1 2 3 4 Row 0: 1585 1660 1735 1810 1885 1: 2875 3010 3145 3280 3415 2: 4165 4360 4555 4750 4945 3: 5455 5710 5965 6220 6475 DTRMM_TEST DTRMM multiplies a triangular matrix A and a rectangular matrix B 1: B = alpha * A * B; 2: B = alpha * A' * B; B = alpha * A * B: Col: 0 1 2 3 4 Row 0: 828 856 884 912 940 1: 970 1000 1030 1060 1090 2: 946 972 998 1024 1050 3: 656 672 688 704 720 B = alpha * A' * B: Col: 0 1 2 3 4 Row 0: 44 48 52 56 60 1: 234 248 262 276 290 2: 670 700 730 760 790 3: 1452 1504 1556 1608 1660 DTRSM_TEST DTRSM solves a linear system involving a triangular matrix A and a rectangular matrix B. 1: Solve A * X = alpha * B; 2: Solve A' * X = alpha * B; 3: Solve X * A = alpha * B; 4: Solve X * A' = alpha * B; X = inv ( A ) * alpha * B: Col: 0 1 2 3 4 Row 0: -7.10938 -6.92708 -6.74479 -6.5625 -6.38021 1: -2.84375 -2.77083 -2.69792 -2.625 -2.55208 2: -1.625 -1.58333 -1.54167 -1.5 -1.45833 3: 10.25 10.5 10.75 11 11.25 X = inv ( A' ) * alpha * B: Col: 0 1 2 3 4 Row 0: 11 12 13 14 15 1: 2.25 2 1.75 1.5 1.25 2: 1.125 1 0.875 0.75 0.625 3: 0.703125 0.625 0.546875 0.46875 0.390625 X = alpha * B * inv ( A ): Col: 0 1 2 3 4 Row 0: 11 -2.25 -1.125 -0.703125 -0.492188 1: 21 -4.75 -2.375 -1.48438 -1.03906 2: 31 -7.25 -3.625 -2.26562 -1.58594 3: 41 -9.75 -4.875 -3.04688 -2.13281 X = alpha * B * inv ( A' ): Col: 0 1 2 3 4 Row 0: 0.820312 0.328125 0.1875 0.125 3 1: 2.46094 0.984375 0.5625 0.375 5 2: 4.10156 1.64062 0.9375 0.625 7 3: 5.74219 2.29688 1.3125 0.875 9 BLAS3_D_TEST Normal end of execution. 19 January 2020 08:48:26 AM