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