htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_gemm.c

   1 #include <gsl/gsl_test.h>
   2 #include <gsl/gsl_ieee_utils.h>
   3 #include <gsl/gsl_math.h>
   4 #include <gsl/gsl_cblas.h>
   5
   6 #include "tests.h"
   7
   8 void
   9 test_gemm (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int order = 101;
  13    int transA = 111;
  14    int transB = 111;
  15    int M = 1;
  16    int N = 2;
  17    int K = 4;
  18    float alpha = 1.0f;
  19    float beta = 0.0f;
  20    float A[] = { 0.199f, 0.237f, 0.456f, 0.377f };
  21    int lda = 4;
  22    float B[] = { 0.842f, -0.734f, 0.323f, -0.957f, -0.303f, -0.873f, -0.871f, -0.819f };
  23    int ldb = 2;
  24    float C[] = { 0.498f, -0.925f };
  25    int ldc = 2;
  26    float C_expected[] = { -0.222426f, -1.07973f };
  27    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  28    {
  29      int i;
  30      for (i = 0; i < 2; i++) {
  31        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1466)");
  32      }
  33    };
  34   };
  35
  36
  37   {
  38    int order = 102;
  39    int transA = 111;
  40    int transB = 111;
  41    int M = 1;
  42    int N = 2;
  43    int K = 4;
  44    float alpha = 1.0f;
  45    float beta = 0.0f;
  46    float A[] = { -0.83f, 0.922f, -0.228f, -0.003f };
  47    int lda = 1;
  48    float B[] = { 0.072f, 0.345f, 0.944f, -0.39f, -0.577f, 0.656f, -0.693f, -0.453f };
  49    int ldb = 4;
  50    float C[] = { 0.583f, 0.522f };
  51    int ldc = 1;
  52    float C_expected[] = { 0.044268f, 1.24311f };
  53    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  54    {
  55      int i;
  56      for (i = 0; i < 2; i++) {
  57        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1467)");
  58      }
  59    };
  60   };
  61
  62
  63   {
  64    int order = 101;
  65    int transA = 111;
  66    int transB = 112;
  67    int M = 1;
  68    int N = 2;
  69    int K = 4;
  70    float alpha = 0.1f;
  71    float beta = 0.1f;
  72    float A[] = { -0.838f, 0.622f, -0.494f, 0.304f };
  73    int lda = 4;
  74    float B[] = { 0.147f, 0.134f, 0.169f, 0.734f, -0.7f, 0.541f, -0.794f, -0.256f };
  75    int ldb = 4;
  76    float C[] = { -0.632f, -0.559f };
  77    int ldc = 2;
  78    float C_expected[] = { -0.0532188f, 0.0678514f };
  79    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  80    {
  81      int i;
  82      for (i = 0; i < 2; i++) {
  83        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1468)");
  84      }
  85    };
  86   };
  87
  88
  89   {
  90    int order = 102;
  91    int transA = 111;
  92    int transB = 112;
  93    int M = 1;
  94    int N = 2;
  95    int K = 4;
  96    float alpha = 0.1f;
  97    float beta = 0.1f;
  98    float A[] = { -0.937f, 0.635f, 0.596f, -0.51f };
  99    int lda = 1;
 100    float B[] = { -0.688f, -0.265f, 0.049f, 0.133f, -0.918f, -0.147f, 0.977f, -0.21f };
 101    int ldb = 2;
 102    float C[] = { 0.844f, 0.999f };
 103    int ldc = 1;
 104    float C_expected[] = { 0.0474373f, 0.135125f };
 105    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 106    {
 107      int i;
 108      for (i = 0; i < 2; i++) {
 109        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1469)");
 110      }
 111    };
 112   };
 113
 114
 115   {
 116    int order = 101;
 117    int transA = 112;
 118    int transB = 111;
 119    int M = 1;
 120    int N = 2;
 121    int K = 4;
 122    float alpha = -0.3f;
 123    float beta = 0.1f;
 124    float A[] = { -0.165f, 0.638f, 0.346f, -0.697f };
 125    int lda = 1;
 126    float B[] = { 0.499f, -0.73f, 0.262f, 0.759f, 0.664f, 0.997f, -0.702f, -0.839f };
 127    int ldb = 2;
 128    float C[] = { 0.17f, 0.425f };
 129    int ldc = 2;
 130    float C_expected[] = { -0.224158f, -0.417831f };
 131    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 132    {
 133      int i;
 134      for (i = 0; i < 2; i++) {
 135        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1470)");
 136      }
 137    };
 138   };
 139
 140
 141   {
 142    int order = 102;
 143    int transA = 112;
 144    int transB = 111;
 145    int M = 1;
 146    int N = 2;
 147    int K = 4;
 148    float alpha = -0.3f;
 149    float beta = 0.1f;
 150    float A[] = { -0.603f, -0.714f, -0.893f, 0.046f };
 151    int lda = 4;
 152    float B[] = { 0.859f, -0.694f, -0.868f, -0.98f, -0.103f, 0.567f, -0.277f, -0.734f };
 153    int ldb = 4;
 154    float C[] = { 0.517f, -0.622f };
 155    int ldc = 1;
 156    float C_expected[] = { -0.160575f, -0.0234604f };
 157    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 158    {
 159      int i;
 160      for (i = 0; i < 2; i++) {
 161        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1471)");
 162      }
 163    };
 164   };
 165
 166
 167   {
 168    int order = 101;
 169    int transA = 112;
 170    int transB = 112;
 171    int M = 1;
 172    int N = 2;
 173    int K = 4;
 174    float alpha = 0.1f;
 175    float beta = 1.0f;
 176    float A[] = { -0.087f, -0.047f, -0.051f, -0.615f };
 177    int lda = 1;
 178    float B[] = { -0.722f, -0.077f, 0.563f, 0.501f, 0.855f, 0.605f, 0.556f, -0.627f };
 179    int ldb = 4;
 180    float C[] = { -0.181f, -0.89f };
 181    int ldc = 2;
 182    float C_expected[] = { -0.208039f, -0.864557f };
 183    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 184    {
 185      int i;
 186      for (i = 0; i < 2; i++) {
 187        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1472)");
 188      }
 189    };
 190   };
 191
 192
 193   {
 194    int order = 102;
 195    int transA = 112;
 196    int transB = 112;
 197    int M = 1;
 198    int N = 2;
 199    int K = 4;
 200    float alpha = 0.1f;
 201    float beta = 1.0f;
 202    float A[] = { -0.753f, -0.074f, -0.247f, -0.19f };
 203    int lda = 4;
 204    float B[] = { 0.061f, 0.743f, 0.22f, -0.682f, 0.733f, 0.417f, 0.772f, 0.665f };
 205    int ldb = 2;
 206    float C[] = { -0.253f, 0.972f };
 207    int ldc = 1;
 208    float C_expected[] = { -0.291994f, 0.898164f };
 209    cblas_sgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 210    {
 211      int i;
 212      for (i = 0; i < 2; i++) {
 213        gsl_test_rel(C[i], C_expected[i], flteps, "sgemm(case 1473)");
 214      }
 215    };
 216   };
 217
 218
 219   {
 220    int order = 101;
 221    int transA = 111;
 222    int transB = 111;
 223    int M = 1;
 224    int N = 2;
 225    int K = 4;
 226    double alpha = 0;
 227    double beta = 0;
 228    double A[] = { 0.017, 0.191, 0.863, -0.97 };
 229    int lda = 4;
 230    double B[] = { -0.207, -0.916, -0.278, 0.403, 0.885, 0.409, -0.772, -0.27 };
 231    int ldb = 2;
 232    double C[] = { -0.274, -0.858 };
 233    int ldc = 2;
 234    double C_expected[] = { 0.0, 0.0 };
 235    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 236    {
 237      int i;
 238      for (i = 0; i < 2; i++) {
 239        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1474)");
 240      }
 241    };
 242   };
 243
 244
 245   {
 246    int order = 102;
 247    int transA = 111;
 248    int transB = 111;
 249    int M = 1;
 250    int N = 2;
 251    int K = 4;
 252    double alpha = 0;
 253    double beta = 0;
 254    double A[] = { 0.571, 0.081, 0.109, 0.988 };
 255    int lda = 1;
 256    double B[] = { -0.048, -0.753, -0.8, -0.89, -0.535, -0.017, -0.018, -0.544 };
 257    int ldb = 4;
 258    double C[] = { -0.876, -0.792 };
 259    int ldc = 1;
 260    double C_expected[] = { 0.0, 0.0 };
 261    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 262    {
 263      int i;
 264      for (i = 0; i < 2; i++) {
 265        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1475)");
 266      }
 267    };
 268   };
 269
 270
 271   {
 272    int order = 101;
 273    int transA = 111;
 274    int transB = 112;
 275    int M = 1;
 276    int N = 2;
 277    int K = 4;
 278    double alpha = -0.3;
 279    double beta = 1;
 280    double A[] = { 0.939, 0.705, 0.977, 0.4 };
 281    int lda = 4;
 282    double B[] = { -0.089, -0.822, 0.937, 0.159, 0.789, -0.413, -0.172, 0.88 };
 283    int ldb = 4;
 284    double C[] = { -0.619, 0.063 };
 285    int ldc = 2;
 286    double C_expected[] = { -0.7137904, -0.1270986 };
 287    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 288    {
 289      int i;
 290      for (i = 0; i < 2; i++) {
 291        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1476)");
 292      }
 293    };
 294   };
 295
 296
 297   {
 298    int order = 102;
 299    int transA = 111;
 300    int transB = 112;
 301    int M = 1;
 302    int N = 2;
 303    int K = 4;
 304    double alpha = -0.3;
 305    double beta = 1;
 306    double A[] = { -0.795, 0.81, 0.388, 0.09 };
 307    int lda = 1;
 308    double B[] = { -0.847, 0.031, -0.938, 0.09, -0.286, -0.478, -0.981, 0.881 };
 309    int ldb = 2;
 310    double C[] = { -0.242, -0.02 };
 311    int ldc = 1;
 312    double C_expected[] = { -0.1562981, -0.0026243 };
 313    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 314    {
 315      int i;
 316      for (i = 0; i < 2; i++) {
 317        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1477)");
 318      }
 319    };
 320   };
 321
 322
 323   {
 324    int order = 101;
 325    int transA = 112;
 326    int transB = 111;
 327    int M = 1;
 328    int N = 2;
 329    int K = 4;
 330    double alpha = -1;
 331    double beta = 0;
 332    double A[] = { -0.556, 0.532, 0.746, 0.673 };
 333    int lda = 1;
 334    double B[] = { -0.525, 0.967, 0.687, -0.024, 0.527, 0.485, 0.109, -0.46 };
 335    int ldb = 2;
 336    double C[] = { -0.495, 0.859 };
 337    int ldc = 2;
 338    double C_expected[] = { -1.123883, 0.49819 };
 339    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 340    {
 341      int i;
 342      for (i = 0; i < 2; i++) {
 343        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1478)");
 344      }
 345    };
 346   };
 347
 348
 349   {
 350    int order = 102;
 351    int transA = 112;
 352    int transB = 111;
 353    int M = 1;
 354    int N = 2;
 355    int K = 4;
 356    double alpha = -1;
 357    double beta = 0;
 358    double A[] = { -0.358, 0.224, -0.941, 0.513 };
 359    int lda = 4;
 360    double B[] = { -0.201, -0.159, -0.586, -0.016, -0.324, 0.411, 0.115, -0.229 };
 361    int ldb = 4;
 362    double C[] = { 0.558, 0.596 };
 363    int ldc = 1;
 364    double C_expected[] = { -0.57956, 0.017636 };
 365    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 366    {
 367      int i;
 368      for (i = 0; i < 2; i++) {
 369        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1479)");
 370      }
 371    };
 372   };
 373
 374
 375   {
 376    int order = 101;
 377    int transA = 112;
 378    int transB = 112;
 379    int M = 1;
 380    int N = 2;
 381    int K = 4;
 382    double alpha = -0.3;
 383    double beta = 1;
 384    double A[] = { -0.586, 0.809, 0.709, -0.524 };
 385    int lda = 1;
 386    double B[] = { 0.768, 0.7, 0.619, -0.478, -0.129, -0.778, -0.432, 0.454 };
 387    int ldb = 4;
 388    double C[] = { 0.042, 0.252 };
 389    int ldc = 2;
 390    double C_expected[] = { -0.1996785, 0.5813976 };
 391    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 392    {
 393      int i;
 394      for (i = 0; i < 2; i++) {
 395        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1480)");
 396      }
 397    };
 398   };
 399
 400
 401   {
 402    int order = 102;
 403    int transA = 112;
 404    int transB = 112;
 405    int M = 1;
 406    int N = 2;
 407    int K = 4;
 408    double alpha = -0.3;
 409    double beta = 1;
 410    double A[] = { -0.164, 0.522, 0.948, -0.624 };
 411    int lda = 4;
 412    double B[] = { -0.142, 0.778, 0.359, 0.622, -0.637, -0.757, -0.282, -0.805 };
 413    int ldb = 2;
 414    double C[] = { -0.09, 0.183 };
 415    int ldc = 1;
 416    double C_expected[] = { -0.0248334, 0.1884672 };
 417    cblas_dgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 418    {
 419      int i;
 420      for (i = 0; i < 2; i++) {
 421        gsl_test_rel(C[i], C_expected[i], dbleps, "dgemm(case 1481)");
 422      }
 423    };
 424   };
 425
 426
 427   {
 428    int order = 101;
 429    int transA = 111;
 430    int transB = 111;
 431    int M = 1;
 432    int N = 2;
 433    int K = 4;
 434    float alpha[2] = {0.0f, 1.0f};
 435    float beta[2] = {0.0f, 0.0f};
 436    float A[] = { -0.082f, -0.281f, -0.096f, 0.913f, 0.974f, -0.706f, -0.773f, 0.522f };
 437    int lda = 4;
 438    float B[] = { 0.745f, -0.664f, 0.352f, -0.733f, 0.304f, -0.555f, -0.493f, -0.089f, 0.188f, 0.631f, 0.235f, 0.152f, -0.299f, -0.731f, -0.686f, -0.332f };
 439    int ldb = 2;
 440    float C[] = { -0.179f, -0.284f, -0.996f, -0.414f };
 441    int ldc = 2;
 442    float C_expected[] = { -1.06679f, 1.47116f, 0.599689f, 0.933532f };
 443    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 444    {
 445      int i;
 446      for (i = 0; i < 2; i++) {
 447        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1482) real");
 448        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1482) imag");
 449      };
 450    };
 451   };
 452
 453
 454   {
 455    int order = 102;
 456    int transA = 111;
 457    int transB = 111;
 458    int M = 1;
 459    int N = 2;
 460    int K = 4;
 461    float alpha[2] = {0.0f, 1.0f};
 462    float beta[2] = {0.0f, 0.0f};
 463    float A[] = { 0.044f, -0.33f, 0.279f, 0.712f, -0.363f, -0.788f, -0.768f, -0.551f };
 464    int lda = 1;
 465    float B[] = { 0.138f, 0.927f, -0.178f, -0.864f, 0.888f, 0.844f, -0.199f, 0.706f, -0.034f, 0.483f, 0.499f, 0.664f, 0.648f, 0.324f, 0.97f, 0.609f };
 466    int ldb = 4;
 467    float C[] = { -0.129f, 0.842f, 0.214f, -0.626f };
 468    int ldc = 1;
 469    float C_expected[] = { 1.81122f, 1.76205f, 1.0574f, -0.564966f };
 470    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 471    {
 472      int i;
 473      for (i = 0; i < 2; i++) {
 474        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1483) real");
 475        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1483) imag");
 476      };
 477    };
 478   };
 479
 480
 481   {
 482    int order = 101;
 483    int transA = 111;
 484    int transB = 112;
 485    int M = 1;
 486    int N = 2;
 487    int K = 4;
 488    float alpha[2] = {0.0f, 0.0f};
 489    float beta[2] = {-1.0f, 0.0f};
 490    float A[] = { 0.812f, -0.471f, 0.241f, 0.795f, 0.439f, 0.131f, -0.636f, 0.531f };
 491    int lda = 4;
 492    float B[] = { 0.062f, 0.807f, 0.873f, 0.372f, 0.239f, 0.804f, 0.537f, -0.954f, -0.396f, 0.838f, 0.081f, 0.15f, 0.489f, -0.438f, 0.165f, 0.429f };
 493    int ldb = 4;
 494    float C[] = { 0.868f, 0.329f, -0.509f, 0.724f };
 495    int ldc = 2;
 496    float C_expected[] = { -0.868f, -0.329f, 0.509f, -0.724f };
 497    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 498    {
 499      int i;
 500      for (i = 0; i < 2; i++) {
 501        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1484) real");
 502        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1484) imag");
 503      };
 504    };
 505   };
 506
 507
 508   {
 509    int order = 102;
 510    int transA = 111;
 511    int transB = 112;
 512    int M = 1;
 513    int N = 2;
 514    int K = 4;
 515    float alpha[2] = {0.0f, 0.0f};
 516    float beta[2] = {-1.0f, 0.0f};
 517    float A[] = { 0.832f, 0.198f, 0.794f, -0.522f, -0.319f, 0.578f, 0.332f, 0.746f };
 518    int lda = 1;
 519    float B[] = { -0.361f, 0.187f, -0.163f, -0.781f, 0.536f, 0.888f, -0.969f, 0.899f, 0.961f, -0.583f, 0.753f, 0.29f, -0.997f, 0.729f, -0.352f, -0.2f };
 520    int ldb = 2;
 521    float C[] = { 0.864f, 0.735f, -0.074f, -0.228f };
 522    int ldc = 1;
 523    float C_expected[] = { -0.864f, -0.735f, 0.074f, 0.228f };
 524    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 525    {
 526      int i;
 527      for (i = 0; i < 2; i++) {
 528        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1485) real");
 529        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1485) imag");
 530      };
 531    };
 532   };
 533
 534
 535   {
 536    int order = 101;
 537    int transA = 111;
 538    int transB = 113;
 539    int M = 1;
 540    int N = 2;
 541    int K = 4;
 542    float alpha[2] = {0.0f, 1.0f};
 543    float beta[2] = {0.0f, 0.1f};
 544    float A[] = { 0.149f, 0.187f, 0.263f, -0.715f, -0.882f, -0.907f, 0.87f, -0.527f };
 545    int lda = 4;
 546    float B[] = { -0.915f, -0.249f, -0.986f, -0.799f, -0.136f, 0.712f, 0.964f, 0.799f, -0.569f, 0.686f, 0.603f, 0.758f, 0.161f, -0.698f, -0.263f, -0.256f };
 547    int ldb = 4;
 548    float C[] = { 0.622f, -0.824f, -0.482f, -0.161f };
 549    int ldc = 2;
 550    float C_expected[] = { -0.246901f, 0.083044f, 1.25556f, 0.009106f };
 551    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 552    {
 553      int i;
 554      for (i = 0; i < 2; i++) {
 555        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1486) real");
 556        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1486) imag");
 557      };
 558    };
 559   };
 560
 561
 562   {
 563    int order = 102;
 564    int transA = 111;
 565    int transB = 113;
 566    int M = 1;
 567    int N = 2;
 568    int K = 4;
 569    float alpha[2] = {0.0f, 1.0f};
 570    float beta[2] = {0.0f, 0.1f};
 571    float A[] = { 0.963f, -0.943f, -0.734f, -0.253f, 0.832f, 0.545f, -0.815f, -0.434f };
 572    int lda = 1;
 573    float B[] = { 0.23f, -0.211f, 0.906f, 0.232f, -0.339f, 0.597f, -0.919f, 0.793f, 0.535f, 0.526f, 0.119f, 0.053f, 0.751f, 0.044f, 0.752f, -0.469f };
 574    int ldb = 2;
 575    float C[] = { 0.483f, -0.266f, -0.224f, -0.692f };
 576    int ldc = 1;
 577    float C_expected[] = { -0.047537f, 0.667177f, 1.02025f, 0.823778f };
 578    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 579    {
 580      int i;
 581      for (i = 0; i < 2; i++) {
 582        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1487) real");
 583        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1487) imag");
 584      };
 585    };
 586   };
 587
 588
 589   {
 590    int order = 101;
 591    int transA = 112;
 592    int transB = 111;
 593    int M = 1;
 594    int N = 2;
 595    int K = 4;
 596    float alpha[2] = {-0.3f, 0.1f};
 597    float beta[2] = {-1.0f, 0.0f};
 598    float A[] = { -0.657f, -0.497f, -0.293f, -0.168f, -0.943f, -0.181f, 0.569f, 0.91f };
 599    int lda = 1;
 600    float B[] = { -0.047f, 0.796f, -0.913f, 0.998f, 0.365f, 0.467f, -0.627f, -0.523f, 0.885f, 0.234f, -0.494f, 0.071f, -0.361f, -0.154f, -0.055f, -0.32f };
 601    int ldb = 2;
 602    float C[] = { 0.956f, 0.268f, 0.152f, 0.717f };
 603    int ldc = 2;
 604    float C_expected[] = { -0.668685f, 0.134477f, -0.715786f, -0.478065f };
 605    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 606    {
 607      int i;
 608      for (i = 0; i < 2; i++) {
 609        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1488) real");
 610        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1488) imag");
 611      };
 612    };
 613   };
 614
 615
 616   {
 617    int order = 102;
 618    int transA = 112;
 619    int transB = 111;
 620    int M = 1;
 621    int N = 2;
 622    int K = 4;
 623    float alpha[2] = {-0.3f, 0.1f};
 624    float beta[2] = {-1.0f, 0.0f};
 625    float A[] = { 0.394f, -0.482f, 0.631f, -0.833f, 0.221f, 0.672f, 0.2f, 0.967f };
 626    int lda = 4;
 627    float B[] = { 0.708f, 0.695f, 0.111f, -0.912f, 0.376f, 0.606f, -0.997f, -0.741f, 0.349f, 0.543f, 0.372f, -0.563f, 0.129f, -0.295f, -0.672f, -0.95f };
 628    int ldb = 4;
 629    float C[] = { 0.436f, 0.752f, 0.074f, 0.209f };
 630    int ldc = 1;
 631    float C_expected[] = { -0.325083f, -0.301952f, -0.283022f, 0.339919f };
 632    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 633    {
 634      int i;
 635      for (i = 0; i < 2; i++) {
 636        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1489) real");
 637        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1489) imag");
 638      };
 639    };
 640   };
 641
 642
 643   {
 644    int order = 101;
 645    int transA = 112;
 646    int transB = 112;
 647    int M = 1;
 648    int N = 2;
 649    int K = 4;
 650    float alpha[2] = {1.0f, 0.0f};
 651    float beta[2] = {-0.3f, 0.1f};
 652    float A[] = { 0.827f, -0.862f, 0.373f, -0.265f, -0.9f, 0.892f, -0.319f, 0.151f };
 653    int lda = 1;
 654    float B[] = { 0.603f, 0.816f, -0.511f, 0.831f, -0.36f, -0.954f, -0.978f, 0.485f, 0.675f, 0.186f, 0.463f, 0.144f, 0.851f, -0.458f, 0.766f, -0.213f };
 655    int ldb = 4;
 656    float C[] = { -0.335f, 0.333f, -0.4f, 0.422f };
 657    int ldc = 2;
 658    float C_expected[] = { 2.7126f, 0.702111f, 0.437661f, 0.691294f };
 659    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 660    {
 661      int i;
 662      for (i = 0; i < 2; i++) {
 663        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1490) real");
 664        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1490) imag");
 665      };
 666    };
 667   };
 668
 669
 670   {
 671    int order = 102;
 672    int transA = 112;
 673    int transB = 112;
 674    int M = 1;
 675    int N = 2;
 676    int K = 4;
 677    float alpha[2] = {1.0f, 0.0f};
 678    float beta[2] = {-0.3f, 0.1f};
 679    float A[] = { 0.966f, 0.476f, -0.013f, -0.655f, 0.773f, -0.543f, -0.231f, -0.353f };
 680    int lda = 4;
 681    float B[] = { -0.684f, 0.144f, 0.018f, -0.77f, -0.688f, 0.909f, -0.094f, -0.938f, -0.757f, 0.574f, -0.479f, 0.473f, 0.0f, 0.064f, -0.168f, 0.858f };
 682    int ldb = 2;
 683    float C[] = { -0.912f, 0.54f, 0.756f, 0.024f };
 684    int ldc = 1;
 685    float C_expected[] = { -0.156236f, 0.839112f, -0.230206f, -0.106256f };
 686    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 687    {
 688      int i;
 689      for (i = 0; i < 2; i++) {
 690        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1491) real");
 691        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1491) imag");
 692      };
 693    };
 694   };
 695
 696
 697   {
 698    int order = 101;
 699    int transA = 112;
 700    int transB = 113;
 701    int M = 1;
 702    int N = 2;
 703    int K = 4;
 704    float alpha[2] = {0.0f, 0.0f};
 705    float beta[2] = {1.0f, 0.0f};
 706    float A[] = { 0.66f, -0.113f, -0.663f, -0.856f, 0.614f, -0.344f, -0.964f, -0.532f };
 707    int lda = 1;
 708    float B[] = { -0.606f, -0.965f, -0.279f, -0.312f, 0.63f, 0.967f, 0.041f, -0.557f, 0.663f, 0.619f, -0.134f, 0.261f, -0.388f, 0.525f, 0.222f, 0.538f };
 709    int ldb = 4;
 710    float C[] = { 0.114f, -0.376f, -0.851f, -0.682f };
 711    int ldc = 2;
 712    float C_expected[] = { 0.114f, -0.376f, -0.851f, -0.682f };
 713    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 714    {
 715      int i;
 716      for (i = 0; i < 2; i++) {
 717        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1492) real");
 718        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1492) imag");
 719      };
 720    };
 721   };
 722
 723
 724   {
 725    int order = 102;
 726    int transA = 112;
 727    int transB = 113;
 728    int M = 1;
 729    int N = 2;
 730    int K = 4;
 731    float alpha[2] = {0.0f, 0.0f};
 732    float beta[2] = {1.0f, 0.0f};
 733    float A[] = { 0.212f, -0.752f, 0.679f, 0.49f, -0.029f, -0.488f, 0.567f, 0.374f };
 734    int lda = 4;
 735    float B[] = { -0.914f, 0.734f, -0.845f, 0.059f, -0.297f, 0.152f, -0.417f, -0.669f, 0.831f, -0.544f, 0.022f, 0.102f, -0.379f, -0.357f, -0.394f, -0.588f };
 736    int ldb = 2;
 737    float C[] = { -0.584f, 0.373f, 0.235f, 0.521f };
 738    int ldc = 1;
 739    float C_expected[] = { -0.584f, 0.373f, 0.235f, 0.521f };
 740    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 741    {
 742      int i;
 743      for (i = 0; i < 2; i++) {
 744        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1493) real");
 745        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1493) imag");
 746      };
 747    };
 748   };
 749
 750
 751   {
 752    int order = 101;
 753    int transA = 113;
 754    int transB = 111;
 755    int M = 1;
 756    int N = 2;
 757    int K = 4;
 758    float alpha[2] = {0.0f, 0.1f};
 759    float beta[2] = {-1.0f, 0.0f};
 760    float A[] = { 0.135f, 0.128f, 0.909f, -0.963f, 0.299f, -0.944f, 0.944f, 0.942f };
 761    int lda = 1;
 762    float B[] = { 0.924f, -0.317f, -0.992f, -0.854f, -0.435f, 0.102f, 0.126f, 0.862f, 0.952f, 0.68f, 0.545f, 0.168f, 0.752f, 0.549f, 0.687f, -0.76f };
 763    int ldb = 2;
 764    float C[] = { -0.369f, -0.33f, 0.849f, -0.632f };
 765    int ldc = 2;
 766    float C_expected[] = { 0.326537f, 0.37603f, -0.86067f, 0.529817f };
 767    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 768    {
 769      int i;
 770      for (i = 0; i < 2; i++) {
 771        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1494) real");
 772        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1494) imag");
 773      };
 774    };
 775   };
 776
 777
 778   {
 779    int order = 102;
 780    int transA = 113;
 781    int transB = 111;
 782    int M = 1;
 783    int N = 2;
 784    int K = 4;
 785    float alpha[2] = {0.0f, 0.1f};
 786    float beta[2] = {-1.0f, 0.0f};
 787    float A[] = { 0.061f, -0.271f, -0.043f, -0.023f, 0.694f, 0.333f, 0.733f, -0.967f };
 788    int lda = 4;
 789    float B[] = { 0.088f, -0.607f, 0.589f, 0.375f, -0.897f, -0.954f, -0.216f, -0.195f, -0.865f, -0.511f, -0.219f, 0.535f, 0.976f, 0.582f, 0.464f, -0.041f };
 790    int ldb = 4;
 791    float C[] = { 0.533f, -0.63f, 0.405f, 0.667f };
 792    int ldc = 1;
 793    float C_expected[] = { -0.459906f, 0.552595f, -0.425391f, -0.533626f };
 794    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 795    {
 796      int i;
 797      for (i = 0; i < 2; i++) {
 798        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1495) real");
 799        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1495) imag");
 800      };
 801    };
 802   };
 803
 804
 805   {
 806    int order = 101;
 807    int transA = 113;
 808    int transB = 112;
 809    int M = 1;
 810    int N = 2;
 811    int K = 4;
 812    float alpha[2] = {0.0f, 0.0f};
 813    float beta[2] = {1.0f, 0.0f};
 814    float A[] = { -0.676f, -0.116f, 0.707f, -0.256f, -0.893f, -0.966f, 0.159f, -0.246f };
 815    int lda = 1;
 816    float B[] = { 0.059f, 0.281f, -0.93f, -0.263f, 0.583f, -0.11f, 0.639f, -0.96f, -0.878f, 0.984f, 0.058f, 0.977f, -0.567f, 0.561f, -0.048f, -0.798f };
 817    int ldb = 4;
 818    float C[] = { 0.362f, -0.808f, 0.428f, -0.112f };
 819    int ldc = 2;
 820    float C_expected[] = { 0.362f, -0.808f, 0.428f, -0.112f };
 821    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 822    {
 823      int i;
 824      for (i = 0; i < 2; i++) {
 825        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1496) real");
 826        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1496) imag");
 827      };
 828    };
 829   };
 830
 831
 832   {
 833    int order = 102;
 834    int transA = 113;
 835    int transB = 112;
 836    int M = 1;
 837    int N = 2;
 838    int K = 4;
 839    float alpha[2] = {0.0f, 0.0f};
 840    float beta[2] = {1.0f, 0.0f};
 841    float A[] = { -0.915f, 0.439f, 0.171f, -0.019f, 0.843f, 0.944f, -0.581f, 0.856f };
 842    int lda = 4;
 843    float B[] = { -0.284f, 0.207f, -0.27f, 0.832f, 0.894f, -0.626f, -0.305f, -0.006f, 0.562f, -0.744f, -0.533f, 0.126f, -0.375f, -0.333f, 0.275f, 0.748f };
 844    int ldb = 2;
 845    float C[] = { -0.763f, -0.829f, 0.708f, -0.613f };
 846    int ldc = 1;
 847    float C_expected[] = { -0.763f, -0.829f, 0.708f, -0.613f };
 848    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 849    {
 850      int i;
 851      for (i = 0; i < 2; i++) {
 852        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1497) real");
 853        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1497) imag");
 854      };
 855    };
 856   };
 857
 858
 859   {
 860    int order = 101;
 861    int transA = 113;
 862    int transB = 113;
 863    int M = 1;
 864    int N = 2;
 865    int K = 4;
 866    float alpha[2] = {0.0f, 0.1f};
 867    float beta[2] = {0.0f, 1.0f};
 868    float A[] = { 0.496f, -0.9f, 0.825f, -0.678f, 0.41f, -0.585f, -0.264f, 0.308f };
 869    int lda = 1;
 870    float B[] = { 0.907f, 0.972f, -0.724f, 0.745f, -0.601f, 0.589f, 0.759f, -0.521f, -0.161f, -0.321f, 0.341f, -0.981f, -0.378f, -0.671f, -0.314f, -0.878f };
 871    int ldb = 4;
 872    float C[] = { -0.293f, 0.07f, 0.087f, -0.542f };
 873    int ldc = 2;
 874    float C_expected[] = { 0.10357f, -0.163927f, 0.444626f, -0.0076744f };
 875    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 876    {
 877      int i;
 878      for (i = 0; i < 2; i++) {
 879        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1498) real");
 880        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1498) imag");
 881      };
 882    };
 883   };
 884
 885
 886   {
 887    int order = 102;
 888    int transA = 113;
 889    int transB = 113;
 890    int M = 1;
 891    int N = 2;
 892    int K = 4;
 893    float alpha[2] = {0.0f, 0.1f};
 894    float beta[2] = {0.0f, 1.0f};
 895    float A[] = { -0.225f, -0.629f, -0.939f, -0.836f, -0.841f, -0.794f, 0.836f, -0.65f };
 896    int lda = 4;
 897    float B[] = { 0.869f, -0.453f, 0.8f, -0.947f, 0.545f, 0.716f, -0.507f, -0.228f, 0.722f, 0.372f, 0.77f, 0.317f, -0.153f, -0.524f, -0.465f, -0.684f };
 898    int ldb = 2;
 899    float C[] = { -0.896f, 0.91f, -0.973f, -0.269f };
 900    int ldc = 1;
 901    float C_expected[] = { -1.18974f, -1.0134f, 0.189027f, -1.14494f };
 902    cblas_cgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 903    {
 904      int i;
 905      for (i = 0; i < 2; i++) {
 906        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "cgemm(case 1499) real");
 907        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "cgemm(case 1499) imag");
 908      };
 909    };
 910   };
 911
 912
 913   {
 914    int order = 101;
 915    int transA = 111;
 916    int transB = 111;
 917    int M = 1;
 918    int N = 2;
 919    int K = 4;
 920    double alpha[2] = {1, 0};
 921    double beta[2] = {-1, 0};
 922    double A[] = { -0.33, 0.457, 0.428, -0.19, 0.86, -0.53, 0.058, -0.942 };
 923    int lda = 4;
 924    double B[] = { 0.434, 0.653, -0.124, 0.191, -0.112, -0.84, -0.72, 0.075, -0.503, -0.109, 0.3, -0.898, 0.489, 0.384, 0.993, -0.804 };
 925    int ldb = 2;
 926    double C[] = { -0.792, -0.155, -0.608, -0.243 };
 927    int ldc = 2;
 928    double C_expected[] = { 0.042563, -0.465908, -0.649991, -1.621116 };
 929    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 930    {
 931      int i;
 932      for (i = 0; i < 2; i++) {
 933        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1500) real");
 934        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1500) imag");
 935      };
 936    };
 937   };
 938
 939
 940   {
 941    int order = 102;
 942    int transA = 111;
 943    int transB = 111;
 944    int M = 1;
 945    int N = 2;
 946    int K = 4;
 947    double alpha[2] = {1, 0};
 948    double beta[2] = {-1, 0};
 949    double A[] = { 0.726, -0.438, -0.23, -0.054, -0.019, 0.902, -0.883, -0.235 };
 950    int lda = 1;
 951    double B[] = { 0.159, -0.18, 0.386, -0.167, 0.971, -0.072, 0.87, -0.839, 0.474, 0.956, -0.235, 0.332, 0.826, -0.056, -0.941, 0.01 };
 952    int ldb = 4;
 953    double C[] = { -0.799, 0.973, -0.549, -0.177 };
 954    int ldc = 1;
 955    double C_expected[] = { -0.181084, 0.257841, 2.251901, 1.558195 };
 956    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 957    {
 958      int i;
 959      for (i = 0; i < 2; i++) {
 960        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1501) real");
 961        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1501) imag");
 962      };
 963    };
 964   };
 965
 966
 967   {
 968    int order = 101;
 969    int transA = 111;
 970    int transB = 112;
 971    int M = 1;
 972    int N = 2;
 973    int K = 4;
 974    double alpha[2] = {0, 0.1};
 975    double beta[2] = {1, 0};
 976    double A[] = { 0.109, 0.892, -0.723, 0.793, 0.109, -0.419, -0.534, 0.448 };
 977    int lda = 4;
 978    double B[] = { -0.875, -0.31, -0.027, 0.067, 0.274, -0.126, -0.548, 0.497, 0.681, 0.388, 0.909, 0.889, 0.982, -0.074, -0.788, 0.233 };
 979    int ldb = 4;
 980    double C[] = { 0.503, 0.067, 0.239, 0.876 };
 981    int ldc = 2;
 982    double C_expected[] = { 0.6553584, 0.0864583, 0.2559136, 0.7518389 };
 983    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 984    {
 985      int i;
 986      for (i = 0; i < 2; i++) {
 987        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1502) real");
 988        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1502) imag");
 989      };
 990    };
 991   };
 992
 993
 994   {
 995    int order = 102;
 996    int transA = 111;
 997    int transB = 112;
 998    int M = 1;
 999    int N = 2;
1000    int K = 4;
1001    double alpha[2] = {0, 0.1};
1002    double beta[2] = {1, 0};
1003    double A[] = { 0.334, 0.192, -0.992, -0.168, 0.154, -0.75, -0.797, -0.76 };
1004    int lda = 1;
1005    double B[] = { -0.82, 0.147, -0.237, 0.68, 0.317, 0.257, -0.406, -0.802, 0.058, 0.012, -0.832, 0.949, -0.263, -0.085, -0.064, 0.492 };
1006    int ldb = 2;
1007    double C[] = { 0.079, -0.602, -0.392, 0.316 };
1008    int ldc = 1;
1009    double C_expected[] = { 0.0980569, -0.6430449, -0.539207, 0.4226848 };
1010    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1011    {
1012      int i;
1013      for (i = 0; i < 2; i++) {
1014        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1503) real");
1015        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1503) imag");
1016      };
1017    };
1018   };
1019
1020
1021   {
1022    int order = 101;
1023    int transA = 111;
1024    int transB = 113;
1025    int M = 1;
1026    int N = 2;
1027    int K = 4;
1028    double alpha[2] = {0, 0};
1029    double beta[2] = {-1, 0};
1030    double A[] = { -0.305, -0.698, -0.072, -0.383, 0.364, -0.656, 0.819, 0.194 };
1031    int lda = 4;
1032    double B[] = { 0.682, 0.498, -0.389, 0.923, -0.853, -0.558, -0.722, -0.085, -0.27, 0.026, -0.107, -0.036, 0.644, -0.327, -0.894, 0.34 };
1033    int ldb = 4;
1034    double C[] = { 0.981, -0.336, -0.377, -0.41 };
1035    int ldc = 2;
1036    double C_expected[] = { -0.981, 0.336, 0.377, 0.41 };
1037    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1038    {
1039      int i;
1040      for (i = 0; i < 2; i++) {
1041        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1504) real");
1042        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1504) imag");
1043      };
1044    };
1045   };
1046
1047
1048   {
1049    int order = 102;
1050    int transA = 111;
1051    int transB = 113;
1052    int M = 1;
1053    int N = 2;
1054    int K = 4;
1055    double alpha[2] = {0, 0};
1056    double beta[2] = {-1, 0};
1057    double A[] = { -0.306, -0.709, -0.196, 0.285, 0.873, -0.802, 0.715, -0.179 };
1058    int lda = 1;
1059    double B[] = { 0.028, 0.109, 0.87, -0.446, 0.735, 0.731, 0.021, -0.186, 0.541, 0.97, -0.333, 0.002, -0.089, -0.01, 0.331, 0.851 };
1060    int ldb = 2;
1061    double C[] = { 0.902, -0.584, -0.695, -0.607 };
1062    int ldc = 1;
1063    double C_expected[] = { -0.902, 0.584, 0.695, 0.607 };
1064    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1065    {
1066      int i;
1067      for (i = 0; i < 2; i++) {
1068        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1505) real");
1069        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1505) imag");
1070      };
1071    };
1072   };
1073
1074
1075   {
1076    int order = 101;
1077    int transA = 112;
1078    int transB = 111;
1079    int M = 1;
1080    int N = 2;
1081    int K = 4;
1082    double alpha[2] = {-1, 0};
1083    double beta[2] = {1, 0};
1084    double A[] = { 0.517, -0.136, 0.72, -0.237, 0.121, -0.66, 0.005, 0.759 };
1085    int lda = 1;
1086    double B[] = { -0.606, 0.049, 0.807, -0.236, -0.258, -0.412, 0.75, -0.659, 0.993, -0.029, -0.968, 0.707, -0.362, -0.005, 0.096, -0.241 };
1087    int ldb = 2;
1088    double C[] = { 0.63, 0.922, 0.025, -0.535 };
1089    int ldc = 2;
1090    double C_expected[] = { 1.117044, 1.983417, -1.276831, -0.447092 };
1091    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1092    {
1093      int i;
1094      for (i = 0; i < 2; i++) {
1095        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1506) real");
1096        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1506) imag");
1097      };
1098    };
1099   };
1100
1101
1102   {
1103    int order = 102;
1104    int transA = 112;
1105    int transB = 111;
1106    int M = 1;
1107    int N = 2;
1108    int K = 4;
1109    double alpha[2] = {-1, 0};
1110    double beta[2] = {1, 0};
1111    double A[] = { 0.064, 0.371, -0.01, -0.262, 0.143, -0.081, 0.1, -0.062 };
1112    int lda = 4;
1113    double B[] = { -0.749, 0.289, -0.239, -0.226, 0.284, 0.668, 0.305, 0.075, -0.36, 0.166, -0.416, 0.234, -0.267, 0.525, 0.116, -0.561 };
1114    int ldb = 4;
1115    double C[] = { 0.671, 0.763, 0.444, -0.246 };
1116    int ldc = 1;
1117    double C_expected[] = { 0.753107, 0.896395, 0.481996, -0.263126 };
1118    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1119    {
1120      int i;
1121      for (i = 0; i < 2; i++) {
1122        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1507) real");
1123        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1507) imag");
1124      };
1125    };
1126   };
1127
1128
1129   {
1130    int order = 101;
1131    int transA = 112;
1132    int transB = 112;
1133    int M = 1;
1134    int N = 2;
1135    int K = 4;
1136    double alpha[2] = {1, 0};
1137    double beta[2] = {-0.3, 0.1};
1138    double A[] = { -0.956, -0.751, 0.671, -0.633, 0.648, -0.042, 0.948, 0.826 };
1139    int lda = 1;
1140    double B[] = { 0.921, 0.506, -0.609, 0.817, -0.686, 0.991, 0.616, -0.482, -0.02, -0.34, 0.559, 0.976, 0.431, 0.385, -0.164, -0.778 };
1141    int ldb = 4;
1142    double C[] = { 0.074, -0.01, 0.165, 0.166 };
1143    int ldc = 2;
1144    double C_expected[] = { 0.166046, 0.491557, 1.473191, -0.033821 };
1145    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1146    {
1147      int i;
1148      for (i = 0; i < 2; i++) {
1149        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1508) real");
1150        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1508) imag");
1151      };
1152    };
1153   };
1154
1155
1156   {
1157    int order = 102;
1158    int transA = 112;
1159    int transB = 112;
1160    int M = 1;
1161    int N = 2;
1162    int K = 4;
1163    double alpha[2] = {1, 0};
1164    double beta[2] = {-0.3, 0.1};
1165    double A[] = { -0.698, -0.062, 0.023, 0.704, 0.443, -0.46, 0.541, 0.296 };
1166    int lda = 4;
1167    double B[] = { 0.787, -0.199, 0.835, -0.276, -0.515, 0.467, -0.76, -0.483, 0.015, -0.394, -0.748, 0.02, 0.573, 0.3, -0.088, -0.238 };
1168    int ldb = 2;
1169    double C[] = { 0.935, -0.655, -0.797, 0.071 };
1170    int ldc = 1;
1171    double C_expected[] = { -1.070679, 0.178755, -0.344714, -0.308137 };
1172    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1173    {
1174      int i;
1175      for (i = 0; i < 2; i++) {
1176        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1509) real");
1177        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1509) imag");
1178      };
1179    };
1180   };
1181
1182
1183   {
1184    int order = 101;
1185    int transA = 112;
1186    int transB = 113;
1187    int M = 1;
1188    int N = 2;
1189    int K = 4;
1190    double alpha[2] = {0, 0.1};
1191    double beta[2] = {-0.3, 0.1};
1192    double A[] = { -0.202, -0.219, 0.741, 0.527, 0.054, 0.16, -0.359, 0.338 };
1193    int lda = 1;
1194    double B[] = { -0.872, 0.995, 0.722, 0.618, -0.27, 0.939, -0.743, 0.547, -0.864, 0.376, -0.997, -0.63, 0.887, -0.454, 0.436, -0.039 };
1195    int ldb = 4;
1196    double C[] = { -0.684, 0.463, -0.386, -0.524 };
1197    int ldc = 2;
1198    double C_expected[] = { 0.1423153, -0.066679, 0.1175618, 0.0012949 };
1199    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1200    {
1201      int i;
1202      for (i = 0; i < 2; i++) {
1203        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1510) real");
1204        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1510) imag");
1205      };
1206    };
1207   };
1208
1209
1210   {
1211    int order = 102;
1212    int transA = 112;
1213    int transB = 113;
1214    int M = 1;
1215    int N = 2;
1216    int K = 4;
1217    double alpha[2] = {0, 0.1};
1218    double beta[2] = {-0.3, 0.1};
1219    double A[] = { -0.855, -0.173, -0.679, 0.824, 0.469, 0.786, 0.757, -0.109 };
1220    int lda = 4;
1221    double B[] = { 0.483, -0.888, -0.757, 0.551, -0.81, 0.23, -0.078, 0.725, -0.592, 0.394, 0.884, 0.802, -0.813, -0.016, -0.853, 0.783 };
1222    int ldb = 2;
1223    double C[] = { 0.181, -0.368, -0.864, -0.784 };
1224    int ldc = 1;
1225    double C_expected[] = { 0.1728438, 0.1183508, 0.2526999, 0.3004174 };
1226    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1227    {
1228      int i;
1229      for (i = 0; i < 2; i++) {
1230        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1511) real");
1231        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1511) imag");
1232      };
1233    };
1234   };
1235
1236
1237   {
1238    int order = 101;
1239    int transA = 113;
1240    int transB = 111;
1241    int M = 1;
1242    int N = 2;
1243    int K = 4;
1244    double alpha[2] = {-1, 0};
1245    double beta[2] = {-0.3, 0.1};
1246    double A[] = { 0.446, -0.65, -0.724, 0.014, 0.792, -0.695, -0.81, -0.358 };
1247    int lda = 1;
1248    double B[] = { -0.08, 0.216, 0.689, 0.699, 0.073, -0.346, 0.821, -0.668, -0.798, 0.869, 0.451, -0.061, -0.41, 0.316, 0.104, -0.514 };
1249    int ldb = 2;
1250    double C[] = { -0.476, 0.211, -0.912, -0.243 };
1251    int ldc = 2;
1252    double C_expected[] = { 1.372475, -0.135616, 0.549353, -1.968747 };
1253    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1254    {
1255      int i;
1256      for (i = 0; i < 2; i++) {
1257        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1512) real");
1258        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1512) imag");
1259      };
1260    };
1261   };
1262
1263
1264   {
1265    int order = 102;
1266    int transA = 113;
1267    int transB = 111;
1268    int M = 1;
1269    int N = 2;
1270    int K = 4;
1271    double alpha[2] = {-1, 0};
1272    double beta[2] = {-0.3, 0.1};
1273    double A[] = { 0.669, 0.046, -0.094, 0.666, 0.23, 0.448, -0.795, -0.142 };
1274    int lda = 4;
1275    double B[] = { 0.037, -0.154, -0.739, 0.905, 0.793, -0.53, -0.34, 0.428, 0.072, -0.263, -0.603, -0.905, 0.681, -0.083, -0.511, -0.337 };
1276    int ldb = 4;
1277    double C[] = { 0.247, 0.575, -0.836, -0.883 };
1278    int ldc = 1;
1279    double C_expected[] = { -0.975939, 0.415528, 0.275533, 0.002716 };
1280    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1281    {
1282      int i;
1283      for (i = 0; i < 2; i++) {
1284        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1513) real");
1285        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1513) imag");
1286      };
1287    };
1288   };
1289
1290
1291   {
1292    int order = 101;
1293    int transA = 113;
1294    int transB = 112;
1295    int M = 1;
1296    int N = 2;
1297    int K = 4;
1298    double alpha[2] = {0, 0};
1299    double beta[2] = {-1, 0};
1300    double A[] = { 0.369, 0.506, 0.217, -0.739, -0.395, 0.16, -0.329, -0.954 };
1301    int lda = 1;
1302    double B[] = { -0.622, -0.945, 0.416, -0.884, 0.797, -0.74, 0.519, -0.789, -0.348, 0.563, -0.398, -0.956, 0.227, 0.84, -0.079, 0.847 };
1303    int ldb = 4;
1304    double C[] = { 0.833, 0.761, 0.074, -0.448 };
1305    int ldc = 2;
1306    double C_expected[] = { -0.833, -0.761, -0.074, 0.448 };
1307    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1308    {
1309      int i;
1310      for (i = 0; i < 2; i++) {
1311        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1514) real");
1312        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1514) imag");
1313      };
1314    };
1315   };
1316
1317
1318   {
1319    int order = 102;
1320    int transA = 113;
1321    int transB = 112;
1322    int M = 1;
1323    int N = 2;
1324    int K = 4;
1325    double alpha[2] = {0, 0};
1326    double beta[2] = {-1, 0};
1327    double A[] = { -0.141, 0.275, 0.717, 0.775, -0.701, -0.689, -0.883, -0.077 };
1328    int lda = 4;
1329    double B[] = { -0.526, -0.437, 0.133, -0.209, -0.83, 0.328, 0.916, -0.337, 0.762, -0.664, -0.566, 0.955, 0.168, 0.488, -0.172, -0.535 };
1330    int ldb = 2;
1331    double C[] = { -0.88, 0.945, 0.416, 0.99 };
1332    int ldc = 1;
1333    double C_expected[] = { 0.88, -0.945, -0.416, -0.99 };
1334    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1335    {
1336      int i;
1337      for (i = 0; i < 2; i++) {
1338        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1515) real");
1339        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1515) imag");
1340      };
1341    };
1342   };
1343
1344
1345   {
1346    int order = 101;
1347    int transA = 113;
1348    int transB = 113;
1349    int M = 1;
1350    int N = 2;
1351    int K = 4;
1352    double alpha[2] = {-0.3, 0.1};
1353    double beta[2] = {0, 0.1};
1354    double A[] = { -0.534, -0.013, -0.258, -0.31, -0.211, -0.883, -0.89, -0.499 };
1355    int lda = 1;
1356    double B[] = { -0.185, -0.798, -0.34, 0.716, 0.035, 0.968, -0.26, 0.784, -0.889, -0.344, -0.685, -0.647, -0.764, 0.03, 0.626, -0.989 };
1357    int ldb = 4;
1358    double C[] = { -0.793, -0.551, 0.182, 0.838 };
1359    int ldc = 2;
1360    double C_expected[] = { -0.5507177, -0.0286821, 0.2222276, 0.5197398 };
1361    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1362    {
1363      int i;
1364      for (i = 0; i < 2; i++) {
1365        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1516) real");
1366        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1516) imag");
1367      };
1368    };
1369   };
1370
1371
1372   {
1373    int order = 102;
1374    int transA = 113;
1375    int transB = 113;
1376    int M = 1;
1377    int N = 2;
1378    int K = 4;
1379    double alpha[2] = {-0.3, 0.1};
1380    double beta[2] = {0, 0.1};
1381    double A[] = { 0.575, -0.128, -0.702, 0.758, 0.383, -0.914, 0.157, 0.368 };
1382    int lda = 4;
1383    double B[] = { 0.572, -0.841, 0.223, -0.334, -0.823, -0.84, 0.671, -0.871, 0.241, 0.927, -0.344, 0.281, -0.034, -0.104, 0.587, -0.329 };
1384    int ldb = 2;
1385    double C[] = { -0.612, 0.167, 0.647, 0.447 };
1386    int ldc = 1;
1387    double C_expected[] = { -0.7876717, 0.0341179, -0.0800018, 0.5717566 };
1388    cblas_zgemm(order, transA, transB, M, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
1389    {
1390      int i;
1391      for (i = 0; i < 2; i++) {
1392        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zgemm(case 1517) real");
1393        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zgemm(case 1517) imag");
1394      };
1395    };
1396   };
1397
1398
1399 }