Blame cblas/test_rotg.c

Packit 67cb25
#include <gsl/gsl_test.h>
Packit 67cb25
#include <gsl/gsl_ieee_utils.h>
Packit 67cb25
#include <gsl/gsl_math.h>
Packit 67cb25
#include <gsl/gsl_cblas.h>
Packit 67cb25
Packit 67cb25
#include "tests.h"
Packit 67cb25
Packit 67cb25
void
Packit 67cb25
test_rotg (void) {
Packit 67cb25
const double flteps = 1e-4, dbleps = 1e-6;
Packit 67cb25
  {
Packit 67cb25
   float a = -1.5f;
Packit 67cb25
   float b = -1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -2.12132034356f;
Packit 67cb25
   float z_expected = 1.41421356237f;
Packit 67cb25
   float c_expected = 0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 166)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 167)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 168)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 169)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.5f;
Packit 67cb25
   float b = -1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.80277563773f;
Packit 67cb25
   float z_expected = 0.554700196225f;
Packit 67cb25
   float c_expected = 0.832050294338f;
Packit 67cb25
   float s_expected = 0.554700196225f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 170)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 171)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 172)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 173)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.5f;
Packit 67cb25
   float b = -0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.50332963784f;
Packit 67cb25
   float z_expected = 0.0665190105238f;
Packit 67cb25
   float c_expected = 0.997785157857f;
Packit 67cb25
   float s_expected = 0.0665190105238f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 174)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 175)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 176)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 177)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.5f;
Packit 67cb25
   float b = 0.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.5f;
Packit 67cb25
   float z_expected = -0.0f;
Packit 67cb25
   float c_expected = 1.0f;
Packit 67cb25
   float s_expected = -0.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 178)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 179)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 180)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 181)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.5f;
Packit 67cb25
   float b = 0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.50332963784f;
Packit 67cb25
   float z_expected = -0.0665190105238f;
Packit 67cb25
   float c_expected = 0.997785157857f;
Packit 67cb25
   float s_expected = -0.0665190105238f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 182)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 183)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 184)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 185)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.5f;
Packit 67cb25
   float b = 1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.80277563773f;
Packit 67cb25
   float z_expected = -0.554700196225f;
Packit 67cb25
   float c_expected = 0.832050294338f;
Packit 67cb25
   float s_expected = -0.554700196225f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 186)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 187)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 188)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 189)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.5f;
Packit 67cb25
   float b = 1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 2.12132034356f;
Packit 67cb25
   float z_expected = -1.41421356237f;
Packit 67cb25
   float c_expected = -0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 190)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 191)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 192)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 193)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.0f;
Packit 67cb25
   float b = -1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.80277563773f;
Packit 67cb25
   float z_expected = 1.80277563773f;
Packit 67cb25
   float c_expected = 0.554700196225f;
Packit 67cb25
   float s_expected = 0.832050294338f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 194)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 195)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 196)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 197)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.0f;
Packit 67cb25
   float b = -1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.41421356237f;
Packit 67cb25
   float z_expected = 1.41421356237f;
Packit 67cb25
   float c_expected = 0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 198)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 199)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 200)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 201)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.0f;
Packit 67cb25
   float b = -0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.00498756211f;
Packit 67cb25
   float z_expected = 0.099503719021f;
Packit 67cb25
   float c_expected = 0.99503719021f;
Packit 67cb25
   float s_expected = 0.099503719021f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 202)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 203)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 204)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 205)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.0f;
Packit 67cb25
   float b = 0.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.0f;
Packit 67cb25
   float z_expected = -0.0f;
Packit 67cb25
   float c_expected = 1.0f;
Packit 67cb25
   float s_expected = -0.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 206)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 207)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 208)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 209)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.0f;
Packit 67cb25
   float b = 0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.00498756211f;
Packit 67cb25
   float z_expected = -0.099503719021f;
Packit 67cb25
   float c_expected = 0.99503719021f;
Packit 67cb25
   float s_expected = -0.099503719021f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 210)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 211)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 212)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 213)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.0f;
Packit 67cb25
   float b = 1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.41421356237f;
