htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_syr.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_syr (void) {
  10 const double flteps = 1e-4, dbleps = 1e-6;
  11   {
  12    int order = 101;
  13    int uplo = 121;
  14    int N = 1;
  15    int lda = 1;
  16    float alpha = 0.1f;
  17    float A[] = { -0.291f };
  18    float X[] = { 0.845f };
  19    int incX = -1;
  20    float A_expected[] = { -0.219597f };
  21    cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
  22    {
  23      int i;
  24      for (i = 0; i < 1; i++) {
  25        gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1402)");
  26      }
  27    };
  28   };
  29
  30
  31   {
  32    int order = 101;
  33    int uplo = 122;
  34    int N = 1;
  35    int lda = 1;
  36    float alpha = 0.1f;
  37    float A[] = { -0.291f };
  38    float X[] = { 0.845f };
  39    int incX = -1;
  40    float A_expected[] = { -0.219597f };
  41    cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
  42    {
  43      int i;
  44      for (i = 0; i < 1; i++) {
  45        gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1403)");
  46      }
  47    };
  48   };
  49
  50
  51   {
  52    int order = 102;
  53    int uplo = 121;
  54    int N = 1;
  55    int lda = 1;
  56    float alpha = 0.1f;
  57    float A[] = { -0.291f };
  58    float X[] = { 0.845f };
  59    int incX = -1;
  60    float A_expected[] = { -0.219597f };
  61    cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
  62    {
  63      int i;
  64      for (i = 0; i < 1; i++) {
  65        gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1404)");
  66      }
  67    };
  68   };
  69
  70
  71   {
  72    int order = 102;
  73    int uplo = 122;
  74    int N = 1;
  75    int lda = 1;
  76    float alpha = 0.1f;
  77    float A[] = { -0.291f };
  78    float X[] = { 0.845f };
  79    int incX = -1;
  80    float A_expected[] = { -0.219597f };
  81    cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
  82    {
  83      int i;
  84      for (i = 0; i < 1; i++) {
  85        gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1405)");
  86      }
  87    };
  88   };
  89
  90
  91   {
  92    int order = 101;
  93    int uplo = 121;
  94    int N = 1;
  95    int lda = 1;
  96    double alpha = -0.3;
  97    double A[] = { -0.65 };
  98    double X[] = { -0.891 };
  99    int incX = -1;
 100    double A_expected[] = { -0.8881643 };
 101    cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
 102    {
 103      int i;
 104      for (i = 0; i < 1; i++) {
 105        gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1406)");
 106      }
 107    };
 108   };
 109
 110
 111   {
 112    int order = 101;
 113    int uplo = 122;
 114    int N = 1;
 115    int lda = 1;
 116    double alpha = -0.3;
 117    double A[] = { -0.65 };
 118    double X[] = { -0.891 };
 119    int incX = -1;
 120    double A_expected[] = { -0.8881643 };
 121    cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
 122    {
 123      int i;
 124      for (i = 0; i < 1; i++) {
 125        gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1407)");
 126      }
 127    };
 128   };
 129
 130
 131   {
 132    int order = 102;
 133    int uplo = 121;
 134    int N = 1;
 135    int lda = 1;
 136    double alpha = -0.3;
 137    double A[] = { -0.65 };
 138    double X[] = { -0.891 };
 139    int incX = -1;
 140    double A_expected[] = { -0.8881643 };
 141    cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
 142    {
 143      int i;
 144      for (i = 0; i < 1; i++) {
 145        gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1408)");
 146      }
 147    };
 148   };
 149
 150
 151   {
 152    int order = 102;
 153    int uplo = 122;
 154    int N = 1;
 155    int lda = 1;
 156    double alpha = -0.3;
 157    double A[] = { -0.65 };
 158    double X[] = { -0.891 };
 159    int incX = -1;
 160    double A_expected[] = { -0.8881643 };
 161    cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
 162    {
 163      int i;
 164      for (i = 0; i < 1; i++) {
 165        gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1409)");
 166      }
 167    };
 168   };
 169
 170
 171 }