test_tpsv.c

该文件为c++的数学函数库!是一个非常有用的编程工具.它含有各种数学函数,为科学计算、工程应用等程序编写提供方便！

   int diag = 132;
   int N = 2;
   double A[] = { 0.355, -0.698, -0.42 };
   double X[] = { -0.702, -0.661 };
   int incX = -1;
   double x_expected[] = { -0.702, -1.150996 };
   cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1337)");
     }
   };
  };


  {
   int order = 102;
   int trans = 112;
   int uplo = 121;
   int diag = 131;
   int N = 2;
   double A[] = { 0.355, -0.698, -0.42 };
   double X[] = { -0.702, -0.661 };
   int incX = -1;
   double x_expected[] = { 4.76584842388, -1.86197183099 };
   cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1338)");
     }
   };
  };


  {
   int order = 102;
   int trans = 112;
   int uplo = 121;
   int diag = 132;
   int N = 2;
   double A[] = { 0.355, -0.698, -0.42 };
   double X[] = { -0.702, -0.661 };
   int incX = -1;
   double x_expected[] = { -1.163378, -0.661 };
   cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1339)");
     }
   };
  };


  {
   int order = 102;
   int trans = 112;
   int uplo = 122;
   int diag = 131;
   int N = 2;
   double A[] = { 0.355, -0.698, -0.42 };
   double X[] = { -0.702, -0.661 };
   int incX = -1;
   double x_expected[] = { 1.67142857143, 1.42438631791 };
   cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1340)");
     }
   };
  };


  {
   int order = 102;
   int trans = 112;
   int uplo = 122;
   int diag = 132;
   int N = 2;
   double A[] = { 0.355, -0.698, -0.42 };
   double X[] = { -0.702, -0.661 };
   int incX = -1;
   double x_expected[] = { -0.702, -1.150996 };
   cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1341)");
     }
   };
  };


  {
   int order = 101;
   int trans = 111;
   int uplo = 121;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1342) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1342) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 111;
   int uplo = 121;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1343) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1343) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 111;
   int uplo = 122;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1344) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1344) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 111;
   int uplo = 122;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1345) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1345) imag");
     };
   };
  };


  {
   int order = 102;
   int trans = 111;
   int uplo = 121;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1346) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1346) imag");
     };
   };
  };


  {
   int order = 102;
   int trans = 111;
   int uplo = 121;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1347) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1347) imag");
     };
   };
  };


  {
   int order = 102;
   int trans = 111;
   int uplo = 122;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1348) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1348) imag");
     };
   };
  };


  {
   int order = 102;
   int trans = 111;
   int uplo = 122;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1349) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1349) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 112;
   int uplo = 121;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1350) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1350) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 112;
   int uplo = 121;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1351) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1351) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 112;
   int uplo = 122;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1352) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1352) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 112;
   int uplo = 122;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1353) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1353) imag");
     };
   };
  };


  {
   int order = 102;
   int trans = 112;
   int uplo = 121;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1354) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1354) imag");
     };
   };
  };


  {
   int order = 102;
   int trans = 112;
   int uplo = 121;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1355) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1355) imag");
     };
   };
  };


  {
   int order = 102;
   int trans = 112;
   int uplo = 122;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1356) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1356) imag");
     };
   };
  };


  {
   int order = 102;
   int trans = 112;
   int uplo = 122;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1357) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1357) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 113;
   int uplo = 121;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { 2.15867f, 1.69498f, 1.69471f, 0.104738f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1358) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1358) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 113;
   int uplo = 121;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -0.613252f, 0.561896f, -0.072f, 0.642f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1359) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1359) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 113;
   int uplo = 122;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { 1.00465f, 0.327432f, 3.44853f, 2.273f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1360) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1360) imag");
     };
   };
  };


  {
   int order = 101;
   int trans = 113;
   int uplo = 122;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -0.133f, 0.998f, -0.806168f, -0.053086f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1361) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1361) imag");
     };
   };
  };


  {
   int order = 102;
   int trans = 113;
   int uplo = 121;
   int diag = 131;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { 2.15867f, 1.69498f, 1.69471f, 0.104738f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1362) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1362) imag");
     };
   };
  };


  {
   int order = 102;
   int trans = 113;
   int uplo = 121;
   int diag = 132;
   int N = 2;
   float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
   float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
   int incX = -1;
   float x_expected[] = { -0.613252f, 0.561896f, -0.072f, 0.642f };
   cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
   {
     int i;
     for (i = 0; i < 2; i++) {
       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1363) real");
       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1363) imag");
     };