htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_asum.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_asum (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int N = 1;
  13    float X[] = { 0.239f };
  14    int incX = -1;
  15    float expected = 0.0f;
  16    float f;
  17    f = cblas_sasum(N, X, incX);
  18    gsl_test_rel(f, expected, flteps, "sasum(case 40)");
  19   };
  20
  21
  22   {
  23    int N = 1;
  24    double X[] = { -0.413 };
  25    int incX = -1;
  26    double expected = 0;
  27    double f;
  28    f = cblas_dasum(N, X, incX);
  29    gsl_test_rel(f, expected, dbleps, "dasum(case 41)");
  30   };
  31
  32
  33   {
  34    int N = 1;
  35    float X[] = { 0.1f, 0.017f };
  36    int incX = -1;
  37    float expected = 0.0f;
  38    float f;
  39    f = cblas_scasum(N, X, incX);
  40    gsl_test_rel(f, expected, flteps, "scasum(case 42)");
  41   };
  42
  43
  44   {
  45    int N = 1;
  46    double X[] = { -0.651, 0.079 };
  47    int incX = -1;
  48    double expected = 0;
  49    double f;
  50    f = cblas_dzasum(N, X, incX);
  51    gsl_test_rel(f, expected, dbleps, "dzasum(case 43)");
  52   };
  53
  54
  55   {
  56    int N = 2;
  57    float X[] = { 0.899f, -0.72f };
  58    int incX = 1;
  59    float expected = 1.619f;
  60    float f;
  61    f = cblas_sasum(N, X, incX);
  62    gsl_test_rel(f, expected, flteps, "sasum(case 44)");
  63   };
  64
  65
  66   {
  67    int N = 2;
  68    double X[] = { 0.271, -0.012 };
  69    int incX = 1;
  70    double expected = 0.283;
  71    double f;
  72    f = cblas_dasum(N, X, incX);
  73    gsl_test_rel(f, expected, dbleps, "dasum(case 45)");
  74   };
  75
  76
  77   {
  78    int N = 2;
  79    float X[] = { -0.567f, -0.645f, 0.098f, 0.256f };
  80    int incX = 1;
  81    float expected = 1.566f;
  82    float f;
  83    f = cblas_scasum(N, X, incX);
  84    gsl_test_rel(f, expected, flteps, "scasum(case 46)");
  85   };
  86
  87
  88   {
  89    int N = 2;
  90    double X[] = { -0.046, -0.671, -0.323, 0.785 };
  91    int incX = 1;
  92    double expected = 1.825;
  93    double f;
  94    f = cblas_dzasum(N, X, incX);
  95    gsl_test_rel(f, expected, dbleps, "dzasum(case 47)");
  96   };
  97
  98
  99   {
 100    int N = 2;
 101    float X[] = { 0.169f, 0.833f };
 102    int incX = -1;
 103    float expected = 0.0f;
 104    float f;
 105    f = cblas_sasum(N, X, incX);
 106    gsl_test_rel(f, expected, flteps, "sasum(case 48)");
 107   };
 108
 109
 110   {
 111    int N = 2;
 112    double X[] = { -0.586, -0.486 };
 113    int incX = -1;
 114    double expected = 0;
 115    double f;
 116    f = cblas_dasum(N, X, incX);
 117    gsl_test_rel(f, expected, dbleps, "dasum(case 49)");
 118   };
 119
 120
 121   {
 122    int N = 2;
 123    float X[] = { -0.314f, -0.318f, -0.835f, -0.807f };
 124    int incX = -1;
 125    float expected = 0.0f;
 126    float f;
 127    f = cblas_scasum(N, X, incX);
 128    gsl_test_rel(f, expected, flteps, "scasum(case 50)");
 129   };
 130
 131
 132   {
 133    int N = 2;
 134    double X[] = { -0.927, 0.152, -0.554, -0.844 };
 135    int incX = -1;
 136    double expected = 0;
 137    double f;
 138    f = cblas_dzasum(N, X, incX);
 139    gsl_test_rel(f, expected, dbleps, "dzasum(case 51)");
 140   };
 141
 142
 143 }