test_syr2k.c

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

#include <gsl/gsl_test.h>
#include <gsl/gsl_ieee_utils.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_cblas.h>

#include "tests.h"

void
test_syr2k (void) {
const double flteps = 1e-4, dbleps = 1e-6;
  {
   int order = 101;
   int uplo = 121;
   int trans = 111;
   int N = 1;
   int K = 2;
   float alpha = 0.1f;
   float beta = 1.0f;
   float A[] = { -0.915f, 0.445f };
   int lda = 2;
   float B[] = { 0.213f, -0.194f };
   int ldb = 2;
   float C[] = { -0.117f };
   int ldc = 1;
   float C_expected[] = { -0.173245f };
   cblas_ssyr2k(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[i], C_expected[i], flteps, "ssyr2k(case 1614)");
     }
   };
  };


  {
   int order = 101;
   int uplo = 122;
   int trans = 111;
   int N = 1;
   int K = 2;
   float alpha = 0.1f;
   float beta = 1.0f;
   float A[] = { 0.089f, -0.889f };
   int lda = 2;
   float B[] = { -0.384f, 0.518f };
   int ldb = 2;
   float C[] = { 0.069f };
   int ldc = 1;
   float C_expected[] = { -0.0299356f };
   cblas_ssyr2k(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[i], C_expected[i], flteps, "ssyr2k(case 1615)");
     }
   };
  };


  {
   int order = 102;
   int uplo = 121;
   int trans = 111;
   int N = 1;
   int K = 2;
   float alpha = 0.1f;
   float beta = 1.0f;
   float A[] = { 0.492f, 0.021f };
   int lda = 1;
   float B[] = { -0.804f, -0.912f };
   int ldb = 1;
   float C[] = { -0.851f };
   int ldc = 1;
   float C_expected[] = { -0.933944f };
   cblas_ssyr2k(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[i], C_expected[i], flteps, "ssyr2k(case 1616)");
     }
   };
  };


  {
   int order = 102;
   int uplo = 122;
   int trans = 111;
   int N = 1;
   int K = 2;
   float alpha = 0.1f;
   float beta = 1.0f;
   float A[] = { -0.376f, 0.689f };
   int lda = 1;
   float B[] = { 0.21f, 0.406f };
   int ldb = 1;
   float C[] = { -0.581f };
   int ldc = 1;
   float C_expected[] = { -0.540845f };
   cblas_ssyr2k(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[i], C_expected[i], flteps, "ssyr2k(case 1617)");
     }
   };
  };


  {
   int order = 101;
   int uplo = 121;
   int trans = 112;
   int N = 1;
   int K = 2;
   float alpha = 1.0f;
   float beta = -0.3f;
   float A[] = { 0.629f, -0.883f };
   int lda = 1;
   float B[] = { -0.165f, 0.02f };
   int ldb = 1;
   float C[] = { 0.236f };
   int ldc = 1;
   float C_expected[] = { -0.31369f };
   cblas_ssyr2k(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[i], C_expected[i], flteps, "ssyr2k(case 1618)");
     }
   };
  };


  {
   int order = 101;
   int uplo = 122;
   int trans = 112;
   int N = 1;
   int K = 2;
   float alpha = 1.0f;
   float beta = -0.3f;
   float A[] = { 0.412f, -0.411f };
   int lda = 1;
   float B[] = { 0.313f, 0.301f };
   int ldb = 1;
   float C[] = { 0.222f };
   int ldc = 1;
   float C_expected[] = { -0.05611f };
   cblas_ssyr2k(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[i], C_expected[i], flteps, "ssyr2k(case 1619)");
     }
   };
  };


  {
   int order = 102;
   int uplo = 121;
   int trans = 112;
   int N = 1;
   int K = 2;
   float alpha = 1.0f;
   float beta = -0.3f;
   float A[] = { -0.02f, 0.593f };
   int lda = 2;
   float B[] = { -0.144f, 0.846f };
   int ldb = 2;
   float C[] = { -0.645f };
   int ldc = 1;
   float C_expected[] = { 1.20262f };
   cblas_ssyr2k(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[i], C_expected[i], flteps, "ssyr2k(case 1620)");
     }
   };
  };


  {
   int order = 102;
   int uplo = 122;
   int trans = 112;
   int N = 1;
   int K = 2;
   float alpha = 1.0f;
   float beta = -0.3f;
   float A[] = { 0.253f, 0.937f };
   int lda = 2;
   float B[] = { 0.24f, -0.27f };
   int ldb = 2;
   float C[] = { 0.128f };
   int ldc = 1;
   float C_expected[] = { -0.42294f };
   cblas_ssyr2k(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[i], C_expected[i], flteps, "ssyr2k(case 1621)");
     }
   };
  };


  {
   int order = 101;
   int uplo = 121;
   int trans = 111;
   int N = 1;
   int K = 2;
   double alpha = 0.1;
   double beta = 0;
   double A[] = { -0.225, 0.857 };
   int lda = 2;
   double B[] = { -0.933, 0.994 };
   int ldb = 2;
   double C[] = { 0.177 };
   int ldc = 1;
   double C_expected[] = { 0.2123566 };
   cblas_dsyr2k(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[i], C_expected[i], dbleps, "dsyr2k(case 1622)");
     }
   };
  };


  {
   int order = 101;
   int uplo = 122;
   int trans = 111;
   int N = 1;
   int K = 2;
   double alpha = 0.1;
   double beta = 0;
   double A[] = { -0.955, 0.112 };
   int lda = 2;
   double B[] = { -0.695, 0.719 };
   int ldb = 2;
   double C[] = { 0.069 };
   int ldc = 1;
   double C_expected[] = { 0.1488506 };
   cblas_dsyr2k(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[i], C_expected[i], dbleps, "dsyr2k(case 1623)");
     }
   };
  };


  {
   int order = 102;
   int uplo = 121;
   int trans = 111;
   int N = 1;
   int K = 2;
   double alpha = 0.1;
   double beta = 0;
   double A[] = { 0.216, 0.911 };
   int lda = 1;
   double B[] = { -0.074, -0.256 };
   int ldb = 1;
   double C[] = { -0.621 };
   int ldc = 1;
   double C_expected[] = { -0.04984 };
   cblas_dsyr2k(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[i], C_expected[i], dbleps, "dsyr2k(case 1624)");
     }
   };
  };


  {
   int order = 102;
   int uplo = 122;
   int trans = 111;
   int N = 1;
   int K = 2;
   double alpha = 0.1;
   double beta = 0;
   double A[] = { -0.343, -0.381 };
   int lda = 1;
   double B[] = { -0.433, -0.087 };
   int ldb = 1;
   double C[] = { -0.889 };
   int ldc = 1;
   double C_expected[] = { 0.0363332 };
   cblas_dsyr2k(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[i], C_expected[i], dbleps, "dsyr2k(case 1625)");
     }
   };
  };


  {
   int order = 101;
   int uplo = 121;
   int trans = 112;
   int N = 1;
   int K = 2;
   double alpha = 1;
   double beta = -1;
   double A[] = { -0.633, 0.219 };
   int lda = 1;
   double B[] = { 0.817, -0.683 };
   int ldb = 1;
   double C[] = { -0.294 };
   int ldc = 1;
   double C_expected[] = { -1.039476 };
   cblas_dsyr2k(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[i], C_expected[i], dbleps, "dsyr2k(case 1626)");
     }
   };
  };


  {
   int order = 101;
   int uplo = 122;
   int trans = 112;
   int N = 1;
   int K = 2;
   double alpha = 1;
   double beta = -1;
   double A[] = { -0.887, -0.43 };
   int lda = 1;
   double B[] = { 0.557, 0.912 };
   int ldb = 1;
   double C[] = { 0.831 };
   int ldc = 1;
   double C_expected[] = { -2.603438 };
   cblas_dsyr2k(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[i], C_expected[i], dbleps, "dsyr2k(case 1627)");
     }
   };
  };


  {
   int order = 102;
   int uplo = 121;
   int trans = 112;
   int N = 1;
   int K = 2;
   double alpha = 1;
   double beta = -1;
   double A[] = { 0.397, -0.173 };
   int lda = 2;
   double B[] = { 0.155, -0.99 };
   int ldb = 2;
   double C[] = { 0.621 };
   int ldc = 1;
   double C_expected[] = { -0.15539 };
   cblas_dsyr2k(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[i], C_expected[i], dbleps, "dsyr2k(case 1628)");
     }
   };
  };


  {
   int order = 102;
   int uplo = 122;
   int trans = 112;
   int N = 1;
   int K = 2;
   double alpha = 1;
   double beta = -1;
   double A[] = { 0.833, -0.52 };
   int lda = 2;
   double B[] = { 0.28, 0.481 };
   int ldb = 2;
   double C[] = { 0.455 };
   int ldc = 1;
   double C_expected[] = { -0.48876 };
   cblas_dsyr2k(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[i], C_expected[i], dbleps, "dsyr2k(case 1629)");
     }
   };
  };


  {
   int order = 101;
   int uplo = 121;
   int trans = 111;