|
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_rot (void) {
|
|
Packit |
67cb25 |
const double flteps = 1e-4, dbleps = 1e-6;
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = 0.0f;
|
|
Packit |
67cb25 |
float s = 0.0f;
|
|
Packit |
67cb25 |
float X[] = { -0.314f };
|
|
Packit |
67cb25 |
int incX = 1;
|
|
Packit |
67cb25 |
float Y[] = { -0.406f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float x_expected[] = { 0.0f };
|
|
Packit |
67cb25 |
float y_expected[] = { 0.0f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 558)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 559)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = 0.866025403784f;
|
|
Packit |
67cb25 |
float s = 0.5f;
|
|
Packit |
67cb25 |
float X[] = { -0.314f };
|
|
Packit |
67cb25 |
int incX = 1;
|
|
Packit |
67cb25 |
float Y[] = { -0.406f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float x_expected[] = { -0.474932f };
|
|
Packit |
67cb25 |
float y_expected[] = { -0.194606f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 560)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 561)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = 0.0f;
|
|
Packit |
67cb25 |
float s = -1.0f;
|
|
Packit |
67cb25 |
float X[] = { -0.314f };
|
|
Packit |
67cb25 |
int incX = 1;
|
|
Packit |
67cb25 |
float Y[] = { -0.406f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float x_expected[] = { 0.406f };
|
|
Packit |
67cb25 |
float y_expected[] = { -0.314f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 562)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 563)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = -1.0f;
|
|
Packit |
67cb25 |
float s = 0.0f;
|
|
Packit |
67cb25 |
float X[] = { -0.314f };
|
|
Packit |
67cb25 |
int incX = 1;
|
|
Packit |
67cb25 |
float Y[] = { -0.406f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float x_expected[] = { 0.314f };
|
|
Packit |
67cb25 |
float y_expected[] = { 0.406f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 564)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 565)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = 0;
|
|
Packit |
67cb25 |
double s = 0;
|
|
Packit |
67cb25 |
double X[] = { -0.493 };
|
|
Packit |
67cb25 |
int incX = 1;
|
|
Packit |
67cb25 |
double Y[] = { -0.014 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double x_expected[] = { 0.0 };
|
|
Packit |
67cb25 |
double y_expected[] = { 0.0 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 566)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 567)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = 0.866025403784;
|
|
Packit |
67cb25 |
double s = 0.5;
|
|
Packit |
67cb25 |
double X[] = { -0.493 };
|
|
Packit |
67cb25 |
int incX = 1;
|
|
Packit |
67cb25 |
double Y[] = { -0.014 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double x_expected[] = { -0.433950524066 };
|
|
Packit |
67cb25 |
double y_expected[] = { 0.234375644347 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 568)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 569)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = 0;
|
|
Packit |
67cb25 |
double s = -1;
|
|
Packit |
67cb25 |
double X[] = { -0.493 };
|
|
Packit |
67cb25 |
int incX = 1;
|
|
Packit |
67cb25 |
double Y[] = { -0.014 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double x_expected[] = { 0.014 };
|
|
Packit |
67cb25 |
double y_expected[] = { -0.493 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 570)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 571)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = -1;
|
|
Packit |
67cb25 |
double s = 0;
|
|
Packit |
67cb25 |
double X[] = { -0.493 };
|
|
Packit |
67cb25 |
int incX = 1;
|
|
Packit |
67cb25 |
double Y[] = { -0.014 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double x_expected[] = { 0.493 };
|
|
Packit |
67cb25 |
double y_expected[] = { 0.014 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 572)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 573)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = 0.0f;
|
|
Packit |
67cb25 |
float s = 0.0f;
|
|
Packit |
67cb25 |
float X[] = { -0.808f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { -0.511f };
|
|
Packit |
67cb25 |
int incY = 1;
|
|
Packit |
67cb25 |
float x_expected[] = { 0.0f };
|
|
Packit |
67cb25 |
float y_expected[] = { 0.0f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 574)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 575)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = 0.866025403784f;
|
|
Packit |
67cb25 |
float s = 0.5f;
|
|
Packit |
67cb25 |
float X[] = { -0.808f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { -0.511f };
|
|
Packit |
67cb25 |
int incY = 1;
|
|
Packit |
67cb25 |
float x_expected[] = { -0.955249f };
|
|
Packit |
67cb25 |
float y_expected[] = { -0.038539f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 576)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 577)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = 0.0f;
|
|
Packit |
67cb25 |
float s = -1.0f;
|
|
Packit |
67cb25 |
float X[] = { -0.808f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { -0.511f };
|
|
Packit |
67cb25 |
int incY = 1;
|
|
Packit |
67cb25 |
float x_expected[] = { 0.511f };
|
|
Packit |
67cb25 |
float y_expected[] = { -0.808f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 578)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 579)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = -1.0f;
|
|
Packit |
67cb25 |
float s = 0.0f;
|
|
Packit |
67cb25 |
float X[] = { -0.808f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { -0.511f };
|
|
Packit |
67cb25 |
int incY = 1;
|
|
Packit |
67cb25 |
float x_expected[] = { 0.808f };
|
|
Packit |
67cb25 |
float y_expected[] = { 0.511f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 580)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 581)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = 0;
|
|
Packit |
67cb25 |
double s = 0;
|
|
Packit |
67cb25 |
double X[] = { -0.176 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { -0.165 };
|
|
Packit |
67cb25 |
int incY = 1;
|
|
Packit |
67cb25 |
double x_expected[] = { 0.0 };
|
|
Packit |
67cb25 |
double y_expected[] = { 0.0 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 582)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 583)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = 0.866025403784;
|
|
Packit |
67cb25 |
double s = 0.5;
|
|
Packit |
67cb25 |
double X[] = { -0.176 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { -0.165 };
|
|
Packit |
67cb25 |
int incY = 1;
|
|
Packit |
67cb25 |
double x_expected[] = { -0.234920471066 };
|
|
Packit |
67cb25 |
double y_expected[] = { -0.0548941916244 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 584)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 585)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = 0;
|
|
Packit |
67cb25 |
double s = -1;
|
|
Packit |
67cb25 |
double X[] = { -0.176 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { -0.165 };
|
|
Packit |
67cb25 |
int incY = 1;
|
|
Packit |
67cb25 |
double x_expected[] = { 0.165 };
|
|
Packit |
67cb25 |
double y_expected[] = { -0.176 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 586)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 587)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = -1;
|
|
Packit |
67cb25 |
double s = 0;
|
|
Packit |
67cb25 |
double X[] = { -0.176 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { -0.165 };
|
|
Packit |
67cb25 |
int incY = 1;
|
|
Packit |
67cb25 |
double x_expected[] = { 0.176 };
|
|
Packit |
67cb25 |
double y_expected[] = { 0.165 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 588)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 589)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = 0.0f;
|
|
Packit |
67cb25 |
float s = 0.0f;
|
|
Packit |
67cb25 |
float X[] = { -0.201f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.087f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float x_expected[] = { 0.0f };
|
|
Packit |
67cb25 |
float y_expected[] = { 0.0f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 590)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 591)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = 0.866025403784f;
|
|
Packit |
67cb25 |
float s = 0.5f;
|
|
Packit |
67cb25 |
float X[] = { -0.201f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.087f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float x_expected[] = { -0.130571f };
|
|
Packit |
67cb25 |
float y_expected[] = { 0.175844f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 592)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 593)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = 0.0f;
|
|
Packit |
67cb25 |
float s = -1.0f;
|
|
Packit |
67cb25 |
float X[] = { -0.201f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.087f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float x_expected[] = { -0.087f };
|
|
Packit |
67cb25 |
float y_expected[] = { -0.201f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 594)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 595)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
float c = -1.0f;
|
|
Packit |
67cb25 |
float s = 0.0f;
|
|
Packit |
67cb25 |
float X[] = { -0.201f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.087f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float x_expected[] = { 0.201f };
|
|
Packit |
67cb25 |
float y_expected[] = { -0.087f };
|
|
Packit |
67cb25 |
cblas_srot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 596)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 597)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = 0;
|
|
Packit |
67cb25 |
double s = 0;
|
|
Packit |
67cb25 |
double X[] = { -0.464 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.7 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double x_expected[] = { 0.0 };
|
|
Packit |
67cb25 |
double y_expected[] = { 0.0 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 598)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 599)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = 0.866025403784;
|
|
Packit |
67cb25 |
double s = 0.5;
|
|
Packit |
67cb25 |
double X[] = { -0.464 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.7 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double x_expected[] = { -0.051835787356 };
|
|
Packit |
67cb25 |
double y_expected[] = { 0.838217782649 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 600)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 601)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = 0;
|
|
Packit |
67cb25 |
double s = -1;
|
|
Packit |
67cb25 |
double X[] = { -0.464 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.7 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double x_expected[] = { -0.7 };
|
|
Packit |
67cb25 |
double y_expected[] = { -0.464 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 602)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 603)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
double c = -1;
|
|
Packit |
67cb25 |
double s = 0;
|
|
Packit |
67cb25 |
double X[] = { -0.464 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.7 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double x_expected[] = { 0.464 };
|
|
Packit |
67cb25 |
double y_expected[] = { -0.7 };
|
|
Packit |
67cb25 |
cblas_drot(N, X, incX, Y, incY, c, s);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 604)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 605)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
}
|