htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_copy.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_copy (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int N = 1;
  13    float X[] = { 0.898f };
  14    int incX = 1;
  15    float Y[] = { 0.699f };
  16    int incY = -1;
  17    float expected[] = { 0.898f };
  18    cblas_scopy(N, X, incX, Y, incY);
  19    {
  20      int i;
  21      for (i = 0; i < 1; i++) {
  22        gsl_test_rel(Y[i], expected[i], flteps, "scopy(case 76)");
  23      }
  24    };
  25   };
  26
  27
  28   {
  29    int N = 1;
  30    double X[] = { 0.002 };
  31    int incX = 1;
  32    double Y[] = { -0.921 };
  33    int incY = -1;
  34    double expected[] = { 0.002 };
  35    cblas_dcopy(N, X, incX, Y, incY);
  36    {
  37      int i;
  38      for (i = 0; i < 1; i++) {
  39        gsl_test_rel(Y[i], expected[i], dbleps, "dcopy(case 77)");
  40      }
  41    };
  42   };
  43
  44
  45   {
  46    int N = 1;
  47    float X[] = { -0.166f, 0.639f };
  48    int incX = 1;
  49    float Y[] = { 0.863f, 0.613f };
  50    int incY = -1;
  51    float expected[] = { -0.166f, 0.639f };
  52    cblas_ccopy(N, X, incX, Y, incY);
  53    {
  54      int i;
  55      for (i = 0; i < 1; i++) {
  56        gsl_test_rel(Y[2*i], expected[2*i], flteps, "ccopy(case 78) real");
  57        gsl_test_rel(Y[2*i+1], expected[2*i+1], flteps, "ccopy(case 78) imag");
  58      };
  59    };
  60   };
  61
  62
  63   {
  64    int N = 1;
  65    double X[] = { 0.315, -0.324 };
  66    int incX = 1;
  67    double Y[] = { -0.312, -0.748 };
  68    int incY = -1;
  69    double expected[] = { 0.315, -0.324 };
  70    cblas_zcopy(N, X, incX, Y, incY);
  71    {
  72      int i;
  73      for (i = 0; i < 1; i++) {
  74        gsl_test_rel(Y[2*i], expected[2*i], dbleps, "zcopy(case 79) real");
  75        gsl_test_rel(Y[2*i+1], expected[2*i+1], dbleps, "zcopy(case 79) imag");
  76      };
  77    };
  78   };
  79
  80
  81   {
  82    int N = 1;
  83    float X[] = { 0.222f };
  84    int incX = -1;
  85    float Y[] = { 0.522f };
  86    int incY = 1;
  87    float expected[] = { 0.222f };
  88    cblas_scopy(N, X, incX, Y, incY);
  89    {
  90      int i;
  91      for (i = 0; i < 1; i++) {
  92        gsl_test_rel(Y[i], expected[i], flteps, "scopy(case 80)");
  93      }
  94    };
  95   };
  96
  97
  98   {
  99    int N = 1;
 100    double X[] = { 0.021 };
 101    int incX = -1;
 102    double Y[] = { 0.898 };
 103    int incY = 1;
 104    double expected[] = { 0.021 };
 105    cblas_dcopy(N, X, incX, Y, incY);
 106    {
 107      int i;
 108      for (i = 0; i < 1; i++) {
 109        gsl_test_rel(Y[i], expected[i], dbleps, "dcopy(case 81)");
 110      }
 111    };
 112   };
 113
 114
 115   {
 116    int N = 1;
 117    float X[] = { 0.376f, 0.229f };
 118    int incX = -1;
 119    float Y[] = { 0.143f, -0.955f };
 120    int incY = 1;
 121    float expected[] = { 0.376f, 0.229f };
 122    cblas_ccopy(N, X, incX, Y, incY);
 123    {
 124      int i;
 125      for (i = 0; i < 1; i++) {
 126        gsl_test_rel(Y[2*i], expected[2*i], flteps, "ccopy(case 82) real");
 127        gsl_test_rel(Y[2*i+1], expected[2*i+1], flteps, "ccopy(case 82) imag");
 128      };
 129    };
 130   };
 131
 132
 133   {
 134    int N = 1;
 135    double X[] = { -0.265, -0.84 };
 136    int incX = -1;
 137    double Y[] = { -0.156, 0.939 };
 138    int incY = 1;
 139    double expected[] = { -0.265, -0.84 };
 140    cblas_zcopy(N, X, incX, Y, incY);
 141    {
 142      int i;
 143      for (i = 0; i < 1; i++) {
 144        gsl_test_rel(Y[2*i], expected[2*i], dbleps, "zcopy(case 83) real");
 145        gsl_test_rel(Y[2*i+1], expected[2*i+1], dbleps, "zcopy(case 83) imag");
 146      };
 147    };
 148   };
 149
 150
 151   {
 152    int N = 1;
 153    float X[] = { 0.074f };
 154    int incX = -1;
 155    float Y[] = { -0.802f };
 156    int incY = -1;
 157    float expected[] = { 0.074f };
 158    cblas_scopy(N, X, incX, Y, incY);
 159    {
 160      int i;
 161      for (i = 0; i < 1; i++) {
 162        gsl_test_rel(Y[i], expected[i], flteps, "scopy(case 84)");
 163      }
 164    };
 165   };
 166
 167
 168   {
 169    int N = 1;
 170    double X[] = { -0.374 };
 171    int incX = -1;
 172    double Y[] = { -0.161 };
 173    int incY = -1;
 174    double expected[] = { -0.374 };
 175    cblas_dcopy(N, X, incX, Y, incY);
 176    {
 177      int i;
 178      for (i = 0; i < 1; i++) {
 179        gsl_test_rel(Y[i], expected[i], dbleps, "dcopy(case 85)");
 180      }
 181    };
 182   };
 183
 184
 185   {
 186    int N = 1;
 187    float X[] = { 0.084f, 0.778f };
 188    int incX = -1;
 189    float Y[] = { 0.31f, -0.797f };
 190    int incY = -1;
 191    float expected[] = { 0.084f, 0.778f };
 192    cblas_ccopy(N, X, incX, Y, incY);
 193    {
 194      int i;
 195      for (i = 0; i < 1; i++) {
 196        gsl_test_rel(Y[2*i], expected[2*i], flteps, "ccopy(case 86) real");
 197        gsl_test_rel(Y[2*i+1], expected[2*i+1], flteps, "ccopy(case 86) imag");
 198      };
 199    };
 200   };
 201
 202
 203   {
 204    int N = 1;
 205    double X[] = { 0.831, -0.282 };
 206    int incX = -1;
 207    double Y[] = { -0.62, 0.32 };
 208    int incY = -1;
 209    double expected[] = { 0.831, -0.282 };
 210    cblas_zcopy(N, X, incX, Y, incY);
 211    {
 212      int i;
 213      for (i = 0; i < 1; i++) {
 214        gsl_test_rel(Y[2*i], expected[2*i], dbleps, "zcopy(case 87) real");
 215        gsl_test_rel(Y[2*i+1], expected[2*i+1], dbleps, "zcopy(case 87) imag");
 216      };
 217    };
 218   };
 219
 220
 221 }