Packit 67cb25
   float z_expected = -1.41421356237f;
Packit 67cb25
   float c_expected = -0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 214)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 215)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 216)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 217)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -1.0f;
Packit 67cb25
   float b = 1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.80277563773f;
Packit 67cb25
   float z_expected = -1.80277563773f;
Packit 67cb25
   float c_expected = -0.554700196225f;
Packit 67cb25
   float s_expected = 0.832050294338f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 218)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 219)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 220)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 221)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -0.1f;
Packit 67cb25
   float b = -1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.50332963784f;
Packit 67cb25
   float z_expected = 15.0332963784f;
Packit 67cb25
   float c_expected = 0.0665190105238f;
Packit 67cb25
   float s_expected = 0.997785157857f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 222)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 223)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 224)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 225)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -0.1f;
Packit 67cb25
   float b = -1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.00498756211f;
Packit 67cb25
   float z_expected = 10.0498756211f;
Packit 67cb25
   float c_expected = 0.099503719021f;
Packit 67cb25
   float s_expected = 0.99503719021f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 226)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 227)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 228)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 229)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -0.1f;
Packit 67cb25
   float b = -0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -0.141421356237f;
Packit 67cb25
   float z_expected = 1.41421356237f;
Packit 67cb25
   float c_expected = 0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 230)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 231)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 232)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 233)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -0.1f;
Packit 67cb25
   float b = 0.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -0.1f;
Packit 67cb25
   float z_expected = -0.0f;
Packit 67cb25
   float c_expected = 1.0f;
Packit 67cb25
   float s_expected = -0.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 234)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 235)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 236)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 237)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -0.1f;
Packit 67cb25
   float b = 0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 0.141421356237f;
Packit 67cb25
   float z_expected = -1.41421356237f;
Packit 67cb25
   float c_expected = -0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 238)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 239)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 240)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 241)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -0.1f;
Packit 67cb25
   float b = 1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.00498756211f;
Packit 67cb25
   float z_expected = -10.0498756211f;
Packit 67cb25
   float c_expected = -0.099503719021f;
Packit 67cb25
   float s_expected = 0.99503719021f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 242)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 243)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 244)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 245)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = -0.1f;
Packit 67cb25
   float b = 1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.50332963784f;
Packit 67cb25
   float z_expected = -15.0332963784f;
Packit 67cb25
   float c_expected = -0.0665190105238f;
Packit 67cb25
   float s_expected = 0.997785157857f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 246)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 247)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 248)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 249)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.0f;
Packit 67cb25
   float b = -1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.5f;
Packit 67cb25
   float z_expected = 1.0f;
Packit 67cb25
   float c_expected = -0.0f;
Packit 67cb25
   float s_expected = 1.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 250)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 251)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 252)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 253)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.0f;
Packit 67cb25
   float b = -1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.0f;
Packit 67cb25
   float z_expected = 1.0f;
Packit 67cb25
   float c_expected = -0.0f;
Packit 67cb25
   float s_expected = 1.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 254)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 255)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 256)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 257)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.0f;
Packit 67cb25
   float b = -0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -0.1f;
Packit 67cb25
   float z_expected = 1.0f;
Packit 67cb25
   float c_expected = -0.0f;
Packit 67cb25
   float s_expected = 1.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 258)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 259)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 260)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 261)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.0f;
Packit 67cb25
   float b = 0.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 0.0f;
Packit 67cb25
   float z_expected = 0.0f;
Packit 67cb25
   float c_expected = 1.0f;
Packit 67cb25
   float s_expected = 0.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 262)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 263)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 264)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 265)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.0f;
Packit 67cb25
   float b = 0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 0.1f;
Packit 67cb25
   float z_expected = 1.0f;
Packit 67cb25
   float c_expected = 0.0f;
Packit 67cb25
   float s_expected = 1.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 266)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 267)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 268)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 269)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.0f;
Packit 67cb25
   float b = 1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.0f;
Packit 67cb25
   float z_expected = 1.0f;
Packit 67cb25
   float c_expected = 0.0f;
Packit 67cb25
   float s_expected = 1.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 270)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 271)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 272)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 273)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.0f;
