diff options
Diffstat (limited to 'gsl-1.9/cblas/test_rot.c')
-rw-r--r-- | gsl-1.9/cblas/test_rot.c | 635 |
1 files changed, 635 insertions, 0 deletions
diff --git a/gsl-1.9/cblas/test_rot.c b/gsl-1.9/cblas/test_rot.c new file mode 100644 index 0000000..d2d6ed3 --- /dev/null +++ b/gsl-1.9/cblas/test_rot.c @@ -0,0 +1,635 @@ +#include <gsl/gsl_test.h> +#include <gsl/gsl_ieee_utils.h> +#include <gsl/gsl_math.h> +#include <gsl/gsl_cblas.h> + +#include "tests.h" + +void +test_rot (void) { +const double flteps = 1e-4, dbleps = 1e-6; + { + int N = 1; + float c = 0.0f; + float s = 0.0f; + float X[] = { -0.314f }; + int incX = 1; + float Y[] = { -0.406f }; + int incY = -1; + float x_expected[] = { 0.0f }; + float y_expected[] = { 0.0f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 558)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 559)"); + } + }; + }; + + + { + int N = 1; + float c = 0.866025403784f; + float s = 0.5f; + float X[] = { -0.314f }; + int incX = 1; + float Y[] = { -0.406f }; + int incY = -1; + float x_expected[] = { -0.474932f }; + float y_expected[] = { -0.194606f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 560)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 561)"); + } + }; + }; + + + { + int N = 1; + float c = 0.0f; + float s = -1.0f; + float X[] = { -0.314f }; + int incX = 1; + float Y[] = { -0.406f }; + int incY = -1; + float x_expected[] = { 0.406f }; + float y_expected[] = { -0.314f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 562)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 563)"); + } + }; + }; + + + { + int N = 1; + float c = -1.0f; + float s = 0.0f; + float X[] = { -0.314f }; + int incX = 1; + float Y[] = { -0.406f }; + int incY = -1; + float x_expected[] = { 0.314f }; + float y_expected[] = { 0.406f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 564)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 565)"); + } + }; + }; + + + { + int N = 1; + double c = 0; + double s = 0; + double X[] = { -0.493 }; + int incX = 1; + double Y[] = { -0.014 }; + int incY = -1; + double x_expected[] = { 0.0 }; + double y_expected[] = { 0.0 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 566)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 567)"); + } + }; + }; + + + { + int N = 1; + double c = 0.866025403784; + double s = 0.5; + double X[] = { -0.493 }; + int incX = 1; + double Y[] = { -0.014 }; + int incY = -1; + double x_expected[] = { -0.433950524066 }; + double y_expected[] = { 0.234375644347 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 568)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 569)"); + } + }; + }; + + + { + int N = 1; + double c = 0; + double s = -1; + double X[] = { -0.493 }; + int incX = 1; + double Y[] = { -0.014 }; + int incY = -1; + double x_expected[] = { 0.014 }; + double y_expected[] = { -0.493 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 570)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 571)"); + } + }; + }; + + + { + int N = 1; + double c = -1; + double s = 0; + double X[] = { -0.493 }; + int incX = 1; + double Y[] = { -0.014 }; + int incY = -1; + double x_expected[] = { 0.493 }; + double y_expected[] = { 0.014 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 572)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 573)"); + } + }; + }; + + + { + int N = 1; + float c = 0.0f; + float s = 0.0f; + float X[] = { -0.808f }; + int incX = -1; + float Y[] = { -0.511f }; + int incY = 1; + float x_expected[] = { 0.0f }; + float y_expected[] = { 0.0f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 574)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 575)"); + } + }; + }; + + + { + int N = 1; + float c = 0.866025403784f; + float s = 0.5f; + float X[] = { -0.808f }; + int incX = -1; + float Y[] = { -0.511f }; + int incY = 1; + float x_expected[] = { -0.955249f }; + float y_expected[] = { -0.038539f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 576)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 577)"); + } + }; + }; + + + { + int N = 1; + float c = 0.0f; + float s = -1.0f; + float X[] = { -0.808f }; + int incX = -1; + float Y[] = { -0.511f }; + int incY = 1; + float x_expected[] = { 0.511f }; + float y_expected[] = { -0.808f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 578)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 579)"); + } + }; + }; + + + { + int N = 1; + float c = -1.0f; + float s = 0.0f; + float X[] = { -0.808f }; + int incX = -1; + float Y[] = { -0.511f }; + int incY = 1; + float x_expected[] = { 0.808f }; + float y_expected[] = { 0.511f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 580)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 581)"); + } + }; + }; + + + { + int N = 1; + double c = 0; + double s = 0; + double X[] = { -0.176 }; + int incX = -1; + double Y[] = { -0.165 }; + int incY = 1; + double x_expected[] = { 0.0 }; + double y_expected[] = { 0.0 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 582)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 583)"); + } + }; + }; + + + { + int N = 1; + double c = 0.866025403784; + double s = 0.5; + double X[] = { -0.176 }; + int incX = -1; + double Y[] = { -0.165 }; + int incY = 1; + double x_expected[] = { -0.234920471066 }; + double y_expected[] = { -0.0548941916244 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 584)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 585)"); + } + }; + }; + + + { + int N = 1; + double c = 0; + double s = -1; + double X[] = { -0.176 }; + int incX = -1; + double Y[] = { -0.165 }; + int incY = 1; + double x_expected[] = { 0.165 }; + double y_expected[] = { -0.176 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 586)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 587)"); + } + }; + }; + + + { + int N = 1; + double c = -1; + double s = 0; + double X[] = { -0.176 }; + int incX = -1; + double Y[] = { -0.165 }; + int incY = 1; + double x_expected[] = { 0.176 }; + double y_expected[] = { 0.165 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 588)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 589)"); + } + }; + }; + + + { + int N = 1; + float c = 0.0f; + float s = 0.0f; + float X[] = { -0.201f }; + int incX = -1; + float Y[] = { 0.087f }; + int incY = -1; + float x_expected[] = { 0.0f }; + float y_expected[] = { 0.0f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 590)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 591)"); + } + }; + }; + + + { + int N = 1; + float c = 0.866025403784f; + float s = 0.5f; + float X[] = { -0.201f }; + int incX = -1; + float Y[] = { 0.087f }; + int incY = -1; + float x_expected[] = { -0.130571f }; + float y_expected[] = { 0.175844f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 592)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 593)"); + } + }; + }; + + + { + int N = 1; + float c = 0.0f; + float s = -1.0f; + float X[] = { -0.201f }; + int incX = -1; + float Y[] = { 0.087f }; + int incY = -1; + float x_expected[] = { -0.087f }; + float y_expected[] = { -0.201f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 594)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 595)"); + } + }; + }; + + + { + int N = 1; + float c = -1.0f; + float s = 0.0f; + float X[] = { -0.201f }; + int incX = -1; + float Y[] = { 0.087f }; + int incY = -1; + float x_expected[] = { 0.201f }; + float y_expected[] = { -0.087f }; + cblas_srot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 596)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 597)"); + } + }; + }; + + + { + int N = 1; + double c = 0; + double s = 0; + double X[] = { -0.464 }; + int incX = -1; + double Y[] = { 0.7 }; + int incY = -1; + double x_expected[] = { 0.0 }; + double y_expected[] = { 0.0 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 598)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 599)"); + } + }; + }; + + + { + int N = 1; + double c = 0.866025403784; + double s = 0.5; + double X[] = { -0.464 }; + int incX = -1; + double Y[] = { 0.7 }; + int incY = -1; + double x_expected[] = { -0.051835787356 }; + double y_expected[] = { 0.838217782649 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 600)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 601)"); + } + }; + }; + + + { + int N = 1; + double c = 0; + double s = -1; + double X[] = { -0.464 }; + int incX = -1; + double Y[] = { 0.7 }; + int incY = -1; + double x_expected[] = { -0.7 }; + double y_expected[] = { -0.464 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 602)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 603)"); + } + }; + }; + + + { + int N = 1; + double c = -1; + double s = 0; + double X[] = { -0.464 }; + int incX = -1; + double Y[] = { 0.7 }; + int incY = -1; + double x_expected[] = { 0.464 }; + double y_expected[] = { -0.7 }; + cblas_drot(N, X, incX, Y, incY, c, s); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 604)"); + } + }; + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 605)"); + } + }; + }; + + +} |