htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_her2.c

   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>
   5
   6 #include "tests.h"
   7
   8 void
   9 test_her2 (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int order = 101;
  13    int uplo = 121;
  14    int N = 1;
  15    int lda = 1;
  16    float alpha[2] = {-1.0f, 0.0f};
  17    float A[] = { -0.821f, 0.954f };
  18    float X[] = { 0.532f, 0.802f };
  19    int incX = -1;
  20    float Y[] = { 0.016f, -0.334f };
  21    int incY = -1;
  22    float A_expected[] = { -0.302288f, 0.0f };
  23    cblas_cher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  24    {
  25      int i;
  26      for (i = 0; i < 1; i++) {
  27        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher2(case 1450) real");
  28        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher2(case 1450) imag");
  29      };
  30    };
  31   };
  32
  33
  34   {
  35    int order = 101;
  36    int uplo = 122;
  37    int N = 1;
  38    int lda = 1;
  39    float alpha[2] = {-1.0f, 0.0f};
  40    float A[] = { -0.821f, 0.954f };
  41    float X[] = { 0.532f, 0.802f };
  42    int incX = -1;
  43    float Y[] = { 0.016f, -0.334f };
  44    int incY = -1;
  45    float A_expected[] = { -0.302288f, 0.0f };
  46    cblas_cher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  47    {
  48      int i;
  49      for (i = 0; i < 1; i++) {
  50        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher2(case 1451) real");
  51        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher2(case 1451) imag");
  52      };
  53    };
  54   };
  55
  56
  57   {
  58    int order = 102;
  59    int uplo = 121;
  60    int N = 1;
  61    int lda = 1;
  62    float alpha[2] = {-1.0f, 0.0f};
  63    float A[] = { -0.821f, 0.954f };
  64    float X[] = { 0.532f, 0.802f };
  65    int incX = -1;
  66    float Y[] = { 0.016f, -0.334f };
  67    int incY = -1;
  68    float A_expected[] = { -0.302288f, 0.0f };
  69    cblas_cher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  70    {
  71      int i;
  72      for (i = 0; i < 1; i++) {
  73        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher2(case 1452) real");
  74        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher2(case 1452) imag");
  75      };
  76    };
  77   };
  78
  79
  80   {
  81    int order = 102;
  82    int uplo = 122;
  83    int N = 1;
  84    int lda = 1;
  85    float alpha[2] = {-1.0f, 0.0f};
  86    float A[] = { -0.821f, 0.954f };
  87    float X[] = { 0.532f, 0.802f };
  88    int incX = -1;
  89    float Y[] = { 0.016f, -0.334f };
  90    int incY = -1;
  91    float A_expected[] = { -0.302288f, 0.0f };
  92    cblas_cher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  93    {
  94      int i;
  95      for (i = 0; i < 1; i++) {
  96        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher2(case 1453) real");
  97        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher2(case 1453) imag");
  98      };
  99    };
 100   };
 101
 102
 103   {
 104    int order = 101;
 105    int uplo = 121;
 106    int N = 1;
 107    int lda = 1;
 108    double alpha[2] = {-0.3, 0.1};
 109    double A[] = { -0.334, 0.286 };
 110    double X[] = { -0.14, -0.135 };
 111    int incX = -1;
 112    double Y[] = { 0.455, 0.358 };
 113    int incY = -1;
 114    double A_expected[] = { -0.264521, 0.0 };
 115    cblas_zher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
 116    {
 117      int i;
 118      for (i = 0; i < 1; i++) {
 119        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher2(case 1454) real");
 120        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher2(case 1454) imag");
 121      };
 122    };
 123   };
 124
 125
 126   {
 127    int order = 101;
 128    int uplo = 122;
 129    int N = 1;
 130    int lda = 1;
 131    double alpha[2] = {-0.3, 0.1};
 132    double A[] = { -0.334, 0.286 };
 133    double X[] = { -0.14, -0.135 };
 134    int incX = -1;
 135    double Y[] = { 0.455, 0.358 };
 136    int incY = -1;
 137    double A_expected[] = { -0.264521, 0.0 };
 138    cblas_zher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
 139    {
 140      int i;
 141      for (i = 0; i < 1; i++) {
 142        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher2(case 1455) real");
 143        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher2(case 1455) imag");
 144      };
 145    };
 146   };
 147
 148
 149   {
 150    int order = 102;
 151    int uplo = 121;
 152    int N = 1;
 153    int lda = 1;
 154    double alpha[2] = {-0.3, 0.1};
 155    double A[] = { -0.334, 0.286 };
 156    double X[] = { -0.14, -0.135 };
 157    int incX = -1;
 158    double Y[] = { 0.455, 0.358 };
 159    int incY = -1;
 160    double A_expected[] = { -0.264521, 0.0 };
 161    cblas_zher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
 162    {
 163      int i;
 164      for (i = 0; i < 1; i++) {
 165        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher2(case 1456) real");
 166        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher2(case 1456) imag");
 167      };
 168    };
 169   };
 170
 171
 172   {
 173    int order = 102;
 174    int uplo = 122;
 175    int N = 1;
 176    int lda = 1;
 177    double alpha[2] = {-0.3, 0.1};
 178    double A[] = { -0.334, 0.286 };
 179    double X[] = { -0.14, -0.135 };
 180    int incX = -1;
 181    double Y[] = { 0.455, 0.358 };
 182    int incY = -1;
 183    double A_expected[] = { -0.264521, 0.0 };
 184    cblas_zher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
 185    {
 186      int i;
 187      for (i = 0; i < 1; i++) {
 188        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher2(case 1457) real");
 189        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher2(case 1457) imag");
 190      };
 191    };
 192   };
 193
 194
 195 }