Packit 67cb25
   float b = 1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.5f;
Packit 67cb25
   float z_expected = 1.0f;
Packit 67cb25
   float c_expected = 0.0f;
Packit 67cb25
   float s_expected = 1.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 274)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 275)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 276)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 277)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.1f;
Packit 67cb25
   float b = -1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.50332963784f;
Packit 67cb25
   float z_expected = -15.0332963784f;
Packit 67cb25
   float c_expected = -0.0665190105238f;
Packit 67cb25
   float s_expected = 0.997785157857f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 278)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 279)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 280)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 281)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.1f;
Packit 67cb25
   float b = -1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.00498756211f;
Packit 67cb25
   float z_expected = -10.0498756211f;
Packit 67cb25
   float c_expected = -0.099503719021f;
Packit 67cb25
   float s_expected = 0.99503719021f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 282)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 283)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 284)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 285)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.1f;
Packit 67cb25
   float b = -0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -0.141421356237f;
Packit 67cb25
   float z_expected = -1.41421356237f;
Packit 67cb25
   float c_expected = -0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 286)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 287)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 288)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 289)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.1f;
Packit 67cb25
   float b = 0.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 0.1f;
Packit 67cb25
   float z_expected = 0.0f;
Packit 67cb25
   float c_expected = 1.0f;
Packit 67cb25
   float s_expected = 0.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 290)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 291)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 292)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 293)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.1f;
Packit 67cb25
   float b = 0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 0.141421356237f;
Packit 67cb25
   float z_expected = 1.41421356237f;
Packit 67cb25
   float c_expected = 0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 294)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 295)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 296)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 297)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.1f;
Packit 67cb25
   float b = 1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.00498756211f;
Packit 67cb25
   float z_expected = 10.0498756211f;
Packit 67cb25
   float c_expected = 0.099503719021f;
Packit 67cb25
   float s_expected = 0.99503719021f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 298)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 299)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 300)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 301)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 0.1f;
Packit 67cb25
   float b = 1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.50332963784f;
Packit 67cb25
   float z_expected = 15.0332963784f;
Packit 67cb25
   float c_expected = 0.0665190105238f;
Packit 67cb25
   float s_expected = 0.997785157857f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 302)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 303)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 304)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 305)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.0f;
Packit 67cb25
   float b = -1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.80277563773f;
Packit 67cb25
   float z_expected = -1.80277563773f;
Packit 67cb25
   float c_expected = -0.554700196225f;
Packit 67cb25
   float s_expected = 0.832050294338f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 306)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 307)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 308)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 309)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.0f;
Packit 67cb25
   float b = -1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -1.41421356237f;
Packit 67cb25
   float z_expected = -1.41421356237f;
Packit 67cb25
   float c_expected = -0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 310)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 311)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 312)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 313)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.0f;
Packit 67cb25
   float b = -0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.00498756211f;
Packit 67cb25
   float z_expected = -0.099503719021f;
Packit 67cb25
   float c_expected = 0.99503719021f;
Packit 67cb25
   float s_expected = -0.099503719021f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 314)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 315)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 316)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 317)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.0f;
Packit 67cb25
   float b = 0.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.0f;
Packit 67cb25
   float z_expected = 0.0f;
Packit 67cb25
   float c_expected = 1.0f;
Packit 67cb25
   float s_expected = 0.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 318)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 319)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 320)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 321)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.0f;
Packit 67cb25
   float b = 0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.00498756211f;
Packit 67cb25
   float z_expected = 0.099503719021f;
Packit 67cb25
   float c_expected = 0.99503719021f;
Packit 67cb25
   float s_expected = 0.099503719021f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 322)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 323)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 324)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 325)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.0f;
Packit 67cb25
   float b = 1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.41421356237f;
Packit 67cb25
   float z_expected = 1.41421356237f;
Packit 67cb25
   float c_expected = 0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 326)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 327)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 328)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 329)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.0f;
Packit 67cb25
   float b = 1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.80277563773f;
Packit 67cb25
   float z_expected = 1.80277563773f;
Packit 67cb25
   float c_expected = 0.554700196225f;
