|
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 |
}
|