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