Packit 67cb25
   float s_expected = 0.832050294338f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 330)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 331)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 332)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 333)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.5f;
Packit 67cb25
   float b = -1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = -2.12132034356f;
Packit 67cb25
   float z_expected = -1.41421356237f;
Packit 67cb25
   float c_expected = -0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 334)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 335)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 336)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 337)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.5f;
Packit 67cb25
   float b = -1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.80277563773f;
Packit 67cb25
   float z_expected = -0.554700196225f;
Packit 67cb25
   float c_expected = 0.832050294338f;
Packit 67cb25
   float s_expected = -0.554700196225f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 338)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 339)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 340)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 341)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.5f;
Packit 67cb25
   float b = -0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.50332963784f;
Packit 67cb25
   float z_expected = -0.0665190105238f;
Packit 67cb25
   float c_expected = 0.997785157857f;
Packit 67cb25
   float s_expected = -0.0665190105238f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 342)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 343)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 344)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 345)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.5f;
Packit 67cb25
   float b = 0.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.5f;
Packit 67cb25
   float z_expected = 0.0f;
Packit 67cb25
   float c_expected = 1.0f;
Packit 67cb25
   float s_expected = 0.0f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 346)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 347)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 348)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 349)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.5f;
Packit 67cb25
   float b = 0.1f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.50332963784f;
Packit 67cb25
   float z_expected = 0.0665190105238f;
Packit 67cb25
   float c_expected = 0.997785157857f;
Packit 67cb25
   float s_expected = 0.0665190105238f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 350)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 351)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 352)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 353)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.5f;
Packit 67cb25
   float b = 1.0f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 1.80277563773f;
Packit 67cb25
   float z_expected = 0.554700196225f;
Packit 67cb25
   float c_expected = 0.832050294338f;
Packit 67cb25
   float s_expected = 0.554700196225f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 354)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 355)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 356)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 357)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   float a = 1.5f;
Packit 67cb25
   float b = 1.5f;
Packit 67cb25
   float c;
Packit 67cb25
   float s;
Packit 67cb25
   float r_expected = 2.12132034356f;
Packit 67cb25
   float z_expected = 1.41421356237f;
Packit 67cb25
   float c_expected = 0.707106781187f;
Packit 67cb25
   float s_expected = 0.707106781187f;
Packit 67cb25
   cblas_srotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, flteps, "srotg(case 358)");
Packit 67cb25
   gsl_test_rel(b, z_expected, flteps, "srotg(case 359)");
Packit 67cb25
   gsl_test_rel(c, c_expected, flteps, "srotg(case 360)");
Packit 67cb25
   gsl_test_rel(s, s_expected, flteps, "srotg(case 361)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1.5;
Packit 67cb25
   double b = -1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -2.12132034356;
Packit 67cb25
   double z_expected = 1.41421356237;
Packit 67cb25
   double c_expected = 0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 362)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 363)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 364)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 365)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1.5;
Packit 67cb25
   double b = -1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.80277563773;
Packit 67cb25
   double z_expected = 0.554700196225;
Packit 67cb25
   double c_expected = 0.832050294338;
Packit 67cb25
   double s_expected = 0.554700196225;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 366)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 367)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 368)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 369)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1.5;
Packit 67cb25
   double b = -0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.50332963784;
Packit 67cb25
   double z_expected = 0.0665190105238;
Packit 67cb25
   double c_expected = 0.997785157857;
Packit 67cb25
   double s_expected = 0.0665190105238;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 370)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 371)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 372)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 373)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1.5;
Packit 67cb25
   double b = 0;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.5;
Packit 67cb25
   double z_expected = -0;
Packit 67cb25
   double c_expected = 1;
Packit 67cb25
   double s_expected = -0;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 374)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 375)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 376)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 377)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1.5;
Packit 67cb25
   double b = 0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.50332963784;
Packit 67cb25
   double z_expected = -0.0665190105238;
Packit 67cb25
   double c_expected = 0.997785157857;
Packit 67cb25
   double s_expected = -0.0665190105238;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 378)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 379)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 380)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 381)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1.5;
Packit 67cb25
   double b = 1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.80277563773;
