htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_spr.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_spr (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int order = 101;
  13    int uplo = 121;
  14    int N = 2;
  15    float alpha = -0.3f;
  16    float Ap[] = { -0.764f, -0.257f, -0.064f };
  17    float X[] = { 0.455f, -0.285f };
  18    int incX = -1;
  19    float Ap_expected[] = { -0.788367f, -0.218097f, -0.126108f };
  20    cblas_sspr(order, uplo, N, alpha, X, incX, Ap);
  21    {
  22      int i;
  23      for (i = 0; i < 3; i++) {
  24        gsl_test_rel(Ap[i], Ap_expected[i], flteps, "sspr(case 1426)");
  25      }
  26    };
  27   };
  28
  29
  30   {
  31    int order = 101;
  32    int uplo = 122;
  33    int N = 2;
  34    float alpha = -0.3f;
  35    float Ap[] = { -0.764f, -0.257f, -0.064f };
  36    float X[] = { 0.455f, -0.285f };
  37    int incX = -1;
  38    float Ap_expected[] = { -0.788367f, -0.218097f, -0.126108f };
  39    cblas_sspr(order, uplo, N, alpha, X, incX, Ap);
  40    {
  41      int i;
  42      for (i = 0; i < 3; i++) {
  43        gsl_test_rel(Ap[i], Ap_expected[i], flteps, "sspr(case 1427)");
  44      }
  45    };
  46   };
  47
  48
  49   {
  50    int order = 102;
  51    int uplo = 121;
  52    int N = 2;
  53    float alpha = -0.3f;
  54    float Ap[] = { -0.764f, -0.257f, -0.064f };
  55    float X[] = { 0.455f, -0.285f };
  56    int incX = -1;
  57    float Ap_expected[] = { -0.788367f, -0.218097f, -0.126108f };
  58    cblas_sspr(order, uplo, N, alpha, X, incX, Ap);
  59    {
  60      int i;
  61      for (i = 0; i < 3; i++) {
  62        gsl_test_rel(Ap[i], Ap_expected[i], flteps, "sspr(case 1428)");
  63      }
  64    };
  65   };
  66
  67
  68   {
  69    int order = 102;
  70    int uplo = 122;
  71    int N = 2;
  72    float alpha = -0.3f;
  73    float Ap[] = { -0.764f, -0.257f, -0.064f };
  74    float X[] = { 0.455f, -0.285f };
  75    int incX = -1;
  76    float Ap_expected[] = { -0.788367f, -0.218097f, -0.126108f };
  77    cblas_sspr(order, uplo, N, alpha, X, incX, Ap);
  78    {
  79      int i;
  80      for (i = 0; i < 3; i++) {
  81        gsl_test_rel(Ap[i], Ap_expected[i], flteps, "sspr(case 1429)");
  82      }
  83    };
  84   };
  85
  86
  87   {
  88    int order = 101;
  89    int uplo = 121;
  90    int N = 2;
  91    double alpha = -1;
  92    double Ap[] = { 0.819, 0.175, -0.809 };
  93    double X[] = { -0.645, -0.222 };
  94    int incX = -1;
  95    double Ap_expected[] = { 0.769716, 0.03181, -1.225025 };
  96    cblas_dspr(order, uplo, N, alpha, X, incX, Ap);
  97    {
  98      int i;
  99      for (i = 0; i < 3; i++) {
 100        gsl_test_rel(Ap[i], Ap_expected[i], dbleps, "dspr(case 1430)");
 101      }
 102    };
 103   };
 104
 105
 106   {
 107    int order = 101;
 108    int uplo = 122;
 109    int N = 2;
 110    double alpha = -1;
 111    double Ap[] = { 0.819, 0.175, -0.809 };
 112    double X[] = { -0.645, -0.222 };
 113    int incX = -1;
 114    double Ap_expected[] = { 0.769716, 0.03181, -1.225025 };
 115    cblas_dspr(order, uplo, N, alpha, X, incX, Ap);
 116    {
 117      int i;
 118      for (i = 0; i < 3; i++) {
 119        gsl_test_rel(Ap[i], Ap_expected[i], dbleps, "dspr(case 1431)");
 120      }
 121    };
 122   };
 123
 124
 125   {
 126    int order = 102;
 127    int uplo = 121;
 128    int N = 2;
 129    double alpha = -1;
 130    double Ap[] = { 0.819, 0.175, -0.809 };
 131    double X[] = { -0.645, -0.222 };
 132    int incX = -1;
 133    double Ap_expected[] = { 0.769716, 0.03181, -1.225025 };
 134    cblas_dspr(order, uplo, N, alpha, X, incX, Ap);
 135    {
 136      int i;
 137      for (i = 0; i < 3; i++) {
 138        gsl_test_rel(Ap[i], Ap_expected[i], dbleps, "dspr(case 1432)");
 139      }
 140    };
 141   };
 142
 143
 144   {
 145    int order = 102;
 146    int uplo = 122;
 147    int N = 2;
 148    double alpha = -1;
 149    double Ap[] = { 0.819, 0.175, -0.809 };
 150    double X[] = { -0.645, -0.222 };
 151    int incX = -1;
 152    double Ap_expected[] = { 0.769716, 0.03181, -1.225025 };
 153    cblas_dspr(order, uplo, N, alpha, X, incX, Ap);
 154    {
 155      int i;
 156      for (i = 0; i < 3; i++) {
 157        gsl_test_rel(Ap[i], Ap_expected[i], dbleps, "dspr(case 1433)");
 158      }
 159    };
 160   };
 161
 162
 163 }