Blame cblas/test_syr.c

Packit 67cb25
#include <gsl/gsl_test.h>
Packit 67cb25
#include <gsl/gsl_ieee_utils.h>
Packit 67cb25
#include <gsl/gsl_math.h>
Packit 67cb25
#include <gsl/gsl_cblas.h>
Packit 67cb25
Packit 67cb25
#include "tests.h"
Packit 67cb25
Packit 67cb25
void
Packit 67cb25
test_syr (void) {
Packit 67cb25
const double flteps = 1e-4, dbleps = 1e-6;
Packit 67cb25
  {
Packit 67cb25
   int order = 101;
Packit 67cb25
   int uplo = 121;
Packit 67cb25
   int N = 1;
Packit 67cb25
   int lda = 1;
Packit 67cb25
   float alpha = 0.1f;
Packit 67cb25
   float A[] = { -0.291f };
Packit 67cb25
   float X[] = { 0.845f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float A_expected[] = { -0.219597f };
Packit 67cb25
   cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1402)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int order = 101;
Packit 67cb25
   int uplo = 122;
Packit 67cb25
   int N = 1;
Packit 67cb25
   int lda = 1;
Packit 67cb25
   float alpha = 0.1f;
Packit 67cb25
   float A[] = { -0.291f };
Packit 67cb25
   float X[] = { 0.845f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float A_expected[] = { -0.219597f };
Packit 67cb25
   cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1403)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int order = 102;
Packit 67cb25
   int uplo = 121;
Packit 67cb25
   int N = 1;
Packit 67cb25
   int lda = 1;
Packit 67cb25
   float alpha = 0.1f;
Packit 67cb25
   float A[] = { -0.291f };
Packit 67cb25
   float X[] = { 0.845f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float A_expected[] = { -0.219597f };
Packit 67cb25
   cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1404)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int order = 102;
Packit 67cb25
   int uplo = 122;
Packit 67cb25
   int N = 1;
Packit 67cb25
   int lda = 1;
Packit 67cb25
   float alpha = 0.1f;
Packit 67cb25
   float A[] = { -0.291f };
Packit 67cb25
   float X[] = { 0.845f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float A_expected[] = { -0.219597f };
Packit 67cb25
   cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1405)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int order = 101;
Packit 67cb25
   int uplo = 121;
Packit 67cb25
   int N = 1;
Packit 67cb25
   int lda = 1;
Packit 67cb25
   double alpha = -0.3;
Packit 67cb25
   double A[] = { -0.65 };
Packit 67cb25
   double X[] = { -0.891 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double A_expected[] = { -0.8881643 };
Packit 67cb25
   cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1406)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int order = 101;
Packit 67cb25
   int uplo = 122;
Packit 67cb25
   int N = 1;
Packit 67cb25
   int lda = 1;
Packit 67cb25
   double alpha = -0.3;
Packit 67cb25
   double A[] = { -0.65 };
Packit 67cb25
   double X[] = { -0.891 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double A_expected[] = { -0.8881643 };
Packit 67cb25
   cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1407)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int order = 102;
Packit 67cb25
   int uplo = 121;
Packit 67cb25
   int N = 1;
Packit 67cb25
   int lda = 1;
Packit 67cb25
   double alpha = -0.3;
Packit 67cb25
   double A[] = { -0.65 };
Packit 67cb25
   double X[] = { -0.891 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double A_expected[] = { -0.8881643 };
Packit 67cb25
   cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1408)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int order = 102;
Packit 67cb25
   int uplo = 122;
Packit 67cb25
   int N = 1;
Packit 67cb25
   int lda = 1;
Packit 67cb25
   double alpha = -0.3;
Packit 67cb25
   double A[] = { -0.65 };
Packit 67cb25
   double X[] = { -0.891 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double A_expected[] = { -0.8881643 };
Packit 67cb25
   cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1409)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
}