Packit 67cb25
   double z_expected = -0.554700196225;
Packit 67cb25
   double c_expected = 0.832050294338;
Packit 67cb25
   double s_expected = -0.554700196225;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 382)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 383)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 384)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 385)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1.5;
Packit 67cb25
   double b = 1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 2.12132034356;
Packit 67cb25
   double z_expected = -1.41421356237;
Packit 67cb25
   double c_expected = -0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 386)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 387)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 388)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 389)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1;
Packit 67cb25
   double b = -1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.80277563773;
Packit 67cb25
   double z_expected = 1.80277563773;
Packit 67cb25
   double c_expected = 0.554700196225;
Packit 67cb25
   double s_expected = 0.832050294338;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 390)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 391)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 392)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 393)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1;
Packit 67cb25
   double b = -1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.41421356237;
Packit 67cb25
   double z_expected = 1.41421356237;
Packit 67cb25
   double c_expected = 0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 394)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 395)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 396)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 397)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1;
Packit 67cb25
   double b = -0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.00498756211;
Packit 67cb25
   double z_expected = 0.099503719021;
Packit 67cb25
   double c_expected = 0.99503719021;
Packit 67cb25
   double s_expected = 0.099503719021;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 398)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 399)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 400)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 401)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1;
Packit 67cb25
   double b = 0;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1;
Packit 67cb25
   double z_expected = -0;
Packit 67cb25
   double c_expected = 1;
Packit 67cb25
   double s_expected = -0;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 402)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 403)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 404)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 405)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1;
Packit 67cb25
   double b = 0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.00498756211;
Packit 67cb25
   double z_expected = -0.099503719021;
Packit 67cb25
   double c_expected = 0.99503719021;
Packit 67cb25
   double s_expected = -0.099503719021;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 406)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 407)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 408)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 409)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1;
Packit 67cb25
   double b = 1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.41421356237;
Packit 67cb25
   double z_expected = -1.41421356237;
Packit 67cb25
   double c_expected = -0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 410)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 411)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 412)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 413)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -1;
Packit 67cb25
   double b = 1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.80277563773;
Packit 67cb25
   double z_expected = -1.80277563773;
Packit 67cb25
   double c_expected = -0.554700196225;
Packit 67cb25
   double s_expected = 0.832050294338;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 414)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 415)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 416)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 417)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -0.1;
Packit 67cb25
   double b = -1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.50332963784;
Packit 67cb25
   double z_expected = 15.0332963784;
Packit 67cb25
   double c_expected = 0.0665190105238;
Packit 67cb25
   double s_expected = 0.997785157857;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 418)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 419)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 420)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 421)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -0.1;
Packit 67cb25
   double b = -1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.00498756211;
Packit 67cb25
   double z_expected = 10.0498756211;
Packit 67cb25
   double c_expected = 0.099503719021;
Packit 67cb25
   double s_expected = 0.99503719021;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 422)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 423)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 424)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 425)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -0.1;
Packit 67cb25
   double b = -0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -0.141421356237;
Packit 67cb25
   double z_expected = 1.41421356237;
Packit 67cb25
   double c_expected = 0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 426)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 427)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 428)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 429)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -0.1;
Packit 67cb25
   double b = 0;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -0.1;
Packit 67cb25
   double z_expected = -0;
Packit 67cb25
   double c_expected = 1;
Packit 67cb25
   double s_expected = -0;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 430)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 431)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 432)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 433)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -0.1;
Packit 67cb25
   double b = 0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 0.141421356237;
Packit 67cb25
   double z_expected = -1.41421356237;
Packit 67cb25
   double c_expected = -0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 434)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 435)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 436)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 437)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -0.1;
Packit 67cb25
   double b = 1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.00498756211;
Packit 67cb25
   double z_expected = -10.0498756211;
Packit 67cb25
   double c_expected = -0.099503719021;
Packit 67cb25
   double s_expected = 0.99503719021;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 438)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 439)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 440)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 441)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = -0.1;
Packit 67cb25
   double b = 1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.50332963784;
Packit 67cb25
   double z_expected = -15.0332963784;
Packit 67cb25
   double c_expected = -0.0665190105238;
