Blame cblas/test_amax.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_amax (void) {
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float X[] = { -0.388f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   int expected = 0;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_isamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "samax(case 52)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double X[] = { 0.247 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   int expected = 0;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_idamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "damax(case 53)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float X[] = { 0.704f, 0.665f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   int expected = 0;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_icamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "camax(case 54)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double X[] = { -0.599, -0.758 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   int expected = 0;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_izamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "zamax(case 55)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float X[] = { 0.909f, 0.037f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   int expected = 0;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_isamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "samax(case 56)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double X[] = { 0.271, -0.426 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   int expected = 1;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_idamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "damax(case 57)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float X[] = { -0.648f, 0.317f, 0.62f, 0.392f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   int expected = 1;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_icamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "camax(case 58)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double X[] = { -0.789, 0.352, 0.562, 0.697 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   int expected = 1;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_izamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "zamax(case 59)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float X[] = { 0.487f, 0.918f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   int expected = 0;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_isamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "samax(case 60)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double X[] = { 0.537, 0.826 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   int expected = 0;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_idamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "damax(case 61)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float X[] = { 0.993f, 0.172f, -0.825f, 0.873f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   int expected = 0;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_icamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "camax(case 62)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double X[] = { 0.913, -0.436, -0.134, 0.129 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   int expected = 0;
Packit 67cb25
   int k;
Packit 67cb25
   k = cblas_izamax(N, X, incX);
Packit 67cb25
   gsl_test_int(k, expected, "zamax(case 63)");
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
}