1 #include <gsl/gsl_test.h>
2 #include <gsl/gsl_ieee_utils.h>
3 #include <gsl/gsl_math.h>
4 #include <gsl/gsl_cblas.h>
10 const double flteps = 1e-4, dbleps = 1e-6;
18 float A[] = { -0.428f };
19 float X[] = { -0.34f };
21 float Y[] = { -0.888f };
23 float y_expected[] = { 1.03352f };
24 cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
27 for (i = 0; i < 1; i++) {
28 gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1054)");
41 float A[] = { -0.428f };
42 float X[] = { -0.34f };
44 float Y[] = { -0.888f };
46 float y_expected[] = { 1.03352f };
47 cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
50 for (i = 0; i < 1; i++) {
51 gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1055)");
64 float A[] = { -0.428f };
65 float X[] = { -0.34f };
67 float Y[] = { -0.888f };
69 float y_expected[] = { 1.03352f };
70 cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
73 for (i = 0; i < 1; i++) {
74 gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1056)");
87 float A[] = { -0.428f };
88 float X[] = { -0.34f };
90 float Y[] = { -0.888f };
92 float y_expected[] = { 1.03352f };
93 cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
96 for (i = 0; i < 1; i++) {
97 gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1057)");
110 float A[] = { -0.428f };
111 float X[] = { -0.34f };
113 float Y[] = { -0.888f };
115 float y_expected[] = { 1.03352f };
116 cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
119 for (i = 0; i < 1; i++) {
120 gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1058)");
133 float A[] = { -0.428f };
134 float X[] = { -0.34f };
136 float Y[] = { -0.888f };
138 float y_expected[] = { 1.03352f };
139 cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
142 for (i = 0; i < 1; i++) {
143 gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1059)");
156 float A[] = { -0.428f };
157 float X[] = { -0.34f };
159 float Y[] = { -0.888f };
161 float y_expected[] = { 1.03352f };
162 cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
165 for (i = 0; i < 1; i++) {
166 gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1060)");
179 float A[] = { -0.428f };
180 float X[] = { -0.34f };
182 float Y[] = { -0.888f };
184 float y_expected[] = { 1.03352f };
185 cblas_ssymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
188 for (i = 0; i < 1; i++) {
189 gsl_test_rel(Y[i], y_expected[i], flteps, "ssymv(case 1061)");
202 double A[] = { 0.544 };
203 double X[] = { -0.601 };
205 double Y[] = { -0.852 };
207 double y_expected[] = { 0.2556 };
208 cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
211 for (i = 0; i < 1; i++) {
212 gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1062)");
225 double A[] = { 0.544 };
226 double X[] = { -0.601 };
228 double Y[] = { -0.852 };
230 double y_expected[] = { 0.2556 };
231 cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
234 for (i = 0; i < 1; i++) {
235 gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1063)");
248 double A[] = { 0.544 };
249 double X[] = { -0.601 };
251 double Y[] = { -0.852 };
253 double y_expected[] = { 0.2556 };
254 cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
257 for (i = 0; i < 1; i++) {
258 gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1064)");
271 double A[] = { 0.544 };
272 double X[] = { -0.601 };
274 double Y[] = { -0.852 };
276 double y_expected[] = { 0.2556 };
277 cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
280 for (i = 0; i < 1; i++) {
281 gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1065)");
294 double A[] = { 0.544 };
295 double X[] = { -0.601 };
297 double Y[] = { -0.852 };
299 double y_expected[] = { 0.2556 };
300 cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
303 for (i = 0; i < 1; i++) {
304 gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1066)");
317 double A[] = { 0.544 };
318 double X[] = { -0.601 };
320 double Y[] = { -0.852 };
322 double y_expected[] = { 0.2556 };
323 cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
326 for (i = 0; i < 1; i++) {
327 gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1067)");
340 double A[] = { 0.544 };
341 double X[] = { -0.601 };
343 double Y[] = { -0.852 };
345 double y_expected[] = { 0.2556 };
346 cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
349 for (i = 0; i < 1; i++) {
350 gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1068)");
363 double A[] = { 0.544 };
364 double X[] = { -0.601 };
366 double Y[] = { -0.852 };
368 double y_expected[] = { 0.2556 };
369 cblas_dsymv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
372 for (i = 0; i < 1; i++) {
373 gsl_test_rel(Y[i], y_expected[i], dbleps, "dsymv(case 1069)");