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