diff options
Diffstat (limited to 'gsl-1.9/cblas/test_symv.c')
-rw-r--r-- | gsl-1.9/cblas/test_symv.c | 379 |
1 files changed, 379 insertions, 0 deletions
diff --git a/gsl-1.9/cblas/test_symv.c b/gsl-1.9/cblas/test_symv.c new file mode 100644 index 0000000..923a0f5 --- /dev/null +++ b/gsl-1.9/cblas/test_symv.c @@ -0,0 +1,379 @@ +#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_symv (void) { +const double flteps = 1e-4, dbleps = 1e-6; + { + int order = 101; + int uplo = 121; + float alpha = 1.0f; + float beta = -1.0f; + int N = 1; + int lda = 1; + float A[] = { -0.428f }; + float X[] = { -0.34f }; + int incX = -1; + float Y[] = { -0.888f }; + int incY = -1; + float y_expected[] = { 1.03352f }; + cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1054)"); + } + }; + }; + + + { + int order = 101; + int uplo = 121; + float alpha = 1.0f; + float beta = -1.0f; + int N = 1; + int lda = 1; + float A[] = { -0.428f }; + float X[] = { -0.34f }; + int incX = -1; + float Y[] = { -0.888f }; + int incY = -1; + float y_expected[] = { 1.03352f }; + cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1055)"); + } + }; + }; + + + { + int order = 101; + int uplo = 122; + float alpha = 1.0f; + float beta = -1.0f; + int N = 1; + int lda = 1; + float A[] = { -0.428f }; + float X[] = { -0.34f }; + int incX = -1; + float Y[] = { -0.888f }; + int incY = -1; + float y_expected[] = { 1.03352f }; + cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1056)"); + } + }; + }; + + + { + int order = 101; + int uplo = 122; + float alpha = 1.0f; + float beta = -1.0f; + int N = 1; + int lda = 1; + float A[] = { -0.428f }; + float X[] = { -0.34f }; + int incX = -1; + float Y[] = { -0.888f }; + int incY = -1; + float y_expected[] = { 1.03352f }; + cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1057)"); + } + }; + }; + + + { + int order = 102; + int uplo = 121; + float alpha = 1.0f; + float beta = -1.0f; + int N = 1; + int lda = 1; + float A[] = { -0.428f }; + float X[] = { -0.34f }; + int incX = -1; + float Y[] = { -0.888f }; + int incY = -1; + float y_expected[] = { 1.03352f }; + cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1058)"); + } + }; + }; + + + { + int order = 102; + int uplo = 121; + float alpha = 1.0f; + float beta = -1.0f; + int N = 1; + int lda = 1; + float A[] = { -0.428f }; + float X[] = { -0.34f }; + int incX = -1; + float Y[] = { -0.888f }; + int incY = -1; + float y_expected[] = { 1.03352f }; + cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1059)"); + } + }; + }; + + + { + int order = 102; + int uplo = 122; + float alpha = 1.0f; + float beta = -1.0f; + int N = 1; + int lda = 1; + float A[] = { -0.428f }; + float X[] = { -0.34f }; + int incX = -1; + float Y[] = { -0.888f }; + int incY = -1; + float y_expected[] = { 1.03352f }; + cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1060)"); + } + }; + }; + + + { + int order = 102; + int uplo = 122; + float alpha = 1.0f; + float beta = -1.0f; + int N = 1; + int lda = 1; + float A[] = { -0.428f }; + float X[] = { -0.34f }; + int incX = -1; + float Y[] = { -0.888f }; + int incY = -1; + float y_expected[] = { 1.03352f }; + cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1061)"); + } + }; + }; + + + { + int order = 101; + int uplo = 121; + double alpha = 0; + double beta = -0.3; + int N = 1; + int lda = 1; + double A[] = { 0.544 }; + double X[] = { -0.601 }; + int incX = -1; + double Y[] = { -0.852 }; + int incY = -1; + double y_expected[] = { 0.2556 }; + cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1062)"); + } + }; + }; + + + { + int order = 101; + int uplo = 121; + double alpha = 0; + double beta = -0.3; + int N = 1; + int lda = 1; + double A[] = { 0.544 }; + double X[] = { -0.601 }; + int incX = -1; + double Y[] = { -0.852 }; + int incY = -1; + double y_expected[] = { 0.2556 }; + cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1063)"); + } + }; + }; + + + { + int order = 101; + int uplo = 122; + double alpha = 0; + double beta = -0.3; + int N = 1; + int lda = 1; + double A[] = { 0.544 }; + double X[] = { -0.601 }; + int incX = -1; + double Y[] = { -0.852 }; + int incY = -1; + double y_expected[] = { 0.2556 }; + cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1064)"); + } + }; + }; + + + { + int order = 101; + int uplo = 122; + double alpha = 0; + double beta = -0.3; + int N = 1; + int lda = 1; + double A[] = { 0.544 }; + double X[] = { -0.601 }; + int incX = -1; + double Y[] = { -0.852 }; + int incY = -1; + double y_expected[] = { 0.2556 }; + cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1065)"); + } + }; + }; + + + { + int order = 102; + int uplo = 121; + double alpha = 0; + double beta = -0.3; + int N = 1; + int lda = 1; + double A[] = { 0.544 }; + double X[] = { -0.601 }; + int incX = -1; + double Y[] = { -0.852 }; + int incY = -1; + double y_expected[] = { 0.2556 }; + cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1066)"); + } + }; + }; + + + { + int order = 102; + int uplo = 121; + double alpha = 0; + double beta = -0.3; + int N = 1; + int lda = 1; + double A[] = { 0.544 }; + double X[] = { -0.601 }; + int incX = -1; + double Y[] = { -0.852 }; + int incY = -1; + double y_expected[] = { 0.2556 }; + cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1067)"); + } + }; + }; + + + { + int order = 102; + int uplo = 122; + double alpha = 0; + double beta = -0.3; + int N = 1; + int lda = 1; + double A[] = { 0.544 }; + double X[] = { -0.601 }; + int incX = -1; + double Y[] = { -0.852 }; + int incY = -1; + double y_expected[] = { 0.2556 }; + cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1068)"); + } + }; + }; + + + { + int order = 102; + int uplo = 122; + double alpha = 0; + double beta = -0.3; + int N = 1; + int lda = 1; + double A[] = { 0.544 }; + double X[] = { -0.601 }; + int incX = -1; + double Y[] = { -0.852 }; + int incY = -1; + double y_expected[] = { 0.2556 }; + cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY); + { + int i; + for (i = 0; i < 1; i++) { + gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1069)"); + } + }; + }; + + +} |