Packit 67cb25
   double s_expected = 0.997785157857;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 442)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 443)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 444)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 445)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0;
Packit 67cb25
   double b = -1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.5;
Packit 67cb25
   double z_expected = 1;
Packit 67cb25
   double c_expected = -0;
Packit 67cb25
   double s_expected = 1;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 446)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 447)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 448)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 449)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0;
Packit 67cb25
   double b = -1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1;
Packit 67cb25
   double z_expected = 1;
Packit 67cb25
   double c_expected = -0;
Packit 67cb25
   double s_expected = 1;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 450)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 451)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 452)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 453)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0;
Packit 67cb25
   double b = -0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -0.1;
Packit 67cb25
   double z_expected = 1;
Packit 67cb25
   double c_expected = -0;
Packit 67cb25
   double s_expected = 1;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 454)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 455)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 456)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 457)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0;
Packit 67cb25
   double b = 0;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 0;
Packit 67cb25
   double z_expected = 0;
Packit 67cb25
   double c_expected = 1;
Packit 67cb25
   double s_expected = 0;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 458)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 459)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 460)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 461)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0;
Packit 67cb25
   double b = 0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 0.1;
Packit 67cb25
   double z_expected = 1;
Packit 67cb25
   double c_expected = 0;
Packit 67cb25
   double s_expected = 1;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 462)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 463)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 464)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 465)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0;
Packit 67cb25
   double b = 1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1;
Packit 67cb25
   double z_expected = 1;
Packit 67cb25
   double c_expected = 0;
Packit 67cb25
   double s_expected = 1;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 466)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 467)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 468)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 469)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0;
Packit 67cb25
   double b = 1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.5;
Packit 67cb25
   double z_expected = 1;
Packit 67cb25
   double c_expected = 0;
Packit 67cb25
   double s_expected = 1;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 470)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 471)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 472)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 473)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0.1;
Packit 67cb25
   double b = -1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.50332963784;
Packit 67cb25
   double z_expected = -15.0332963784;
Packit 67cb25
   double c_expected = -0.0665190105238;
Packit 67cb25
   double s_expected = 0.997785157857;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 474)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 475)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 476)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 477)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0.1;
Packit 67cb25
   double b = -1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.00498756211;
Packit 67cb25
   double z_expected = -10.0498756211;
Packit 67cb25
   double c_expected = -0.099503719021;
Packit 67cb25
   double s_expected = 0.99503719021;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 478)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 479)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 480)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 481)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0.1;
Packit 67cb25
   double b = -0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -0.141421356237;
Packit 67cb25
   double z_expected = -1.41421356237;
Packit 67cb25
   double c_expected = -0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 482)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 483)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 484)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 485)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0.1;
Packit 67cb25
   double b = 0;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 0.1;
Packit 67cb25
   double z_expected = 0;
Packit 67cb25
   double c_expected = 1;
Packit 67cb25
   double s_expected = 0;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 486)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 487)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 488)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 489)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0.1;
Packit 67cb25
   double b = 0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 0.141421356237;
Packit 67cb25
   double z_expected = 1.41421356237;
Packit 67cb25
   double c_expected = 0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 490)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 491)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 492)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 493)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0.1;
Packit 67cb25
   double b = 1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.00498756211;
Packit 67cb25
   double z_expected = 10.0498756211;
Packit 67cb25
   double c_expected = 0.099503719021;
Packit 67cb25
   double s_expected = 0.99503719021;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 494)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 495)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 496)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 497)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 0.1;
Packit 67cb25
   double b = 1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.50332963784;
Packit 67cb25
   double z_expected = 15.0332963784;
Packit 67cb25
   double c_expected = 0.0665190105238;
Packit 67cb25
   double s_expected = 0.997785157857;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 498)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 499)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 500)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 501)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1;
Packit 67cb25
   double b = -1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.80277563773;
Packit 67cb25
   double z_expected = -1.80277563773;
Packit 67cb25
   double c_expected = -0.554700196225;
Packit 67cb25
   double s_expected = 0.832050294338;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 502)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 503)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 504)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 505)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1;
Packit 67cb25
   double b = -1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -1.41421356237;
