htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_syr2k.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_syr2k (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int order = 101;
  13    int uplo = 121;
  14    int trans = 111;
  15    int N = 1;
  16    int K = 2;
  17    float alpha = 0.1f;
  18    float beta = 1.0f;
  19    float A[] = { -0.915f, 0.445f };
  20    int lda = 2;
  21    float B[] = { 0.213f, -0.194f };
  22    int ldb = 2;
  23    float C[] = { -0.117f };
  24    int ldc = 1;
  25    float C_expected[] = { -0.173245f };
  26    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  27    {
  28      int i;
  29      for (i = 0; i < 1; i++) {
  30        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1614)");
  31      }
  32    };
  33   };
  34
  35
  36   {
  37    int order = 101;
  38    int uplo = 122;
  39    int trans = 111;
  40    int N = 1;
  41    int K = 2;
  42    float alpha = 0.1f;
  43    float beta = 1.0f;
  44    float A[] = { 0.089f, -0.889f };
  45    int lda = 2;
  46    float B[] = { -0.384f, 0.518f };
  47    int ldb = 2;
  48    float C[] = { 0.069f };
  49    int ldc = 1;
  50    float C_expected[] = { -0.0299356f };
  51    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  52    {
  53      int i;
  54      for (i = 0; i < 1; i++) {
  55        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1615)");
  56      }
  57    };
  58   };
  59
  60
  61   {
  62    int order = 102;
  63    int uplo = 121;
  64    int trans = 111;
  65    int N = 1;
  66    int K = 2;
  67    float alpha = 0.1f;
  68    float beta = 1.0f;
  69    float A[] = { 0.492f, 0.021f };
  70    int lda = 1;
  71    float B[] = { -0.804f, -0.912f };
  72    int ldb = 1;
  73    float C[] = { -0.851f };
  74    int ldc = 1;
  75    float C_expected[] = { -0.933944f };
  76    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
  77    {
  78      int i;
  79      for (i = 0; i < 1; i++) {
  80        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1616)");
  81      }
  82    };
  83   };
  84
  85
  86   {
  87    int order = 102;
  88    int uplo = 122;
  89    int trans = 111;
  90    int N = 1;
  91    int K = 2;
  92    float alpha = 0.1f;
  93    float beta = 1.0f;
  94    float A[] = { -0.376f, 0.689f };
  95    int lda = 1;
  96    float B[] = { 0.21f, 0.406f };
  97    int ldb = 1;
  98    float C[] = { -0.581f };
  99    int ldc = 1;
 100    float C_expected[] = { -0.540845f };
 101    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 102    {
 103      int i;
 104      for (i = 0; i < 1; i++) {
 105        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1617)");
 106      }
 107    };
 108   };
 109
 110
 111   {
 112    int order = 101;
 113    int uplo = 121;
 114    int trans = 112;
 115    int N = 1;
 116    int K = 2;
 117    float alpha = 1.0f;
 118    float beta = -0.3f;
 119    float A[] = { 0.629f, -0.883f };
 120    int lda = 1;
 121    float B[] = { -0.165f, 0.02f };
 122    int ldb = 1;
 123    float C[] = { 0.236f };
 124    int ldc = 1;
 125    float C_expected[] = { -0.31369f };
 126    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 127    {
 128      int i;
 129      for (i = 0; i < 1; i++) {
 130        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1618)");
 131      }
 132    };
 133   };
 134
 135
 136   {
 137    int order = 101;
 138    int uplo = 122;
 139    int trans = 112;
 140    int N = 1;
 141    int K = 2;
 142    float alpha = 1.0f;
 143    float beta = -0.3f;
 144    float A[] = { 0.412f, -0.411f };
 145    int lda = 1;
 146    float B[] = { 0.313f, 0.301f };
 147    int ldb = 1;
 148    float C[] = { 0.222f };
 149    int ldc = 1;
 150    float C_expected[] = { -0.05611f };
 151    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 152    {
 153      int i;
 154      for (i = 0; i < 1; i++) {
 155        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1619)");
 156      }
 157    };
 158   };
 159
 160
 161   {
 162    int order = 102;
 163    int uplo = 121;
 164    int trans = 112;
 165    int N = 1;
 166    int K = 2;
 167    float alpha = 1.0f;
 168    float beta = -0.3f;
 169    float A[] = { -0.02f, 0.593f };
 170    int lda = 2;
 171    float B[] = { -0.144f, 0.846f };
 172    int ldb = 2;
 173    float C[] = { -0.645f };
 174    int ldc = 1;
 175    float C_expected[] = { 1.20262f };
 176    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 177    {
 178      int i;
 179      for (i = 0; i < 1; i++) {
 180        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1620)");
 181      }
 182    };
 183   };
 184
 185
 186   {
 187    int order = 102;
 188    int uplo = 122;
 189    int trans = 112;
 190    int N = 1;
 191    int K = 2;
 192    float alpha = 1.0f;
 193    float beta = -0.3f;
 194    float A[] = { 0.253f, 0.937f };
 195    int lda = 2;
 196    float B[] = { 0.24f, -0.27f };
 197    int ldb = 2;
 198    float C[] = { 0.128f };
 199    int ldc = 1;
 200    float C_expected[] = { -0.42294f };
 201    cblas_ssyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 202    {
 203      int i;
 204      for (i = 0; i < 1; i++) {
 205        gsl_test_rel(C[i], C_expected[i], flteps, "ssyr2k(case 1621)");
 206      }
 207    };
 208   };
 209
 210
 211   {
 212    int order = 101;
 213    int uplo = 121;
 214    int trans = 111;
 215    int N = 1;
 216    int K = 2;
 217    double alpha = 0.1;
 218    double beta = 0;
 219    double A[] = { -0.225, 0.857 };
 220    int lda = 2;
 221    double B[] = { -0.933, 0.994 };
 222    int ldb = 2;
 223    double C[] = { 0.177 };
 224    int ldc = 1;
 225    double C_expected[] = { 0.2123566 };
 226    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 227    {
 228      int i;
 229      for (i = 0; i < 1; i++) {
 230        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1622)");
 231      }
 232    };
 233   };
 234
 235
 236   {
 237    int order = 101;
 238    int uplo = 122;
 239    int trans = 111;
 240    int N = 1;
 241    int K = 2;
 242    double alpha = 0.1;
 243    double beta = 0;
 244    double A[] = { -0.955, 0.112 };
 245    int lda = 2;
 246    double B[] = { -0.695, 0.719 };
 247    int ldb = 2;
 248    double C[] = { 0.069 };
 249    int ldc = 1;
 250    double C_expected[] = { 0.1488506 };
 251    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 252    {
 253      int i;
 254      for (i = 0; i < 1; i++) {
 255        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1623)");
 256      }
 257    };
 258   };
 259
 260
 261   {
 262    int order = 102;
 263    int uplo = 121;
 264    int trans = 111;
 265    int N = 1;
 266    int K = 2;
 267    double alpha = 0.1;
 268    double beta = 0;
 269    double A[] = { 0.216, 0.911 };
 270    int lda = 1;
 271    double B[] = { -0.074, -0.256 };
 272    int ldb = 1;
 273    double C[] = { -0.621 };
 274    int ldc = 1;
 275    double C_expected[] = { -0.04984 };
 276    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 277    {
 278      int i;
 279      for (i = 0; i < 1; i++) {
 280        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1624)");
 281      }
 282    };
 283   };
 284
 285
 286   {
 287    int order = 102;
 288    int uplo = 122;
 289    int trans = 111;
 290    int N = 1;
 291    int K = 2;
 292    double alpha = 0.1;
 293    double beta = 0;
 294    double A[] = { -0.343, -0.381 };
 295    int lda = 1;
 296    double B[] = { -0.433, -0.087 };
 297    int ldb = 1;
 298    double C[] = { -0.889 };
 299    int ldc = 1;
 300    double C_expected[] = { 0.0363332 };
 301    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 302    {
 303      int i;
 304      for (i = 0; i < 1; i++) {
 305        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1625)");
 306      }
 307    };
 308   };
 309
 310
 311   {
 312    int order = 101;
 313    int uplo = 121;
 314    int trans = 112;
 315    int N = 1;
 316    int K = 2;
 317    double alpha = 1;
 318    double beta = -1;
 319    double A[] = { -0.633, 0.219 };
 320    int lda = 1;
 321    double B[] = { 0.817, -0.683 };
 322    int ldb = 1;
 323    double C[] = { -0.294 };
 324    int ldc = 1;
 325    double C_expected[] = { -1.039476 };
 326    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 327    {
 328      int i;
 329      for (i = 0; i < 1; i++) {
 330        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1626)");
 331      }
 332    };
 333   };
 334
 335
 336   {
 337    int order = 101;
 338    int uplo = 122;
 339    int trans = 112;
 340    int N = 1;
 341    int K = 2;
 342    double alpha = 1;
 343    double beta = -1;
 344    double A[] = { -0.887, -0.43 };
 345    int lda = 1;
 346    double B[] = { 0.557, 0.912 };
 347    int ldb = 1;
 348    double C[] = { 0.831 };
 349    int ldc = 1;
 350    double C_expected[] = { -2.603438 };
 351    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 352    {
 353      int i;
 354      for (i = 0; i < 1; i++) {
 355        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1627)");
 356      }
 357    };
 358   };
 359
 360
 361   {
 362    int order = 102;
 363    int uplo = 121;
 364    int trans = 112;
 365    int N = 1;
 366    int K = 2;
 367    double alpha = 1;
 368    double beta = -1;
 369    double A[] = { 0.397, -0.173 };
 370    int lda = 2;
 371    double B[] = { 0.155, -0.99 };
 372    int ldb = 2;
 373    double C[] = { 0.621 };
 374    int ldc = 1;
 375    double C_expected[] = { -0.15539 };
 376    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 377    {
 378      int i;
 379      for (i = 0; i < 1; i++) {
 380        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1628)");
 381      }
 382    };
 383   };
 384
 385
 386   {
 387    int order = 102;
 388    int uplo = 122;
 389    int trans = 112;
 390    int N = 1;
 391    int K = 2;
 392    double alpha = 1;
 393    double beta = -1;
 394    double A[] = { 0.833, -0.52 };
 395    int lda = 2;
 396    double B[] = { 0.28, 0.481 };
 397    int ldb = 2;
 398    double C[] = { 0.455 };
 399    int ldc = 1;
 400    double C_expected[] = { -0.48876 };
 401    cblas_dsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 402    {
 403      int i;
 404      for (i = 0; i < 1; i++) {
 405        gsl_test_rel(C[i], C_expected[i], dbleps, "dsyr2k(case 1629)");
 406      }
 407    };
 408   };
 409
 410
 411   {
 412    int order = 101;
 413    int uplo = 121;
 414    int trans = 111;
 415    int N = 1;
 416    int K = 2;
 417    float alpha[2] = {0.0f, 0.1f};
 418    float beta[2] = {0.0f, 0.0f};
 419    float A[] = { -0.248f, -0.037f, -0.124f, 0.998f };
 420    int lda = 2;
 421    float B[] = { -0.608f, -0.115f, -0.718f, -0.551f };
 422    int ldb = 2;
 423    float C[] = { 0.187f, -0.329f };
 424    int ldc = 1;
 425    float C_expected[] = { 0.119445f, 0.157092f };
 426    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 427    {
 428      int i;
 429      for (i = 0; i < 1; i++) {
 430        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1630) real");
 431        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1630) imag");
 432      };
 433    };
 434   };
 435
 436
 437   {
 438    int order = 101;
 439    int uplo = 122;
 440    int trans = 111;
 441    int N = 1;
 442    int K = 2;
 443    float alpha[2] = {0.0f, 0.1f};
 444    float beta[2] = {0.0f, 0.0f};
 445    float A[] = { 0.068f, 0.751f, -0.449f, -0.598f };
 446    int lda = 2;
 447    float B[] = { 0.616f, 0.805f, -0.635f, 0.773f };
 448    int ldb = 2;
 449    float C[] = { -0.287f, 0.917f };
 450    int ldc = 1;
 451    float C_expected[] = { -0.110002f, 0.0369404f };
 452    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 453    {
 454      int i;
 455      for (i = 0; i < 1; i++) {
 456        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1631) real");
 457        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1631) imag");
 458      };
 459    };
 460   };
 461
 462
 463   {
 464    int order = 102;
 465    int uplo = 121;
 466    int trans = 111;
 467    int N = 1;
 468    int K = 2;
 469    float alpha[2] = {0.0f, 0.1f};
 470    float beta[2] = {0.0f, 0.0f};
 471    float A[] = { -0.396f, -0.603f, -0.131f, -0.288f };
 472    int lda = 1;
 473    float B[] = { -0.64f, -0.444f, -0.085f, 0.936f };
 474    int ldb = 1;
 475    float C[] = { 0.375f, -0.434f };
 476    int ldc = 1;
 477    float C_expected[] = { -0.0927216f, 0.0532822f };
 478    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 479    {
 480      int i;
 481      for (i = 0; i < 1; i++) {
 482        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1632) real");
 483        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1632) imag");
 484      };
 485    };
 486   };
 487
 488
 489   {
 490    int order = 102;
 491    int uplo = 122;
 492    int trans = 111;
 493    int N = 1;
 494    int K = 2;
 495    float alpha[2] = {0.0f, 0.1f};
 496    float beta[2] = {0.0f, 0.0f};
 497    float A[] = { 0.655f, 0.16f, 0.45f, -0.747f };
 498    int lda = 1;
 499    float B[] = { 0.923f, 0.432f, -0.986f, 0.259f };
 500    int ldb = 1;
 501    float C[] = { 0.752f, 0.576f };
 502    int ldc = 1;
 503    float C_expected[] = { -0.256746f, 0.0570436f };
 504    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 505    {
 506      int i;
 507      for (i = 0; i < 1; i++) {
 508        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1633) real");
 509        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1633) imag");
 510      };
 511    };
 512   };
 513
 514
 515   {
 516    int order = 101;
 517    int uplo = 121;
 518    int trans = 112;
 519    int N = 1;
 520    int K = 2;
 521    float alpha[2] = {0.0f, 0.1f};
 522    float beta[2] = {0.0f, 0.0f};
 523    float A[] = { -0.765f, 0.487f, 0.7f, 0.768f };
 524    int lda = 1;
 525    float B[] = { -0.529f, 0.056f, -0.584f, 0.928f };
 526    int ldb = 1;
 527    float C[] = { -0.426f, 0.836f };
 528    int ldc = 1;
 529    float C_expected[] = { 0.019875f, -0.148818f };
 530    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 531    {
 532      int i;
 533      for (i = 0; i < 1; i++) {
 534        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1634) real");
 535        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1634) imag");
 536      };
 537    };
 538   };
 539
 540
 541   {
 542    int order = 101;
 543    int uplo = 122;
 544    int trans = 112;
 545    int N = 1;
 546    int K = 2;
 547    float alpha[2] = {0.0f, 0.1f};
 548    float beta[2] = {0.0f, 0.0f};
 549    float A[] = { 0.25f, 0.489f, 0.8f, -0.642f };
 550    int lda = 1;
 551    float B[] = { -0.732f, -0.856f, -0.654f, 0.591f };
 552    int ldb = 1;
 553    float C[] = { -0.101f, 0.322f };
 554    int ldc = 1;
 555    float C_expected[] = { -0.064144f, 0.0183612f };
 556    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 557    {
 558      int i;
 559      for (i = 0; i < 1; i++) {
 560        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1635) real");
 561        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1635) imag");
 562      };
 563    };
 564   };
 565
 566
 567   {
 568    int order = 102;
 569    int uplo = 121;
 570    int trans = 112;
 571    int N = 1;
 572    int K = 2;
 573    float alpha[2] = {0.0f, 0.1f};
 574    float beta[2] = {0.0f, 0.0f};
 575    float A[] = { -0.579f, -0.971f, 0.521f, -0.824f };
 576    int lda = 2;
 577    float B[] = { -0.227f, 0.907f, 0.457f, -0.274f };
 578    int ldb = 2;
 579    float C[] = { 0.21f, -0.718f };
 580    int ldc = 1;
 581    float C_expected[] = { 0.164812f, 0.20489f };
 582    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 583    {
 584      int i;
 585      for (i = 0; i < 1; i++) {
 586        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1636) real");
 587        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1636) imag");
 588      };
 589    };
 590   };
 591
 592
 593   {
 594    int order = 102;
 595    int uplo = 122;
 596    int trans = 112;
 597    int N = 1;
 598    int K = 2;
 599    float alpha[2] = {0.0f, 0.1f};
 600    float beta[2] = {0.0f, 0.0f};
 601    float A[] = { -0.83f, -0.512f, -0.667f, -0.436f };
 602    int lda = 2;
 603    float B[] = { -0.443f, 0.82f, -0.259f, -0.618f };
 604    int ldb = 2;
 605    float C[] = { 0.583f, 0.668f };
 606    int ldc = 1;
 607    float C_expected[] = { -0.0142692f, 0.138167f };
 608    cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 609    {
 610      int i;
 611      for (i = 0; i < 1; i++) {
 612        gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1637) real");
 613        gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1637) imag");
 614      };
 615    };
 616   };
 617
 618
 619   {
 620    int order = 101;
 621    int uplo = 121;
 622    int trans = 111;
 623    int N = 1;
 624    int K = 2;
 625    double alpha[2] = {0, 0};
 626    double beta[2] = {-0.3, 0.1};
 627    double A[] = { -0.315, 0.03, 0.281, 0.175 };
 628    int lda = 2;
 629    double B[] = { -0.832, -0.964, 0.291, 0.476 };
 630    int ldb = 2;
 631    double C[] = { -0.341, 0.743 };
 632    int ldc = 1;
 633    double C_expected[] = { 0.028, -0.257 };
 634    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 635    {
 636      int i;
 637      for (i = 0; i < 1; i++) {
 638        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1638) real");
 639        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1638) imag");
 640      };
 641    };
 642   };
 643
 644
 645   {
 646    int order = 101;
 647    int uplo = 122;
 648    int trans = 111;
 649    int N = 1;
 650    int K = 2;
 651    double alpha[2] = {0, 0};
 652    double beta[2] = {-0.3, 0.1};
 653    double A[] = { -0.159, -0.489, -0.11, 0.611 };
 654    int lda = 2;
 655    double B[] = { -0.285, -0.048, -0.673, -0.492 };
 656    int ldb = 2;
 657    double C[] = { 0.496, -0.626 };
 658    int ldc = 1;
 659    double C_expected[] = { -0.0862, 0.2374 };
 660    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 661    {
 662      int i;
 663      for (i = 0; i < 1; i++) {
 664        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1639) real");
 665        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1639) imag");
 666      };
 667    };
 668   };
 669
 670
 671   {
 672    int order = 102;
 673    int uplo = 121;
 674    int trans = 111;
 675    int N = 1;
 676    int K = 2;
 677    double alpha[2] = {0, 0};
 678    double beta[2] = {-0.3, 0.1};
 679    double A[] = { 0.796, 0.872, -0.919, 0.748 };
 680    int lda = 1;
 681    double B[] = { -0.945, 0.915, -0.252, -0.276 };
 682    int ldb = 1;
 683    double C[] = { 0.07, -0.957 };
 684    int ldc = 1;
 685    double C_expected[] = { 0.0747, 0.2941 };
 686    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 687    {
 688      int i;
 689      for (i = 0; i < 1; i++) {
 690        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1640) real");
 691        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1640) imag");
 692      };
 693    };
 694   };
 695
 696
 697   {
 698    int order = 102;
 699    int uplo = 122;
 700    int trans = 111;
 701    int N = 1;
 702    int K = 2;
 703    double alpha[2] = {0, 0};
 704    double beta[2] = {-0.3, 0.1};
 705    double A[] = { 0.984, 0.526, 0.284, 0.806 };
 706    int lda = 1;
 707    double B[] = { -0.509, -0.178, 0.188, -0.221 };
 708    int ldb = 1;
 709    double C[] = { -0.388, 0.795 };
 710    int ldc = 1;
 711    double C_expected[] = { 0.0369, -0.2773 };
 712    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 713    {
 714      int i;
 715      for (i = 0; i < 1; i++) {
 716        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1641) real");
 717        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1641) imag");
 718      };
 719    };
 720   };
 721
 722
 723   {
 724    int order = 101;
 725    int uplo = 121;
 726    int trans = 112;
 727    int N = 1;
 728    int K = 2;
 729    double alpha[2] = {1, 0};
 730    double beta[2] = {0, 0.1};
 731    double A[] = { 0.628, 0.846, -0.645, 0.032 };
 732    int lda = 1;
 733    double B[] = { 0.545, -0.54, 0.493, -0.035 };
 734    int ldb = 1;
 735    double C[] = { -0.16, -0.06 };
 736    int ldc = 1;
 737    double C_expected[] = { 0.97047, 0.304602 };
 738    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 739    {
 740      int i;
 741      for (i = 0; i < 1; i++) {
 742        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1642) real");
 743        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1642) imag");
 744      };
 745    };
 746   };
 747
 748
 749   {
 750    int order = 101;
 751    int uplo = 122;
 752    int trans = 112;
 753    int N = 1;
 754    int K = 2;
 755    double alpha[2] = {1, 0};
 756    double beta[2] = {0, 0.1};
 757    double A[] = { -0.556, -0.946, 0.177, -0.859 };
 758    int lda = 1;
 759    double B[] = { 0.423, -0.91, 0.736, -0.251 };
 760    int ldb = 1;
 761    double C[] = { -0.478, 0.519 };
 762    int ldc = 1;
 763    double C_expected[] = { -2.41467, -1.189498 };
 764    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 765    {
 766      int i;
 767      for (i = 0; i < 1; i++) {
 768        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1643) real");
 769        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1643) imag");
 770      };
 771    };
 772   };
 773
 774
 775   {
 776    int order = 102;
 777    int uplo = 121;
 778    int trans = 112;
 779    int N = 1;
 780    int K = 2;
 781    double alpha[2] = {1, 0};
 782    double beta[2] = {0, 0.1};
 783    double A[] = { -0.582, 0.09, -0.176, 0.784 };
 784    int lda = 2;
 785    double B[] = { 0.687, -0.859, 0.945, 0.756 };
 786    int ldb = 2;
 787    double C[] = { -0.663, -0.186 };
 788    int ldc = 1;
 789    double C_expected[] = { -2.144496, 2.272884 };
 790    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 791    {
 792      int i;
 793      for (i = 0; i < 1; i++) {
 794        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1644) real");
 795        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1644) imag");
 796      };
 797    };
 798   };
 799
 800
 801   {
 802    int order = 102;
 803    int uplo = 122;
 804    int trans = 112;
 805    int N = 1;
 806    int K = 2;
 807    double alpha[2] = {1, 0};
 808    double beta[2] = {0, 0.1};
 809    double A[] = { 0.231, -0.452, -0.112, -0.837 };
 810    int lda = 2;
 811    double B[] = { -0.258, 0.464, -0.224, 0.893 };
 812    int ldb = 2;
 813    double C[] = { -0.448, 0.046 };
 814    int ldc = 1;
 815    double C_expected[] = { 1.840718, 0.577744 };
 816    cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
 817    {
 818      int i;
 819      for (i = 0; i < 1; i++) {
 820        gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1645) real");
 821        gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1645) imag");
 822      };
 823    };
 824   };
 825
 826
 827 }