htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_tbmv.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_tbmv (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int order = 101;
  13    int trans = 111;
  14    int uplo = 121;
  15    int diag = 131;
  16    int N = 3;
  17    int K = 1;
  18    int lda = 3;
  19    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  20    float X[] = { -0.089f, -0.688f, -0.203f };
  21    int incX = -1;
  22    float x_expected[] = { 0.017088f, 0.315595f, 0.243875f };
  23    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  24    {
  25      int i;
  26      for (i = 0; i < 3; i++) {
  27        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 894)");
  28      }
  29    };
  30   };
  31
  32
  33   {
  34    int order = 101;
  35    int trans = 111;
  36    int uplo = 121;
  37    int diag = 132;
  38    int N = 3;
  39    int K = 1;
  40    int lda = 3;
  41    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  42    float X[] = { -0.089f, -0.688f, -0.203f };
  43    int incX = -1;
  44    float x_expected[] = { -0.089f, -0.721909f, 0.129992f };
  45    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  46    {
  47      int i;
  48      for (i = 0; i < 3; i++) {
  49        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 895)");
  50      }
  51    };
  52   };
  53
  54
  55   {
  56    int order = 101;
  57    int trans = 111;
  58    int uplo = 122;
  59    int diag = 131;
  60    int N = 3;
  61    int K = 1;
  62    int lda = 3;
  63    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  64    float X[] = { -0.089f, -0.688f, -0.203f };
  65    int incX = -1;
  66    float x_expected[] = { 0.156927f, -0.159004f, 0.098252f };
  67    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  68    {
  69      int i;
  70      for (i = 0; i < 3; i++) {
  71        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 896)");
  72      }
  73    };
  74   };
  75
  76
  77   {
  78    int order = 101;
  79    int trans = 111;
  80    int uplo = 122;
  81    int diag = 132;
  82    int N = 3;
  83    int K = 1;
  84    int lda = 3;
  85    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  86    float X[] = { -0.089f, -0.688f, -0.203f };
  87    int incX = -1;
  88    float x_expected[] = { 0.043096f, -0.584876f, -0.203f };
  89    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  90    {
  91      int i;
  92      for (i = 0; i < 3; i++) {
  93        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 897)");
  94      }
  95    };
  96   };
  97
  98
  99   {
 100    int order = 102;
 101    int trans = 111;
 102    int uplo = 121;
 103    int diag = 131;
 104    int N = 3;
 105    int K = 1;
 106    int lda = 3;
 107    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
 108    float X[] = { -0.089f, -0.688f, -0.203f };
 109    int incX = -1;
 110    float x_expected[] = { 0.024831f, -0.24504f, 0.447756f };
 111    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 112    {
 113      int i;
 114      for (i = 0; i < 3; i++) {
 115        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 898)");
 116      }
 117    };
 118   };
 119
 120
 121   {
 122    int order = 102;
 123    int trans = 111;
 124    int uplo = 121;
 125    int diag = 132;
 126    int N = 3;
 127    int K = 1;
 128    int lda = 3;
 129    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
 130    float X[] = { -0.089f, -0.688f, -0.203f };
 131    int incX = -1;
 132    float x_expected[] = { -0.089f, -0.670912f, 0.146504f };
 133    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 134    {
 135      int i;
 136      for (i = 0; i < 3; i++) {
 137        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 899)");
 138      }
 139    };
 140   };
 141
 142
 143   {
 144    int order = 102;
 145    int trans = 111;
 146    int uplo = 122;
 147    int diag = 131;
 148    int N = 3;
 149    int K = 1;
 150    int lda = 3;
 151    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
 152    float X[] = { -0.089f, -0.688f, -0.203f };
 153    int incX = -1;
 154    float x_expected[] = { -0.24504f, 0.447756f, -0.089117f };
 155    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 156    {
 157      int i;
 158      for (i = 0; i < 3; i++) {
 159        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 900)");
 160      }
 161    };
 162   };
 163
 164
 165   {
 166    int order = 102;
 167    int trans = 111;
 168    int uplo = 122;
 169    int diag = 132;
 170    int N = 3;
 171    int K = 1;
 172    int lda = 3;
 173    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
 174    float X[] = { -0.089f, -0.688f, -0.203f };
 175    int incX = -1;
 176    float x_expected[] = { -0.351128f, -0.589748f, -0.203f };
 177    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 178    {
 179      int i;
 180      for (i = 0; i < 3; i++) {
 181        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 901)");
 182      }
 183    };
 184   };
 185
 186
 187   {
 188    int order = 101;
 189    int trans = 112;
 190    int uplo = 121;
 191    int diag = 131;
 192    int N = 3;
 193    int K = 1;
 194    int lda = 3;
 195    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
 196    float X[] = { 0.023f, -0.501f, -0.562f };
 197    int incX = -1;
 198    float x_expected[] = { 0.156047f, 0.189418f, -0.52828f };
 199    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 200    {
 201      int i;
 202      for (i = 0; i < 3; i++) {
 203        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 902)");
 204      }
 205    };
 206   };
 207
 208
 209   {
 210    int order = 101;
 211    int trans = 112;
 212    int uplo = 121;
 213    int diag = 132;
 214    int N = 3;
 215    int K = 1;
 216    int lda = 3;
 217    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
 218    float X[] = { 0.023f, -0.501f, -0.562f };
 219    int incX = -1;
 220    float x_expected[] = { 0.194342f, -0.449858f, -0.562f };
 221    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 222    {
 223      int i;
 224      for (i = 0; i < 3; i++) {
 225        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 903)");
 226      }
 227    };
 228   };
 229
 230
 231   {
 232    int order = 101;
 233    int trans = 112;
 234    int uplo = 122;
 235    int diag = 131;
 236    int N = 3;
 237    int K = 1;
 238    int lda = 3;
 239    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
 240    float X[] = { 0.023f, -0.501f, -0.562f };
 241    int incX = -1;
 242    float x_expected[] = { -0.0046f, 0.156047f, 0.189418f };
 243    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 244    {
 245      int i;
 246      for (i = 0; i < 3; i++) {
 247        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 904)");
 248      }
 249    };
 250   };
 251
 252
 253   {
 254    int order = 101;
 255    int trans = 112;
 256    int uplo = 122;
 257    int diag = 132;
 258    int N = 3;
 259    int K = 1;
 260    int lda = 3;
 261    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
 262    float X[] = { 0.023f, -0.501f, -0.562f };
 263    int incX = -1;
 264    float x_expected[] = { 0.023f, -0.516295f, -0.423724f };
 265    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 266    {
 267      int i;
 268      for (i = 0; i < 3; i++) {
 269        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 905)");
 270      }
 271    };
 272   };
 273
 274
 275   {
 276    int order = 102;
 277    int trans = 112;
 278    int uplo = 121;
 279    int diag = 131;
 280    int N = 3;
 281    int K = 1;
 282    int lda = 3;
 283    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
 284    float X[] = { 0.023f, -0.501f, -0.562f };
 285    int incX = -1;
 286    float x_expected[] = { 0.328565f, 0.326454f, 0.051142f };
 287    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 288    {
 289      int i;
 290      for (i = 0; i < 3; i++) {
 291        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 906)");
 292      }
 293    };
 294   };
 295
 296
 297   {
 298    int order = 102;
 299    int trans = 112;
 300    int uplo = 121;
 301    int diag = 132;
 302    int N = 3;
 303    int K = 1;
 304    int lda = 3;
 305    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
 306    float X[] = { 0.023f, -0.501f, -0.562f };
 307    int incX = -1;
 308    float x_expected[] = { 0.356165f, -0.345888f, -0.562f };
 309    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 310    {
 311      int i;
 312      for (i = 0; i < 3; i++) {
 313        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 907)");
 314      }
 315    };
 316   };
 317
 318
 319   {
 320    int order = 102;
 321    int trans = 112;
 322    int uplo = 122;
 323    int diag = 131;
 324    int N = 3;
 325    int K = 1;
 326    int lda = 3;
 327    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
 328    float X[] = { 0.023f, -0.501f, -0.562f };
 329    int incX = -1;
 330    float x_expected[] = { -0.015295f, 0.13041f, -0.482689f };
 331    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 332    {
 333      int i;
 334      for (i = 0; i < 3; i++) {
 335        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 908)");
 336      }
 337    };
 338   };
 339
 340
 341   {
 342    int order = 102;
 343    int trans = 112;
 344    int uplo = 122;
 345    int diag = 132;
 346    int N = 3;
 347    int K = 1;
 348    int lda = 3;
 349    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
 350    float X[] = { 0.023f, -0.501f, -0.562f };
 351    int incX = -1;
 352    float x_expected[] = { 0.023f, -0.508866f, -0.516409f };
 353    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 354    {
 355      int i;
 356      for (i = 0; i < 3; i++) {
 357        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 909)");
 358      }
 359    };
 360   };
 361
 362
 363   {
 364    int order = 101;
 365    int trans = 111;
 366    int uplo = 121;
 367    int diag = 131;
 368    int N = 3;
 369    int K = 1;
 370    int lda = 3;
 371    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
 372    double X[] = { -0.77, -0.818, 0.337 };
 373    int incX = -1;
 374    double x_expected[] = { 0.50204, 0.563918, -0.590448 };
 375    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 376    {
 377      int i;
 378      for (i = 0; i < 3; i++) {
 379        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 910)");
 380      }
 381    };
 382   };
 383
 384
 385   {
 386    int order = 101;
 387    int trans = 111;
 388    int uplo = 121;
 389    int diag = 132;
 390    int N = 3;
 391    int K = 1;
 392    int lda = 3;
 393    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
 394    double X[] = { -0.77, -0.818, 0.337 };
 395    int incX = -1;
 396    double x_expected[] = { -0.77, -0.95429, -0.44419 };
 397    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 398    {
 399      int i;
 400      for (i = 0; i < 3; i++) {
 401        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 911)");
 402      }
 403    };
 404   };
 405
 406
 407   {
 408    int order = 101;
 409    int trans = 111;
 410    int uplo = 122;
 411    int diag = 131;
 412    int N = 3;
 413    int K = 1;
 414    int lda = 3;
 415    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
 416    double X[] = { -0.77, -0.818, 0.337 };
 417    int incX = -1;
 418    double x_expected[] = { 1.214016, -0.433258, 0.321835 };
 419    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 420    {
 421      int i;
 422      for (i = 0; i < 3; i++) {
 423        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 912)");
 424      }
 425    };
 426   };
 427
 428
 429   {
 430    int order = 101;
 431    int trans = 111;
 432    int uplo = 122;
 433    int diag = 132;
 434    int N = 3;
 435    int K = 1;
 436    int lda = 3;
 437    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
 438    double X[] = { -0.77, -0.818, 0.337 };
 439    int incX = -1;
 440    double x_expected[] = { -0.236664, -1.106472, 0.337 };
 441    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 442    {
 443      int i;
 444      for (i = 0; i < 3; i++) {
 445        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 913)");
 446      }
 447    };
 448   };
 449
 450
 451   {
 452    int order = 102;
 453    int trans = 111;
 454    int uplo = 121;
 455    int diag = 131;
 456    int N = 3;
 457    int K = 1;
 458    int lda = 3;
 459    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
 460    double X[] = { -0.77, -0.818, 0.337 };
 461    int incX = -1;
 462    double x_expected[] = { 0.68068, 0.357254, 1.022043 };
 463    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 464    {
 465      int i;
 466      for (i = 0; i < 3; i++) {
 467        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 914)");
 468      }
 469    };
 470   };
 471
 472
 473   {
 474    int order = 102;
 475    int trans = 111;
 476    int uplo = 121;
 477    int diag = 132;
 478    int N = 3;
 479    int K = 1;
 480    int lda = 3;
 481    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
 482    double X[] = { -0.77, -0.818, 0.337 };
 483    int incX = -1;
 484    double x_expected[] = { -0.77, -0.31596, 1.037208 };
 485    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 486    {
 487      int i;
 488      for (i = 0; i < 3; i++) {
 489        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 915)");
 490      }
 491    };
 492   };
 493
 494
 495   {
 496    int order = 102;
 497    int trans = 111;
 498    int uplo = 122;
 499    int diag = 131;
 500    int N = 3;
 501    int K = 1;
 502    int lda = 3;
 503    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
 504    double X[] = { -0.77, -0.818, 0.337 };
 505    int incX = -1;
 506    double x_expected[] = { 0.357254, 1.022043, 0.190742 };
 507    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 508    {
 509      int i;
 510      for (i = 0; i < 3; i++) {
 511        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 916)");
 512      }
 513    };
 514   };
 515
 516
 517   {
 518    int order = 102;
 519    int trans = 111;
 520    int uplo = 122;
 521    int diag = 132;
 522    int N = 3;
 523    int K = 1;
 524    int lda = 3;
 525    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
 526    double X[] = { -0.77, -0.818, 0.337 };
 527    int incX = -1;
 528    double x_expected[] = { -0.914786, -0.496165, 0.337 };
 529    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 530    {
 531      int i;
 532      for (i = 0; i < 3; i++) {
 533        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 917)");
 534      }
 535    };
 536   };
 537
 538
 539   {
 540    int order = 101;
 541    int trans = 112;
 542    int uplo = 121;
 543    int diag = 131;
 544    int N = 3;
 545    int K = 1;
 546    int lda = 3;
 547    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
 548    double X[] = { -0.851, 0.481, 0.155 };
 549    int incX = -1;
 550    double x_expected[] = { 0.610833, -0.293243, 0.02914 };
 551    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 552    {
 553      int i;
 554      for (i = 0; i < 3; i++) {
 555        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 918)");
 556      }
 557    };
 558   };
 559
 560
 561   {
 562    int order = 101;
 563    int trans = 112;
 564    int uplo = 121;
 565    int diag = 132;
 566    int N = 3;
 567    int K = 1;
 568    int lda = 3;
 569    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
 570    double X[] = { -0.851, 0.481, 0.155 };
 571    int incX = -1;
 572    double x_expected[] = { -0.635031, 0.574, 0.155 };
 573    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 574    {
 575      int i;
 576      for (i = 0; i < 3; i++) {
 577        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 919)");
 578      }
 579    };
 580   };
 581
 582
 583   {
 584    int order = 101;
 585    int trans = 112;
 586    int uplo = 122;
 587    int diag = 131;
 588    int N = 3;
 589    int K = 1;
 590    int lda = 3;
 591    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
 592    double X[] = { -0.851, 0.481, 0.155 };
 593    int incX = -1;
 594    double x_expected[] = { 0.024679, 0.610833, -0.293243 };
 595    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 596    {
 597      int i;
 598      for (i = 0; i < 3; i++) {
 599        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 920)");
 600      }
 601    };
 602   };
 603
 604
 605   {
 606    int order = 101;
 607    int trans = 112;
 608    int uplo = 122;
 609    int diag = 132;
 610    int N = 3;
 611    int K = 1;
 612    int lda = 3;
 613    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
 614    double X[] = { -0.851, 0.481, 0.155 };
 615    int incX = -1;
 616    double x_expected[] = { -0.851, 0.875864, -0.231243 };
 617    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 618    {
 619      int i;
 620      for (i = 0; i < 3; i++) {
 621        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 921)");
 622      }
 623    };
 624   };
 625
 626
 627   {
 628    int order = 102;
 629    int trans = 112;
 630    int uplo = 121;
 631    int diag = 131;
 632    int N = 3;
 633    int K = 1;
 634    int lda = 3;
 635    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
 636    double X[] = { -0.851, 0.481, 0.155 };
 637    int incX = -1;
 638    double x_expected[] = { -0.198505, 0.091504, 0.093 };
 639    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 640    {
 641      int i;
 642      for (i = 0; i < 3; i++) {
 643        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 922)");
 644      }
 645    };
 646   };
 647
 648
 649   {
 650    int order = 102;
 651    int trans = 112;
 652    int uplo = 121;
 653    int diag = 132;
 654    int N = 3;
 655    int K = 1;
 656    int lda = 3;
 657    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
 658    double X[] = { -0.851, 0.481, 0.155 };
 659    int incX = -1;
 660    double x_expected[] = { -1.074184, 0.356535, 0.155 };
 661    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 662    {
 663      int i;
 664      for (i = 0; i < 3; i++) {
 665        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 923)");
 666      }
 667    };
 668   };
 669
 670
 671   {
 672    int order = 102;
 673    int trans = 112;
 674    int uplo = 122;
 675    int diag = 131;
 676    int N = 3;
 677    int K = 1;
 678    int lda = 3;
 679    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
 680    double X[] = { -0.851, 0.481, 0.155 };
 681    int incX = -1;
 682    double x_expected[] = { 0.394864, -0.768342, 0.31774 };
 683    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 684    {
 685      int i;
 686      for (i = 0; i < 3; i++) {
 687        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 924)");
 688      }
 689    };
 690   };
 691
 692
 693   {
 694    int order = 102;
 695    int trans = 112;
 696    int uplo = 122;
 697    int diag = 132;
 698    int N = 3;
 699    int K = 1;
 700    int lda = 3;
 701    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
 702    double X[] = { -0.851, 0.481, 0.155 };
 703    int incX = -1;
 704    double x_expected[] = { -0.851, 0.098901, 0.4436 };
 705    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 706    {
 707      int i;
 708      for (i = 0; i < 3; i++) {
 709        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 925)");
 710      }
 711    };
 712   };
 713
 714
 715   {
 716    int order = 101;
 717    int trans = 111;
 718    int uplo = 121;
 719    int diag = 131;
 720    int N = 3;
 721    int K = 1;
 722    int lda = 3;
 723    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
 724    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
 725    int incX = -1;
 726    float x_expected[] = { -0.113114f, -0.051704f, -0.403567f, -0.288349f, -0.223936f, 0.841145f };
 727    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 728    {
 729      int i;
 730      for (i = 0; i < 3; i++) {
 731        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 926) real");
 732        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 926) imag");
 733      };
 734    };
 735   };
 736
 737
 738   {
 739    int order = 101;
 740    int trans = 111;
 741    int uplo = 121;
 742    int diag = 132;
 743    int N = 3;
 744    int K = 1;
 745    int lda = 3;
 746    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
 747    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
 748    int incX = -1;
 749    float x_expected[] = { -0.46f, 0.069f, -0.14027f, -0.23208f, -0.537722f, 0.841425f };
 750    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 751    {
 752      int i;
 753      for (i = 0; i < 3; i++) {
 754        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 927) real");
 755        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 927) imag");
 756      };
 757    };
 758   };
 759
 760
 761   {
 762    int order = 101;
 763    int trans = 111;
 764    int uplo = 122;
 765    int diag = 131;
 766    int N = 3;
 767    int K = 1;
 768    int lda = 3;
 769    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
 770    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
 771    int incX = -1;
 772    float x_expected[] = { -0.099689f, 0.487805f, 0.353793f, 0.325411f, -0.225658f, -0.776023f };
 773    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 774    {
 775      int i;
 776      for (i = 0; i < 3; i++) {
 777        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 928) real");
 778        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 928) imag");
 779      };
 780    };
 781   };
 782
 783
 784   {
 785    int order = 101;
 786    int trans = 111;
 787    int uplo = 122;
 788    int diag = 132;
 789    int N = 3;
 790    int K = 1;
 791    int lda = 3;
 792    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
 793    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
 794    int incX = -1;
 795    float x_expected[] = { -0.39057f, 0.113296f, 0.388863f, 0.131011f, -0.236f, 0.605f };
 796    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 797    {
 798      int i;
 799      for (i = 0; i < 3; i++) {
 800        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 929) real");
 801        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 929) imag");
 802      };
 803    };
 804   };
 805
 806
 807   {
 808    int order = 102;
 809    int trans = 111;
 810    int uplo = 121;
 811    int diag = 131;
 812    int N = 3;
 813    int K = 1;
 814    int lda = 3;
 815    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
 816    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
 817    int incX = -1;
 818    float x_expected[] = { -0.169119f, 0.443509f, 0.159816f, 0.139696f, -0.180955f, -0.835292f };
 819    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 820    {
 821      int i;
 822      for (i = 0; i < 3; i++) {
 823        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 930) real");
 824        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 930) imag");
 825      };
 826    };
 827   };
 828
 829
 830   {
 831    int order = 102;
 832    int trans = 111;
 833    int uplo = 121;
 834    int diag = 132;
 835    int N = 3;
 836    int K = 1;
 837    int lda = 3;
 838    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
 839    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
 840    int incX = -1;
 841    float x_expected[] = { -0.46f, 0.069f, 0.194886f, -0.054704f, -0.191297f, 0.545731f };
 842    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 843    {
 844      int i;
 845      for (i = 0; i < 3; i++) {
 846        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 931) real");
 847        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 931) imag");
 848      };
 849    };
 850   };
 851
 852
 853   {
 854    int order = 102;
 855    int trans = 111;
 856    int uplo = 122;
 857    int diag = 131;
 858    int N = 3;
 859    int K = 1;
 860    int lda = 3;
 861    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
 862    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
 863    int incX = -1;
 864    float x_expected[] = { 0.159816f, 0.139696f, -0.180955f, -0.835292f, 0.077786f, 0.60472f };
 865    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 866    {
 867      int i;
 868      for (i = 0; i < 3; i++) {
 869        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 932) real");
 870        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 932) imag");
 871      };
 872    };
 873   };
 874
 875
 876   {
 877    int order = 102;
 878    int trans = 111;
 879    int uplo = 122;
 880    int diag = 132;
 881    int N = 3;
 882    int K = 1;
 883    int lda = 3;
 884    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
 885    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
 886    int incX = -1;
 887    float x_expected[] = { -0.18707f, 0.2604f, 0.082342f, -0.779023f, -0.236f, 0.605f };
 888    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 889    {
 890      int i;
 891      for (i = 0; i < 3; i++) {
 892        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 933) real");
 893        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 933) imag");
 894      };
 895    };
 896   };
 897
 898
 899   {
 900    int order = 101;
 901    int trans = 112;
 902    int uplo = 121;
 903    int diag = 131;
 904    int N = 3;
 905    int K = 1;
 906    int lda = 3;
 907    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
 908    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
 909    int incX = -1;
 910    float x_expected[] = { 0.647885f, 0.621535f, -0.104407f, 0.05309f, 0.732704f, 0.055982f };
 911    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 912    {
 913      int i;
 914      for (i = 0; i < 3; i++) {
 915        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 934) real");
 916        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 934) imag");
 917      };
 918    };
 919   };
 920
 921
 922   {
 923    int order = 101;
 924    int trans = 112;
 925    int uplo = 121;
 926    int diag = 132;
 927    int N = 3;
 928    int K = 1;
 929    int lda = 3;
 930    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
 931    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
 932    int incX = -1;
 933    float x_expected[] = { 1.2955f, 0.190774f, -0.247934f, 0.982616f, -0.894f, -0.116f };
 934    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 935    {
 936      int i;
 937      for (i = 0; i < 3; i++) {
 938        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 935) real");
 939        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 935) imag");
 940      };
 941    };
 942   };
 943
 944
 945   {
 946    int order = 101;
 947    int trans = 112;
 948    int uplo = 122;
 949    int diag = 131;
 950    int N = 3;
 951    int K = 1;
 952    int lda = 3;
 953    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
 954    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
 955    int incX = -1;
 956    float x_expected[] = { 0.096482f, -0.071661f, 0.647885f, 0.621535f, -0.104407f, 0.05309f };
 957    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 958    {
 959      int i;
 960      for (i = 0; i < 3; i++) {
 961        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 936) real");
 962        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 936) imag");
 963      };
 964    };
 965   };
 966
 967
 968   {
 969    int order = 101;
 970    int trans = 112;
 971    int uplo = 122;
 972    int diag = 132;
 973    int N = 3;
 974    int K = 1;
 975    int lda = 3;
 976    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
 977    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
 978    int incX = -1;
 979    float x_expected[] = { 0.411f, -0.308f, -1.14861f, 0.933761f, -1.66247f, -0.234526f };
 980    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
 981    {
 982      int i;
 983      for (i = 0; i < 3; i++) {
 984        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 937) real");
 985        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 937) imag");
 986      };
 987    };
 988   };
 989
 990
 991   {
 992    int order = 102;
 993    int trans = 112;
 994    int uplo = 121;
 995    int diag = 131;
 996    int N = 3;
 997    int K = 1;
 998    int lda = 3;
 999    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
1000    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
1001    int incX = -1;
1002    float x_expected[] = { 0.632361f, -0.409373f, 0.578489f, 0.012724f, 0.664066f, 0.171616f };
1003    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1004    {
1005      int i;
1006      for (i = 0; i < 3; i++) {
1007        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 938) real");
1008        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 938) imag");
1009      };
1010    };
1011   };
1012
1013
1014   {
1015    int order = 102;
1016    int trans = 112;
1017    int uplo = 121;
1018    int diag = 132;
1019    int N = 3;
1020    int K = 1;
1021    int lda = 3;
1022    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
1023    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
1024    int incX = -1;
1025    float x_expected[] = { 0.946879f, -0.645712f, -1.21801f, 0.32495f, -0.894f, -0.116f };
1026    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1027    {
1028      int i;
1029      for (i = 0; i < 3; i++) {
1030        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 939) real");
1031        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 939) imag");
1032      };
1033    };
1034   };
1035
1036
1037   {
1038    int order = 102;
1039    int trans = 112;
1040    int uplo = 122;
1041    int diag = 131;
1042    int N = 3;
1043    int K = 1;
1044    int lda = 3;
1045    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
1046    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
1047    int incX = -1;
1048    float x_expected[] = { -0.236612f, 0.122761f, -1.12184f, -0.358823f, 1.4975f, -0.470595f };
1049    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1050    {
1051      int i;
1052      for (i = 0; i < 3; i++) {
1053        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 940) real");
1054        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 940) imag");
1055      };
1056    };
1057   };
1058
1059
1060   {
1061    int order = 102;
1062    int trans = 112;
1063    int uplo = 122;
1064    int diag = 132;
1065    int N = 3;
1066    int K = 1;
1067    int lda = 3;
1068    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
1069    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
1070    int incX = -1;
1071    float x_expected[] = { 0.411f, -0.308f, -1.26537f, 0.570703f, -0.129206f, -0.642577f };
1072    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1073    {
1074      int i;
1075      for (i = 0; i < 3; i++) {
1076        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 941) real");
1077        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 941) imag");
1078      };
1079    };
1080   };
1081
1082
1083   {
1084    int order = 101;
1085    int trans = 113;
1086    int uplo = 121;
1087    int diag = 131;
1088    int N = 3;
1089    int K = 1;
1090    int lda = 3;
1091    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
1092    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
1093    int incX = -1;
1094    float x_expected[] = { 0.413357f, 0.178267f, -0.114618f, -1.35595f, -0.513288f, 0.611332f };
1095    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1096    {
1097      int i;
1098      for (i = 0; i < 3; i++) {
1099        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 942) real");
1100        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 942) imag");
1101      };
1102    };
1103   };
1104
1105
1106   {
1107    int order = 101;
1108    int trans = 113;
1109    int uplo = 121;
1110    int diag = 132;
1111    int N = 3;
1112    int K = 1;
1113    int lda = 3;
1114    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
1115    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
1116    int incX = -1;
1117    float x_expected[] = { 0.368428f, 0.071217f, -0.954366f, -0.390486f, 0.694f, -0.954f };
1118    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1119    {
1120      int i;
1121      for (i = 0; i < 3; i++) {
1122        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 943) real");
1123        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 943) imag");
1124      };
1125    };
1126   };
1127
1128
1129   {
1130    int order = 101;
1131    int trans = 113;
1132    int uplo = 122;
1133    int diag = 131;
1134    int N = 3;
1135    int K = 1;
1136    int lda = 3;
1137    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
1138    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
1139    int incX = -1;
1140    float x_expected[] = { -0.084786f, -0.059464f, 0.413357f, 0.178267f, -0.114618f, -1.35595f };
1141    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1142    {
1143      int i;
1144      for (i = 0; i < 3; i++) {
1145        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 944) real");
1146        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 944) imag");
1147      };
1148    };
1149   };
1150
1151
1152   {
1153    int order = 101;
1154    int trans = 113;
1155    int uplo = 122;
1156    int diag = 132;
1157    int N = 3;
1158    int K = 1;
1159    int lda = 3;
1160    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
1161    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
1162    int incX = -1;
1163    float x_expected[] = { 0.065f, -0.082f, -0.636071f, 0.80005f, 0.787748f, -1.14446f };
1164    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1165    {
1166      int i;
1167      for (i = 0; i < 3; i++) {
1168        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 945) real");
1169        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 945) imag");
1170      };
1171    };
1172   };
1173
1174
1175   {
1176    int order = 102;
1177    int trans = 113;
1178    int uplo = 121;
1179    int diag = 131;
1180    int N = 3;
1181    int K = 1;
1182    int lda = 3;
1183    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
1184    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
1185    int incX = -1;
1186    float x_expected[] = { -1.18498f, -0.424201f, 0.230196f, 0.374209f, -0.208366f, -1.16549f };
1187    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1188    {
1189      int i;
1190      for (i = 0; i < 3; i++) {
1191        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 946) real");
1192        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 946) imag");
1193      };
1194    };
1195   };
1196
1197
1198   {
1199    int order = 102;
1200    int trans = 113;
1201    int uplo = 121;
1202    int diag = 132;
1203    int N = 3;
1204    int K = 1;
1205    int lda = 3;
1206    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
1207    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
1208    int incX = -1;
1209    float x_expected[] = { -1.03519f, -0.446737f, -0.819232f, 0.995992f, 0.694f, -0.954f };
1210    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1211    {
1212      int i;
1213      for (i = 0; i < 3; i++) {
1214        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 947) real");
1215        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 947) imag");
1216      };
1217    };
1218   };
1219
1220
1221   {
1222    int order = 102;
1223    int trans = 113;
1224    int uplo = 122;
1225    int diag = 131;
1226    int N = 3;
1227    int K = 1;
1228    int lda = 3;
1229    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
1230    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
1231    int incX = -1;
1232    float x_expected[] = { 0.109929f, 0.02505f, 0.062939f, -0.202464f, -0.470658f, 1.69006f };
1233    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1234    {
1235      int i;
1236      for (i = 0; i < 3; i++) {
1237        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 948) real");
1238        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 948) imag");
1239      };
1240    };
1241   };
1242
1243
1244   {
1245    int order = 102;
1246    int trans = 113;
1247    int uplo = 122;
1248    int diag = 132;
1249    int N = 3;
1250    int K = 1;
1251    int lda = 3;
1252    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
1253    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
1254    int incX = -1;
1255    float x_expected[] = { 0.065f, -0.082f, -0.776809f, 0.762996f, 0.73663f, 0.124729f };
1256    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1257    {
1258      int i;
1259      for (i = 0; i < 3; i++) {
1260        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 949) real");
1261        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 949) imag");
1262      };
1263    };
1264   };
1265
1266
1267   {
1268    int order = 101;
1269    int trans = 111;
1270    int uplo = 121;
1271    int diag = 131;
1272    int N = 3;
1273    int K = 1;
1274    int lda = 3;
1275    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
1276    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
1277    int incX = -1;
1278    double x_expected[] = { -0.010019, -0.1678, -0.042017, -1.112094, 0.010004, -0.480427 };
1279    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1280    {
1281      int i;
1282      for (i = 0; i < 3; i++) {
1283        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 950) real");
1284        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 950) imag");
1285      };
1286    };
1287   };
1288
1289
1290   {
1291    int order = 101;
1292    int trans = 111;
1293    int uplo = 121;
1294    int diag = 132;
1295    int N = 3;
1296    int K = 1;
1297    int lda = 3;
1298    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
1299    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
1300    int incX = -1;
1301    double x_expected[] = { 0.064, 0.169, -0.80842, -0.715637, -0.829924, -0.212971 };
1302    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1303    {
1304      int i;
1305      for (i = 0; i < 3; i++) {
1306        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 951) real");
1307        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 951) imag");
1308      };
1309    };
1310   };
1311
1312
1313   {
1314    int order = 101;
1315    int trans = 111;
1316    int uplo = 122;
1317    int diag = 131;
1318    int N = 3;
1319    int K = 1;
1320    int lda = 3;
1321    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
1322    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
1323    int incX = -1;
1324    double x_expected[] = { 0.634014, 0.796937, -0.585538, -0.895375, -0.125887, 0.010019 };
1325    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1326    {
1327      int i;
1328      for (i = 0; i < 3; i++) {
1329        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 952) real");
1330        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 952) imag");
1331      };
1332    };
1333   };
1334
1335
1336   {
1337    int order = 101;
1338    int trans = 111;
1339    int uplo = 122;
1340    int diag = 132;
1341    int N = 3;
1342    int K = 1;
1343    int lda = 3;
1344    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
1345    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
1346    int incX = -1;
1347    double x_expected[] = { 0.567497, 1.085122, -1.217792, -1.322566, -0.641, -0.103 };
1348    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1349    {
1350      int i;
1351      for (i = 0; i < 3; i++) {
1352        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 953) real");
1353        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 953) imag");
1354      };
1355    };
1356   };
1357
1358
1359   {
1360    int order = 102;
1361    int trans = 111;
1362    int uplo = 121;
1363    int diag = 131;
1364    int N = 3;
1365    int K = 1;
1366    int lda = 3;
1367    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
1368    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
1369    int incX = -1;
1370    double x_expected[] = { 0.130517, -0.119185, -0.187765, -0.519609, -0.169484, -1.165438 };
1371    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1372    {
1373      int i;
1374      for (i = 0; i < 3; i++) {
1375        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 954) real");
1376        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 954) imag");
1377      };
1378    };
1379   };
1380
1381
1382   {
1383    int order = 102;
1384    int trans = 111;
1385    int uplo = 121;
1386    int diag = 132;
1387    int N = 3;
1388    int K = 1;
1389    int lda = 3;
1390    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
1391    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
1392    int incX = -1;
1393    double x_expected[] = { 0.064, 0.169, -0.820019, -0.9468, -0.684597, -1.278457 };
1394    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1395    {
1396      int i;
1397      for (i = 0; i < 3; i++) {
1398        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 955) real");
1399        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 955) imag");
1400      };
1401    };
1402   };
1403
1404
1405   {
1406    int order = 102;
1407    int trans = 111;
1408    int uplo = 122;
1409    int diag = 131;
1410    int N = 3;
1411    int K = 1;
1412    int lda = 3;
1413    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
1414    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
1415    int incX = -1;
1416    double x_expected[] = { -0.187765, -0.519609, -0.169484, -1.165438, 0.198928, -0.370456 };
1417    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1418    {
1419      int i;
1420      for (i = 0; i < 3; i++) {
1421        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 956) real");
1422        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 956) imag");
1423      };
1424    };
1425   };
1426
1427
1428   {
1429    int order = 102;
1430    int trans = 111;
1431    int uplo = 122;
1432    int diag = 132;
1433    int N = 3;
1434    int K = 1;
1435    int lda = 3;
1436    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
1437    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
1438    int incX = -1;
1439    double x_expected[] = { -0.113746, -0.182809, -0.935887, -0.768981, -0.641, -0.103 };
1440    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1441    {
1442      int i;
1443      for (i = 0; i < 3; i++) {
1444        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 957) real");
1445        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 957) imag");
1446      };
1447    };
1448   };
1449
1450
1451   {
1452    int order = 101;
1453    int trans = 112;
1454    int uplo = 121;
1455    int diag = 131;
1456    int N = 3;
1457    int K = 1;
1458    int lda = 3;
1459    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
1460    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
1461    int incX = -1;
1462    double x_expected[] = { -0.436746, 0.963714, -1.087615, -0.018695, 0.30063, 0.12958 };
1463    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1464    {
1465      int i;
1466      for (i = 0; i < 3; i++) {
1467        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 958) real");
1468        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 958) imag");
1469      };
1470    };
1471   };
1472
1473
1474   {
1475    int order = 101;
1476    int trans = 112;
1477    int uplo = 121;
1478    int diag = 132;
1479    int N = 3;
1480    int K = 1;
1481    int lda = 3;
1482    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
1483    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
1484    int incX = -1;
1485    double x_expected[] = { 0.895682, 1.407174, 0.2408, -0.14282, -0.649, 0.188 };
1486    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1487    {
1488      int i;
1489      for (i = 0; i < 3; i++) {
1490        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 959) real");
1491        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 959) imag");
1492      };
1493    };
1494   };
1495
1496
1497   {
1498    int order = 101;
1499    int trans = 112;
1500    int uplo = 122;
1501    int diag = 131;
1502    int N = 3;
1503    int K = 1;
1504    int lda = 3;
1505    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
1506    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
1507    int incX = -1;
1508    double x_expected[] = { 0.785744, -0.3966, -0.436746, 0.963714, -1.087615, -0.018695 };
1509    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1510    {
1511      int i;
1512      for (i = 0; i < 3; i++) {
1513        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 960) real");
1514        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 960) imag");
1515      };
1516    };
1517   };
1518
1519
1520   {
1521    int order = 101;
1522    int trans = 112;
1523    int uplo = 122;
1524    int diag = 132;
1525    int N = 3;
1526    int K = 1;
1527    int lda = 3;
1528    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
1529    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
1530    int incX = -1;
1531    double x_expected[] = { 0.884, 0.636, 0.472572, 0.47454, -1.056415, 0.594125 };
1532    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1533    {
1534      int i;
1535      for (i = 0; i < 3; i++) {
1536        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 961) real");
1537        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 961) imag");
1538      };
1539    };
1540   };
1541
1542
1543   {
1544    int order = 102;
1545    int trans = 112;
1546    int uplo = 121;
1547    int diag = 131;
1548    int N = 3;
1549    int K = 1;
1550    int lda = 3;
1551    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
1552    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
1553    int incX = -1;
1554    double x_expected[] = { 0.464705, -0.108078, 0.094975, 0.376323, -0.6802, -0.42482 };
1555    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1556    {
1557      int i;
1558      for (i = 0; i < 3; i++) {
1559        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 962) real");
1560        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 962) imag");
1561      };
1562    };
1563   };
1564
1565
1566   {
1567    int order = 102;
1568    int trans = 112;
1569    int uplo = 121;
1570    int diag = 132;
1571    int N = 3;
1572    int K = 1;
1573    int lda = 3;
1574    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
1575    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
1576    int incX = -1;
1577    double x_expected[] = { 0.562961, 0.924522, 1.004293, -0.112851, -0.649, 0.188 };
1578    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1579    {
1580      int i;
1581      for (i = 0; i < 3; i++) {
1582        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 963) real");
1583        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 963) imag");
1584      };
1585    };
1586   };
1587
1588
1589   {
1590    int order = 102;
1591    int trans = 112;
1592    int uplo = 122;
1593    int diag = 131;
1594    int N = 3;
1595    int K = 1;
1596    int lda = 3;
1597    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
1598    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
1599    int incX = -1;
1600    double x_expected[] = { -0.448428, 0.19254, -0.674583, 1.236189, 0.780774, 1.167088 };
1601    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1602    {
1603      int i;
1604      for (i = 0; i < 3; i++) {
1605        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 964) real");
1606        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 964) imag");
1607      };
1608    };
1609   };
1610
1611
1612   {
1613    int order = 102;
1614    int trans = 112;
1615    int uplo = 122;
1616    int diag = 132;
1617    int N = 3;
1618    int K = 1;
1619    int lda = 3;
1620    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
1621    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
1622    int incX = -1;
1623    double x_expected[] = { 0.884, 0.636, 0.653832, 1.112064, -0.168856, 1.225508 };
1624    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1625    {
1626      int i;
1627      for (i = 0; i < 3; i++) {
1628        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 965) real");
1629        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 965) imag");
1630      };
1631    };
1632   };
1633
1634
1635   {
1636    int order = 101;
1637    int trans = 113;
1638    int uplo = 121;
1639    int diag = 131;
1640    int N = 3;
1641    int K = 1;
1642    int lda = 3;
1643    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
1644    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
1645    int incX = -1;
1646    double x_expected[] = { -0.59515, 0.077106, -0.27658, -0.637356, 0.407252, -0.308844 };
1647    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1648    {
1649      int i;
1650      for (i = 0; i < 3; i++) {
1651        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 966) real");
1652        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 966) imag");
1653      };
1654    };
1655   };
1656
1657
1658   {
1659    int order = 101;
1660    int trans = 113;
1661    int uplo = 121;
1662    int diag = 132;
1663    int N = 3;
1664    int K = 1;
1665    int lda = 3;
1666    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
1667    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
1668    int incX = -1;
1669    double x_expected[] = { -1.46131, 0.537642, 0.624614, 0.762252, 0.326, 0.428 };
1670    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1671    {
1672      int i;
1673      for (i = 0; i < 3; i++) {
1674        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 967) real");
1675        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 967) imag");
1676      };
1677    };
1678   };
1679
1680
1681   {
1682    int order = 101;
1683    int trans = 113;
1684    int uplo = 122;
1685    int diag = 131;
1686    int N = 3;
1687    int K = 1;
1688    int lda = 3;
1689    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
1690    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
1691    int incX = -1;
1692    double x_expected[] = { -0.536274, 0.421806, -0.59515, 0.077106, -0.27658, -0.637356 };
1693    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1694    {
1695      int i;
1696      for (i = 0; i < 3; i++) {
1697        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 968) real");
1698        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 968) imag");
1699      };
1700    };
1701   };
1702
1703
1704   {
1705    int order = 101;
1706    int trans = 113;
1707    int uplo = 122;
1708    int diag = 132;
1709    int N = 3;
1710    int K = 1;
1711    int lda = 3;
1712    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
1713    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
1714    int incX = -1;
1715    double x_expected[] = { -0.591, -0.084, 0.98216, 0.400464, 0.131806, -0.026608 };
1716    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1717    {
1718      int i;
1719      for (i = 0; i < 3; i++) {
1720        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 969) real");
1721        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 969) imag");
1722      };
1723    };
1724   };
1725
1726
1727   {
1728    int order = 102;
1729    int trans = 113;
1730    int uplo = 121;
1731    int diag = 131;
1732    int N = 3;
1733    int K = 1;
1734    int lda = 3;
1735    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
1736    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
1737    int incX = -1;
1738    double x_expected[] = { -1.68293, 0.796222, -0.96062, 0.415172, -0.082386, -0.182748 };
1739    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1740    {
1741      int i;
1742      for (i = 0; i < 3; i++) {
1743        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 970) real");
1744        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 970) imag");
1745      };
1746    };
1747   };
1748
1749
1750   {
1751    int order = 102;
1752    int trans = 113;
1753    int uplo = 121;
1754    int diag = 132;
1755    int N = 3;
1756    int K = 1;
1757    int lda = 3;
1758    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
1759    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
1760    int incX = -1;
1761    double x_expected[] = { -1.737656, 0.290416, 0.61669, 0.73853, 0.326, 0.428 };
1762    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1763    {
1764      int i;
1765      for (i = 0; i < 3; i++) {
1766        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 971) real");
1767        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 971) imag");
1768      };
1769    };
1770   };
1771
1772
1773   {
1774    int order = 102;
1775    int trans = 113;
1776    int uplo = 122;
1777    int diag = 131;
1778    int N = 3;
1779    int K = 1;
1780    int lda = 3;
1781    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
1782    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
1783    int incX = -1;
1784    double x_expected[] = { 0.27516, -0.544536, -0.10627, -0.988374, 0.229991, -0.711267 };
1785    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1786    {
1787      int i;
1788      for (i = 0; i < 3; i++) {
1789        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 972) real");
1790        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 972) imag");
1791      };
1792    };
1793   };
1794
1795
1796   {
1797    int order = 102;
1798    int trans = 113;
1799    int uplo = 122;
1800    int diag = 132;
1801    int N = 3;
1802    int K = 1;
1803    int lda = 3;
1804    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
1805    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
1806    int incX = -1;
1807    double x_expected[] = { -0.591, -0.084, 0.794924, 0.411234, 0.148739, 0.025577 };
1808    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1809    {
1810      int i;
1811      for (i = 0; i < 3; i++) {
1812        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 973) real");
1813        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 973) imag");
1814      };
1815    };
1816   };
1817
1818
1819 }