Packit 67cb25
   double z_expected = -1.41421356237;
Packit 67cb25
   double c_expected = -0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 506)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 507)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 508)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 509)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1;
Packit 67cb25
   double b = -0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.00498756211;
Packit 67cb25
   double z_expected = -0.099503719021;
Packit 67cb25
   double c_expected = 0.99503719021;
Packit 67cb25
   double s_expected = -0.099503719021;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 510)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 511)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 512)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 513)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1;
Packit 67cb25
   double b = 0;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1;
Packit 67cb25
   double z_expected = 0;
Packit 67cb25
   double c_expected = 1;
Packit 67cb25
   double s_expected = 0;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 514)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 515)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 516)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 517)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1;
Packit 67cb25
   double b = 0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.00498756211;
Packit 67cb25
   double z_expected = 0.099503719021;
Packit 67cb25
   double c_expected = 0.99503719021;
Packit 67cb25
   double s_expected = 0.099503719021;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 518)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 519)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 520)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 521)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1;
Packit 67cb25
   double b = 1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.41421356237;
Packit 67cb25
   double z_expected = 1.41421356237;
Packit 67cb25
   double c_expected = 0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 522)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 523)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 524)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 525)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1;
Packit 67cb25
   double b = 1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.80277563773;
Packit 67cb25
   double z_expected = 1.80277563773;
Packit 67cb25
   double c_expected = 0.554700196225;
Packit 67cb25
   double s_expected = 0.832050294338;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 526)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 527)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 528)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 529)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1.5;
Packit 67cb25
   double b = -1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = -2.12132034356;
Packit 67cb25
   double z_expected = -1.41421356237;
Packit 67cb25
   double c_expected = -0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 530)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 531)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 532)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 533)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1.5;
Packit 67cb25
   double b = -1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.80277563773;
Packit 67cb25
   double z_expected = -0.554700196225;
Packit 67cb25
   double c_expected = 0.832050294338;
Packit 67cb25
   double s_expected = -0.554700196225;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 534)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 535)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 536)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 537)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1.5;
Packit 67cb25
   double b = -0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.50332963784;
Packit 67cb25
   double z_expected = -0.0665190105238;
Packit 67cb25
   double c_expected = 0.997785157857;
Packit 67cb25
   double s_expected = -0.0665190105238;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 538)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 539)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 540)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 541)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1.5;
Packit 67cb25
   double b = 0;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.5;
Packit 67cb25
   double z_expected = 0;
Packit 67cb25
   double c_expected = 1;
Packit 67cb25
   double s_expected = 0;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 542)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 543)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 544)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 545)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1.5;
Packit 67cb25
   double b = 0.1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.50332963784;
Packit 67cb25
   double z_expected = 0.0665190105238;
Packit 67cb25
   double c_expected = 0.997785157857;
Packit 67cb25
   double s_expected = 0.0665190105238;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 546)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 547)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 548)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 549)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1.5;
Packit 67cb25
   double b = 1;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 1.80277563773;
Packit 67cb25
   double z_expected = 0.554700196225;
Packit 67cb25
   double c_expected = 0.832050294338;
Packit 67cb25
   double s_expected = 0.554700196225;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 550)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 551)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 552)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 553)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   double a = 1.5;
Packit 67cb25
   double b = 1.5;
Packit 67cb25
   double c;
Packit 67cb25
   double s;
Packit 67cb25
   double r_expected = 2.12132034356;
Packit 67cb25
   double z_expected = 1.41421356237;
Packit 67cb25
   double c_expected = 0.707106781187;
Packit 67cb25
   double s_expected = 0.707106781187;
Packit 67cb25
   cblas_drotg(&a, &b, &c, &s);
Packit 67cb25
   gsl_test_rel(a, r_expected, dbleps, "drotg(case 554)");
Packit 67cb25
   gsl_test_rel(b, z_expected, dbleps, "drotg(case 555)");
Packit 67cb25
   gsl_test_rel(c, c_expected, dbleps, "drotg(case 556)");
Packit 67cb25
   gsl_test_rel(s, s_expected, dbleps, "drotg(case 557)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
}