test_syr2k.c

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

   int N = 1;
   int K = 2;
   float alpha[2] = {0.0f, 0.1f};
   float beta[2] = {0.0f, 0.0f};
   float A[] = { -0.248f, -0.037f, -0.124f, 0.998f };
   int lda = 2;
   float B[] = { -0.608f, -0.115f, -0.718f, -0.551f };
   int ldb = 2;
   float C[] = { 0.187f, -0.329f };
   int ldc = 1;
   float C_expected[] = { 0.119445f, 0.157092f };
   cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1630) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1630) imag");
     };
   };
  };


  {
   int order = 101;
   int uplo = 122;
   int trans = 111;
   int N = 1;
   int K = 2;
   float alpha[2] = {0.0f, 0.1f};
   float beta[2] = {0.0f, 0.0f};
   float A[] = { 0.068f, 0.751f, -0.449f, -0.598f };
   int lda = 2;
   float B[] = { 0.616f, 0.805f, -0.635f, 0.773f };
   int ldb = 2;
   float C[] = { -0.287f, 0.917f };
   int ldc = 1;
   float C_expected[] = { -0.110002f, 0.0369404f };
   cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1631) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1631) imag");
     };
   };
  };


  {
   int order = 102;
   int uplo = 121;
   int trans = 111;
   int N = 1;
   int K = 2;
   float alpha[2] = {0.0f, 0.1f};
   float beta[2] = {0.0f, 0.0f};
   float A[] = { -0.396f, -0.603f, -0.131f, -0.288f };
   int lda = 1;
   float B[] = { -0.64f, -0.444f, -0.085f, 0.936f };
   int ldb = 1;
   float C[] = { 0.375f, -0.434f };
   int ldc = 1;
   float C_expected[] = { -0.0927216f, 0.0532822f };
   cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1632) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1632) imag");
     };
   };
  };


  {
   int order = 102;
   int uplo = 122;
   int trans = 111;
   int N = 1;
   int K = 2;
   float alpha[2] = {0.0f, 0.1f};
   float beta[2] = {0.0f, 0.0f};
   float A[] = { 0.655f, 0.16f, 0.45f, -0.747f };
   int lda = 1;
   float B[] = { 0.923f, 0.432f, -0.986f, 0.259f };
   int ldb = 1;
   float C[] = { 0.752f, 0.576f };
   int ldc = 1;
   float C_expected[] = { -0.256746f, 0.0570436f };
   cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1633) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1633) imag");
     };
   };
  };


  {
   int order = 101;
   int uplo = 121;
   int trans = 112;
   int N = 1;
   int K = 2;
   float alpha[2] = {0.0f, 0.1f};
   float beta[2] = {0.0f, 0.0f};
   float A[] = { -0.765f, 0.487f, 0.7f, 0.768f };
   int lda = 1;
   float B[] = { -0.529f, 0.056f, -0.584f, 0.928f };
   int ldb = 1;
   float C[] = { -0.426f, 0.836f };
   int ldc = 1;
   float C_expected[] = { 0.019875f, -0.148818f };
   cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1634) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1634) imag");
     };
   };
  };


  {
   int order = 101;
   int uplo = 122;
   int trans = 112;
   int N = 1;
   int K = 2;
   float alpha[2] = {0.0f, 0.1f};
   float beta[2] = {0.0f, 0.0f};
   float A[] = { 0.25f, 0.489f, 0.8f, -0.642f };
   int lda = 1;
   float B[] = { -0.732f, -0.856f, -0.654f, 0.591f };
   int ldb = 1;
   float C[] = { -0.101f, 0.322f };
   int ldc = 1;
   float C_expected[] = { -0.064144f, 0.0183612f };
   cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1635) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1635) imag");
     };
   };
  };


  {
   int order = 102;
   int uplo = 121;
   int trans = 112;
   int N = 1;
   int K = 2;
   float alpha[2] = {0.0f, 0.1f};
   float beta[2] = {0.0f, 0.0f};
   float A[] = { -0.579f, -0.971f, 0.521f, -0.824f };
   int lda = 2;
   float B[] = { -0.227f, 0.907f, 0.457f, -0.274f };
   int ldb = 2;
   float C[] = { 0.21f, -0.718f };
   int ldc = 1;
   float C_expected[] = { 0.164812f, 0.20489f };
   cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1636) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1636) imag");
     };
   };
  };


  {
   int order = 102;
   int uplo = 122;
   int trans = 112;
   int N = 1;
   int K = 2;
   float alpha[2] = {0.0f, 0.1f};
   float beta[2] = {0.0f, 0.0f};
   float A[] = { -0.83f, -0.512f, -0.667f, -0.436f };
   int lda = 2;
   float B[] = { -0.443f, 0.82f, -0.259f, -0.618f };
   int ldb = 2;
   float C[] = { 0.583f, 0.668f };
   int ldc = 1;
   float C_expected[] = { -0.0142692f, 0.138167f };
   cblas_csyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], flteps, "csyr2k(case 1637) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], flteps, "csyr2k(case 1637) imag");
     };
   };
  };


  {
   int order = 101;
   int uplo = 121;
   int trans = 111;
   int N = 1;
   int K = 2;
   double alpha[2] = {0, 0};
   double beta[2] = {-0.3, 0.1};
   double A[] = { -0.315, 0.03, 0.281, 0.175 };
   int lda = 2;
   double B[] = { -0.832, -0.964, 0.291, 0.476 };
   int ldb = 2;
   double C[] = { -0.341, 0.743 };
   int ldc = 1;
   double C_expected[] = { 0.028, -0.257 };
   cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1638) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1638) imag");
     };
   };
  };


  {
   int order = 101;
   int uplo = 122;
   int trans = 111;
   int N = 1;
   int K = 2;
   double alpha[2] = {0, 0};
   double beta[2] = {-0.3, 0.1};
   double A[] = { -0.159, -0.489, -0.11, 0.611 };
   int lda = 2;
   double B[] = { -0.285, -0.048, -0.673, -0.492 };
   int ldb = 2;
   double C[] = { 0.496, -0.626 };
   int ldc = 1;
   double C_expected[] = { -0.0862, 0.2374 };
   cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1639) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1639) imag");
     };
   };
  };


  {
   int order = 102;
   int uplo = 121;
   int trans = 111;
   int N = 1;
   int K = 2;
   double alpha[2] = {0, 0};
   double beta[2] = {-0.3, 0.1};
   double A[] = { 0.796, 0.872, -0.919, 0.748 };
   int lda = 1;
   double B[] = { -0.945, 0.915, -0.252, -0.276 };
   int ldb = 1;
   double C[] = { 0.07, -0.957 };
   int ldc = 1;
   double C_expected[] = { 0.0747, 0.2941 };
   cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1640) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1640) imag");
     };
   };
  };


  {
   int order = 102;
   int uplo = 122;
   int trans = 111;
   int N = 1;
   int K = 2;
   double alpha[2] = {0, 0};
   double beta[2] = {-0.3, 0.1};
   double A[] = { 0.984, 0.526, 0.284, 0.806 };
   int lda = 1;
   double B[] = { -0.509, -0.178, 0.188, -0.221 };
   int ldb = 1;
   double C[] = { -0.388, 0.795 };
   int ldc = 1;
   double C_expected[] = { 0.0369, -0.2773 };
   cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1641) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1641) imag");
     };
   };
  };


  {
   int order = 101;
   int uplo = 121;
   int trans = 112;
   int N = 1;
   int K = 2;
   double alpha[2] = {1, 0};
   double beta[2] = {0, 0.1};
   double A[] = { 0.628, 0.846, -0.645, 0.032 };
   int lda = 1;
   double B[] = { 0.545, -0.54, 0.493, -0.035 };
   int ldb = 1;
   double C[] = { -0.16, -0.06 };
   int ldc = 1;
   double C_expected[] = { 0.97047, 0.304602 };
   cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1642) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1642) imag");
     };
   };
  };


  {
   int order = 101;
   int uplo = 122;
   int trans = 112;
   int N = 1;
   int K = 2;
   double alpha[2] = {1, 0};
   double beta[2] = {0, 0.1};
   double A[] = { -0.556, -0.946, 0.177, -0.859 };
   int lda = 1;
   double B[] = { 0.423, -0.91, 0.736, -0.251 };
   int ldb = 1;
   double C[] = { -0.478, 0.519 };
   int ldc = 1;
   double C_expected[] = { -2.41467, -1.189498 };
   cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1643) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1643) imag");
     };
   };
  };


  {
   int order = 102;
   int uplo = 121;
   int trans = 112;
   int N = 1;
   int K = 2;
   double alpha[2] = {1, 0};
   double beta[2] = {0, 0.1};
   double A[] = { -0.582, 0.09, -0.176, 0.784 };
   int lda = 2;
   double B[] = { 0.687, -0.859, 0.945, 0.756 };
   int ldb = 2;
   double C[] = { -0.663, -0.186 };
   int ldc = 1;
   double C_expected[] = { -2.144496, 2.272884 };
   cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1644) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1644) imag");
     };
   };
  };


  {
   int order = 102;
   int uplo = 122;
   int trans = 112;
   int N = 1;
   int K = 2;
   double alpha[2] = {1, 0};
   double beta[2] = {0, 0.1};
   double A[] = { 0.231, -0.452, -0.112, -0.837 };
   int lda = 2;
   double B[] = { -0.258, 0.464, -0.224, 0.893 };
   int ldb = 2;
   double C[] = { -0.448, 0.046 };
   int ldc = 1;
   double C_expected[] = { 1.840718, 0.577744 };
   cblas_zsyr2k(order, uplo, trans, N, K, alpha, A, lda, B, ldb, beta, C, ldc);
   {
     int i;
     for (i = 0; i < 1; i++) {
       gsl_test_rel(C[2*i], C_expected[2*i], dbleps, "zsyr2k(case 1645) real");
       gsl_test_rel(C[2*i+1], C_expected[2*i+1], dbleps, "zsyr2k(case 1645) imag");
     };
   };
  };


}