| #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_scal (void) { |
| const double flteps = 1e-4, dbleps = 1e-6; |
| { |
| int N = 1; |
| float alpha = 0.0f; |
| float X[] = { 0.651f }; |
| int incX = -1; |
| float expected[] = { 0.651f }; |
| cblas_sscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[i], expected[i], flteps, "sscal(case 112)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| float alpha = 0.1f; |
| float X[] = { 0.651f }; |
| int incX = -1; |
| float expected[] = { 0.651f }; |
| cblas_sscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[i], expected[i], flteps, "sscal(case 113)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| float alpha = 1.0f; |
| float X[] = { 0.651f }; |
| int incX = -1; |
| float expected[] = { 0.651f }; |
| cblas_sscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[i], expected[i], flteps, "sscal(case 114)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| double alpha = 0; |
| double X[] = { 0.686 }; |
| int incX = -1; |
| double expected[] = { 0.686 }; |
| cblas_dscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 115)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| double alpha = 0.1; |
| double X[] = { 0.686 }; |
| int incX = -1; |
| double expected[] = { 0.686 }; |
| cblas_dscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 116)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| double alpha = 1; |
| double X[] = { 0.686 }; |
| int incX = -1; |
| double expected[] = { 0.686 }; |
| cblas_dscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 117)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| float alpha[2] = {0.0f, 0.0f}; |
| float X[] = { 0.986f, -0.775f }; |
| int incX = -1; |
| float expected[] = { 0.986f, -0.775f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 118) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 118) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| float alpha[2] = {0.1f, 0.0f}; |
| float X[] = { 0.986f, -0.775f }; |
| int incX = -1; |
| float expected[] = { 0.986f, -0.775f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 119) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 119) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| float alpha[2] = {1.0f, 0.0f}; |
| float X[] = { 0.986f, -0.775f }; |
| int incX = -1; |
| float expected[] = { 0.986f, -0.775f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 120) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 120) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| float alpha[2] = {0.0f, 0.1f}; |
| float X[] = { 0.986f, -0.775f }; |
| int incX = -1; |
| float expected[] = { 0.986f, -0.775f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 121) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 121) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| float alpha[2] = {0.1f, 0.2f}; |
| float X[] = { 0.986f, -0.775f }; |
| int incX = -1; |
| float expected[] = { 0.986f, -0.775f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 122) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 122) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| float alpha[2] = {1.0f, 0.3f}; |
| float X[] = { 0.986f, -0.775f }; |
| int incX = -1; |
| float expected[] = { 0.986f, -0.775f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 123) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 123) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| double alpha[2] = {0, 0}; |
| double X[] = { 0.454, -0.478 }; |
| int incX = -1; |
| double expected[] = { 0.454, -0.478 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 124) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 124) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| double alpha[2] = {0.1, 0}; |
| double X[] = { 0.454, -0.478 }; |
| int incX = -1; |
| double expected[] = { 0.454, -0.478 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 125) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 125) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| double alpha[2] = {1, 0}; |
| double X[] = { 0.454, -0.478 }; |
| int incX = -1; |
| double expected[] = { 0.454, -0.478 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 126) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 126) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| double alpha[2] = {0, 0.1}; |
| double X[] = { 0.454, -0.478 }; |
| int incX = -1; |
| double expected[] = { 0.454, -0.478 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 127) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 127) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| double alpha[2] = {0.1, 0.2}; |
| double X[] = { 0.454, -0.478 }; |
| int incX = -1; |
| double expected[] = { 0.454, -0.478 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 128) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 128) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 1; |
| double alpha[2] = {1, 0.3}; |
| double X[] = { 0.454, -0.478 }; |
| int incX = -1; |
| double expected[] = { 0.454, -0.478 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 1; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 129) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 129) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha = 0.0f; |
| float X[] = { 0.389f, -0.236f }; |
| int incX = 1; |
| float expected[] = { 0.0f, -0.0f }; |
| cblas_sscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], flteps, "sscal(case 130)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha = 0.1f; |
| float X[] = { 0.389f, -0.236f }; |
| int incX = 1; |
| float expected[] = { 0.0389f, -0.0236f }; |
| cblas_sscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], flteps, "sscal(case 131)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha = 1.0f; |
| float X[] = { 0.389f, -0.236f }; |
| int incX = 1; |
| float expected[] = { 0.389f, -0.236f }; |
| cblas_sscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], flteps, "sscal(case 132)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha = 0; |
| double X[] = { -0.429, -0.183 }; |
| int incX = 1; |
| double expected[] = { -0.0, -0.0 }; |
| cblas_dscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 133)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha = 0.1; |
| double X[] = { -0.429, -0.183 }; |
| int incX = 1; |
| double expected[] = { -0.0429, -0.0183 }; |
| cblas_dscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 134)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha = 1; |
| double X[] = { -0.429, -0.183 }; |
| int incX = 1; |
| double expected[] = { -0.429, -0.183 }; |
| cblas_dscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 135)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {0.0f, 0.0f}; |
| float X[] = { -0.603f, 0.239f, 0.339f, -0.58f }; |
| int incX = 1; |
| float expected[] = { -0.0f, 0.0f, 0.0f, 0.0f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 136) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 136) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {0.1f, 0.0f}; |
| float X[] = { -0.603f, 0.239f, 0.339f, -0.58f }; |
| int incX = 1; |
| float expected[] = { -0.0603f, 0.0239f, 0.0339f, -0.058f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 137) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 137) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {1.0f, 0.0f}; |
| float X[] = { -0.603f, 0.239f, 0.339f, -0.58f }; |
| int incX = 1; |
| float expected[] = { -0.603f, 0.239f, 0.339f, -0.58f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 138) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 138) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {0.0f, 0.1f}; |
| float X[] = { -0.603f, 0.239f, 0.339f, -0.58f }; |
| int incX = 1; |
| float expected[] = { -0.0239f, -0.0603f, 0.058f, 0.0339f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 139) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 139) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {0.1f, 0.2f}; |
| float X[] = { -0.603f, 0.239f, 0.339f, -0.58f }; |
| int incX = 1; |
| float expected[] = { -0.1081f, -0.0967f, 0.1499f, 0.0098f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 140) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 140) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {1.0f, 0.3f}; |
| float X[] = { -0.603f, 0.239f, 0.339f, -0.58f }; |
| int incX = 1; |
| float expected[] = { -0.6747f, 0.0581f, 0.513f, -0.4783f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 141) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 141) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {0, 0}; |
| double X[] = { -0.956, 0.613, 0.443, 0.503 }; |
| int incX = 1; |
| double expected[] = { -0.0, 0.0, 0.0, 0.0 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 142) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 142) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {0.1, 0}; |
| double X[] = { -0.956, 0.613, 0.443, 0.503 }; |
| int incX = 1; |
| double expected[] = { -0.0956, 0.0613, 0.0443, 0.0503 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 143) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 143) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {1, 0}; |
| double X[] = { -0.956, 0.613, 0.443, 0.503 }; |
| int incX = 1; |
| double expected[] = { -0.956, 0.613, 0.443, 0.503 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 144) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 144) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {0, 0.1}; |
| double X[] = { -0.956, 0.613, 0.443, 0.503 }; |
| int incX = 1; |
| double expected[] = { -0.0613, -0.0956, -0.0503, 0.0443 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 145) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 145) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {0.1, 0.2}; |
| double X[] = { -0.956, 0.613, 0.443, 0.503 }; |
| int incX = 1; |
| double expected[] = { -0.2182, -0.1299, -0.0563, 0.1389 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 146) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 146) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {1, 0.3}; |
| double X[] = { -0.956, 0.613, 0.443, 0.503 }; |
| int incX = 1; |
| double expected[] = { -1.1399, 0.3262, 0.2921, 0.6359 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 147) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 147) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha = 0.0f; |
| float X[] = { 0.629f, -0.419f }; |
| int incX = -1; |
| float expected[] = { 0.629f, -0.419f }; |
| cblas_sscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], flteps, "sscal(case 148)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha = 0.1f; |
| float X[] = { 0.629f, -0.419f }; |
| int incX = -1; |
| float expected[] = { 0.629f, -0.419f }; |
| cblas_sscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], flteps, "sscal(case 149)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha = 1.0f; |
| float X[] = { 0.629f, -0.419f }; |
| int incX = -1; |
| float expected[] = { 0.629f, -0.419f }; |
| cblas_sscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], flteps, "sscal(case 150)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha = 0; |
| double X[] = { 0.398, -0.656 }; |
| int incX = -1; |
| double expected[] = { 0.398, -0.656 }; |
| cblas_dscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 151)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha = 0.1; |
| double X[] = { 0.398, -0.656 }; |
| int incX = -1; |
| double expected[] = { 0.398, -0.656 }; |
| cblas_dscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 152)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha = 1; |
| double X[] = { 0.398, -0.656 }; |
| int incX = -1; |
| double expected[] = { 0.398, -0.656 }; |
| cblas_dscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 153)"); |
| } |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {0.0f, 0.0f}; |
| float X[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| int incX = -1; |
| float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 154) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 154) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {0.1f, 0.0f}; |
| float X[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| int incX = -1; |
| float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 155) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 155) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {1.0f, 0.0f}; |
| float X[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| int incX = -1; |
| float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 156) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 156) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {0.0f, 0.1f}; |
| float X[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| int incX = -1; |
| float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 157) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 157) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {0.1f, 0.2f}; |
| float X[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| int incX = -1; |
| float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 158) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 158) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| float alpha[2] = {1.0f, 0.3f}; |
| float X[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| int incX = -1; |
| float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f }; |
| cblas_cscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 159) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 159) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {0, 0}; |
| double X[] = { 0.521, -0.811, 0.556, -0.147 }; |
| int incX = -1; |
| double expected[] = { 0.521, -0.811, 0.556, -0.147 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 160) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 160) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {0.1, 0}; |
| double X[] = { 0.521, -0.811, 0.556, -0.147 }; |
| int incX = -1; |
| double expected[] = { 0.521, -0.811, 0.556, -0.147 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 161) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 161) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {1, 0}; |
| double X[] = { 0.521, -0.811, 0.556, -0.147 }; |
| int incX = -1; |
| double expected[] = { 0.521, -0.811, 0.556, -0.147 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 162) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 162) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {0, 0.1}; |
| double X[] = { 0.521, -0.811, 0.556, -0.147 }; |
| int incX = -1; |
| double expected[] = { 0.521, -0.811, 0.556, -0.147 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 163) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 163) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {0.1, 0.2}; |
| double X[] = { 0.521, -0.811, 0.556, -0.147 }; |
| int incX = -1; |
| double expected[] = { 0.521, -0.811, 0.556, -0.147 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 164) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 164) imag"); |
| }; |
| }; |
| }; |
| |
| |
| { |
| int N = 2; |
| double alpha[2] = {1, 0.3}; |
| double X[] = { 0.521, -0.811, 0.556, -0.147 }; |
| int incX = -1; |
| double expected[] = { 0.521, -0.811, 0.556, -0.147 }; |
| cblas_zscal(N, alpha, X, incX); |
| { |
| int i; |
| for (i = 0; i < 2; i++) { |
| gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 165) real"); |
| gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 165) imag"); |
| }; |
| }; |
| }; |
| |
| |
| } |