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.987f };
19 float X[] = { -0.138f };
21 float x_expected[] = { -0.136206f };
22 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
25 for (i = 0; i < 1; i++) {
26 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 814)");
39 float A[] = { 0.987f };
40 float X[] = { -0.138f };
42 float x_expected[] = { -0.138f };
43 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
46 for (i = 0; i < 1; i++) {
47 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 815)");
60 float A[] = { 0.987f };
61 float X[] = { -0.138f };
63 float x_expected[] = { -0.136206f };
64 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
67 for (i = 0; i < 1; i++) {
68 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 816)");
81 float A[] = { 0.987f };
82 float X[] = { -0.138f };
84 float x_expected[] = { -0.138f };
85 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
88 for (i = 0; i < 1; i++) {
89 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 817)");
102 float A[] = { 0.987f };
103 float X[] = { -0.138f };
105 float x_expected[] = { -0.136206f };
106 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
109 for (i = 0; i < 1; i++) {
110 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 818)");
123 float A[] = { 0.987f };
124 float X[] = { -0.138f };
126 float x_expected[] = { -0.138f };
127 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
130 for (i = 0; i < 1; i++) {
131 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 819)");
144 float A[] = { 0.987f };
145 float X[] = { -0.138f };
147 float x_expected[] = { -0.136206f };
148 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
151 for (i = 0; i < 1; i++) {
152 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 820)");
165 float A[] = { 0.987f };
166 float X[] = { -0.138f };
168 float x_expected[] = { -0.138f };
169 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
172 for (i = 0; i < 1; i++) {
173 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 821)");
186 float A[] = { -0.329f };
187 float X[] = { 0.463f };
189 float x_expected[] = { -0.152327f };
190 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
193 for (i = 0; i < 1; i++) {
194 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 822)");
207 float A[] = { -0.329f };
208 float X[] = { 0.463f };
210 float x_expected[] = { 0.463f };
211 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
214 for (i = 0; i < 1; i++) {
215 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 823)");
228 float A[] = { -0.329f };
229 float X[] = { 0.463f };
231 float x_expected[] = { -0.152327f };
232 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
235 for (i = 0; i < 1; i++) {
236 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 824)");
249 float A[] = { -0.329f };
250 float X[] = { 0.463f };
252 float x_expected[] = { 0.463f };
253 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
256 for (i = 0; i < 1; i++) {
257 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 825)");
270 float A[] = { -0.329f };
271 float X[] = { 0.463f };
273 float x_expected[] = { -0.152327f };
274 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
277 for (i = 0; i < 1; i++) {
278 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 826)");
291 float A[] = { -0.329f };
292 float X[] = { 0.463f };
294 float x_expected[] = { 0.463f };
295 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
298 for (i = 0; i < 1; i++) {
299 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 827)");
312 float A[] = { -0.329f };
313 float X[] = { 0.463f };
315 float x_expected[] = { -0.152327f };
316 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
319 for (i = 0; i < 1; i++) {
320 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 828)");
333 float A[] = { -0.329f };
334 float X[] = { 0.463f };
336 float x_expected[] = { 0.463f };
337 cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);
340 for (i = 0; i < 1; i++) {
341 gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 829)");
354 double A[] = { -0.429 };
355 double X[] = { -0.899 };
357 double x_expected[] = { 0.385671 };
358 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
361 for (i = 0; i < 1; i++) {
362 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 830)");
375 double A[] = { -0.429 };
376 double X[] = { -0.899 };
378 double x_expected[] = { -0.899 };
379 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
382 for (i = 0; i < 1; i++) {
383 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 831)");
396 double A[] = { -0.429 };
397 double X[] = { -0.899 };
399 double x_expected[] = { 0.385671 };
400 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
403 for (i = 0; i < 1; i++) {
404 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 832)");
417 double A[] = { -0.429 };
418 double X[] = { -0.899 };
420 double x_expected[] = { -0.899 };
421 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
424 for (i = 0; i < 1; i++) {
425 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 833)");
438 double A[] = { -0.429 };
439 double X[] = { -0.899 };
441 double x_expected[] = { 0.385671 };
442 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
445 for (i = 0; i < 1; i++) {
446 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 834)");
459 double A[] = { -0.429 };
460 double X[] = { -0.899 };
462 double x_expected[] = { -0.899 };
463 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
466 for (i = 0; i < 1; i++) {
467 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 835)");
480 double A[] = { -0.429 };
481 double X[] = { -0.899 };
483 double x_expected[] = { 0.385671 };
484 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
487 for (i = 0; i < 1; i++) {
488 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 836)");
501 double A[] = { -0.429 };
502 double X[] = { -0.899 };
504 double x_expected[] = { -0.899 };
505 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
508 for (i = 0; i < 1; i++) {
509 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 837)");
522 double A[] = { 0.842 };
523 double X[] = { 0.192 };
525 double x_expected[] = { 0.161664 };
526 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
529 for (i = 0; i < 1; i++) {
530 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 838)");
543 double A[] = { 0.842 };
544 double X[] = { 0.192 };
546 double x_expected[] = { 0.192 };
547 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
550 for (i = 0; i < 1; i++) {
551 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 839)");
564 double A[] = { 0.842 };
565 double X[] = { 0.192 };
567 double x_expected[] = { 0.161664 };
568 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
571 for (i = 0; i < 1; i++) {
572 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 840)");
585 double A[] = { 0.842 };
586 double X[] = { 0.192 };
588 double x_expected[] = { 0.192 };
589 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
592 for (i = 0; i < 1; i++) {
593 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 841)");
606 double A[] = { 0.842 };
607 double X[] = { 0.192 };
609 double x_expected[] = { 0.161664 };
610 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
613 for (i = 0; i < 1; i++) {
614 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 842)");
627 double A[] = { 0.842 };
628 double X[] = { 0.192 };
630 double x_expected[] = { 0.192 };
631 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
634 for (i = 0; i < 1; i++) {
635 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 843)");
648 double A[] = { 0.842 };
649 double X[] = { 0.192 };
651 double x_expected[] = { 0.161664 };
652 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
655 for (i = 0; i < 1; i++) {
656 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 844)");
669 double A[] = { 0.842 };
670 double X[] = { 0.192 };
672 double x_expected[] = { 0.192 };
673 cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);
676 for (i = 0; i < 1; i++) {
677 gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 845)");
690 float A[] = { -0.162f, -0.108f };
691 float X[] = { 0.542f, 0.461f };
693 float x_expected[] = { -0.038016f, -0.133218f };
694 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
697 for (i = 0; i < 1; i++) {
698 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 846) real");
699 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 846) imag");
712 float A[] = { -0.162f, -0.108f };
713 float X[] = { 0.542f, 0.461f };
715 float x_expected[] = { 0.542f, 0.461f };
716 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
719 for (i = 0; i < 1; i++) {
720 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 847) real");
721 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 847) imag");
734 float A[] = { -0.162f, -0.108f };
735 float X[] = { 0.542f, 0.461f };
737 float x_expected[] = { -0.038016f, -0.133218f };
738 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
741 for (i = 0; i < 1; i++) {
742 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 848) real");
743 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 848) imag");
756 float A[] = { -0.162f, -0.108f };
757 float X[] = { 0.542f, 0.461f };
759 float x_expected[] = { 0.542f, 0.461f };
760 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
763 for (i = 0; i < 1; i++) {
764 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 849) real");
765 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 849) imag");
778 float A[] = { -0.162f, -0.108f };
779 float X[] = { 0.542f, 0.461f };
781 float x_expected[] = { -0.038016f, -0.133218f };
782 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
785 for (i = 0; i < 1; i++) {
786 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 850) real");
787 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 850) imag");
800 float A[] = { -0.162f, -0.108f };
801 float X[] = { 0.542f, 0.461f };
803 float x_expected[] = { 0.542f, 0.461f };
804 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
807 for (i = 0; i < 1; i++) {
808 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 851) real");
809 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 851) imag");
822 float A[] = { -0.162f, -0.108f };
823 float X[] = { 0.542f, 0.461f };
825 float x_expected[] = { -0.038016f, -0.133218f };
826 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
829 for (i = 0; i < 1; i++) {
830 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 852) real");
831 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 852) imag");
844 float A[] = { -0.162f, -0.108f };
845 float X[] = { 0.542f, 0.461f };
847 float x_expected[] = { 0.542f, 0.461f };
848 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
851 for (i = 0; i < 1; i++) {
852 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 853) real");
853 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 853) imag");
866 float A[] = { 0.547f, 0.583f };
867 float X[] = { -0.302f, 0.434f };
869 float x_expected[] = { -0.418216f, 0.061332f };
870 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
873 for (i = 0; i < 1; i++) {
874 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 854) real");
875 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 854) imag");
888 float A[] = { 0.547f, 0.583f };
889 float X[] = { -0.302f, 0.434f };
891 float x_expected[] = { -0.302f, 0.434f };
892 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
895 for (i = 0; i < 1; i++) {
896 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 855) real");
897 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 855) imag");
910 float A[] = { 0.547f, 0.583f };
911 float X[] = { -0.302f, 0.434f };
913 float x_expected[] = { -0.418216f, 0.061332f };
914 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
917 for (i = 0; i < 1; i++) {
918 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 856) real");
919 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 856) imag");
932 float A[] = { 0.547f, 0.583f };
933 float X[] = { -0.302f, 0.434f };
935 float x_expected[] = { -0.302f, 0.434f };
936 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
939 for (i = 0; i < 1; i++) {
940 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 857) real");
941 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 857) imag");
954 float A[] = { 0.547f, 0.583f };
955 float X[] = { -0.302f, 0.434f };
957 float x_expected[] = { -0.418216f, 0.061332f };
958 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
961 for (i = 0; i < 1; i++) {
962 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 858) real");
963 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 858) imag");
976 float A[] = { 0.547f, 0.583f };
977 float X[] = { -0.302f, 0.434f };
979 float x_expected[] = { -0.302f, 0.434f };
980 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
983 for (i = 0; i < 1; i++) {
984 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 859) real");
985 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 859) imag");
998 float A[] = { 0.547f, 0.583f };
999 float X[] = { -0.302f, 0.434f };
1001 float x_expected[] = { -0.418216f, 0.061332f };
1002 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1005 for (i = 0; i < 1; i++) {
1006 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 860) real");
1007 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 860) imag");
1020 float A[] = { 0.547f, 0.583f };
1021 float X[] = { -0.302f, 0.434f };
1023 float x_expected[] = { -0.302f, 0.434f };
1024 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1027 for (i = 0; i < 1; i++) {
1028 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 861) real");
1029 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 861) imag");
1042 float A[] = { 0.216f, 0.192f };
1043 float X[] = { -0.564f, -0.297f };
1045 float x_expected[] = { -0.178848f, 0.044136f };
1046 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1049 for (i = 0; i < 1; i++) {
1050 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 862) real");
1051 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 862) imag");
1064 float A[] = { 0.216f, 0.192f };
1065 float X[] = { -0.564f, -0.297f };
1067 float x_expected[] = { -0.564f, -0.297f };
1068 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1071 for (i = 0; i < 1; i++) {
1072 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 863) real");
1073 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 863) imag");
1086 float A[] = { 0.216f, 0.192f };
1087 float X[] = { -0.564f, -0.297f };
1089 float x_expected[] = { -0.178848f, 0.044136f };
1090 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1093 for (i = 0; i < 1; i++) {
1094 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 864) real");
1095 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 864) imag");
1108 float A[] = { 0.216f, 0.192f };
1109 float X[] = { -0.564f, -0.297f };
1111 float x_expected[] = { -0.564f, -0.297f };
1112 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1115 for (i = 0; i < 1; i++) {
1116 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 865) real");
1117 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 865) imag");
1130 float A[] = { 0.216f, 0.192f };
1131 float X[] = { -0.564f, -0.297f };
1133 float x_expected[] = { -0.178848f, 0.044136f };
1134 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1137 for (i = 0; i < 1; i++) {
1138 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 866) real");
1139 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 866) imag");
1152 float A[] = { 0.216f, 0.192f };
1153 float X[] = { -0.564f, -0.297f };
1155 float x_expected[] = { -0.564f, -0.297f };
1156 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1159 for (i = 0; i < 1; i++) {
1160 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 867) real");
1161 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 867) imag");
1174 float A[] = { 0.216f, 0.192f };
1175 float X[] = { -0.564f, -0.297f };
1177 float x_expected[] = { -0.178848f, 0.044136f };
1178 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1181 for (i = 0; i < 1; i++) {
1182 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 868) real");
1183 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 868) imag");
1196 float A[] = { 0.216f, 0.192f };
1197 float X[] = { -0.564f, -0.297f };
1199 float x_expected[] = { -0.564f, -0.297f };
1200 cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1203 for (i = 0; i < 1; i++) {
1204 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 869) real");
1205 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 869) imag");
1218 double A[] = { 0.693, -0.22 };
1219 double X[] = { -0.101, 0.889 };
1221 double x_expected[] = { 0.125587, 0.638297 };
1222 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1225 for (i = 0; i < 1; i++) {
1226 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 870) real");
1227 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 870) imag");
1240 double A[] = { 0.693, -0.22 };
1241 double X[] = { -0.101, 0.889 };
1243 double x_expected[] = { -0.101, 0.889 };
1244 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1247 for (i = 0; i < 1; i++) {
1248 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 871) real");
1249 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 871) imag");
1262 double A[] = { 0.693, -0.22 };
1263 double X[] = { -0.101, 0.889 };
1265 double x_expected[] = { 0.125587, 0.638297 };
1266 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1269 for (i = 0; i < 1; i++) {
1270 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 872) real");
1271 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 872) imag");
1284 double A[] = { 0.693, -0.22 };
1285 double X[] = { -0.101, 0.889 };
1287 double x_expected[] = { -0.101, 0.889 };
1288 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1291 for (i = 0; i < 1; i++) {
1292 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 873) real");
1293 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 873) imag");
1306 double A[] = { 0.693, -0.22 };
1307 double X[] = { -0.101, 0.889 };
1309 double x_expected[] = { 0.125587, 0.638297 };
1310 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1313 for (i = 0; i < 1; i++) {
1314 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 874) real");
1315 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 874) imag");
1328 double A[] = { 0.693, -0.22 };
1329 double X[] = { -0.101, 0.889 };
1331 double x_expected[] = { -0.101, 0.889 };
1332 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1335 for (i = 0; i < 1; i++) {
1336 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 875) real");
1337 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 875) imag");
1350 double A[] = { 0.693, -0.22 };
1351 double X[] = { -0.101, 0.889 };
1353 double x_expected[] = { 0.125587, 0.638297 };
1354 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1357 for (i = 0; i < 1; i++) {
1358 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 876) real");
1359 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 876) imag");
1372 double A[] = { 0.693, -0.22 };
1373 double X[] = { -0.101, 0.889 };
1375 double x_expected[] = { -0.101, 0.889 };
1376 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1379 for (i = 0; i < 1; i++) {
1380 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 877) real");
1381 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 877) imag");
1394 double A[] = { -0.216, -0.623 };
1395 double X[] = { 0.048, 0.293 };
1397 double x_expected[] = { 0.172171, -0.093192 };
1398 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1401 for (i = 0; i < 1; i++) {
1402 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 878) real");
1403 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 878) imag");
1416 double A[] = { -0.216, -0.623 };
1417 double X[] = { 0.048, 0.293 };
1419 double x_expected[] = { 0.048, 0.293 };
1420 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1423 for (i = 0; i < 1; i++) {
1424 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 879) real");
1425 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 879) imag");
1438 double A[] = { -0.216, -0.623 };
1439 double X[] = { 0.048, 0.293 };
1441 double x_expected[] = { 0.172171, -0.093192 };
1442 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1445 for (i = 0; i < 1; i++) {
1446 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 880) real");
1447 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 880) imag");
1460 double A[] = { -0.216, -0.623 };
1461 double X[] = { 0.048, 0.293 };
1463 double x_expected[] = { 0.048, 0.293 };
1464 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1467 for (i = 0; i < 1; i++) {
1468 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 881) real");
1469 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 881) imag");
1482 double A[] = { -0.216, -0.623 };
1483 double X[] = { 0.048, 0.293 };
1485 double x_expected[] = { 0.172171, -0.093192 };
1486 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1489 for (i = 0; i < 1; i++) {
1490 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 882) real");
1491 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 882) imag");
1504 double A[] = { -0.216, -0.623 };
1505 double X[] = { 0.048, 0.293 };
1507 double x_expected[] = { 0.048, 0.293 };
1508 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1511 for (i = 0; i < 1; i++) {
1512 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 883) real");
1513 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 883) imag");
1526 double A[] = { -0.216, -0.623 };
1527 double X[] = { 0.048, 0.293 };
1529 double x_expected[] = { 0.172171, -0.093192 };
1530 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1533 for (i = 0; i < 1; i++) {
1534 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 884) real");
1535 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 884) imag");
1548 double A[] = { -0.216, -0.623 };
1549 double X[] = { 0.048, 0.293 };
1551 double x_expected[] = { 0.048, 0.293 };
1552 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1555 for (i = 0; i < 1; i++) {
1556 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 885) real");
1557 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 885) imag");
1570 double A[] = { -0.345, -0.851 };
1571 double X[] = { -0.708, 0.298 };
1573 double x_expected[] = { -0.009338, -0.705318 };
1574 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1577 for (i = 0; i < 1; i++) {
1578 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 886) real");
1579 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 886) imag");
1592 double A[] = { -0.345, -0.851 };
1593 double X[] = { -0.708, 0.298 };
1595 double x_expected[] = { -0.708, 0.298 };
1596 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1599 for (i = 0; i < 1; i++) {
1600 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 887) real");
1601 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 887) imag");
1614 double A[] = { -0.345, -0.851 };
1615 double X[] = { -0.708, 0.298 };
1617 double x_expected[] = { -0.009338, -0.705318 };
1618 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1621 for (i = 0; i < 1; i++) {
1622 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 888) real");
1623 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 888) imag");
1636 double A[] = { -0.345, -0.851 };
1637 double X[] = { -0.708, 0.298 };
1639 double x_expected[] = { -0.708, 0.298 };
1640 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1643 for (i = 0; i < 1; i++) {
1644 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 889) real");
1645 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 889) imag");
1658 double A[] = { -0.345, -0.851 };
1659 double X[] = { -0.708, 0.298 };
1661 double x_expected[] = { -0.009338, -0.705318 };
1662 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1665 for (i = 0; i < 1; i++) {
1666 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 890) real");
1667 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 890) imag");
1680 double A[] = { -0.345, -0.851 };
1681 double X[] = { -0.708, 0.298 };
1683 double x_expected[] = { -0.708, 0.298 };
1684 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1687 for (i = 0; i < 1; i++) {
1688 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 891) real");
1689 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 891) imag");
1702 double A[] = { -0.345, -0.851 };
1703 double X[] = { -0.708, 0.298 };
1705 double x_expected[] = { -0.009338, -0.705318 };
1706 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1709 for (i = 0; i < 1; i++) {
1710 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 892) real");
1711 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 892) imag");
1724 double A[] = { -0.345, -0.851 };
1725 double X[] = { -0.708, 0.298 };
1727 double x_expected[] = { -0.708, 0.298 };
1728 cblas_ztrmv(order, uplo, trans, diag, N, A, lda, X, incX);
1731 for (i = 0; i < 1; i++) {
1732 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztrmv(case 893) real");
1733 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztrmv(case 893) imag");