htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_tpsv.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_tpsv (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 = 2;
  17    float A[] = { -0.381f, 0.53f, 0.451f };
  18    float X[] = { 0.144f, 0.032f };
  19    int incX = -1;
  20    float x_expected[] = { 0.31929f, 0.360168f };
  21    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  22    {
  23      int i;
  24      for (i = 0; i < 2; i++) {
  25        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1310)");
  26      }
  27    };
  28   };
  29
  30
  31   {
  32    int order = 101;
  33    int trans = 111;
  34    int uplo = 121;
  35    int diag = 132;
  36    int N = 2;
  37    float A[] = { -0.381f, 0.53f, 0.451f };
  38    float X[] = { 0.144f, 0.032f };
  39    int incX = -1;
  40    float x_expected[] = { 0.144f, -0.04432f };
  41    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  42    {
  43      int i;
  44      for (i = 0; i < 2; i++) {
  45        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1311)");
  46      }
  47    };
  48   };
  49
  50
  51   {
  52    int order = 101;
  53    int trans = 111;
  54    int uplo = 122;
  55    int diag = 131;
  56    int N = 2;
  57    float A[] = { -0.381f, 0.53f, 0.451f };
  58    float X[] = { 0.144f, 0.032f };
  59    int incX = -1;
  60    float x_expected[] = { 0.417992f, -0.0839895f };
  61    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  62    {
  63      int i;
  64      for (i = 0; i < 2; i++) {
  65        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1312)");
  66      }
  67    };
  68   };
  69
  70
  71   {
  72    int order = 101;
  73    int trans = 111;
  74    int uplo = 122;
  75    int diag = 132;
  76    int N = 2;
  77    float A[] = { -0.381f, 0.53f, 0.451f };
  78    float X[] = { 0.144f, 0.032f };
  79    int incX = -1;
  80    float x_expected[] = { 0.12704f, 0.032f };
  81    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  82    {
  83      int i;
  84      for (i = 0; i < 2; i++) {
  85        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1313)");
  86      }
  87    };
  88   };
  89
  90
  91   {
  92    int order = 102;
  93    int trans = 111;
  94    int uplo = 121;
  95    int diag = 131;
  96    int N = 2;
  97    float A[] = { -0.381f, 0.53f, 0.451f };
  98    float X[] = { 0.144f, 0.032f };
  99    int incX = -1;
 100    float x_expected[] = { 0.31929f, 0.360168f };
 101    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 102    {
 103      int i;
 104      for (i = 0; i < 2; i++) {
 105        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1314)");
 106      }
 107    };
 108   };
 109
 110
 111   {
 112    int order = 102;
 113    int trans = 111;
 114    int uplo = 121;
 115    int diag = 132;
 116    int N = 2;
 117    float A[] = { -0.381f, 0.53f, 0.451f };
 118    float X[] = { 0.144f, 0.032f };
 119    int incX = -1;
 120    float x_expected[] = { 0.144f, -0.04432f };
 121    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 122    {
 123      int i;
 124      for (i = 0; i < 2; i++) {
 125        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1315)");
 126      }
 127    };
 128   };
 129
 130
 131   {
 132    int order = 102;
 133    int trans = 111;
 134    int uplo = 122;
 135    int diag = 131;
 136    int N = 2;
 137    float A[] = { -0.381f, 0.53f, 0.451f };
 138    float X[] = { 0.144f, 0.032f };
 139    int incX = -1;
 140    float x_expected[] = { 0.417992f, -0.0839895f };
 141    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 142    {
 143      int i;
 144      for (i = 0; i < 2; i++) {
 145        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1316)");
 146      }
 147    };
 148   };
 149
 150
 151   {
 152    int order = 102;
 153    int trans = 111;
 154    int uplo = 122;
 155    int diag = 132;
 156    int N = 2;
 157    float A[] = { -0.381f, 0.53f, 0.451f };
 158    float X[] = { 0.144f, 0.032f };
 159    int incX = -1;
 160    float x_expected[] = { 0.12704f, 0.032f };
 161    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 162    {
 163      int i;
 164      for (i = 0; i < 2; i++) {
 165        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1317)");
 166      }
 167    };
 168   };
 169
 170
 171   {
 172    int order = 101;
 173    int trans = 112;
 174    int uplo = 121;
 175    int diag = 131;
 176    int N = 2;
 177    float A[] = { -0.381f, 0.53f, 0.451f };
 178    float X[] = { 0.144f, 0.032f };
 179    int incX = -1;
 180    float x_expected[] = { 0.417992f, -0.0839895f };
 181    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 182    {
 183      int i;
 184      for (i = 0; i < 2; i++) {
 185        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1318)");
 186      }
 187    };
 188   };
 189
 190
 191   {
 192    int order = 101;
 193    int trans = 112;
 194    int uplo = 121;
 195    int diag = 132;
 196    int N = 2;
 197    float A[] = { -0.381f, 0.53f, 0.451f };
 198    float X[] = { 0.144f, 0.032f };
 199    int incX = -1;
 200    float x_expected[] = { 0.12704f, 0.032f };
 201    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 202    {
 203      int i;
 204      for (i = 0; i < 2; i++) {
 205        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1319)");
 206      }
 207    };
 208   };
 209
 210
 211   {
 212    int order = 101;
 213    int trans = 112;
 214    int uplo = 122;
 215    int diag = 131;
 216    int N = 2;
 217    float A[] = { -0.381f, 0.53f, 0.451f };
 218    float X[] = { 0.144f, 0.032f };
 219    int incX = -1;
 220    float x_expected[] = { 0.31929f, 0.360168f };
 221    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 222    {
 223      int i;
 224      for (i = 0; i < 2; i++) {
 225        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1320)");
 226      }
 227    };
 228   };
 229
 230
 231   {
 232    int order = 101;
 233    int trans = 112;
 234    int uplo = 122;
 235    int diag = 132;
 236    int N = 2;
 237    float A[] = { -0.381f, 0.53f, 0.451f };
 238    float X[] = { 0.144f, 0.032f };
 239    int incX = -1;
 240    float x_expected[] = { 0.144f, -0.04432f };
 241    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 242    {
 243      int i;
 244      for (i = 0; i < 2; i++) {
 245        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1321)");
 246      }
 247    };
 248   };
 249
 250
 251   {
 252    int order = 102;
 253    int trans = 112;
 254    int uplo = 121;
 255    int diag = 131;
 256    int N = 2;
 257    float A[] = { -0.381f, 0.53f, 0.451f };
 258    float X[] = { 0.144f, 0.032f };
 259    int incX = -1;
 260    float x_expected[] = { 0.417992f, -0.0839895f };
 261    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 262    {
 263      int i;
 264      for (i = 0; i < 2; i++) {
 265        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1322)");
 266      }
 267    };
 268   };
 269
 270
 271   {
 272    int order = 102;
 273    int trans = 112;
 274    int uplo = 121;
 275    int diag = 132;
 276    int N = 2;
 277    float A[] = { -0.381f, 0.53f, 0.451f };
 278    float X[] = { 0.144f, 0.032f };
 279    int incX = -1;
 280    float x_expected[] = { 0.12704f, 0.032f };
 281    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 282    {
 283      int i;
 284      for (i = 0; i < 2; i++) {
 285        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1323)");
 286      }
 287    };
 288   };
 289
 290
 291   {
 292    int order = 102;
 293    int trans = 112;
 294    int uplo = 122;
 295    int diag = 131;
 296    int N = 2;
 297    float A[] = { -0.381f, 0.53f, 0.451f };
 298    float X[] = { 0.144f, 0.032f };
 299    int incX = -1;
 300    float x_expected[] = { 0.31929f, 0.360168f };
 301    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 302    {
 303      int i;
 304      for (i = 0; i < 2; i++) {
 305        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1324)");
 306      }
 307    };
 308   };
 309
 310
 311   {
 312    int order = 102;
 313    int trans = 112;
 314    int uplo = 122;
 315    int diag = 132;
 316    int N = 2;
 317    float A[] = { -0.381f, 0.53f, 0.451f };
 318    float X[] = { 0.144f, 0.032f };
 319    int incX = -1;
 320    float x_expected[] = { 0.144f, -0.04432f };
 321    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
 322    {
 323      int i;
 324      for (i = 0; i < 2; i++) {
 325        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1325)");
 326      }
 327    };
 328   };
 329
 330
 331   {
 332    int order = 101;
 333    int trans = 111;
 334    int uplo = 121;
 335    int diag = 131;
 336    int N = 2;
 337    double A[] = { 0.355, -0.698, -0.42 };
 338    double X[] = { -0.702, -0.661 };
 339    int incX = -1;
 340    double x_expected[] = { 1.67142857143, 1.42438631791 };
 341    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 342    {
 343      int i;
 344      for (i = 0; i < 2; i++) {
 345        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1326)");
 346      }
 347    };
 348   };
 349
 350
 351   {
 352    int order = 101;
 353    int trans = 111;
 354    int uplo = 121;
 355    int diag = 132;
 356    int N = 2;
 357    double A[] = { 0.355, -0.698, -0.42 };
 358    double X[] = { -0.702, -0.661 };
 359    int incX = -1;
 360    double x_expected[] = { -0.702, -1.150996 };
 361    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 362    {
 363      int i;
 364      for (i = 0; i < 2; i++) {
 365        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1327)");
 366      }
 367    };
 368   };
 369
 370
 371   {
 372    int order = 101;
 373    int trans = 111;
 374    int uplo = 122;
 375    int diag = 131;
 376    int N = 2;
 377    double A[] = { 0.355, -0.698, -0.42 };
 378    double X[] = { -0.702, -0.661 };
 379    int incX = -1;
 380    double x_expected[] = { 4.76584842388, -1.86197183099 };
 381    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 382    {
 383      int i;
 384      for (i = 0; i < 2; i++) {
 385        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1328)");
 386      }
 387    };
 388   };
 389
 390
 391   {
 392    int order = 101;
 393    int trans = 111;
 394    int uplo = 122;
 395    int diag = 132;
 396    int N = 2;
 397    double A[] = { 0.355, -0.698, -0.42 };
 398    double X[] = { -0.702, -0.661 };
 399    int incX = -1;
 400    double x_expected[] = { -1.163378, -0.661 };
 401    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 402    {
 403      int i;
 404      for (i = 0; i < 2; i++) {
 405        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1329)");
 406      }
 407    };
 408   };
 409
 410
 411   {
 412    int order = 102;
 413    int trans = 111;
 414    int uplo = 121;
 415    int diag = 131;
 416    int N = 2;
 417    double A[] = { 0.355, -0.698, -0.42 };
 418    double X[] = { -0.702, -0.661 };
 419    int incX = -1;
 420    double x_expected[] = { 1.67142857143, 1.42438631791 };
 421    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 422    {
 423      int i;
 424      for (i = 0; i < 2; i++) {
 425        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1330)");
 426      }
 427    };
 428   };
 429
 430
 431   {
 432    int order = 102;
 433    int trans = 111;
 434    int uplo = 121;
 435    int diag = 132;
 436    int N = 2;
 437    double A[] = { 0.355, -0.698, -0.42 };
 438    double X[] = { -0.702, -0.661 };
 439    int incX = -1;
 440    double x_expected[] = { -0.702, -1.150996 };
 441    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 442    {
 443      int i;
 444      for (i = 0; i < 2; i++) {
 445        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1331)");
 446      }
 447    };
 448   };
 449
 450
 451   {
 452    int order = 102;
 453    int trans = 111;
 454    int uplo = 122;
 455    int diag = 131;
 456    int N = 2;
 457    double A[] = { 0.355, -0.698, -0.42 };
 458    double X[] = { -0.702, -0.661 };
 459    int incX = -1;
 460    double x_expected[] = { 4.76584842388, -1.86197183099 };
 461    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 462    {
 463      int i;
 464      for (i = 0; i < 2; i++) {
 465        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1332)");
 466      }
 467    };
 468   };
 469
 470
 471   {
 472    int order = 102;
 473    int trans = 111;
 474    int uplo = 122;
 475    int diag = 132;
 476    int N = 2;
 477    double A[] = { 0.355, -0.698, -0.42 };
 478    double X[] = { -0.702, -0.661 };
 479    int incX = -1;
 480    double x_expected[] = { -1.163378, -0.661 };
 481    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 482    {
 483      int i;
 484      for (i = 0; i < 2; i++) {
 485        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1333)");
 486      }
 487    };
 488   };
 489
 490
 491   {
 492    int order = 101;
 493    int trans = 112;
 494    int uplo = 121;
 495    int diag = 131;
 496    int N = 2;
 497    double A[] = { 0.355, -0.698, -0.42 };
 498    double X[] = { -0.702, -0.661 };
 499    int incX = -1;
 500    double x_expected[] = { 4.76584842388, -1.86197183099 };
 501    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 502    {
 503      int i;
 504      for (i = 0; i < 2; i++) {
 505        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1334)");
 506      }
 507    };
 508   };
 509
 510
 511   {
 512    int order = 101;
 513    int trans = 112;
 514    int uplo = 121;
 515    int diag = 132;
 516    int N = 2;
 517    double A[] = { 0.355, -0.698, -0.42 };
 518    double X[] = { -0.702, -0.661 };
 519    int incX = -1;
 520    double x_expected[] = { -1.163378, -0.661 };
 521    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 522    {
 523      int i;
 524      for (i = 0; i < 2; i++) {
 525        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1335)");
 526      }
 527    };
 528   };
 529
 530
 531   {
 532    int order = 101;
 533    int trans = 112;
 534    int uplo = 122;
 535    int diag = 131;
 536    int N = 2;
 537    double A[] = { 0.355, -0.698, -0.42 };
 538    double X[] = { -0.702, -0.661 };
 539    int incX = -1;
 540    double x_expected[] = { 1.67142857143, 1.42438631791 };
 541    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 542    {
 543      int i;
 544      for (i = 0; i < 2; i++) {
 545        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1336)");
 546      }
 547    };
 548   };
 549
 550
 551   {
 552    int order = 101;
 553    int trans = 112;
 554    int uplo = 122;
 555    int diag = 132;
 556    int N = 2;
 557    double A[] = { 0.355, -0.698, -0.42 };
 558    double X[] = { -0.702, -0.661 };
 559    int incX = -1;
 560    double x_expected[] = { -0.702, -1.150996 };
 561    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 562    {
 563      int i;
 564      for (i = 0; i < 2; i++) {
 565        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1337)");
 566      }
 567    };
 568   };
 569
 570
 571   {
 572    int order = 102;
 573    int trans = 112;
 574    int uplo = 121;
 575    int diag = 131;
 576    int N = 2;
 577    double A[] = { 0.355, -0.698, -0.42 };
 578    double X[] = { -0.702, -0.661 };
 579    int incX = -1;
 580    double x_expected[] = { 4.76584842388, -1.86197183099 };
 581    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 582    {
 583      int i;
 584      for (i = 0; i < 2; i++) {
 585        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1338)");
 586      }
 587    };
 588   };
 589
 590
 591   {
 592    int order = 102;
 593    int trans = 112;
 594    int uplo = 121;
 595    int diag = 132;
 596    int N = 2;
 597    double A[] = { 0.355, -0.698, -0.42 };
 598    double X[] = { -0.702, -0.661 };
 599    int incX = -1;
 600    double x_expected[] = { -1.163378, -0.661 };
 601    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 602    {
 603      int i;
 604      for (i = 0; i < 2; i++) {
 605        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1339)");
 606      }
 607    };
 608   };
 609
 610
 611   {
 612    int order = 102;
 613    int trans = 112;
 614    int uplo = 122;
 615    int diag = 131;
 616    int N = 2;
 617    double A[] = { 0.355, -0.698, -0.42 };
 618    double X[] = { -0.702, -0.661 };
 619    int incX = -1;
 620    double x_expected[] = { 1.67142857143, 1.42438631791 };
 621    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 622    {
 623      int i;
 624      for (i = 0; i < 2; i++) {
 625        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1340)");
 626      }
 627    };
 628   };
 629
 630
 631   {
 632    int order = 102;
 633    int trans = 112;
 634    int uplo = 122;
 635    int diag = 132;
 636    int N = 2;
 637    double A[] = { 0.355, -0.698, -0.42 };
 638    double X[] = { -0.702, -0.661 };
 639    int incX = -1;
 640    double x_expected[] = { -0.702, -1.150996 };
 641    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
 642    {
 643      int i;
 644      for (i = 0; i < 2; i++) {
 645        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1341)");
 646      }
 647    };
 648   };
 649
 650
 651   {
 652    int order = 101;
 653    int trans = 111;
 654    int uplo = 121;
 655    int diag = 131;
 656    int N = 2;
 657    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 658    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 659    int incX = -1;
 660    float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
 661    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 662    {
 663      int i;
 664      for (i = 0; i < 2; i++) {
 665        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1342) real");
 666        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1342) imag");
 667      };
 668    };
 669   };
 670
 671
 672   {
 673    int order = 101;
 674    int trans = 111;
 675    int uplo = 121;
 676    int diag = 132;
 677    int N = 2;
 678    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 679    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 680    int incX = -1;
 681    float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
 682    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 683    {
 684      int i;
 685      for (i = 0; i < 2; i++) {
 686        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1343) real");
 687        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1343) imag");
 688      };
 689    };
 690   };
 691
 692
 693   {
 694    int order = 101;
 695    int trans = 111;
 696    int uplo = 122;
 697    int diag = 131;
 698    int N = 2;
 699    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 700    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 701    int incX = -1;
 702    float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
 703    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 704    {
 705      int i;
 706      for (i = 0; i < 2; i++) {
 707        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1344) real");
 708        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1344) imag");
 709      };
 710    };
 711   };
 712
 713
 714   {
 715    int order = 101;
 716    int trans = 111;
 717    int uplo = 122;
 718    int diag = 132;
 719    int N = 2;
 720    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 721    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 722    int incX = -1;
 723    float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
 724    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 725    {
 726      int i;
 727      for (i = 0; i < 2; i++) {
 728        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1345) real");
 729        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1345) imag");
 730      };
 731    };
 732   };
 733
 734
 735   {
 736    int order = 102;
 737    int trans = 111;
 738    int uplo = 121;
 739    int diag = 131;
 740    int N = 2;
 741    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 742    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 743    int incX = -1;
 744    float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
 745    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 746    {
 747      int i;
 748      for (i = 0; i < 2; i++) {
 749        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1346) real");
 750        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1346) imag");
 751      };
 752    };
 753   };
 754
 755
 756   {
 757    int order = 102;
 758    int trans = 111;
 759    int uplo = 121;
 760    int diag = 132;
 761    int N = 2;
 762    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 763    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 764    int incX = -1;
 765    float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
 766    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 767    {
 768      int i;
 769      for (i = 0; i < 2; i++) {
 770        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1347) real");
 771        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1347) imag");
 772      };
 773    };
 774   };
 775
 776
 777   {
 778    int order = 102;
 779    int trans = 111;
 780    int uplo = 122;
 781    int diag = 131;
 782    int N = 2;
 783    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 784    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 785    int incX = -1;
 786    float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
 787    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 788    {
 789      int i;
 790      for (i = 0; i < 2; i++) {
 791        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1348) real");
 792        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1348) imag");
 793      };
 794    };
 795   };
 796
 797
 798   {
 799    int order = 102;
 800    int trans = 111;
 801    int uplo = 122;
 802    int diag = 132;
 803    int N = 2;
 804    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 805    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 806    int incX = -1;
 807    float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
 808    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 809    {
 810      int i;
 811      for (i = 0; i < 2; i++) {
 812        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1349) real");
 813        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1349) imag");
 814      };
 815    };
 816   };
 817
 818
 819   {
 820    int order = 101;
 821    int trans = 112;
 822    int uplo = 121;
 823    int diag = 131;
 824    int N = 2;
 825    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 826    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 827    int incX = -1;
 828    float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
 829    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 830    {
 831      int i;
 832      for (i = 0; i < 2; i++) {
 833        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1350) real");
 834        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1350) imag");
 835      };
 836    };
 837   };
 838
 839
 840   {
 841    int order = 101;
 842    int trans = 112;
 843    int uplo = 121;
 844    int diag = 132;
 845    int N = 2;
 846    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 847    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 848    int incX = -1;
 849    float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
 850    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 851    {
 852      int i;
 853      for (i = 0; i < 2; i++) {
 854        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1351) real");
 855        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1351) imag");
 856      };
 857    };
 858   };
 859
 860
 861   {
 862    int order = 101;
 863    int trans = 112;
 864    int uplo = 122;
 865    int diag = 131;
 866    int N = 2;
 867    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 868    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 869    int incX = -1;
 870    float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
 871    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 872    {
 873      int i;
 874      for (i = 0; i < 2; i++) {
 875        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1352) real");
 876        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1352) imag");
 877      };
 878    };
 879   };
 880
 881
 882   {
 883    int order = 101;
 884    int trans = 112;
 885    int uplo = 122;
 886    int diag = 132;
 887    int N = 2;
 888    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 889    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 890    int incX = -1;
 891    float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
 892    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 893    {
 894      int i;
 895      for (i = 0; i < 2; i++) {
 896        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1353) real");
 897        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1353) imag");
 898      };
 899    };
 900   };
 901
 902
 903   {
 904    int order = 102;
 905    int trans = 112;
 906    int uplo = 121;
 907    int diag = 131;
 908    int N = 2;
 909    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 910    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 911    int incX = -1;
 912    float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
 913    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 914    {
 915      int i;
 916      for (i = 0; i < 2; i++) {
 917        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1354) real");
 918        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1354) imag");
 919      };
 920    };
 921   };
 922
 923
 924   {
 925    int order = 102;
 926    int trans = 112;
 927    int uplo = 121;
 928    int diag = 132;
 929    int N = 2;
 930    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 931    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 932    int incX = -1;
 933    float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
 934    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 935    {
 936      int i;
 937      for (i = 0; i < 2; i++) {
 938        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1355) real");
 939        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1355) imag");
 940      };
 941    };
 942   };
 943
 944
 945   {
 946    int order = 102;
 947    int trans = 112;
 948    int uplo = 122;
 949    int diag = 131;
 950    int N = 2;
 951    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 952    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 953    int incX = -1;
 954    float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
 955    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 956    {
 957      int i;
 958      for (i = 0; i < 2; i++) {
 959        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1356) real");
 960        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1356) imag");
 961      };
 962    };
 963   };
 964
 965
 966   {
 967    int order = 102;
 968    int trans = 112;
 969    int uplo = 122;
 970    int diag = 132;
 971    int N = 2;
 972    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 973    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 974    int incX = -1;
 975    float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
 976    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 977    {
 978      int i;
 979      for (i = 0; i < 2; i++) {
 980        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1357) real");
 981        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1357) imag");
 982      };
 983    };
 984   };
 985
 986
 987   {
 988    int order = 101;
 989    int trans = 113;
 990    int uplo = 121;
 991    int diag = 131;
 992    int N = 2;
 993    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
 994    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
 995    int incX = -1;
 996    float x_expected[] = { 2.15867f, 1.69498f, 1.69471f, 0.104738f };
 997    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
 998    {
 999      int i;
1000      for (i = 0; i < 2; i++) {
1001        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1358) real");
1002        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1358) imag");
1003      };
1004    };
1005   };
1006
1007
1008   {
1009    int order = 101;
1010    int trans = 113;
1011    int uplo = 121;
1012    int diag = 132;
1013    int N = 2;
1014    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
1015    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
1016    int incX = -1;
1017    float x_expected[] = { -0.613252f, 0.561896f, -0.072f, 0.642f };
1018    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
1019    {
1020      int i;
1021      for (i = 0; i < 2; i++) {
1022        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1359) real");
1023        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1359) imag");
1024      };
1025    };
1026   };
1027
1028
1029   {
1030    int order = 101;
1031    int trans = 113;
1032    int uplo = 122;
1033    int diag = 131;
1034    int N = 2;
1035    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
1036    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
1037    int incX = -1;
1038    float x_expected[] = { 1.00465f, 0.327432f, 3.44853f, 2.273f };
1039    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
1040    {
1041      int i;
1042      for (i = 0; i < 2; i++) {
1043        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1360) real");
1044        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1360) imag");
1045      };
1046    };
1047   };
1048
1049
1050   {
1051    int order = 101;
1052    int trans = 113;
1053    int uplo = 122;
1054    int diag = 132;
1055    int N = 2;
1056    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
1057    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
1058    int incX = -1;
1059    float x_expected[] = { -0.133f, 0.998f, -0.806168f, -0.053086f };
1060    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
1061    {
1062      int i;
1063      for (i = 0; i < 2; i++) {
1064        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1361) real");
1065        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1361) imag");
1066      };
1067    };
1068   };
1069
1070
1071   {
1072    int order = 102;
1073    int trans = 113;
1074    int uplo = 121;
1075    int diag = 131;
1076    int N = 2;
1077    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
1078    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
1079    int incX = -1;
1080    float x_expected[] = { 2.15867f, 1.69498f, 1.69471f, 0.104738f };
1081    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
1082    {
1083      int i;
1084      for (i = 0; i < 2; i++) {
1085        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1362) real");
1086        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1362) imag");
1087      };
1088    };
1089   };
1090
1091
1092   {
1093    int order = 102;
1094    int trans = 113;
1095    int uplo = 121;
1096    int diag = 132;
1097    int N = 2;
1098    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
1099    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
1100    int incX = -1;
1101    float x_expected[] = { -0.613252f, 0.561896f, -0.072f, 0.642f };
1102    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
1103    {
1104      int i;
1105      for (i = 0; i < 2; i++) {
1106        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1363) real");
1107        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1363) imag");
1108      };
1109    };
1110   };
1111
1112
1113   {
1114    int order = 102;
1115    int trans = 113;
1116    int uplo = 122;
1117    int diag = 131;
1118    int N = 2;
1119    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
1120    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
1121    int incX = -1;
1122    float x_expected[] = { 1.00465f, 0.327432f, 3.44853f, 2.273f };
1123    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
1124    {
1125      int i;
1126      for (i = 0; i < 2; i++) {
1127        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1364) real");
1128        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1364) imag");
1129      };
1130    };
1131   };
1132
1133
1134   {
1135    int order = 102;
1136    int trans = 113;
1137    int uplo = 122;
1138    int diag = 132;
1139    int N = 2;
1140    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
1141    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
1142    int incX = -1;
1143    float x_expected[] = { -0.133f, 0.998f, -0.806168f, -0.053086f };
1144    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
1145    {
1146      int i;
1147      for (i = 0; i < 2; i++) {
1148        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1365) real");
1149        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1365) imag");
1150      };
1151    };
1152   };
1153
1154
1155   {
1156    int order = 101;
1157    int trans = 111;
1158    int uplo = 121;
1159    int diag = 131;
1160    int N = 2;
1161    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1162    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1163    int incX = -1;
1164    double x_expected[] = { 0.430509772467, -0.0927067365535, -0.611144484555, 0.999982608216 };
1165    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1166    {
1167      int i;
1168      for (i = 0; i < 2; i++) {
1169        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1366) real");
1170        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1366) imag");
1171      };
1172    };
1173   };
1174
1175
1176   {
1177    int order = 101;
1178    int trans = 111;
1179    int uplo = 121;
1180    int diag = 132;
1181    int N = 2;
1182    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1183    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1184    int incX = -1;
1185    double x_expected[] = { 0.344, -0.143, -0.795928, -0.523879 };
1186    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1187    {
1188      int i;
1189      for (i = 0; i < 2; i++) {
1190        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1367) real");
1191        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1367) imag");
1192      };
1193    };
1194   };
1195
1196
1197   {
1198    int order = 101;
1199    int trans = 111;
1200    int uplo = 122;
1201    int diag = 131;
1202    int N = 2;
1203    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1204    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1205    int incX = -1;
1206    double x_expected[] = { 1.00136463678, -1.84591534629, -1.12140892769, 0.696784840869 };
1207    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1208    {
1209      int i;
1210      for (i = 0; i < 2; i++) {
1211        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1368) real");
1212        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1368) imag");
1213      };
1214    };
1215   };
1216
1217
1218   {
1219    int order = 101;
1220    int trans = 111;
1221    int uplo = 122;
1222    int diag = 132;
1223    int N = 2;
1224    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1225    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1226    int incX = -1;
1227    double x_expected[] = { 1.707508, -0.042521, -0.668, -0.945 };
1228    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1229    {
1230      int i;
1231      for (i = 0; i < 2; i++) {
1232        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1369) real");
1233        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1369) imag");
1234      };
1235    };
1236   };
1237
1238
1239   {
1240    int order = 102;
1241    int trans = 111;
1242    int uplo = 121;
1243    int diag = 131;
1244    int N = 2;
1245    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1246    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1247    int incX = -1;
1248    double x_expected[] = { 0.430509772467, -0.0927067365535, -0.611144484555, 0.999982608216 };
1249    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1250    {
1251      int i;
1252      for (i = 0; i < 2; i++) {
1253        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1370) real");
1254        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1370) imag");
1255      };
1256    };
1257   };
1258
1259
1260   {
1261    int order = 102;
1262    int trans = 111;
1263    int uplo = 121;
1264    int diag = 132;
1265    int N = 2;
1266    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1267    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1268    int incX = -1;
1269    double x_expected[] = { 0.344, -0.143, -0.795928, -0.523879 };
1270    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1271    {
1272      int i;
1273      for (i = 0; i < 2; i++) {
1274        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1371) real");
1275        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1371) imag");
1276      };
1277    };
1278   };
1279
1280
1281   {
1282    int order = 102;
1283    int trans = 111;
1284    int uplo = 122;
1285    int diag = 131;
1286    int N = 2;
1287    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1288    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1289    int incX = -1;
1290    double x_expected[] = { 1.00136463678, -1.84591534629, -1.12140892769, 0.696784840869 };
1291    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1292    {
1293      int i;
1294      for (i = 0; i < 2; i++) {
1295        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1372) real");
1296        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1372) imag");
1297      };
1298    };
1299   };
1300
1301
1302   {
1303    int order = 102;
1304    int trans = 111;
1305    int uplo = 122;
1306    int diag = 132;
1307    int N = 2;
1308    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1309    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1310    int incX = -1;
1311    double x_expected[] = { 1.707508, -0.042521, -0.668, -0.945 };
1312    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1313    {
1314      int i;
1315      for (i = 0; i < 2; i++) {
1316        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1373) real");
1317        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1373) imag");
1318      };
1319    };
1320   };
1321
1322
1323   {
1324    int order = 101;
1325    int trans = 112;
1326    int uplo = 121;
1327    int diag = 131;
1328    int N = 2;
1329    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1330    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1331    int incX = -1;
1332    double x_expected[] = { 1.00136463678, -1.84591534629, -1.12140892769, 0.696784840869 };
1333    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1334    {
1335      int i;
1336      for (i = 0; i < 2; i++) {
1337        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1374) real");
1338        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1374) imag");
1339      };
1340    };
1341   };
1342
1343
1344   {
1345    int order = 101;
1346    int trans = 112;
1347    int uplo = 121;
1348    int diag = 132;
1349    int N = 2;
1350    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1351    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1352    int incX = -1;
1353    double x_expected[] = { 1.707508, -0.042521, -0.668, -0.945 };
1354    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1355    {
1356      int i;
1357      for (i = 0; i < 2; i++) {
1358        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1375) real");
1359        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1375) imag");
1360      };
1361    };
1362   };
1363
1364
1365   {
1366    int order = 101;
1367    int trans = 112;
1368    int uplo = 122;
1369    int diag = 131;
1370    int N = 2;
1371    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1372    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1373    int incX = -1;
1374    double x_expected[] = { 0.430509772467, -0.0927067365535, -0.611144484555, 0.999982608216 };
1375    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1376    {
1377      int i;
1378      for (i = 0; i < 2; i++) {
1379        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1376) real");
1380        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1376) imag");
1381      };
1382    };
1383   };
1384
1385
1386   {
1387    int order = 101;
1388    int trans = 112;
1389    int uplo = 122;
1390    int diag = 132;
1391    int N = 2;
1392    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1393    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1394    int incX = -1;
1395    double x_expected[] = { 0.344, -0.143, -0.795928, -0.523879 };
1396    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1397    {
1398      int i;
1399      for (i = 0; i < 2; i++) {
1400        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1377) real");
1401        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1377) imag");
1402      };
1403    };
1404   };
1405
1406
1407   {
1408    int order = 102;
1409    int trans = 112;
1410    int uplo = 121;
1411    int diag = 131;
1412    int N = 2;
1413    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1414    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1415    int incX = -1;
1416    double x_expected[] = { 1.00136463678, -1.84591534629, -1.12140892769, 0.696784840869 };
1417    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1418    {
1419      int i;
1420      for (i = 0; i < 2; i++) {
1421        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1378) real");
1422        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1378) imag");
1423      };
1424    };
1425   };
1426
1427
1428   {
1429    int order = 102;
1430    int trans = 112;
1431    int uplo = 121;
1432    int diag = 132;
1433    int N = 2;
1434    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1435    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1436    int incX = -1;
1437    double x_expected[] = { 1.707508, -0.042521, -0.668, -0.945 };
1438    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1439    {
1440      int i;
1441      for (i = 0; i < 2; i++) {
1442        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1379) real");
1443        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1379) imag");
1444      };
1445    };
1446   };
1447
1448
1449   {
1450    int order = 102;
1451    int trans = 112;
1452    int uplo = 122;
1453    int diag = 131;
1454    int N = 2;
1455    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1456    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1457    int incX = -1;
1458    double x_expected[] = { 0.430509772467, -0.0927067365535, -0.611144484555, 0.999982608216 };
1459    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1460    {
1461      int i;
1462      for (i = 0; i < 2; i++) {
1463        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1380) real");
1464        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1380) imag");
1465      };
1466    };
1467   };
1468
1469
1470   {
1471    int order = 102;
1472    int trans = 112;
1473    int uplo = 122;
1474    int diag = 132;
1475    int N = 2;
1476    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1477    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1478    int incX = -1;
1479    double x_expected[] = { 0.344, -0.143, -0.795928, -0.523879 };
1480    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1481    {
1482      int i;
1483      for (i = 0; i < 2; i++) {
1484        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1381) real");
1485        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1381) imag");
1486      };
1487    };
1488   };
1489
1490
1491   {
1492    int order = 101;
1493    int trans = 113;
1494    int uplo = 121;
1495    int diag = 131;
1496    int N = 2;
1497    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1498    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1499    int incX = -1;
1500    double x_expected[] = { -1.47384781823, -0.286556198408, 1.03098932879, -0.824698794397 };
1501    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1502    {
1503      int i;
1504      for (i = 0; i < 2; i++) {
1505        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1382) real");
1506        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1382) imag");
1507      };
1508    };
1509   };
1510
1511
1512   {
1513    int order = 101;
1514    int trans = 113;
1515    int uplo = 121;
1516    int diag = 132;
1517    int N = 2;
1518    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1519    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1520    int incX = -1;
1521    double x_expected[] = { -0.016172, 1.175911, -0.668, -0.945 };
1522    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1523    {
1524      int i;
1525      for (i = 0; i < 2; i++) {
1526        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1383) real");
1527        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1383) imag");
1528      };
1529    };
1530   };
1531
1532
1533   {
1534    int order = 101;
1535    int trans = 113;
1536    int uplo = 122;
1537    int diag = 131;
1538    int N = 2;
1539    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1540    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1541    int incX = -1;
1542    double x_expected[] = { 0.369363905801, -0.239798891331, 1.1759505739, -1.40027235656 };
1543    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1544    {
1545      int i;
1546      for (i = 0; i < 2; i++) {
1547        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1384) real");
1548        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1384) imag");
1549      };
1550    };
1551   };
1552
1553
1554   {
1555    int order = 101;
1556    int trans = 113;
1557    int uplo = 122;
1558    int diag = 132;
1559    int N = 2;
1560    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1561    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1562    int incX = -1;
1563    double x_expected[] = { 0.344, -0.143, -1.05676, -1.151335 };
1564    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1565    {
1566      int i;
1567      for (i = 0; i < 2; i++) {
1568        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1385) real");
1569        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1385) imag");
1570      };
1571    };
1572   };
1573
1574
1575   {
1576    int order = 102;
1577    int trans = 113;
1578    int uplo = 121;
1579    int diag = 131;
1580    int N = 2;
1581    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1582    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1583    int incX = -1;
1584    double x_expected[] = { -1.47384781823, -0.286556198408, 1.03098932879, -0.824698794397 };
1585    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1586    {
1587      int i;
1588      for (i = 0; i < 2; i++) {
1589        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1386) real");
1590        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1386) imag");
1591      };
1592    };
1593   };
1594
1595
1596   {
1597    int order = 102;
1598    int trans = 113;
1599    int uplo = 121;
1600    int diag = 132;
1601    int N = 2;
1602    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1603    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1604    int incX = -1;
1605    double x_expected[] = { -0.016172, 1.175911, -0.668, -0.945 };
1606    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1607    {
1608      int i;
1609      for (i = 0; i < 2; i++) {
1610        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1387) real");
1611        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1387) imag");
1612      };
1613    };
1614   };
1615
1616
1617   {
1618    int order = 102;
1619    int trans = 113;
1620    int uplo = 122;
1621    int diag = 131;
1622    int N = 2;
1623    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1624    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1625    int incX = -1;
1626    double x_expected[] = { 0.369363905801, -0.239798891331, 1.1759505739, -1.40027235656 };
1627    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1628    {
1629      int i;
1630      for (i = 0; i < 2; i++) {
1631        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1388) real");
1632        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1388) imag");
1633      };
1634    };
1635   };
1636
1637
1638   {
1639    int order = 102;
1640    int trans = 113;
1641    int uplo = 122;
1642    int diag = 132;
1643    int N = 2;
1644    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
1645    double X[] = { 0.344, -0.143, -0.668, -0.945 };
1646    int incX = -1;
1647    double x_expected[] = { 0.344, -0.143, -1.05676, -1.151335 };
1648    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
1649    {
1650      int i;
1651      for (i = 0; i < 2; i++) {
1652        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1389) real");
1653        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1389) imag");
1654      };
1655    };
1656   };
1657
1658
1659 }