htswanalysis/MACS/lib/gsl/gsl-1.11/cblas/test_her.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_her (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 = 1.0f;
  17    float A[] = { 0.188f, 0.856f };
  18    float X[] = { -0.832f, -0.151f };
  19    int incX = -1;
  20    float A_expected[] = { 0.903025f, 0.0f };
  21    cblas_cher(order, uplo, N, alpha, X, incX, A, lda);
  22    {
  23      int i;
  24      for (i = 0; i < 1; i++) {
  25        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher(case 1410) real");
  26        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher(case 1410) imag");
  27      };
  28    };
  29   };
  30
  31
  32   {
  33    int order = 101;
  34    int uplo = 122;
  35    int N = 1;
  36    int lda = 1;
  37    float alpha = 1.0f;
  38    float A[] = { 0.188f, 0.856f };
  39    float X[] = { -0.832f, -0.151f };
  40    int incX = -1;
  41    float A_expected[] = { 0.903025f, 0.0f };
  42    cblas_cher(order, uplo, N, alpha, X, incX, A, lda);
  43    {
  44      int i;
  45      for (i = 0; i < 1; i++) {
  46        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher(case 1411) real");
  47        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher(case 1411) imag");
  48      };
  49    };
  50   };
  51
  52
  53   {
  54    int order = 102;
  55    int uplo = 121;
  56    int N = 1;
  57    int lda = 1;
  58    float alpha = 1.0f;
  59    float A[] = { 0.188f, 0.856f };
  60    float X[] = { -0.832f, -0.151f };
  61    int incX = -1;
  62    float A_expected[] = { 0.903025f, 0.0f };
  63    cblas_cher(order, uplo, N, alpha, X, incX, A, lda);
  64    {
  65      int i;
  66      for (i = 0; i < 1; i++) {
  67        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher(case 1412) real");
  68        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher(case 1412) imag");
  69      };
  70    };
  71   };
  72
  73
  74   {
  75    int order = 102;
  76    int uplo = 122;
  77    int N = 1;
  78    int lda = 1;
  79    float alpha = 1.0f;
  80    float A[] = { 0.188f, 0.856f };
  81    float X[] = { -0.832f, -0.151f };
  82    int incX = -1;
  83    float A_expected[] = { 0.903025f, 0.0f };
  84    cblas_cher(order, uplo, N, alpha, X, incX, A, lda);
  85    {
  86      int i;
  87      for (i = 0; i < 1; i++) {
  88        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher(case 1413) real");
  89        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher(case 1413) imag");
  90      };
  91    };
  92   };
  93
  94
  95   {
  96    int order = 101;
  97    int uplo = 121;
  98    int N = 1;
  99    int lda = 1;
 100    double alpha = 0.1;
 101    double A[] = { 0.257, 0.326 };
 102    double X[] = { 0.319, -0.009 };
 103    int incX = -1;
 104    double A_expected[] = { 0.2671842, 0.0 };
 105    cblas_zher(order, uplo, N, alpha, X, incX, A, lda);
 106    {
 107      int i;
 108      for (i = 0; i < 1; i++) {
 109        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher(case 1414) real");
 110        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher(case 1414) imag");
 111      };
 112    };
 113   };
 114
 115
 116   {
 117    int order = 101;
 118    int uplo = 122;
 119    int N = 1;
 120    int lda = 1;
 121    double alpha = 0.1;
 122    double A[] = { 0.257, 0.326 };
 123    double X[] = { 0.319, -0.009 };
 124    int incX = -1;
 125    double A_expected[] = { 0.2671842, 0.0 };
 126    cblas_zher(order, uplo, N, alpha, X, incX, A, lda);
 127    {
 128      int i;
 129      for (i = 0; i < 1; i++) {
 130        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher(case 1415) real");
 131        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher(case 1415) imag");
 132      };
 133    };
 134   };
 135
 136
 137   {
 138    int order = 102;
 139    int uplo = 121;
 140    int N = 1;
 141    int lda = 1;
 142    double alpha = 0.1;
 143    double A[] = { 0.257, 0.326 };
 144    double X[] = { 0.319, -0.009 };
 145    int incX = -1;
 146    double A_expected[] = { 0.2671842, 0.0 };
 147    cblas_zher(order, uplo, N, alpha, X, incX, A, lda);
 148    {
 149      int i;
 150      for (i = 0; i < 1; i++) {
 151        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher(case 1416) real");
 152        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher(case 1416) imag");
 153      };
 154    };
 155   };
 156
 157
 158   {
 159    int order = 102;
 160    int uplo = 122;
 161    int N = 1;
 162    int lda = 1;
 163    double alpha = 0.1;
 164    double A[] = { 0.257, 0.326 };
 165    double X[] = { 0.319, -0.009 };
 166    int incX = -1;
 167    double A_expected[] = { 0.2671842, 0.0 };
 168    cblas_zher(order, uplo, N, alpha, X, incX, A, lda);
 169    {
 170      int i;
 171      for (i = 0; i < 1; i++) {
 172        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher(case 1417) real");
 173        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher(case 1417) imag");
 174      };
 175    };
 176   };
 177
 178
 179 }