#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_trmv (void) {

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

  {

   int order = 101;

   int trans = 111;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { 0.987f };

   float X[] = { -0.138f };

   int incX = -1;

   float x_expected[] = { -0.136206f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 814)");

     }

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 121;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { 0.987f };

   float X[] = { -0.138f };

   int incX = -1;

   float x_expected[] = { -0.138f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 815)");

     }

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 122;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { 0.987f };

   float X[] = { -0.138f };

   int incX = -1;

   float x_expected[] = { -0.136206f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 816)");

     }

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 122;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { 0.987f };

   float X[] = { -0.138f };

   int incX = -1;

   float x_expected[] = { -0.138f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 817)");

     }

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { 0.987f };

   float X[] = { -0.138f };

   int incX = -1;

   float x_expected[] = { -0.136206f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 818)");

     }

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 121;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { 0.987f };

   float X[] = { -0.138f };

   int incX = -1;

   float x_expected[] = { -0.138f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 819)");

     }

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 122;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { 0.987f };

   float X[] = { -0.138f };

   int incX = -1;

   float x_expected[] = { -0.136206f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 820)");

     }

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 122;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { 0.987f };

   float X[] = { -0.138f };

   int incX = -1;

   float x_expected[] = { -0.138f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 821)");

     }

   };

  };

  {

   int order = 101;

   int trans = 112;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { -0.329f };

   float X[] = { 0.463f };

   int incX = -1;

   float x_expected[] = { -0.152327f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 822)");

     }

   };

  };

  {

   int order = 101;

   int trans = 112;

   int uplo = 121;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { -0.329f };

   float X[] = { 0.463f };

   int incX = -1;

   float x_expected[] = { 0.463f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 823)");

     }

   };

  };

  {

   int order = 101;

   int trans = 112;

   int uplo = 122;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { -0.329f };

   float X[] = { 0.463f };

   int incX = -1;

   float x_expected[] = { -0.152327f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 824)");

     }

   };

  };

  {

   int order = 101;

   int trans = 112;

   int uplo = 122;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { -0.329f };

   float X[] = { 0.463f };

   int incX = -1;

   float x_expected[] = { 0.463f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 825)");

     }

   };

  };

  {

   int order = 102;

   int trans = 112;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { -0.329f };

   float X[] = { 0.463f };

   int incX = -1;

   float x_expected[] = { -0.152327f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 826)");

     }

   };

  };

  {

   int order = 102;

   int trans = 112;

   int uplo = 121;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { -0.329f };

   float X[] = { 0.463f };

   int incX = -1;

   float x_expected[] = { 0.463f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 827)");

     }

   };

  };

  {

   int order = 102;

   int trans = 112;

   int uplo = 122;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { -0.329f };

   float X[] = { 0.463f };

   int incX = -1;

   float x_expected[] = { -0.152327f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 828)");

     }

   };

  };

  {

   int order = 102;

   int trans = 112;

   int uplo = 122;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { -0.329f };

   float X[] = { 0.463f };

   int incX = -1;

   float x_expected[] = { 0.463f };

   cblas_strmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], flteps, "strmv(case 829)");

     }

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   double A[] = { -0.429 };

   double X[] = { -0.899 };

   int incX = -1;

   double x_expected[] = { 0.385671 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 830)");

     }

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 121;

   int diag = 132;

   int N = 1;

   int lda = 1;

   double A[] = { -0.429 };

   double X[] = { -0.899 };

   int incX = -1;

   double x_expected[] = { -0.899 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 831)");

     }

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 122;

   int diag = 131;

   int N = 1;

   int lda = 1;

   double A[] = { -0.429 };

   double X[] = { -0.899 };

   int incX = -1;

   double x_expected[] = { 0.385671 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 832)");

     }

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 122;

   int diag = 132;

   int N = 1;

   int lda = 1;

   double A[] = { -0.429 };

   double X[] = { -0.899 };

   int incX = -1;

   double x_expected[] = { -0.899 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 833)");

     }

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   double A[] = { -0.429 };

   double X[] = { -0.899 };

   int incX = -1;

   double x_expected[] = { 0.385671 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 834)");

     }

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 121;

   int diag = 132;

   int N = 1;

   int lda = 1;

   double A[] = { -0.429 };

   double X[] = { -0.899 };

   int incX = -1;

   double x_expected[] = { -0.899 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 835)");

     }

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 122;

   int diag = 131;

   int N = 1;

   int lda = 1;

   double A[] = { -0.429 };

   double X[] = { -0.899 };

   int incX = -1;

   double x_expected[] = { 0.385671 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 836)");

     }

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 122;

   int diag = 132;

   int N = 1;

   int lda = 1;

   double A[] = { -0.429 };

   double X[] = { -0.899 };

   int incX = -1;

   double x_expected[] = { -0.899 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 837)");

     }

   };

  };

  {

   int order = 101;

   int trans = 112;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   double A[] = { 0.842 };

   double X[] = { 0.192 };

   int incX = -1;

   double x_expected[] = { 0.161664 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 838)");

     }

   };

  };

  {

   int order = 101;

   int trans = 112;

   int uplo = 121;

   int diag = 132;

   int N = 1;

   int lda = 1;

   double A[] = { 0.842 };

   double X[] = { 0.192 };

   int incX = -1;

   double x_expected[] = { 0.192 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 839)");

     }

   };

  };

  {

   int order = 101;

   int trans = 112;

   int uplo = 122;

   int diag = 131;

   int N = 1;

   int lda = 1;

   double A[] = { 0.842 };

   double X[] = { 0.192 };

   int incX = -1;

   double x_expected[] = { 0.161664 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 840)");

     }

   };

  };

  {

   int order = 101;

   int trans = 112;

   int uplo = 122;

   int diag = 132;

   int N = 1;

   int lda = 1;

   double A[] = { 0.842 };

   double X[] = { 0.192 };

   int incX = -1;

   double x_expected[] = { 0.192 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 841)");

     }

   };

  };

  {

   int order = 102;

   int trans = 112;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   double A[] = { 0.842 };

   double X[] = { 0.192 };

   int incX = -1;

   double x_expected[] = { 0.161664 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 842)");

     }

   };

  };

  {

   int order = 102;

   int trans = 112;

   int uplo = 121;

   int diag = 132;

   int N = 1;

   int lda = 1;

   double A[] = { 0.842 };

   double X[] = { 0.192 };

   int incX = -1;

   double x_expected[] = { 0.192 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 843)");

     }

   };

  };

  {

   int order = 102;

   int trans = 112;

   int uplo = 122;

   int diag = 131;

   int N = 1;

   int lda = 1;

   double A[] = { 0.842 };

   double X[] = { 0.192 };

   int incX = -1;

   double x_expected[] = { 0.161664 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 844)");

     }

   };

  };

  {

   int order = 102;

   int trans = 112;

   int uplo = 122;

   int diag = 132;

   int N = 1;

   int lda = 1;

   double A[] = { 0.842 };

   double X[] = { 0.192 };

   int incX = -1;

   double x_expected[] = { 0.192 };

   cblas_dtrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[i], x_expected[i], dbleps, "dtrmv(case 845)");

     }

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { -0.162f, -0.108f };

   float X[] = { 0.542f, 0.461f };

   int incX = -1;

   float x_expected[] = { -0.038016f, -0.133218f };

   cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 846) real");

       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 846) imag");

     };

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 121;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { -0.162f, -0.108f };

   float X[] = { 0.542f, 0.461f };

   int incX = -1;

   float x_expected[] = { 0.542f, 0.461f };

   cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 847) real");

       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 847) imag");

     };

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 122;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { -0.162f, -0.108f };

   float X[] = { 0.542f, 0.461f };

   int incX = -1;

   float x_expected[] = { -0.038016f, -0.133218f };

   cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 848) real");

       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 848) imag");

     };

   };

  };

  {

   int order = 101;

   int trans = 111;

   int uplo = 122;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { -0.162f, -0.108f };

   float X[] = { 0.542f, 0.461f };

   int incX = -1;

   float x_expected[] = { 0.542f, 0.461f };

   cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 849) real");

       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 849) imag");

     };

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { -0.162f, -0.108f };

   float X[] = { 0.542f, 0.461f };

   int incX = -1;

   float x_expected[] = { -0.038016f, -0.133218f };

   cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 850) real");

       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 850) imag");

     };

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 121;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { -0.162f, -0.108f };

   float X[] = { 0.542f, 0.461f };

   int incX = -1;

   float x_expected[] = { 0.542f, 0.461f };

   cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 851) real");

       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 851) imag");

     };

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 122;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { -0.162f, -0.108f };

   float X[] = { 0.542f, 0.461f };

   int incX = -1;

   float x_expected[] = { -0.038016f, -0.133218f };

   cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 852) real");

       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 852) imag");

     };

   };

  };

  {

   int order = 102;

   int trans = 111;

   int uplo = 122;

   int diag = 132;

   int N = 1;

   int lda = 1;

   float A[] = { -0.162f, -0.108f };

   float X[] = { 0.542f, 0.461f };

   int incX = -1;

   float x_expected[] = { 0.542f, 0.461f };

   cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 853) real");

       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 853) imag");

     };

   };

  };

  {

   int order = 101;

   int trans = 112;

   int uplo = 121;

   int diag = 131;

   int N = 1;

   int lda = 1;

   float A[] = { 0.547f, 0.583f };

   float X[] = { -0.302f, 0.434f };

   int incX = -1;

   float x_expected[] = { -0.418216f, 0.061332f };

   cblas_ctrmv(order, uplo, trans, diag, N, A, lda, X, incX);

   {

     int i;

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

       gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctrmv(case 854) real");

       gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctrmv(case 854) imag");

     };

   };

  };