htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_hemv.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_hemv (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int order = 101;
  13    int uplo = 121;
  14    float alpha[2] = {1.0f, 0.0f};
  15    float beta[2] = {-0.3f, 0.1f};
  16    int N = 1;
  17    int lda = 1;
  18    float A[] = { -0.434f, 0.837f };
  19    float X[] = { 0.209f, -0.935f };
  20    int incX = -1;
  21    float Y[] = { 0.346f, -0.412f };
  22    int incY = -1;
  23    float y_expected[] = { -0.153306f, 0.56399f };
  24    cblas_chemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
  25    {
  26      int i;
  27      for (i = 0; i < 1; i++) {
  28        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chemv(case 1070) real");
  29        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chemv(case 1070) imag");
  30      };
  31    };
  32   };
  33
  34
  35   {
  36    int order = 101;
  37    int uplo = 121;
  38    float alpha[2] = {1.0f, 0.0f};
  39    float beta[2] = {-0.3f, 0.1f};
  40    int N = 1;
  41    int lda = 1;
  42    float A[] = { -0.434f, 0.837f };
  43    float X[] = { 0.209f, -0.935f };
  44    int incX = -1;
  45    float Y[] = { 0.346f, -0.412f };
  46    int incY = -1;
  47    float y_expected[] = { -0.153306f, 0.56399f };
  48    cblas_chemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
  49    {
  50      int i;
  51      for (i = 0; i < 1; i++) {
  52        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chemv(case 1071) real");
  53        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chemv(case 1071) imag");
  54      };
  55    };
  56   };
  57
  58
  59   {
  60    int order = 101;
  61    int uplo = 122;
  62    float alpha[2] = {1.0f, 0.0f};
  63    float beta[2] = {-0.3f, 0.1f};
  64    int N = 1;
  65    int lda = 1;
  66    float A[] = { -0.434f, 0.837f };
  67    float X[] = { 0.209f, -0.935f };
  68    int incX = -1;
  69    float Y[] = { 0.346f, -0.412f };
  70    int incY = -1;
  71    float y_expected[] = { -0.153306f, 0.56399f };
  72    cblas_chemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
  73    {
  74      int i;
  75      for (i = 0; i < 1; i++) {
  76        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chemv(case 1072) real");
  77        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chemv(case 1072) imag");
  78      };
  79    };
  80   };
  81
  82
  83   {
  84    int order = 101;
  85    int uplo = 122;
  86    float alpha[2] = {1.0f, 0.0f};
  87    float beta[2] = {-0.3f, 0.1f};
  88    int N = 1;
  89    int lda = 1;
  90    float A[] = { -0.434f, 0.837f };
  91    float X[] = { 0.209f, -0.935f };
  92    int incX = -1;
  93    float Y[] = { 0.346f, -0.412f };
  94    int incY = -1;
  95    float y_expected[] = { -0.153306f, 0.56399f };
  96    cblas_chemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
  97    {
  98      int i;
  99      for (i = 0; i < 1; i++) {
 100        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chemv(case 1073) real");
 101        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chemv(case 1073) imag");
 102      };
 103    };
 104   };
 105
 106
 107   {
 108    int order = 102;
 109    int uplo = 121;
 110    float alpha[2] = {1.0f, 0.0f};
 111    float beta[2] = {-0.3f, 0.1f};
 112    int N = 1;
 113    int lda = 1;
 114    float A[] = { -0.434f, 0.837f };
 115    float X[] = { 0.209f, -0.935f };
 116    int incX = -1;
 117    float Y[] = { 0.346f, -0.412f };
 118    int incY = -1;
 119    float y_expected[] = { -0.153306f, 0.56399f };
 120    cblas_chemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 121    {
 122      int i;
 123      for (i = 0; i < 1; i++) {
 124        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chemv(case 1074) real");
 125        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chemv(case 1074) imag");
 126      };
 127    };
 128   };
 129
 130
 131   {
 132    int order = 102;
 133    int uplo = 121;
 134    float alpha[2] = {1.0f, 0.0f};
 135    float beta[2] = {-0.3f, 0.1f};
 136    int N = 1;
 137    int lda = 1;
 138    float A[] = { -0.434f, 0.837f };
 139    float X[] = { 0.209f, -0.935f };
 140    int incX = -1;
 141    float Y[] = { 0.346f, -0.412f };
 142    int incY = -1;
 143    float y_expected[] = { -0.153306f, 0.56399f };
 144    cblas_chemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 145    {
 146      int i;
 147      for (i = 0; i < 1; i++) {
 148        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chemv(case 1075) real");
 149        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chemv(case 1075) imag");
 150      };
 151    };
 152   };
 153
 154
 155   {
 156    int order = 102;
 157    int uplo = 122;
 158    float alpha[2] = {1.0f, 0.0f};
 159    float beta[2] = {-0.3f, 0.1f};
 160    int N = 1;
 161    int lda = 1;
 162    float A[] = { -0.434f, 0.837f };
 163    float X[] = { 0.209f, -0.935f };
 164    int incX = -1;
 165    float Y[] = { 0.346f, -0.412f };
 166    int incY = -1;
 167    float y_expected[] = { -0.153306f, 0.56399f };
 168    cblas_chemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 169    {
 170      int i;
 171      for (i = 0; i < 1; i++) {
 172        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chemv(case 1076) real");
 173        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chemv(case 1076) imag");
 174      };
 175    };
 176   };
 177
 178
 179   {
 180    int order = 102;
 181    int uplo = 122;
 182    float alpha[2] = {1.0f, 0.0f};
 183    float beta[2] = {-0.3f, 0.1f};
 184    int N = 1;
 185    int lda = 1;
 186    float A[] = { -0.434f, 0.837f };
 187    float X[] = { 0.209f, -0.935f };
 188    int incX = -1;
 189    float Y[] = { 0.346f, -0.412f };
 190    int incY = -1;
 191    float y_expected[] = { -0.153306f, 0.56399f };
 192    cblas_chemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 193    {
 194      int i;
 195      for (i = 0; i < 1; i++) {
 196        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chemv(case 1077) real");
 197        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chemv(case 1077) imag");
 198      };
 199    };
 200   };
 201
 202
 203   {
 204    int order = 101;
 205    int uplo = 121;
 206    double alpha[2] = {0, 0};
 207    double beta[2] = {1, 0};
 208    int N = 1;
 209    int lda = 1;
 210    double A[] = { 0.036, -0.966 };
 211    double X[] = { -0.695, 0.886 };
 212    int incX = -1;
 213    double Y[] = { 0.486, 0.629 };
 214    int incY = -1;
 215    double y_expected[] = { 0.486, 0.629 };
 216    cblas_zhemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 217    {
 218      int i;
 219      for (i = 0; i < 1; i++) {
 220        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhemv(case 1078) real");
 221        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhemv(case 1078) imag");
 222      };
 223    };
 224   };
 225
 226
 227   {
 228    int order = 101;
 229    int uplo = 121;
 230    double alpha[2] = {0, 0};
 231    double beta[2] = {1, 0};
 232    int N = 1;
 233    int lda = 1;
 234    double A[] = { 0.036, -0.966 };
 235    double X[] = { -0.695, 0.886 };
 236    int incX = -1;
 237    double Y[] = { 0.486, 0.629 };
 238    int incY = -1;
 239    double y_expected[] = { 0.486, 0.629 };
 240    cblas_zhemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 241    {
 242      int i;
 243      for (i = 0; i < 1; i++) {
 244        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhemv(case 1079) real");
 245        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhemv(case 1079) imag");
 246      };
 247    };
 248   };
 249
 250
 251   {
 252    int order = 101;
 253    int uplo = 122;
 254    double alpha[2] = {0, 0};
 255    double beta[2] = {1, 0};
 256    int N = 1;
 257    int lda = 1;
 258    double A[] = { 0.036, -0.966 };
 259    double X[] = { -0.695, 0.886 };
 260    int incX = -1;
 261    double Y[] = { 0.486, 0.629 };
 262    int incY = -1;
 263    double y_expected[] = { 0.486, 0.629 };
 264    cblas_zhemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 265    {
 266      int i;
 267      for (i = 0; i < 1; i++) {
 268        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhemv(case 1080) real");
 269        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhemv(case 1080) imag");
 270      };
 271    };
 272   };
 273
 274
 275   {
 276    int order = 101;
 277    int uplo = 122;
 278    double alpha[2] = {0, 0};
 279    double beta[2] = {1, 0};
 280    int N = 1;
 281    int lda = 1;
 282    double A[] = { 0.036, -0.966 };
 283    double X[] = { -0.695, 0.886 };
 284    int incX = -1;
 285    double Y[] = { 0.486, 0.629 };
 286    int incY = -1;
 287    double y_expected[] = { 0.486, 0.629 };
 288    cblas_zhemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 289    {
 290      int i;
 291      for (i = 0; i < 1; i++) {
 292        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhemv(case 1081) real");
 293        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhemv(case 1081) imag");
 294      };
 295    };
 296   };
 297
 298
 299   {
 300    int order = 102;
 301    int uplo = 121;
 302    double alpha[2] = {0, 0};
 303    double beta[2] = {1, 0};
 304    int N = 1;
 305    int lda = 1;
 306    double A[] = { 0.036, -0.966 };
 307    double X[] = { -0.695, 0.886 };
 308    int incX = -1;
 309    double Y[] = { 0.486, 0.629 };
 310    int incY = -1;
 311    double y_expected[] = { 0.486, 0.629 };
 312    cblas_zhemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 313    {
 314      int i;
 315      for (i = 0; i < 1; i++) {
 316        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhemv(case 1082) real");
 317        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhemv(case 1082) imag");
 318      };
 319    };
 320   };
 321
 322
 323   {
 324    int order = 102;
 325    int uplo = 121;
 326    double alpha[2] = {0, 0};
 327    double beta[2] = {1, 0};
 328    int N = 1;
 329    int lda = 1;
 330    double A[] = { 0.036, -0.966 };
 331    double X[] = { -0.695, 0.886 };
 332    int incX = -1;
 333    double Y[] = { 0.486, 0.629 };
 334    int incY = -1;
 335    double y_expected[] = { 0.486, 0.629 };
 336    cblas_zhemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 337    {
 338      int i;
 339      for (i = 0; i < 1; i++) {
 340        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhemv(case 1083) real");
 341        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhemv(case 1083) imag");
 342      };
 343    };
 344   };
 345
 346
 347   {
 348    int order = 102;
 349    int uplo = 122;
 350    double alpha[2] = {0, 0};
 351    double beta[2] = {1, 0};
 352    int N = 1;
 353    int lda = 1;
 354    double A[] = { 0.036, -0.966 };
 355    double X[] = { -0.695, 0.886 };
 356    int incX = -1;
 357    double Y[] = { 0.486, 0.629 };
 358    int incY = -1;
 359    double y_expected[] = { 0.486, 0.629 };
 360    cblas_zhemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 361    {
 362      int i;
 363      for (i = 0; i < 1; i++) {
 364        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhemv(case 1084) real");
 365        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhemv(case 1084) imag");
 366      };
 367    };
 368   };
 369
 370
 371   {
 372    int order = 102;
 373    int uplo = 122;
 374    double alpha[2] = {0, 0};
 375    double beta[2] = {1, 0};
 376    int N = 1;
 377    int lda = 1;
 378    double A[] = { 0.036, -0.966 };
 379    double X[] = { -0.695, 0.886 };
 380    int incX = -1;
 381    double Y[] = { 0.486, 0.629 };
 382    int incY = -1;
 383    double y_expected[] = { 0.486, 0.629 };
 384    cblas_zhemv(order, uplo, N, alpha, A, lda, X, incX, beta, Y, incY);
 385    {
 386      int i;
 387      for (i = 0; i < 1; i++) {
 388        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhemv(case 1085) real");
 389        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhemv(case 1085) imag");
 390      };
 391    };
 392   };
 393
 394
 395 }