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;
17 float A[] = { -0.587f, 0.14f, 0.841f };
18 float X[] = { -0.213f, 0.885f };
20 float x_expected[] = { -0.179133f, -0.549315f };
21 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
24 for (i = 0; i < 2; i++) {
25 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 974)");
37 float A[] = { -0.587f, 0.14f, 0.841f };
38 float X[] = { -0.213f, 0.885f };
40 float x_expected[] = { -0.213f, 0.85518f };
41 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
44 for (i = 0; i < 2; i++) {
45 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 975)");
57 float A[] = { -0.587f, 0.14f, 0.841f };
58 float X[] = { -0.213f, 0.885f };
60 float x_expected[] = { -0.055233f, -0.519495f };
61 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
64 for (i = 0; i < 2; i++) {
65 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 976)");
77 float A[] = { -0.587f, 0.14f, 0.841f };
78 float X[] = { -0.213f, 0.885f };
80 float x_expected[] = { -0.0891f, 0.885f };
81 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
84 for (i = 0; i < 2; i++) {
85 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 977)");
97 float A[] = { -0.587f, 0.14f, 0.841f };
98 float X[] = { -0.213f, 0.885f };
100 float x_expected[] = { -0.179133f, -0.549315f };
101 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
104 for (i = 0; i < 2; i++) {
105 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 978)");
117 float A[] = { -0.587f, 0.14f, 0.841f };
118 float X[] = { -0.213f, 0.885f };
120 float x_expected[] = { -0.213f, 0.85518f };
121 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
124 for (i = 0; i < 2; i++) {
125 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 979)");
137 float A[] = { -0.587f, 0.14f, 0.841f };
138 float X[] = { -0.213f, 0.885f };
140 float x_expected[] = { -0.055233f, -0.519495f };
141 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
144 for (i = 0; i < 2; i++) {
145 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 980)");
157 float A[] = { -0.587f, 0.14f, 0.841f };
158 float X[] = { -0.213f, 0.885f };
160 float x_expected[] = { -0.0891f, 0.885f };
161 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
164 for (i = 0; i < 2; i++) {
165 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 981)");
177 float A[] = { -0.765f, 0.968f, -0.956f };
178 float X[] = { 0.243f, -0.274f };
180 float x_expected[] = { -0.49754f, 0.20961f };
181 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
184 for (i = 0; i < 2; i++) {
185 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 982)");
197 float A[] = { -0.765f, 0.968f, -0.956f };
198 float X[] = { 0.243f, -0.274f };
200 float x_expected[] = { -0.022232f, -0.274f };
201 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
204 for (i = 0; i < 2; i++) {
205 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 983)");
217 float A[] = { -0.765f, 0.968f, -0.956f };
218 float X[] = { 0.243f, -0.274f };
220 float x_expected[] = { -0.232308f, 0.444834f };
221 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
224 for (i = 0; i < 2; i++) {
225 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 984)");
237 float A[] = { -0.765f, 0.968f, -0.956f };
238 float X[] = { 0.243f, -0.274f };
240 float x_expected[] = { 0.243f, -0.038776f };
241 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
244 for (i = 0; i < 2; i++) {
245 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 985)");
257 float A[] = { -0.765f, 0.968f, -0.956f };
258 float X[] = { 0.243f, -0.274f };
260 float x_expected[] = { -0.49754f, 0.20961f };
261 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
264 for (i = 0; i < 2; i++) {
265 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 986)");
277 float A[] = { -0.765f, 0.968f, -0.956f };
278 float X[] = { 0.243f, -0.274f };
280 float x_expected[] = { -0.022232f, -0.274f };
281 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
284 for (i = 0; i < 2; i++) {
285 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 987)");
297 float A[] = { -0.765f, 0.968f, -0.956f };
298 float X[] = { 0.243f, -0.274f };
300 float x_expected[] = { -0.232308f, 0.444834f };
301 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
304 for (i = 0; i < 2; i++) {
305 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 988)");
317 float A[] = { -0.765f, 0.968f, -0.956f };
318 float X[] = { 0.243f, -0.274f };
320 float x_expected[] = { 0.243f, -0.038776f };
321 cblas_stpmv(order, uplo, trans, diag, N, A, X, incX);
324 for (i = 0; i < 2; i++) {
325 gsl_test_rel(X[i], x_expected[i], flteps, "stpmv(case 989)");
337 double A[] = { 0.393, -0.221, 0.356 };
338 double X[] = { -0.062, -0.221 };
340 double x_expected[] = { -0.022072, -0.073151 };
341 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
344 for (i = 0; i < 2; i++) {
345 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 990)");
357 double A[] = { 0.393, -0.221, 0.356 };
358 double X[] = { -0.062, -0.221 };
360 double x_expected[] = { -0.062, -0.207298 };
361 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
364 for (i = 0; i < 2; i++) {
365 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 991)");
377 double A[] = { 0.393, -0.221, 0.356 };
378 double X[] = { -0.062, -0.221 };
380 double x_expected[] = { 0.026769, -0.086853 };
381 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
384 for (i = 0; i < 2; i++) {
385 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 992)");
397 double A[] = { 0.393, -0.221, 0.356 };
398 double X[] = { -0.062, -0.221 };
400 double x_expected[] = { -0.013159, -0.221 };
401 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
404 for (i = 0; i < 2; i++) {
405 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 993)");
417 double A[] = { 0.393, -0.221, 0.356 };
418 double X[] = { -0.062, -0.221 };
420 double x_expected[] = { -0.022072, -0.073151 };
421 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
424 for (i = 0; i < 2; i++) {
425 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 994)");
437 double A[] = { 0.393, -0.221, 0.356 };
438 double X[] = { -0.062, -0.221 };
440 double x_expected[] = { -0.062, -0.207298 };
441 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
444 for (i = 0; i < 2; i++) {
445 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 995)");
457 double A[] = { 0.393, -0.221, 0.356 };
458 double X[] = { -0.062, -0.221 };
460 double x_expected[] = { 0.026769, -0.086853 };
461 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
464 for (i = 0; i < 2; i++) {
465 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 996)");
477 double A[] = { 0.393, -0.221, 0.356 };
478 double X[] = { -0.062, -0.221 };
480 double x_expected[] = { -0.013159, -0.221 };
481 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
484 for (i = 0; i < 2; i++) {
485 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 997)");
497 double A[] = { 0.694, 0.501, 0.019 };
498 double X[] = { -0.928, 0.365 };
500 double x_expected[] = { 0.165233, 0.25331 };
501 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
504 for (i = 0; i < 2; i++) {
505 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 998)");
517 double A[] = { 0.694, 0.501, 0.019 };
518 double X[] = { -0.928, 0.365 };
520 double x_expected[] = { -0.745135, 0.365 };
521 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
524 for (i = 0; i < 2; i++) {
525 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 999)");
537 double A[] = { 0.694, 0.501, 0.019 };
538 double X[] = { -0.928, 0.365 };
540 double x_expected[] = { -0.017632, -0.211618 };
541 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
544 for (i = 0; i < 2; i++) {
545 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 1000)");
557 double A[] = { 0.694, 0.501, 0.019 };
558 double X[] = { -0.928, 0.365 };
560 double x_expected[] = { -0.928, -0.099928 };
561 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
564 for (i = 0; i < 2; i++) {
565 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 1001)");
577 double A[] = { 0.694, 0.501, 0.019 };
578 double X[] = { -0.928, 0.365 };
580 double x_expected[] = { 0.165233, 0.25331 };
581 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
584 for (i = 0; i < 2; i++) {
585 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 1002)");
597 double A[] = { 0.694, 0.501, 0.019 };
598 double X[] = { -0.928, 0.365 };
600 double x_expected[] = { -0.745135, 0.365 };
601 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
604 for (i = 0; i < 2; i++) {
605 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 1003)");
617 double A[] = { 0.694, 0.501, 0.019 };
618 double X[] = { -0.928, 0.365 };
620 double x_expected[] = { -0.017632, -0.211618 };
621 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
624 for (i = 0; i < 2; i++) {
625 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 1004)");
637 double A[] = { 0.694, 0.501, 0.019 };
638 double X[] = { -0.928, 0.365 };
640 double x_expected[] = { -0.928, -0.099928 };
641 cblas_dtpmv(order, uplo, trans, diag, N, A, X, incX);
644 for (i = 0; i < 2; i++) {
645 gsl_test_rel(X[i], x_expected[i], dbleps, "dtpmv(case 1005)");
657 float A[] = { 0.362f, -0.849f, -0.612f, -0.718f, 0.503f, -0.923f };
658 float X[] = { 0.904f, 0.461f, -0.367f, 0.153f };
660 float x_expected[] = { 0.880215f, -0.602509f, -0.225207f, -0.564235f };
661 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
664 for (i = 0; i < 2; i++) {
665 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1006) real");
666 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1006) imag");
678 float A[] = { 0.362f, -0.849f, -0.612f, -0.718f, 0.503f, -0.923f };
679 float X[] = { 0.904f, 0.461f, -0.367f, 0.153f };
681 float x_expected[] = { 0.904f, 0.461f, -0.58925f, -0.778204f };
682 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
685 for (i = 0; i < 2; i++) {
686 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1007) real");
687 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1007) imag");
699 float A[] = { 0.362f, -0.849f, -0.612f, -0.718f, 0.503f, -0.923f };
700 float X[] = { 0.904f, 0.461f, -0.367f, 0.153f };
702 float x_expected[] = { 1.21467f, -0.432639f, -0.002957f, 0.366969f };
703 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
706 for (i = 0; i < 2; i++) {
707 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1008) real");
708 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1008) imag");
720 float A[] = { 0.362f, -0.849f, -0.612f, -0.718f, 0.503f, -0.923f };
721 float X[] = { 0.904f, 0.461f, -0.367f, 0.153f };
723 float x_expected[] = { 1.23846f, 0.63087f, -0.367f, 0.153f };
724 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
727 for (i = 0; i < 2; i++) {
728 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1009) real");
729 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1009) imag");
741 float A[] = { 0.362f, -0.849f, -0.612f, -0.718f, 0.503f, -0.923f };
742 float X[] = { 0.904f, 0.461f, -0.367f, 0.153f };
744 float x_expected[] = { 0.880215f, -0.602509f, -0.225207f, -0.564235f };
745 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
748 for (i = 0; i < 2; i++) {
749 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1010) real");
750 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1010) imag");
762 float A[] = { 0.362f, -0.849f, -0.612f, -0.718f, 0.503f, -0.923f };
763 float X[] = { 0.904f, 0.461f, -0.367f, 0.153f };
765 float x_expected[] = { 0.904f, 0.461f, -0.58925f, -0.778204f };
766 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
769 for (i = 0; i < 2; i++) {
770 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1011) real");
771 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1011) imag");
783 float A[] = { 0.362f, -0.849f, -0.612f, -0.718f, 0.503f, -0.923f };
784 float X[] = { 0.904f, 0.461f, -0.367f, 0.153f };
786 float x_expected[] = { 1.21467f, -0.432639f, -0.002957f, 0.366969f };
787 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
790 for (i = 0; i < 2; i++) {
791 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1012) real");
792 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1012) imag");
804 float A[] = { 0.362f, -0.849f, -0.612f, -0.718f, 0.503f, -0.923f };
805 float X[] = { 0.904f, 0.461f, -0.367f, 0.153f };
807 float x_expected[] = { 1.23846f, 0.63087f, -0.367f, 0.153f };
808 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
811 for (i = 0; i < 2; i++) {
812 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1013) real");
813 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1013) imag");
825 float A[] = { -0.876f, -0.697f, -0.519f, -0.223f, 0.526f, -0.077f };
826 float X[] = { 0.338f, -0.807f, 0.444f, -0.748f };
828 float x_expected[] = { -0.281591f, -0.161308f, -0.9103f, 0.34578f };
829 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
832 for (i = 0; i < 2; i++) {
833 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1014) real");
834 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1014) imag");
846 float A[] = { -0.876f, -0.697f, -0.519f, -0.223f, 0.526f, -0.077f };
847 float X[] = { 0.338f, -0.807f, 0.444f, -0.748f };
849 float x_expected[] = { -0.05924f, -0.5178f, 0.444f, -0.748f };
850 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
853 for (i = 0; i < 2; i++) {
854 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1015) real");
855 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1015) imag");
867 float A[] = { -0.876f, -0.697f, -0.519f, -0.223f, 0.526f, -0.077f };
868 float X[] = { 0.338f, -0.807f, 0.444f, -0.748f };
870 float x_expected[] = { 0.115649f, -0.450508f, -1.26568f, 0.689239f };
871 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
874 for (i = 0; i < 2; i++) {
875 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1016) real");
876 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1016) imag");
888 float A[] = { -0.876f, -0.697f, -0.519f, -0.223f, 0.526f, -0.077f };
889 float X[] = { 0.338f, -0.807f, 0.444f, -0.748f };
891 float x_expected[] = { 0.338f, -0.807f, 0.088617f, -0.404541f };
892 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
895 for (i = 0; i < 2; i++) {
896 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1017) real");
897 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1017) imag");
909 float A[] = { -0.876f, -0.697f, -0.519f, -0.223f, 0.526f, -0.077f };
910 float X[] = { 0.338f, -0.807f, 0.444f, -0.748f };
912 float x_expected[] = { -0.281591f, -0.161308f, -0.9103f, 0.34578f };
913 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
916 for (i = 0; i < 2; i++) {
917 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1018) real");
918 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1018) imag");
930 float A[] = { -0.876f, -0.697f, -0.519f, -0.223f, 0.526f, -0.077f };
931 float X[] = { 0.338f, -0.807f, 0.444f, -0.748f };
933 float x_expected[] = { -0.05924f, -0.5178f, 0.444f, -0.748f };
934 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
937 for (i = 0; i < 2; i++) {
938 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1019) real");
939 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1019) imag");
951 float A[] = { -0.876f, -0.697f, -0.519f, -0.223f, 0.526f, -0.077f };
952 float X[] = { 0.338f, -0.807f, 0.444f, -0.748f };
954 float x_expected[] = { 0.115649f, -0.450508f, -1.26568f, 0.689239f };
955 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
958 for (i = 0; i < 2; i++) {
959 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1020) real");
960 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1020) imag");
972 float A[] = { -0.876f, -0.697f, -0.519f, -0.223f, 0.526f, -0.077f };
973 float X[] = { 0.338f, -0.807f, 0.444f, -0.748f };
975 float x_expected[] = { 0.338f, -0.807f, 0.088617f, -0.404541f };
976 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
979 for (i = 0; i < 2; i++) {
980 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1021) real");
981 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1021) imag");
993 float A[] = { 0.869f, -0.091f, -0.859f, 0.008f, -0.921f, -0.321f };
994 float X[] = { -0.122f, -0.364f, 0.602f, -0.96f };
996 float x_expected[] = { -0.295592f, 1.11591f, 0.610498f, -0.779458f };
997 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
1000 for (i = 0; i < 2; i++) {
1001 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1022) real");
1002 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1022) imag");
1014 float A[] = { 0.869f, -0.091f, -0.859f, 0.008f, -0.921f, -0.321f };
1015 float X[] = { -0.122f, -0.364f, 0.602f, -0.96f };
1017 float x_expected[] = { -0.646798f, 0.455824f, 0.602f, -0.96f };
1018 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
1021 for (i = 0; i < 2; i++) {
1022 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1023) real");
1023 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1023) imag");
1035 float A[] = { 0.869f, -0.091f, -0.859f, 0.008f, -0.921f, -0.321f };
1036 float X[] = { -0.122f, -0.364f, 0.602f, -0.96f };
1038 float x_expected[] = { 0.229206f, 0.296082f, 0.712384f, -0.465806f };
1039 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
1042 for (i = 0; i < 2; i++) {
1043 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1024) real");
1044 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1024) imag");
1056 float A[] = { 0.869f, -0.091f, -0.859f, 0.008f, -0.921f, -0.321f };
1057 float X[] = { -0.122f, -0.364f, 0.602f, -0.96f };
1059 float x_expected[] = { -0.122f, -0.364f, 0.703886f, -0.646348f };
1060 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
1063 for (i = 0; i < 2; i++) {
1064 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1025) real");
1065 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1025) imag");
1077 float A[] = { 0.869f, -0.091f, -0.859f, 0.008f, -0.921f, -0.321f };
1078 float X[] = { -0.122f, -0.364f, 0.602f, -0.96f };
1080 float x_expected[] = { -0.295592f, 1.11591f, 0.610498f, -0.779458f };
1081 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
1084 for (i = 0; i < 2; i++) {
1085 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1026) real");
1086 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1026) imag");
1098 float A[] = { 0.869f, -0.091f, -0.859f, 0.008f, -0.921f, -0.321f };
1099 float X[] = { -0.122f, -0.364f, 0.602f, -0.96f };
1101 float x_expected[] = { -0.646798f, 0.455824f, 0.602f, -0.96f };
1102 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
1105 for (i = 0; i < 2; i++) {
1106 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1027) real");
1107 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1027) imag");
1119 float A[] = { 0.869f, -0.091f, -0.859f, 0.008f, -0.921f, -0.321f };
1120 float X[] = { -0.122f, -0.364f, 0.602f, -0.96f };
1122 float x_expected[] = { 0.229206f, 0.296082f, 0.712384f, -0.465806f };
1123 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
1126 for (i = 0; i < 2; i++) {
1127 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1028) real");
1128 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1028) imag");
1140 float A[] = { 0.869f, -0.091f, -0.859f, 0.008f, -0.921f, -0.321f };
1141 float X[] = { -0.122f, -0.364f, 0.602f, -0.96f };
1143 float x_expected[] = { -0.122f, -0.364f, 0.703886f, -0.646348f };
1144 cblas_ctpmv(order, uplo, trans, diag, N, A, X, incX);
1147 for (i = 0; i < 2; i++) {
1148 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpmv(case 1029) real");
1149 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpmv(case 1029) imag");
1161 double A[] = { 0.254, 0.263, -0.271, -0.595, -0.182, -0.672 };
1162 double X[] = { -0.042, -0.705, -0.255, -0.854 };
1164 double x_expected[] = { -0.466116, 0.156534, -0.248261, -0.067936 };
1165 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1168 for (i = 0; i < 2; i++) {
1169 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1030) real");
1170 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1030) imag");
1182 double A[] = { 0.254, 0.263, -0.271, -0.595, -0.182, -0.672 };
1183 double X[] = { -0.042, -0.705, -0.255, -0.854 };
1185 double x_expected[] = { -0.042, -0.705, -0.663093, -0.637955 };
1186 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1189 for (i = 0; i < 2; i++) {
1190 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1031) real");
1191 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1031) imag");
1203 double A[] = { 0.254, 0.263, -0.271, -0.595, -0.182, -0.672 };
1204 double X[] = { -0.042, -0.705, -0.255, -0.854 };
1206 double x_expected[] = { -0.905141, 0.539693, 0.159832, -0.283981 };
1207 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1210 for (i = 0; i < 2; i++) {
1211 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1032) real");
1212 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1032) imag");
1224 double A[] = { 0.254, 0.263, -0.271, -0.595, -0.182, -0.672 };
1225 double X[] = { -0.042, -0.705, -0.255, -0.854 };
1227 double x_expected[] = { -0.481025, -0.321841, -0.255, -0.854 };
1228 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1231 for (i = 0; i < 2; i++) {
1232 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1033) real");
1233 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1033) imag");
1245 double A[] = { 0.254, 0.263, -0.271, -0.595, -0.182, -0.672 };
1246 double X[] = { -0.042, -0.705, -0.255, -0.854 };
1248 double x_expected[] = { -0.466116, 0.156534, -0.248261, -0.067936 };
1249 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1252 for (i = 0; i < 2; i++) {
1253 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1034) real");
1254 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1034) imag");
1266 double A[] = { 0.254, 0.263, -0.271, -0.595, -0.182, -0.672 };
1267 double X[] = { -0.042, -0.705, -0.255, -0.854 };
1269 double x_expected[] = { -0.042, -0.705, -0.663093, -0.637955 };
1270 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1273 for (i = 0; i < 2; i++) {
1274 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1035) real");
1275 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1035) imag");
1287 double A[] = { 0.254, 0.263, -0.271, -0.595, -0.182, -0.672 };
1288 double X[] = { -0.042, -0.705, -0.255, -0.854 };
1290 double x_expected[] = { -0.905141, 0.539693, 0.159832, -0.283981 };
1291 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1294 for (i = 0; i < 2; i++) {
1295 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1036) real");
1296 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1036) imag");
1308 double A[] = { 0.254, 0.263, -0.271, -0.595, -0.182, -0.672 };
1309 double X[] = { -0.042, -0.705, -0.255, -0.854 };
1311 double x_expected[] = { -0.481025, -0.321841, -0.255, -0.854 };
1312 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1315 for (i = 0; i < 2; i++) {
1316 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1037) real");
1317 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1037) imag");
1329 double A[] = { 0.421, -0.407, -0.595, -0.387, 0.884, -0.498 };
1330 double X[] = { -0.008, 0.904, -0.689, -0.679 };
1332 double x_expected[] = { 0.590302, 1.473768, -0.566422, -0.005436 };
1333 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1336 for (i = 0; i < 2; i++) {
1337 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1038) real");
1338 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1038) imag");
1350 double A[] = { 0.421, -0.407, -0.595, -0.387, 0.884, -0.498 };
1351 double X[] = { -0.008, 0.904, -0.689, -0.679 };
1353 double x_expected[] = { 0.139182, 1.574648, -0.689, -0.679 };
1354 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1357 for (i = 0; i < 2; i++) {
1358 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1039) real");
1359 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1039) imag");
1371 double A[] = { 0.421, -0.407, -0.595, -0.387, 0.884, -0.498 };
1372 double X[] = { -0.008, 0.904, -0.689, -0.679 };
1374 double x_expected[] = { 0.44312, 0.80312, -0.211814, -0.54022 };
1375 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1378 for (i = 0; i < 2; i++) {
1379 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1040) real");
1380 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1040) imag");
1392 double A[] = { 0.421, -0.407, -0.595, -0.387, 0.884, -0.498 };
1393 double X[] = { -0.008, 0.904, -0.689, -0.679 };
1395 double x_expected[] = { -0.008, 0.904, -0.334392, -1.213784 };
1396 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1399 for (i = 0; i < 2; i++) {
1400 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1041) real");
1401 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1041) imag");
1413 double A[] = { 0.421, -0.407, -0.595, -0.387, 0.884, -0.498 };
1414 double X[] = { -0.008, 0.904, -0.689, -0.679 };
1416 double x_expected[] = { 0.590302, 1.473768, -0.566422, -0.005436 };
1417 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1420 for (i = 0; i < 2; i++) {
1421 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1042) real");
1422 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1042) imag");
1434 double A[] = { 0.421, -0.407, -0.595, -0.387, 0.884, -0.498 };
1435 double X[] = { -0.008, 0.904, -0.689, -0.679 };
1437 double x_expected[] = { 0.139182, 1.574648, -0.689, -0.679 };
1438 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1441 for (i = 0; i < 2; i++) {
1442 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1043) real");
1443 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1043) imag");
1455 double A[] = { 0.421, -0.407, -0.595, -0.387, 0.884, -0.498 };
1456 double X[] = { -0.008, 0.904, -0.689, -0.679 };
1458 double x_expected[] = { 0.44312, 0.80312, -0.211814, -0.54022 };
1459 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1462 for (i = 0; i < 2; i++) {
1463 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1044) real");
1464 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1044) imag");
1476 double A[] = { 0.421, -0.407, -0.595, -0.387, 0.884, -0.498 };
1477 double X[] = { -0.008, 0.904, -0.689, -0.679 };
1479 double x_expected[] = { -0.008, 0.904, -0.334392, -1.213784 };
1480 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1483 for (i = 0; i < 2; i++) {
1484 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1045) real");
1485 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1045) imag");
1497 double A[] = { -0.743, -0.078, 0.77, 0.505, 0.157, -0.986 };
1498 double X[] = { -0.641, 0.565, -0.406, -0.948 };
1500 double x_expected[] = { -1.449087, -1.068251, 0.375602, 0.672696 };
1501 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1504 for (i = 0; i < 2; i++) {
1505 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1046) real");
1506 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1046) imag");
1518 double A[] = { -0.743, -0.078, 0.77, 0.505, 0.157, -0.986 };
1519 double X[] = { -0.641, 0.565, -0.406, -0.948 };
1521 double x_expected[] = { -1.43236, 0.04007, -0.406, -0.948 };
1522 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1525 for (i = 0; i < 2; i++) {
1526 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1047) real");
1527 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1047) imag");
1539 double A[] = { -0.743, -0.078, 0.77, 0.505, 0.157, -0.986 };
1540 double X[] = { -0.641, 0.565, -0.406, -0.948 };
1542 double x_expected[] = { -0.657727, -0.543321, 0.167357, 1.431451 };
1543 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1546 for (i = 0; i < 2; i++) {
1547 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1048) real");
1548 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1048) imag");
1560 double A[] = { -0.743, -0.078, 0.77, 0.505, 0.157, -0.986 };
1561 double X[] = { -0.641, 0.565, -0.406, -0.948 };
1563 double x_expected[] = { -0.641, 0.565, -0.614245, -0.189245 };
1564 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1567 for (i = 0; i < 2; i++) {
1568 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1049) real");
1569 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1049) imag");
1581 double A[] = { -0.743, -0.078, 0.77, 0.505, 0.157, -0.986 };
1582 double X[] = { -0.641, 0.565, -0.406, -0.948 };
1584 double x_expected[] = { -1.449087, -1.068251, 0.375602, 0.672696 };
1585 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1588 for (i = 0; i < 2; i++) {
1589 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1050) real");
1590 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1050) imag");
1602 double A[] = { -0.743, -0.078, 0.77, 0.505, 0.157, -0.986 };
1603 double X[] = { -0.641, 0.565, -0.406, -0.948 };
1605 double x_expected[] = { -1.43236, 0.04007, -0.406, -0.948 };
1606 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1609 for (i = 0; i < 2; i++) {
1610 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1051) real");
1611 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1051) imag");
1623 double A[] = { -0.743, -0.078, 0.77, 0.505, 0.157, -0.986 };
1624 double X[] = { -0.641, 0.565, -0.406, -0.948 };
1626 double x_expected[] = { -0.657727, -0.543321, 0.167357, 1.431451 };
1627 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1630 for (i = 0; i < 2; i++) {
1631 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1052) real");
1632 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1052) imag");
1644 double A[] = { -0.743, -0.078, 0.77, 0.505, 0.157, -0.986 };
1645 double X[] = { -0.641, 0.565, -0.406, -0.948 };
1647 double x_expected[] = { -0.641, 0.565, -0.614245, -0.189245 };
1648 cblas_ztpmv(order, uplo, trans, diag, N, A, X, incX);
1651 for (i = 0; i < 2; i++) {
1652 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpmv(case 1053) real");
1653 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpmv(case 1053) imag");