htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_swap.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_swap (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int N = 1;
  13    float X[] = { 0.539f };
  14    int incX = 1;
  15    float Y[] = { -0.262f };
  16    int incY = -1;
  17    float expected1[] = { -0.262f };
  18    float expected2[] = { 0.539f };
  19    cblas_sswap(N, X, incX, Y, incY);
  20    {
  21      int i;
  22      for (i = 0; i < 1; i++) {
  23        gsl_test_rel(X[i], expected1[i], flteps, "sswap(case 88)");
  24      }
  25    };
  26    {
  27      int i;
  28      for (i = 0; i < 1; i++) {
  29        gsl_test_rel(Y[i], expected2[i], flteps, "sswap(case 89)");
  30      }
  31    };
  32   };
  33
  34
  35   {
  36    int N = 1;
  37    double X[] = { 0.906 };
  38    int incX = 1;
  39    double Y[] = { 0.373 };
  40    int incY = -1;
  41    double expected1[] = { 0.373 };
  42    double expected2[] = { 0.906 };
  43    cblas_dswap(N, X, incX, Y, incY);
  44    {
  45      int i;
  46      for (i = 0; i < 1; i++) {
  47        gsl_test_rel(X[i], expected1[i], dbleps, "dswap(case 90)");
  48      }
  49    };
  50    {
  51      int i;
  52      for (i = 0; i < 1; i++) {
  53        gsl_test_rel(Y[i], expected2[i], dbleps, "dswap(case 91)");
  54      }
  55    };
  56   };
  57
  58
  59   {
  60    int N = 1;
  61    float X[] = { -0.316f, -0.529f };
  62    int incX = 1;
  63    float Y[] = { -0.313f, 0.363f };
  64    int incY = -1;
  65    float expected1[] = { -0.313f, 0.363f };
  66    float expected2[] = { -0.316f, -0.529f };
  67    cblas_cswap(N, X, incX, Y, incY);
  68    {
  69      int i;
  70      for (i = 0; i < 1; i++) {
  71        gsl_test_rel(X[2*i], expected1[2*i], flteps, "cswap(case 92) real");
  72        gsl_test_rel(X[2*i+1], expected1[2*i+1], flteps, "cswap(case 92) imag");
  73      };
  74    };
  75    {
  76      int i;
  77      for (i = 0; i < 1; i++) {
  78        gsl_test_rel(Y[2*i], expected2[2*i], flteps, "cswap(case 93) real");
  79        gsl_test_rel(Y[2*i+1], expected2[2*i+1], flteps, "cswap(case 93) imag");
  80      };
  81    };
  82   };
  83
  84
  85   {
  86    int N = 1;
  87    double X[] = { 0.512, -0.89 };
  88    int incX = 1;
  89    double Y[] = { -0.225, -0.511 };
  90    int incY = -1;
  91    double expected1[] = { -0.225, -0.511 };
  92    double expected2[] = { 0.512, -0.89 };
  93    cblas_zswap(N, X, incX, Y, incY);
  94    {
  95      int i;
  96      for (i = 0; i < 1; i++) {
  97        gsl_test_rel(X[2*i], expected1[2*i], dbleps, "zswap(case 94) real");
  98        gsl_test_rel(X[2*i+1], expected1[2*i+1], dbleps, "zswap(case 94) imag");
  99      };
 100    };
 101    {
 102      int i;
 103      for (i = 0; i < 1; i++) {
 104        gsl_test_rel(Y[2*i], expected2[2*i], dbleps, "zswap(case 95) real");
 105        gsl_test_rel(Y[2*i+1], expected2[2*i+1], dbleps, "zswap(case 95) imag");
 106      };
 107    };
 108   };
 109
 110
 111   {
 112    int N = 1;
 113    float X[] = { 0.336f };
 114    int incX = -1;
 115    float Y[] = { -0.431f };
 116    int incY = 1;
 117    float expected1[] = { -0.431f };
 118    float expected2[] = { 0.336f };
 119    cblas_sswap(N, X, incX, Y, incY);
 120    {
 121      int i;
 122      for (i = 0; i < 1; i++) {
 123        gsl_test_rel(X[i], expected1[i], flteps, "sswap(case 96)");
 124      }
 125    };
 126    {
 127      int i;
 128      for (i = 0; i < 1; i++) {
 129        gsl_test_rel(Y[i], expected2[i], flteps, "sswap(case 97)");
 130      }
 131    };
 132   };
 133
 134
 135   {
 136    int N = 1;
 137    double X[] = { 0.764 };
 138    int incX = -1;
 139    double Y[] = { -0.293 };
 140    int incY = 1;
 141    double expected1[] = { -0.293 };
 142    double expected2[] = { 0.764 };
 143    cblas_dswap(N, X, incX, Y, incY);
 144    {
 145      int i;
 146      for (i = 0; i < 1; i++) {
 147        gsl_test_rel(X[i], expected1[i], dbleps, "dswap(case 98)");
 148      }
 149    };
 150    {
 151      int i;
 152      for (i = 0; i < 1; i++) {
 153        gsl_test_rel(Y[i], expected2[i], dbleps, "dswap(case 99)");
 154      }
 155    };
 156   };
 157
 158
 159   {
 160    int N = 1;
 161    float X[] = { -0.239f, 0.361f };
 162    int incX = -1;
 163    float Y[] = { 0.149f, 0.347f };
 164    int incY = 1;
 165    float expected1[] = { 0.149f, 0.347f };
 166    float expected2[] = { -0.239f, 0.361f };
 167    cblas_cswap(N, X, incX, Y, incY);
 168    {
 169      int i;
 170      for (i = 0; i < 1; i++) {
 171        gsl_test_rel(X[2*i], expected1[2*i], flteps, "cswap(case 100) real");
 172        gsl_test_rel(X[2*i+1], expected1[2*i+1], flteps, "cswap(case 100) imag");
 173      };
 174    };
 175    {
 176      int i;
 177      for (i = 0; i < 1; i++) {
 178        gsl_test_rel(Y[2*i], expected2[2*i], flteps, "cswap(case 101) real");
 179        gsl_test_rel(Y[2*i+1], expected2[2*i+1], flteps, "cswap(case 101) imag");
 180      };
 181    };
 182   };
 183
 184
 185   {
 186    int N = 1;
 187    double X[] = { -0.171, -0.936 };
 188    int incX = -1;
 189    double Y[] = { 0.495, -0.835 };
 190    int incY = 1;
 191    double expected1[] = { 0.495, -0.835 };
 192    double expected2[] = { -0.171, -0.936 };
 193    cblas_zswap(N, X, incX, Y, incY);
 194    {
 195      int i;
 196      for (i = 0; i < 1; i++) {
 197        gsl_test_rel(X[2*i], expected1[2*i], dbleps, "zswap(case 102) real");
 198        gsl_test_rel(X[2*i+1], expected1[2*i+1], dbleps, "zswap(case 102) imag");
 199      };
 200    };
 201    {
 202      int i;
 203      for (i = 0; i < 1; i++) {
 204        gsl_test_rel(Y[2*i], expected2[2*i], dbleps, "zswap(case 103) real");
 205        gsl_test_rel(Y[2*i+1], expected2[2*i+1], dbleps, "zswap(case 103) imag");
 206      };
 207    };
 208   };
 209
 210
 211   {
 212    int N = 1;
 213    float X[] = { -0.405f };
 214    int incX = -1;
 215    float Y[] = { -0.213f };
 216    int incY = -1;
 217    float expected1[] = { -0.213f };
 218    float expected2[] = { -0.405f };
 219    cblas_sswap(N, X, incX, Y, incY);
 220    {
 221      int i;
 222      for (i = 0; i < 1; i++) {
 223        gsl_test_rel(X[i], expected1[i], flteps, "sswap(case 104)");
 224      }
 225    };
 226    {
 227      int i;
 228      for (i = 0; i < 1; i++) {
 229        gsl_test_rel(Y[i], expected2[i], flteps, "sswap(case 105)");
 230      }
 231    };
 232   };
 233
 234
 235   {
 236    int N = 1;
 237    double X[] = { -0.761 };
 238    int incX = -1;
 239    double Y[] = { -0.585 };
 240    int incY = -1;
 241    double expected1[] = { -0.585 };
 242    double expected2[] = { -0.761 };
 243    cblas_dswap(N, X, incX, Y, incY);
 244    {
 245      int i;
 246      for (i = 0; i < 1; i++) {
 247        gsl_test_rel(X[i], expected1[i], dbleps, "dswap(case 106)");
 248      }
 249    };
 250    {
 251      int i;
 252      for (i = 0; i < 1; i++) {
 253        gsl_test_rel(Y[i], expected2[i], dbleps, "dswap(case 107)");
 254      }
 255    };
 256   };
 257
 258
 259   {
 260    int N = 1;
 261    float X[] = { 0.853f, 0.146f };
 262    int incX = -1;
 263    float Y[] = { 0.009f, -0.178f };
 264    int incY = -1;
 265    float expected1[] = { 0.009f, -0.178f };
 266    float expected2[] = { 0.853f, 0.146f };
 267    cblas_cswap(N, X, incX, Y, incY);
 268    {
 269      int i;
 270      for (i = 0; i < 1; i++) {
 271        gsl_test_rel(X[2*i], expected1[2*i], flteps, "cswap(case 108) real");
 272        gsl_test_rel(X[2*i+1], expected1[2*i+1], flteps, "cswap(case 108) imag");
 273      };
 274    };
 275    {
 276      int i;
 277      for (i = 0; i < 1; i++) {
 278        gsl_test_rel(Y[2*i], expected2[2*i], flteps, "cswap(case 109) real");
 279        gsl_test_rel(Y[2*i+1], expected2[2*i+1], flteps, "cswap(case 109) imag");
 280      };
 281    };
 282   };
 283
 284
 285   {
 286    int N = 1;
 287    double X[] = { -0.228, 0.386 };
 288    int incX = -1;
 289    double Y[] = { 0.988, -0.084 };
 290    int incY = -1;
 291    double expected1[] = { 0.988, -0.084 };
 292    double expected2[] = { -0.228, 0.386 };
 293    cblas_zswap(N, X, incX, Y, incY);
 294    {
 295      int i;
 296      for (i = 0; i < 1; i++) {
 297        gsl_test_rel(X[2*i], expected1[2*i], dbleps, "zswap(case 110) real");
 298        gsl_test_rel(X[2*i+1], expected1[2*i+1], dbleps, "zswap(case 110) imag");
 299      };
 300    };
 301    {
 302      int i;
 303      for (i = 0; i < 1; i++) {
 304        gsl_test_rel(Y[2*i], expected2[2*i], dbleps, "zswap(case 111) real");
 305        gsl_test_rel(Y[2*i+1], expected2[2*i+1], dbleps, "zswap(case 111) imag");
 306      };
 307    };
 308   };
 309
 310
 311 }