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;
19 float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
20 float X[] = { -0.305f, 0.61f, -0.831f };
22 float x_expected[] = { -0.354651f, -2.40855f, 0.481076f };
23 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
26 for (i = 0; i < 3; i++) {
27 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1230)");
41 float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
42 float X[] = { -0.305f, 0.61f, -0.831f };
44 float x_expected[] = { -0.305f, 0.84973f, -1.00859f };
45 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
48 for (i = 0; i < 3; i++) {
49 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1231)");
63 float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
64 float X[] = { -0.305f, 0.61f, -0.831f };
66 float x_expected[] = { -2.71619f, -1.09055f, -3.97608f };
67 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
70 for (i = 0; i < 3; i++) {
71 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1232)");
85 float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
86 float X[] = { -0.305f, 0.61f, -0.831f };
88 float x_expected[] = { -0.56589f, 0.303361f, -0.831f };
89 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
92 for (i = 0; i < 3; i++) {
93 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1233)");
107 float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
108 float X[] = { -0.305f, 0.61f, -0.831f };
110 float x_expected[] = { 1.30901f, -0.656172f, -5.13458f };
111 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
114 for (i = 0; i < 3; i++) {
115 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1234)");
129 float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
130 float X[] = { -0.305f, 0.61f, -0.831f };
132 float x_expected[] = { -0.305f, 0.8723f, -0.509121f };
133 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
136 for (i = 0; i < 3; i++) {
137 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1235)");
151 float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
152 float X[] = { -0.305f, 0.61f, -0.831f };
154 float x_expected[] = { 0.524539f, -0.961964f, 1.22026f };
155 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
158 for (i = 0; i < 3; i++) {
159 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1236)");
173 float A[] = { -0.681f, 0.209f, 0.436f, -0.369f, 0.786f, -0.84f, 0.86f, -0.233f, 0.734f };
174 float X[] = { -0.305f, 0.61f, -0.831f };
176 float x_expected[] = { -0.920972f, 0.783679f, -0.831f };
177 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
180 for (i = 0; i < 3; i++) {
181 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1237)");
195 float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
196 float X[] = { 0.144f, 0.635f, 0.116f };
198 float x_expected[] = { 16.8676f, 17.3503f, 5.27273f };
199 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
202 for (i = 0; i < 3; i++) {
203 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1238)");
217 float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
218 float X[] = { 0.144f, 0.635f, 0.116f };
220 float x_expected[] = { 0.209676f, 0.54278f, 0.116f };
221 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
224 for (i = 0; i < 3; i++) {
225 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1239)");
239 float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
240 float X[] = { 0.144f, 0.635f, 0.116f };
242 float x_expected[] = { 0.212077f, -5.01482f, -1.14722f };
243 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
246 for (i = 0; i < 3; i++) {
247 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1240)");
261 float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
262 float X[] = { 0.144f, 0.635f, 0.116f };
264 float x_expected[] = { 0.144f, 0.615848f, 0.242249f };
265 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
268 for (i = 0; i < 3; i++) {
269 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1241)");
283 float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
284 float X[] = { 0.144f, 0.635f, 0.116f };
286 float x_expected[] = { 1.28844f, -5.49514f, 0.145912f };
287 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
290 for (i = 0; i < 3; i++) {
291 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1242)");
305 float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
306 float X[] = { 0.144f, 0.635f, 0.116f };
308 float x_expected[] = { 0.0563823f, 0.65878f, 0.116f };
309 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
312 for (i = 0; i < 3; i++) {
313 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1243)");
327 float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
328 float X[] = { 0.144f, 0.635f, 0.116f };
330 float x_expected[] = { 1.08271f, -3.73662f, 140.301f };
331 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
334 for (i = 0; i < 3; i++) {
335 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1244)");
349 float A[] = { 0.022f, 0.795f, -0.389f, -0.205f, -0.121f, 0.323f, 0.133f, 0.679f, 0.742f };
350 float X[] = { 0.144f, 0.635f, 0.116f };
352 float x_expected[] = { 0.144f, 0.652424f, -0.402677f };
353 cblas_stbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
356 for (i = 0; i < 3; i++) {
357 gsl_test_rel(X[i], x_expected[i], flteps, "stbsv(case 1245)");
371 double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
372 double X[] = { 0.332, 0.588, 0.252 };
374 double x_expected[] = { -0.967930029155, 0.138412575592, 0.506166027443 };
375 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
378 for (i = 0; i < 3; i++) {
379 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1246)");
393 double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
394 double X[] = { 0.332, 0.588, 0.252 };
396 double x_expected[] = { 0.332, 0.819736, 0.615143048 };
397 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
400 for (i = 0; i < 3; i++) {
401 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1247)");
415 double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
416 double X[] = { 0.332, 0.588, 0.252 };
418 double x_expected[] = { -0.364842154056, -0.326531140246, -0.568848758465 };
419 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
422 for (i = 0; i < 3; i++) {
423 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1248)");
437 double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
438 double X[] = { 0.332, 0.588, 0.252 };
440 double x_expected[] = { 0.588397988, 0.747516, 0.252 };
441 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
444 for (i = 0; i < 3; i++) {
445 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1249)");
459 double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
460 double X[] = { 0.332, 0.588, 0.252 };
462 double x_expected[] = { -0.550580431177, -0.571849444278, 0.248263427151 };
463 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
466 for (i = 0; i < 3; i++) {
467 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1250)");
481 double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
482 double X[] = { 0.332, 0.588, 0.252 };
484 double x_expected[] = { 0.332, 0.701876, 0.696287508 };
485 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
488 for (i = 0; i < 3; i++) {
489 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1251)");
503 double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
504 double X[] = { 0.332, 0.588, 0.252 };
506 double x_expected[] = { 1.50217883761, -1.21382140588, 0.407108239095 };
507 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
510 for (i = 0; i < 3; i++) {
511 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1252)");
525 double A[] = { 0.619, -0.443, 0.957, -0.633, -0.698, 0.783, -0.343, -0.603, 0.735 };
526 double X[] = { 0.332, 0.588, 0.252 };
528 double x_expected[] = { 0.820345928, 0.699636, 0.252 };
529 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
532 for (i = 0; i < 3; i++) {
533 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1253)");
547 double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
548 double X[] = { 0.58, 0.558, -0.54 };
550 double x_expected[] = { 18.994209959, 20.323927329, 2.7135678392 };
551 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
554 for (i = 0; i < 3; i++) {
555 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1254)");
569 double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
570 double X[] = { 0.58, 0.558, -0.54 };
572 double x_expected[] = { 1.06925836, 0.72162, -0.54 };
573 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
576 for (i = 0; i < 3; i++) {
577 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1255)");
591 double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
592 double X[] = { 0.58, 0.558, -0.54 };
594 double x_expected[] = { -3.27683615819, -4.47682615869, -1.97425326753 };
595 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
598 for (i = 0; i < 3; i++) {
599 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1256)");
613 double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
614 double X[] = { 0.58, 0.558, -0.54 };
616 double x_expected[] = { 0.58, 0.11952, -0.53844624 };
617 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
620 for (i = 0; i < 3; i++) {
621 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1257)");
635 double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
636 double X[] = { 0.58, 0.558, -0.54 };
638 double x_expected[] = { -6.6461072986, -0.788837290809, -1.78217821782 };
639 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
642 for (i = 0; i < 3; i++) {
643 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1258)");
657 double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
658 double X[] = { 0.58, 0.558, -0.54 };
660 double x_expected[] = { 0.16345912, 0.55098, -0.54 };
661 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
664 for (i = 0; i < 3; i++) {
665 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1259)");
679 double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
680 double X[] = { 0.58, 0.558, -0.54 };
682 double x_expected[] = { 0.767195767196, -82.9352869353, -123.564783625 };
683 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
686 for (i = 0; i < 3; i++) {
687 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1260)");
701 double A[] = { -0.199, 0.303, -0.705, -0.013, -0.678, 0.547, 0.756, -0.177, -0.079 };
702 double X[] = { 0.58, 0.558, -0.54 };
704 double x_expected[] = { 0.58, 0.95124, -0.82822572 };
705 cblas_dtbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
708 for (i = 0; i < 3; i++) {
709 gsl_test_rel(X[i], x_expected[i], dbleps, "dtbsv(case 1261)");
723 float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
724 float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
726 float x_expected[] = { 1.28871f, 0.289887f, 1.76043f, 1.27481f, 1.56506f, -2.35181f };
727 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
730 for (i = 0; i < 3; i++) {
731 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1262) real");
732 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1262) imag");
746 float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
747 float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
749 float x_expected[] = { 0.11f, 0.787f, -1.04259f, 0.18935f, 0.228474f, -0.564917f };
750 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
753 for (i = 0; i < 3; i++) {
754 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1263) real");
755 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1263) imag");
769 float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
770 float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
772 float x_expected[] = { -0.0906249f, 3.09442f, -1.60036f, 1.28475f, -0.582941f, 0.0383898f };
773 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
776 for (i = 0; i < 3; i++) {
777 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1264) real");
778 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1264) imag");
792 float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
793 float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
795 float x_expected[] = { 1.05233f, 0.79657f, -0.566883f, 1.46031f, -0.437f, 0.592f };
796 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
799 for (i = 0; i < 3; i++) {
800 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1265) real");
801 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1265) imag");
815 float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
816 float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
818 float x_expected[] = { -0.735844f, 1.11782f, -0.28244f, 1.16117f, -0.66707f, 0.938302f };
819 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
822 for (i = 0; i < 3; i++) {
823 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1266) real");
824 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1266) imag");
838 float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
839 float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
841 float x_expected[] = { 0.11f, 0.787f, -0.406239f, 0.580226f, -0.171935f, 1.2125f };
842 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
845 for (i = 0; i < 3; i++) {
846 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1267) real");
847 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1267) imag");
861 float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
862 float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
864 float x_expected[] = { 1.70081f, 2.20477f, 1.32753f, -0.522112f, 0.0223652f, -0.62248f };
865 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
868 for (i = 0; i < 3; i++) {
869 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1268) real");
870 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1268) imag");
884 float A[] = { -0.975f, -0.667f, 0.813f, -0.962f, -0.961f, 0.226f, -0.503f, 0.809f, 0.81f, -0.162f, -0.027f, -0.044f, 0.212f, 0.563f, 0.446f, -0.392f, 0.798f, -0.07f };
885 float X[] = { 0.11f, 0.787f, -0.826f, 0.809f, -0.437f, 0.592f };
887 float x_expected[] = { 0.967596f, 0.693563f, -1.04022f, -0.09269f, -0.437f, 0.592f };
888 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
891 for (i = 0; i < 3; i++) {
892 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1269) real");
893 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1269) imag");
907 float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
908 float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
910 float x_expected[] = { -1.11985f, 0.801655f, 0.273814f, -1.09438f, -0.52531f, 0.166748f };
911 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
914 for (i = 0; i < 3; i++) {
915 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1270) real");
916 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1270) imag");
930 float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
931 float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
933 float x_expected[] = { 0.266087f, 0.618557f, 0.031897f, -0.914419f, -0.134f, 0.179f };
934 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
937 for (i = 0; i < 3; i++) {
938 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1271) real");
939 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1271) imag");
953 float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
954 float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
956 float x_expected[] = { -0.762749f, -0.016292f, 1.59299f, 0.158751f, -4.75603f, -1.78591f };
957 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
960 for (i = 0; i < 3; i++) {
961 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1272) real");
962 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1272) imag");
976 float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
977 float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
979 float x_expected[] = { -0.509f, 0.608f, -0.332731f, -1.24444f, 0.262904f, 1.21961f };
980 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
983 for (i = 0; i < 3; i++) {
984 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1273) real");
985 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1273) imag");
999 float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
1000 float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
1002 float x_expected[] = { -1.76046f, 0.0455463f, 1.38348f, 0.700097f, -0.669451f, 0.321896f };
1003 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1006 for (i = 0; i < 3; i++) {
1007 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1274) real");
1008 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1274) imag");
1022 float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
1023 float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
1025 float x_expected[] = { 0.151523f, 0.78611f, 0.120309f, -1.01387f, -0.134f, 0.179f };
1026 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1029 for (i = 0; i < 3; i++) {
1030 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1275) real");
1031 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1275) imag");
1045 float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
1046 float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
1048 float x_expected[] = { -1.00779f, -0.620278f, 0.81164f, -1.90759f, -1.32022f, 1.48356f };
1049 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1052 for (i = 0; i < 3; i++) {
1053 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1276) real");
1054 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1276) imag");
1068 float A[] = { 0.33f, -0.236f, 0.267f, -0.139f, 0.25f, 0.509f, 0.86f, -0.089f, -0.018f, -0.847f, 0.424f, -0.573f, 0.097f, -0.663f, 0.65f, -0.811f, 0.283f, 0.032f };
1069 float X[] = { -0.509f, 0.608f, 0.021f, -0.848f, -0.134f, 0.179f };
1071 float x_expected[] = { -0.509f, 0.608f, -0.503138f, -1.26818f, 0.176615f, 0.447668f };
1072 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1075 for (i = 0; i < 3; i++) {
1076 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1277) real");
1077 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1277) imag");
1091 float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
1092 float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
1094 float x_expected[] = { -0.613838f, -1.13321f, -1.34847f, 0.0432903f, 0.0879552f, -0.479334f };
1095 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1098 for (i = 0; i < 3; i++) {
1099 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1278) real");
1100 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1278) imag");
1114 float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
1115 float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
1117 float x_expected[] = { 0.76323f, -1.23595f, 0.943058f, -0.618694f, 0.296f, 0.034f };
1118 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1121 for (i = 0; i < 3; i++) {
1122 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1279) real");
1123 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1279) imag");
1137 float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
1138 float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
1140 float x_expected[] = { -1.15557f, -2.50103f, -3.85402f, -1.04833f, 0.414582f, 5.91218f };
1141 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1144 for (i = 0; i < 3; i++) {
1145 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1280) real");
1146 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1280) imag");
1160 float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
1161 float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
1163 float x_expected[] = { -0.037f, -0.599f, 1.39953f, -0.064424f, 1.0801f, -0.481747f };
1164 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1167 for (i = 0; i < 3; i++) {
1168 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1281) real");
1169 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1281) imag");
1183 float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
1184 float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
1186 float x_expected[] = { -3.0802f, -9.09377f, -1.05845f, 0.99239f, 0.259763f, -0.687744f };
1187 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1190 for (i = 0; i < 3; i++) {
1191 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1282) real");
1192 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1282) imag");
1206 float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
1207 float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
1209 float x_expected[] = { -0.513897f, 0.632031f, 1.14112f, -0.580648f, 0.296f, 0.034f };
1210 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1213 for (i = 0; i < 3; i++) {
1214 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1283) real");
1215 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1283) imag");
1229 float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
1230 float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
1232 float x_expected[] = { 0.360899f, -0.456643f, -2.31803f, 0.257877f, 1.56928f, -0.922115f };
1233 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1236 for (i = 0; i < 3; i++) {
1237 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1284) real");
1238 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1284) imag");
1252 float A[] = { 0.041f, -0.61f, 0.099f, -0.393f, 0.357f, -0.984f, -0.576f, -0.342f, -0.903f, -0.083f, -0.157f, -0.694f, 0.768f, 0.688f, 0.203f, -0.079f, 0.298f, -0.424f };
1253 float X[] = { -0.037f, -0.599f, 0.959f, -0.499f, 0.296f, 0.034f };
1255 float x_expected[] = { -0.037f, -0.599f, 0.875872f, -1.03683f, -0.198184f, -0.207572f };
1256 cblas_ctbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1259 for (i = 0; i < 3; i++) {
1260 gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbsv(case 1285) real");
1261 gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbsv(case 1285) imag");
1275 double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
1276 double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
1278 double x_expected[] = { 0.0490338308139, -0.158433417494, 0.261604043488, 1.28058846321, 1.77633350191, -1.07039599422 };
1279 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1282 for (i = 0; i < 3; i++) {
1283 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1286) real");
1284 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1286) imag");
1298 double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
1299 double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
1301 double x_expected[] = { -0.123, 0.122, 0.96534, 0.346049, 1.067212328, 0.445330131 };
1302 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1305 for (i = 0; i < 3; i++) {
1306 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1287) real");
1307 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1287) imag");
1321 double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
1322 double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
1324 double x_expected[] = { 72.7437666278, 10.4206532927, -4.34946941374, -14.8012581742, 2.01859491883, -1.53922125931 };
1325 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1328 for (i = 0; i < 3; i++) {
1329 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1288) real");
1330 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1288) imag");
1344 double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
1345 double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
1347 double x_expected[] = { -0.464775024, 0.662224708, -0.0457, 0.610264, 0.942, 0.98 };
1348 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1351 for (i = 0; i < 3; i++) {
1352 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1289) real");
1353 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1289) imag");
1367 double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
1368 double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
1370 double x_expected[] = { -0.591747295323, -0.534096923761, -4.60251824353, 1.70172936273, -4.94687072873, -3.32536493524 };
1371 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1374 for (i = 0; i < 3; i++) {
1375 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1290) real");
1376 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1290) imag");
1390 double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
1391 double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
1393 double x_expected[] = { -0.123, 0.122, 0.807692, 0.373091, 0.384974988, 1.400879194 };
1394 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1397 for (i = 0; i < 3; i++) {
1398 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1291) real");
1399 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1291) imag");
1413 double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
1414 double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
1416 double x_expected[] = { -0.129998778267, -0.116630230861, 0.993340886904, 0.530739563688, 1.55891621291, -0.284019181928 };
1417 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1420 for (i = 0; i < 3; i++) {
1421 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1292) real");
1422 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1292) imag");
1436 double A[] = { 0.474, 0.715, 0.061, 0.532, 0.004, -0.318, 0.37, -0.692, -0.166, 0.039, -0.946, 0.857, -0.922, -0.491, 0.012, -0.217, -0.674, -0.429 };
1437 double X[] = { -0.123, 0.122, 0.981, 0.321, 0.942, 0.98 };
1439 double x_expected[] = { 0.107496032, 0.025821594, 1.444898, -0.239924, 0.942, 0.98 };
1440 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1443 for (i = 0; i < 3; i++) {
1444 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1293) real");
1445 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1293) imag");
1459 double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
1460 double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
1462 double x_expected[] = { -0.825842176606, 0.212941473892, -0.548817434511, -0.703261551538, 0.0746069436827, 0.425751789407 };
1463 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1466 for (i = 0; i < 3; i++) {
1467 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1294) real");
1468 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1294) imag");
1482 double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
1483 double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
1485 double x_expected[] = { -0.619710352, 0.018225936, 1.211252, 0.891864, 0.293, -0.434 };
1486 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1489 for (i = 0; i < 3; i++) {
1490 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1295) real");
1491 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1295) imag");
1505 double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
1506 double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
1508 double x_expected[] = { 0.203289119964, 1.58288482537, -1.7720160159, 0.479463518178, -0.511241930019, -1.79333888299 };
1509 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1512 for (i = 0; i < 3; i++) {
1513 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1296) real");
1514 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1296) imag");
1528 double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
1529 double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
1531 double x_expected[] = { -0.373, 0.566, 0.618602, -0.084689, 0.887531803, -0.570220771 };
1532 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1535 for (i = 0; i < 3; i++) {
1536 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1297) real");
1537 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1297) imag");
1551 double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
1552 double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
1554 double x_expected[] = { 1.72799012007, 13.4612400765, 4.46126528205, -0.0212528722047, 0.627282377919, 0.302760084926 };
1555 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1558 for (i = 0; i < 3; i++) {
1559 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1298) real");
1560 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1298) imag");
1574 double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
1575 double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
1577 double x_expected[] = { -1.280839615, 1.560525655, 1.167331, 0.179227, 0.293, -0.434 };
1578 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1581 for (i = 0; i < 3; i++) {
1582 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1299) real");
1583 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1299) imag");
1597 double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
1598 double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
1600 double x_expected[] = { -0.594503951847, 0.00287302167266, -1.08185265666, -0.859860374254, 0.0331027077244, 1.28233265933 };
1601 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1604 for (i = 0; i < 3; i++) {
1605 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1300) real");
1606 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1300) imag");
1620 double A[] = { -0.872, -0.841, 0.108, -0.744, 0.231, -0.513, -0.973, 0.087, 0.348, 0.196, 0.447, 0.307, 0.632, -0.949, 0.322, 0.277, 0.282, 0.831 };
1621 double X[] = { -0.373, 0.566, 0.92, 0.627, 0.293, -0.434 };
1623 double x_expected[] = { -0.373, 0.566, 1.16074, 0.50314, -0.20669608, 0.37525144 };
1624 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1627 for (i = 0; i < 3; i++) {
1628 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1301) real");
1629 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1301) imag");
1643 double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
1644 double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
1646 double x_expected[] = { 0.0654496252357, 0.224007771015, -0.752486084395, -0.554870892947, -0.587163401057, 0.166737652215 };
1647 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1650 for (i = 0; i < 3; i++) {
1651 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1302) real");
1652 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1302) imag");
1666 double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
1667 double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
1669 double x_expected[] = { -0.595558802, -1.147174647, 0.589506, -0.500919, -0.126, 0.459 };
1670 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1673 for (i = 0; i < 3; i++) {
1674 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1303) real");
1675 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1303) imag");
1689 double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
1690 double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
1692 double x_expected[] = { 3.39346077201, 0.652889512141, -2.33602680355, -2.7859245153, -5.04672104102, -0.334110541026 };
1693 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1696 for (i = 0; i < 3; i++) {
1697 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1304) real");
1698 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1304) imag");
1712 double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
1713 double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
1715 double x_expected[] = { 0.028, -0.804, -0.109456, -0.217192, -0.41110804, 0.41693792 };
1716 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1719 for (i = 0; i < 3; i++) {
1720 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1305) real");
1721 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1305) imag");
1735 double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
1736 double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
1738 double x_expected[] = { 7.16970224467, -0.772071373678, 0.833386981173, -0.673826630129, -0.26524050899, 0.465327628365 };
1739 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1742 for (i = 0; i < 3; i++) {
1743 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1306) real");
1744 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1306) imag");
1758 double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
1759 double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
1761 double x_expected[] = { 0.471459157, -1.566755859, 0.940839, 0.357132, -0.126, 0.459 };
1762 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1765 for (i = 0; i < 3; i++) {
1766 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1307) real");
1767 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1307) imag");
1781 double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
1782 double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
1784 double x_expected[] = { -0.909961830373, 0.118063054039, -0.0169425582229, -1.00055409731, -1.37205489923, 0.994032418785 };
1785 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1788 for (i = 0; i < 3; i++) {
1789 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1308) real");
1790 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1308) imag");
1804 double A[] = { 0.404, 0.667, 0.861, 0.22, 0.298, -0.858, -0.682, -0.969, 0.327, -0.86, 0.125, 0.606, -0.143, -0.865, -0.036, 0.23, -0.776, 0.079 };
1805 double X[] = { 0.028, -0.804, 0.582, -0.078, -0.126, 0.459 };
1807 double x_expected[] = { 0.028, -0.804, -0.118596, 0.160828, -0.059271004, 0.294435972 };
1808 cblas_ztbsv(order, uplo, trans, diag, N, K, A, lda, X, incX);
1811 for (i = 0; i < 3; i++) {
1812 gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbsv(case 1309) real");
1813 gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbsv(case 1309) imag");