#include <gsl/gsl_test.h>

#include <gsl/gsl_ieee_utils.h>

#include <gsl/gsl_math.h>

#include <gsl/gsl_cblas.h>

#include "tests.h"

void

test_spmv (void) {

const double flteps = 1e-4, dbleps = 1e-6;

  {

   int order = 101;

   int uplo = 121;

   float alpha = 0.1f;

   float beta = -0.3f;

   int N = 2;

   float A[] = { -0.174f, 0.878f, 0.478f };

   float X[] = { 0.503f, 0.313f };

   int incX = -1;

   float Y[] = { -0.565f, -0.109f };

   int incY = -1;

   float y_expected[] = { 0.221025f, 0.0714172f };

   cblas_sspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], flteps, "sspmv(case 1134)");

     }

   };

  };

  {

   int order = 101;

   int uplo = 121;

   float alpha = 0.1f;

   float beta = -0.3f;

   int N = 2;

   float A[] = { -0.174f, 0.878f, 0.478f };

   float X[] = { 0.503f, 0.313f };

   int incX = -1;

   float Y[] = { -0.565f, -0.109f };

   int incY = -1;

   float y_expected[] = { 0.221025f, 0.0714172f };

   cblas_sspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], flteps, "sspmv(case 1135)");

     }

   };

  };

  {

   int order = 101;

   int uplo = 122;

   float alpha = 0.1f;

   float beta = -0.3f;

   int N = 2;

   float A[] = { -0.174f, 0.878f, 0.478f };

   float X[] = { 0.503f, 0.313f };

   int incX = -1;

   float Y[] = { -0.565f, -0.109f };

   int incY = -1;

   float y_expected[] = { 0.221025f, 0.0714172f };

   cblas_sspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], flteps, "sspmv(case 1136)");

     }

   };

  };

  {

   int order = 101;

   int uplo = 122;

   float alpha = 0.1f;

   float beta = -0.3f;

   int N = 2;

   float A[] = { -0.174f, 0.878f, 0.478f };

   float X[] = { 0.503f, 0.313f };

   int incX = -1;

   float Y[] = { -0.565f, -0.109f };

   int incY = -1;

   float y_expected[] = { 0.221025f, 0.0714172f };

   cblas_sspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], flteps, "sspmv(case 1137)");

     }

   };

  };

  {

   int order = 102;

   int uplo = 121;

   float alpha = 0.1f;

   float beta = -0.3f;

   int N = 2;

   float A[] = { -0.174f, 0.878f, 0.478f };

   float X[] = { 0.503f, 0.313f };

   int incX = -1;

   float Y[] = { -0.565f, -0.109f };

   int incY = -1;

   float y_expected[] = { 0.221025f, 0.0714172f };

   cblas_sspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], flteps, "sspmv(case 1138)");

     }

   };

  };

  {

   int order = 102;

   int uplo = 121;

   float alpha = 0.1f;

   float beta = -0.3f;

   int N = 2;

   float A[] = { -0.174f, 0.878f, 0.478f };

   float X[] = { 0.503f, 0.313f };

   int incX = -1;

   float Y[] = { -0.565f, -0.109f };

   int incY = -1;

   float y_expected[] = { 0.221025f, 0.0714172f };

   cblas_sspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], flteps, "sspmv(case 1139)");

     }

   };

  };

  {

   int order = 102;

   int uplo = 122;

   float alpha = 0.1f;

   float beta = -0.3f;

   int N = 2;

   float A[] = { -0.174f, 0.878f, 0.478f };

   float X[] = { 0.503f, 0.313f };

   int incX = -1;

   float Y[] = { -0.565f, -0.109f };

   int incY = -1;

   float y_expected[] = { 0.221025f, 0.0714172f };

   cblas_sspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], flteps, "sspmv(case 1140)");

     }

   };

  };

  {

   int order = 102;

   int uplo = 122;

   float alpha = 0.1f;

   float beta = -0.3f;

   int N = 2;

   float A[] = { -0.174f, 0.878f, 0.478f };

   float X[] = { 0.503f, 0.313f };

   int incX = -1;

   float Y[] = { -0.565f, -0.109f };

   int incY = -1;

   float y_expected[] = { 0.221025f, 0.0714172f };

   cblas_sspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], flteps, "sspmv(case 1141)");

     }

   };

  };

  {

   int order = 101;

   int uplo = 121;

   double alpha = -1;

   double beta = 0.1;

   int N = 2;

   double A[] = { -0.181, -0.071, -0.038 };

   double X[] = { -0.015, 0.132 };

   int incX = -1;

   double Y[] = { -0.449, -0.219 };

   int incY = -1;

   double y_expected[] = { -0.036098, 9.27e-04 };

   cblas_dspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], dbleps, "dspmv(case 1142)");

     }

   };

  };

  {

   int order = 101;

   int uplo = 121;

   double alpha = -1;

   double beta = 0.1;

   int N = 2;

   double A[] = { -0.181, -0.071, -0.038 };

   double X[] = { -0.015, 0.132 };

   int incX = -1;

   double Y[] = { -0.449, -0.219 };

   int incY = -1;

   double y_expected[] = { -0.036098, 9.27e-04 };

   cblas_dspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], dbleps, "dspmv(case 1143)");

     }

   };

  };

  {

   int order = 101;

   int uplo = 122;

   double alpha = -1;

   double beta = 0.1;

   int N = 2;

   double A[] = { -0.181, -0.071, -0.038 };

   double X[] = { -0.015, 0.132 };

   int incX = -1;

   double Y[] = { -0.449, -0.219 };

   int incY = -1;

   double y_expected[] = { -0.036098, 9.27e-04 };

   cblas_dspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], dbleps, "dspmv(case 1144)");

     }

   };

  };

  {

   int order = 101;

   int uplo = 122;

   double alpha = -1;

   double beta = 0.1;

   int N = 2;

   double A[] = { -0.181, -0.071, -0.038 };

   double X[] = { -0.015, 0.132 };

   int incX = -1;

   double Y[] = { -0.449, -0.219 };

   int incY = -1;

   double y_expected[] = { -0.036098, 9.27e-04 };

   cblas_dspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], dbleps, "dspmv(case 1145)");

     }

   };

  };

  {

   int order = 102;

   int uplo = 121;

   double alpha = -1;

   double beta = 0.1;

   int N = 2;

   double A[] = { -0.181, -0.071, -0.038 };

   double X[] = { -0.015, 0.132 };

   int incX = -1;

   double Y[] = { -0.449, -0.219 };

   int incY = -1;

   double y_expected[] = { -0.036098, 9.27e-04 };

   cblas_dspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], dbleps, "dspmv(case 1146)");

     }

   };

  };

  {

   int order = 102;

   int uplo = 121;

   double alpha = -1;

   double beta = 0.1;

   int N = 2;

   double A[] = { -0.181, -0.071, -0.038 };

   double X[] = { -0.015, 0.132 };

   int incX = -1;

   double Y[] = { -0.449, -0.219 };

   int incY = -1;

   double y_expected[] = { -0.036098, 9.27e-04 };

   cblas_dspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], dbleps, "dspmv(case 1147)");

     }

   };

  };

  {

   int order = 102;

   int uplo = 122;

   double alpha = -1;

   double beta = 0.1;

   int N = 2;

   double A[] = { -0.181, -0.071, -0.038 };

   double X[] = { -0.015, 0.132 };

   int incX = -1;

   double Y[] = { -0.449, -0.219 };

   int incY = -1;

   double y_expected[] = { -0.036098, 9.27e-04 };

   cblas_dspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], dbleps, "dspmv(case 1148)");

     }

   };

  };

  {

   int order = 102;

   int uplo = 122;

   double alpha = -1;

   double beta = 0.1;

   int N = 2;

   double A[] = { -0.181, -0.071, -0.038 };

   double X[] = { -0.015, 0.132 };

   int incX = -1;

   double Y[] = { -0.449, -0.219 };

   int incY = -1;

   double y_expected[] = { -0.036098, 9.27e-04 };

   cblas_dspmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);

   {

     int i;

     for (i = 0; i < 2; i++) {

       gsl_test_rel(Y[i], y_expected[i], dbleps, "dspmv(case 1149)");

     }

   };

  };

}