htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_rot.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_rot (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int N = 1;
  13    float c = 0.0f;
  14    float s = 0.0f;
  15    float X[] = { -0.314f };
  16    int incX = 1;
  17    float Y[] = { -0.406f };
  18    int incY = -1;
  19    float x_expected[] = { 0.0f };
  20    float y_expected[] = { 0.0f };
  21    cblas_srot(N, X, incX, Y, incY, c, s);
  22    {
  23      int i;
  24      for (i = 0; i < 1; i++) {
  25        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 558)");
  26      }
  27    };
  28    {
  29      int i;
  30      for (i = 0; i < 1; i++) {
  31        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 559)");
  32      }
  33    };
  34   };
  35
  36
  37   {
  38    int N = 1;
  39    float c = 0.866025403784f;
  40    float s = 0.5f;
  41    float X[] = { -0.314f };
  42    int incX = 1;
  43    float Y[] = { -0.406f };
  44    int incY = -1;
  45    float x_expected[] = { -0.474932f };
  46    float y_expected[] = { -0.194606f };
  47    cblas_srot(N, X, incX, Y, incY, c, s);
  48    {
  49      int i;
  50      for (i = 0; i < 1; i++) {
  51        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 560)");
  52      }
  53    };
  54    {
  55      int i;
  56      for (i = 0; i < 1; i++) {
  57        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 561)");
  58      }
  59    };
  60   };
  61
  62
  63   {
  64    int N = 1;
  65    float c = 0.0f;
  66    float s = -1.0f;
  67    float X[] = { -0.314f };
  68    int incX = 1;
  69    float Y[] = { -0.406f };
  70    int incY = -1;
  71    float x_expected[] = { 0.406f };
  72    float y_expected[] = { -0.314f };
  73    cblas_srot(N, X, incX, Y, incY, c, s);
  74    {
  75      int i;
  76      for (i = 0; i < 1; i++) {
  77        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 562)");
  78      }
  79    };
  80    {
  81      int i;
  82      for (i = 0; i < 1; i++) {
  83        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 563)");
  84      }
  85    };
  86   };
  87
  88
  89   {
  90    int N = 1;
  91    float c = -1.0f;
  92    float s = 0.0f;
  93    float X[] = { -0.314f };
  94    int incX = 1;
  95    float Y[] = { -0.406f };
  96    int incY = -1;
  97    float x_expected[] = { 0.314f };
  98    float y_expected[] = { 0.406f };
  99    cblas_srot(N, X, incX, Y, incY, c, s);
 100    {
 101      int i;
 102      for (i = 0; i < 1; i++) {
 103        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 564)");
 104      }
 105    };
 106    {
 107      int i;
 108      for (i = 0; i < 1; i++) {
 109        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 565)");
 110      }
 111    };
 112   };
 113
 114
 115   {
 116    int N = 1;
 117    double c = 0;
 118    double s = 0;
 119    double X[] = { -0.493 };
 120    int incX = 1;
 121    double Y[] = { -0.014 };
 122    int incY = -1;
 123    double x_expected[] = { 0.0 };
 124    double y_expected[] = { 0.0 };
 125    cblas_drot(N, X, incX, Y, incY, c, s);
 126    {
 127      int i;
 128      for (i = 0; i < 1; i++) {
 129        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 566)");
 130      }
 131    };
 132    {
 133      int i;
 134      for (i = 0; i < 1; i++) {
 135        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 567)");
 136      }
 137    };
 138   };
 139
 140
 141   {
 142    int N = 1;
 143    double c = 0.866025403784;
 144    double s = 0.5;
 145    double X[] = { -0.493 };
 146    int incX = 1;
 147    double Y[] = { -0.014 };
 148    int incY = -1;
 149    double x_expected[] = { -0.433950524066 };
 150    double y_expected[] = { 0.234375644347 };
 151    cblas_drot(N, X, incX, Y, incY, c, s);
 152    {
 153      int i;
 154      for (i = 0; i < 1; i++) {
 155        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 568)");
 156      }
 157    };
 158    {
 159      int i;
 160      for (i = 0; i < 1; i++) {
 161        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 569)");
 162      }
 163    };
 164   };
 165
 166
 167   {
 168    int N = 1;
 169    double c = 0;
 170    double s = -1;
 171    double X[] = { -0.493 };
 172    int incX = 1;
 173    double Y[] = { -0.014 };
 174    int incY = -1;
 175    double x_expected[] = { 0.014 };
 176    double y_expected[] = { -0.493 };
 177    cblas_drot(N, X, incX, Y, incY, c, s);
 178    {
 179      int i;
 180      for (i = 0; i < 1; i++) {
 181        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 570)");
 182      }
 183    };
 184    {
 185      int i;
 186      for (i = 0; i < 1; i++) {
 187        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 571)");
 188      }
 189    };
 190   };
 191
 192
 193   {
 194    int N = 1;
 195    double c = -1;
 196    double s = 0;
 197    double X[] = { -0.493 };
 198    int incX = 1;
 199    double Y[] = { -0.014 };
 200    int incY = -1;
 201    double x_expected[] = { 0.493 };
 202    double y_expected[] = { 0.014 };
 203    cblas_drot(N, X, incX, Y, incY, c, s);
 204    {
 205      int i;
 206      for (i = 0; i < 1; i++) {
 207        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 572)");
 208      }
 209    };
 210    {
 211      int i;
 212      for (i = 0; i < 1; i++) {
 213        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 573)");
 214      }
 215    };
 216   };
 217
 218
 219   {
 220    int N = 1;
 221    float c = 0.0f;
 222    float s = 0.0f;
 223    float X[] = { -0.808f };
 224    int incX = -1;
 225    float Y[] = { -0.511f };
 226    int incY = 1;
 227    float x_expected[] = { 0.0f };
 228    float y_expected[] = { 0.0f };
 229    cblas_srot(N, X, incX, Y, incY, c, s);
 230    {
 231      int i;
 232      for (i = 0; i < 1; i++) {
 233        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 574)");
 234      }
 235    };
 236    {
 237      int i;
 238      for (i = 0; i < 1; i++) {
 239        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 575)");
 240      }
 241    };
 242   };
 243
 244
 245   {
 246    int N = 1;
 247    float c = 0.866025403784f;
 248    float s = 0.5f;
 249    float X[] = { -0.808f };
 250    int incX = -1;
 251    float Y[] = { -0.511f };
 252    int incY = 1;
 253    float x_expected[] = { -0.955249f };
 254    float y_expected[] = { -0.038539f };
 255    cblas_srot(N, X, incX, Y, incY, c, s);
 256    {
 257      int i;
 258      for (i = 0; i < 1; i++) {
 259        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 576)");
 260      }
 261    };
 262    {
 263      int i;
 264      for (i = 0; i < 1; i++) {
 265        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 577)");
 266      }
 267    };
 268   };
 269
 270
 271   {
 272    int N = 1;
 273    float c = 0.0f;
 274    float s = -1.0f;
 275    float X[] = { -0.808f };
 276    int incX = -1;
 277    float Y[] = { -0.511f };
 278    int incY = 1;
 279    float x_expected[] = { 0.511f };
 280    float y_expected[] = { -0.808f };
 281    cblas_srot(N, X, incX, Y, incY, c, s);
 282    {
 283      int i;
 284      for (i = 0; i < 1; i++) {
 285        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 578)");
 286      }
 287    };
 288    {
 289      int i;
 290      for (i = 0; i < 1; i++) {
 291        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 579)");
 292      }
 293    };
 294   };
 295
 296
 297   {
 298    int N = 1;
 299    float c = -1.0f;
 300    float s = 0.0f;
 301    float X[] = { -0.808f };
 302    int incX = -1;
 303    float Y[] = { -0.511f };
 304    int incY = 1;
 305    float x_expected[] = { 0.808f };
 306    float y_expected[] = { 0.511f };
 307    cblas_srot(N, X, incX, Y, incY, c, s);
 308    {
 309      int i;
 310      for (i = 0; i < 1; i++) {
 311        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 580)");
 312      }
 313    };
 314    {
 315      int i;
 316      for (i = 0; i < 1; i++) {
 317        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 581)");
 318      }
 319    };
 320   };
 321
 322
 323   {
 324    int N = 1;
 325    double c = 0;
 326    double s = 0;
 327    double X[] = { -0.176 };
 328    int incX = -1;
 329    double Y[] = { -0.165 };
 330    int incY = 1;
 331    double x_expected[] = { 0.0 };
 332    double y_expected[] = { 0.0 };
 333    cblas_drot(N, X, incX, Y, incY, c, s);
 334    {
 335      int i;
 336      for (i = 0; i < 1; i++) {
 337        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 582)");
 338      }
 339    };
 340    {
 341      int i;
 342      for (i = 0; i < 1; i++) {
 343        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 583)");
 344      }
 345    };
 346   };
 347
 348
 349   {
 350    int N = 1;
 351    double c = 0.866025403784;
 352    double s = 0.5;
 353    double X[] = { -0.176 };
 354    int incX = -1;
 355    double Y[] = { -0.165 };
 356    int incY = 1;
 357    double x_expected[] = { -0.234920471066 };
 358    double y_expected[] = { -0.0548941916244 };
 359    cblas_drot(N, X, incX, Y, incY, c, s);
 360    {
 361      int i;
 362      for (i = 0; i < 1; i++) {
 363        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 584)");
 364      }
 365    };
 366    {
 367      int i;
 368      for (i = 0; i < 1; i++) {
 369        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 585)");
 370      }
 371    };
 372   };
 373
 374
 375   {
 376    int N = 1;
 377    double c = 0;
 378    double s = -1;
 379    double X[] = { -0.176 };
 380    int incX = -1;
 381    double Y[] = { -0.165 };
 382    int incY = 1;
 383    double x_expected[] = { 0.165 };
 384    double y_expected[] = { -0.176 };
 385    cblas_drot(N, X, incX, Y, incY, c, s);
 386    {
 387      int i;
 388      for (i = 0; i < 1; i++) {
 389        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 586)");
 390      }
 391    };
 392    {
 393      int i;
 394      for (i = 0; i < 1; i++) {
 395        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 587)");
 396      }
 397    };
 398   };
 399
 400
 401   {
 402    int N = 1;
 403    double c = -1;
 404    double s = 0;
 405    double X[] = { -0.176 };
 406    int incX = -1;
 407    double Y[] = { -0.165 };
 408    int incY = 1;
 409    double x_expected[] = { 0.176 };
 410    double y_expected[] = { 0.165 };
 411    cblas_drot(N, X, incX, Y, incY, c, s);
 412    {
 413      int i;
 414      for (i = 0; i < 1; i++) {
 415        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 588)");
 416      }
 417    };
 418    {
 419      int i;
 420      for (i = 0; i < 1; i++) {
 421        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 589)");
 422      }
 423    };
 424   };
 425
 426
 427   {
 428    int N = 1;
 429    float c = 0.0f;
 430    float s = 0.0f;
 431    float X[] = { -0.201f };
 432    int incX = -1;
 433    float Y[] = { 0.087f };
 434    int incY = -1;
 435    float x_expected[] = { 0.0f };
 436    float y_expected[] = { 0.0f };
 437    cblas_srot(N, X, incX, Y, incY, c, s);
 438    {
 439      int i;
 440      for (i = 0; i < 1; i++) {
 441        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 590)");
 442      }
 443    };
 444    {
 445      int i;
 446      for (i = 0; i < 1; i++) {
 447        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 591)");
 448      }
 449    };
 450   };
 451
 452
 453   {
 454    int N = 1;
 455    float c = 0.866025403784f;
 456    float s = 0.5f;
 457    float X[] = { -0.201f };
 458    int incX = -1;
 459    float Y[] = { 0.087f };
 460    int incY = -1;
 461    float x_expected[] = { -0.130571f };
 462    float y_expected[] = { 0.175844f };
 463    cblas_srot(N, X, incX, Y, incY, c, s);
 464    {
 465      int i;
 466      for (i = 0; i < 1; i++) {
 467        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 592)");
 468      }
 469    };
 470    {
 471      int i;
 472      for (i = 0; i < 1; i++) {
 473        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 593)");
 474      }
 475    };
 476   };
 477
 478
 479   {
 480    int N = 1;
 481    float c = 0.0f;
 482    float s = -1.0f;
 483    float X[] = { -0.201f };
 484    int incX = -1;
 485    float Y[] = { 0.087f };
 486    int incY = -1;
 487    float x_expected[] = { -0.087f };
 488    float y_expected[] = { -0.201f };
 489    cblas_srot(N, X, incX, Y, incY, c, s);
 490    {
 491      int i;
 492      for (i = 0; i < 1; i++) {
 493        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 594)");
 494      }
 495    };
 496    {
 497      int i;
 498      for (i = 0; i < 1; i++) {
 499        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 595)");
 500      }
 501    };
 502   };
 503
 504
 505   {
 506    int N = 1;
 507    float c = -1.0f;
 508    float s = 0.0f;
 509    float X[] = { -0.201f };
 510    int incX = -1;
 511    float Y[] = { 0.087f };
 512    int incY = -1;
 513    float x_expected[] = { 0.201f };
 514    float y_expected[] = { -0.087f };
 515    cblas_srot(N, X, incX, Y, incY, c, s);
 516    {
 517      int i;
 518      for (i = 0; i < 1; i++) {
 519        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 596)");
 520      }
 521    };
 522    {
 523      int i;
 524      for (i = 0; i < 1; i++) {
 525        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 597)");
 526      }
 527    };
 528   };
 529
 530
 531   {
 532    int N = 1;
 533    double c = 0;
 534    double s = 0;
 535    double X[] = { -0.464 };
 536    int incX = -1;
 537    double Y[] = { 0.7 };
 538    int incY = -1;
 539    double x_expected[] = { 0.0 };
 540    double y_expected[] = { 0.0 };
 541    cblas_drot(N, X, incX, Y, incY, c, s);
 542    {
 543      int i;
 544      for (i = 0; i < 1; i++) {
 545        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 598)");
 546      }
 547    };
 548    {
 549      int i;
 550      for (i = 0; i < 1; i++) {
 551        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 599)");
 552      }
 553    };
 554   };
 555
 556
 557   {
 558    int N = 1;
 559    double c = 0.866025403784;
 560    double s = 0.5;
 561    double X[] = { -0.464 };
 562    int incX = -1;
 563    double Y[] = { 0.7 };
 564    int incY = -1;
 565    double x_expected[] = { -0.051835787356 };
 566    double y_expected[] = { 0.838217782649 };
 567    cblas_drot(N, X, incX, Y, incY, c, s);
 568    {
 569      int i;
 570      for (i = 0; i < 1; i++) {
 571        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 600)");
 572      }
 573    };
 574    {
 575      int i;
 576      for (i = 0; i < 1; i++) {
 577        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 601)");
 578      }
 579    };
 580   };
 581
 582
 583   {
 584    int N = 1;
 585    double c = 0;
 586    double s = -1;
 587    double X[] = { -0.464 };
 588    int incX = -1;
 589    double Y[] = { 0.7 };
 590    int incY = -1;
 591    double x_expected[] = { -0.7 };
 592    double y_expected[] = { -0.464 };
 593    cblas_drot(N, X, incX, Y, incY, c, s);
 594    {
 595      int i;
 596      for (i = 0; i < 1; i++) {
 597        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 602)");
 598      }
 599    };
 600    {
 601      int i;
 602      for (i = 0; i < 1; i++) {
 603        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 603)");
 604      }
 605    };
 606   };
 607
 608
 609   {
 610    int N = 1;
 611    double c = -1;
 612    double s = 0;
 613    double X[] = { -0.464 };
 614    int incX = -1;
 615    double Y[] = { 0.7 };
 616    int incY = -1;
 617    double x_expected[] = { 0.464 };
 618    double y_expected[] = { -0.7 };
 619    cblas_drot(N, X, incX, Y, incY, c, s);
 620    {
 621      int i;
 622      for (i = 0; i < 1; i++) {
 623        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 604)");
 624      }
 625    };
 626    {
 627      int i;
 628      for (i = 0; i < 1; i++) {
 629        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 605)");
 630      }
 631    };
 632   };
 633
 634
 635 }