Blame cblas/test_scal.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_scal (void) {
Packit 67cb25
const double flteps = 1e-4, dbleps = 1e-6;
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float alpha = 0.0f;
Packit 67cb25
   float X[] = { 0.651f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.651f };
Packit 67cb25
   cblas_sscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], flteps, "sscal(case 112)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float alpha = 0.1f;
Packit 67cb25
   float X[] = { 0.651f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.651f };
Packit 67cb25
   cblas_sscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], flteps, "sscal(case 113)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float alpha = 1.0f;
Packit 67cb25
   float X[] = { 0.651f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.651f };
Packit 67cb25
   cblas_sscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], flteps, "sscal(case 114)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double alpha = 0;
Packit 67cb25
   double X[] = { 0.686 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.686 };
Packit 67cb25
   cblas_dscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 115)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double alpha = 0.1;
Packit 67cb25
   double X[] = { 0.686 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.686 };
Packit 67cb25
   cblas_dscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 116)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double alpha = 1;
Packit 67cb25
   double X[] = { 0.686 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.686 };
Packit 67cb25
   cblas_dscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 117)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float alpha[2] = {0.0f, 0.0f};
Packit 67cb25
   float X[] = { 0.986f, -0.775f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.986f, -0.775f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 118) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 118) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float alpha[2] = {0.1f, 0.0f};
Packit 67cb25
   float X[] = { 0.986f, -0.775f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.986f, -0.775f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 119) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 119) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float alpha[2] = {1.0f, 0.0f};
Packit 67cb25
   float X[] = { 0.986f, -0.775f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.986f, -0.775f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 120) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 120) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float alpha[2] = {0.0f, 0.1f};
Packit 67cb25
   float X[] = { 0.986f, -0.775f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.986f, -0.775f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 121) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 121) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float alpha[2] = {0.1f, 0.2f};
Packit 67cb25
   float X[] = { 0.986f, -0.775f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.986f, -0.775f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 122) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 122) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   float alpha[2] = {1.0f, 0.3f};
Packit 67cb25
   float X[] = { 0.986f, -0.775f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.986f, -0.775f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 123) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 123) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double alpha[2] = {0, 0};
Packit 67cb25
   double X[] = { 0.454, -0.478 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.454, -0.478 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 124) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 124) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double alpha[2] = {0.1, 0};
Packit 67cb25
   double X[] = { 0.454, -0.478 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.454, -0.478 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 125) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 125) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double alpha[2] = {1, 0};
Packit 67cb25
   double X[] = { 0.454, -0.478 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.454, -0.478 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 126) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 126) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double alpha[2] = {0, 0.1};
Packit 67cb25
   double X[] = { 0.454, -0.478 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.454, -0.478 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 127) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 127) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double alpha[2] = {0.1, 0.2};
Packit 67cb25
   double X[] = { 0.454, -0.478 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.454, -0.478 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 128) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 128) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 1;
Packit 67cb25
   double alpha[2] = {1, 0.3};
Packit 67cb25
   double X[] = { 0.454, -0.478 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.454, -0.478 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 1; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 129) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 129) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha = 0.0f;
Packit 67cb25
   float X[] = { 0.389f, -0.236f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   float expected[] = { 0.0f, -0.0f };
Packit 67cb25
   cblas_sscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], flteps, "sscal(case 130)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha = 0.1f;
Packit 67cb25
   float X[] = { 0.389f, -0.236f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   float expected[] = { 0.0389f, -0.0236f };
Packit 67cb25
   cblas_sscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], flteps, "sscal(case 131)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha = 1.0f;
Packit 67cb25
   float X[] = { 0.389f, -0.236f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   float expected[] = { 0.389f, -0.236f };
Packit 67cb25
   cblas_sscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], flteps, "sscal(case 132)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha = 0;
Packit 67cb25
   double X[] = { -0.429, -0.183 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   double expected[] = { -0.0, -0.0 };
Packit 67cb25
   cblas_dscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 133)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha = 0.1;
Packit 67cb25
   double X[] = { -0.429, -0.183 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   double expected[] = { -0.0429, -0.0183 };
Packit 67cb25
   cblas_dscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 134)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha = 1;
Packit 67cb25
   double X[] = { -0.429, -0.183 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   double expected[] = { -0.429, -0.183 };
Packit 67cb25
   cblas_dscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 135)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {0.0f, 0.0f};
Packit 67cb25
   float X[] = { -0.603f, 0.239f, 0.339f, -0.58f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   float expected[] = { -0.0f, 0.0f, 0.0f, 0.0f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 136) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 136) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {0.1f, 0.0f};
Packit 67cb25
   float X[] = { -0.603f, 0.239f, 0.339f, -0.58f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   float expected[] = { -0.0603f, 0.0239f, 0.0339f, -0.058f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 137) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 137) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {1.0f, 0.0f};
Packit 67cb25
   float X[] = { -0.603f, 0.239f, 0.339f, -0.58f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   float expected[] = { -0.603f, 0.239f, 0.339f, -0.58f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 138) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 138) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {0.0f, 0.1f};
Packit 67cb25
   float X[] = { -0.603f, 0.239f, 0.339f, -0.58f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   float expected[] = { -0.0239f, -0.0603f, 0.058f, 0.0339f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 139) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 139) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {0.1f, 0.2f};
Packit 67cb25
   float X[] = { -0.603f, 0.239f, 0.339f, -0.58f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   float expected[] = { -0.1081f, -0.0967f, 0.1499f, 0.0098f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 140) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 140) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {1.0f, 0.3f};
Packit 67cb25
   float X[] = { -0.603f, 0.239f, 0.339f, -0.58f };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   float expected[] = { -0.6747f, 0.0581f, 0.513f, -0.4783f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 141) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 141) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {0, 0};
Packit 67cb25
   double X[] = { -0.956, 0.613, 0.443, 0.503 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   double expected[] = { -0.0, 0.0, 0.0, 0.0 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 142) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 142) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {0.1, 0};
Packit 67cb25
   double X[] = { -0.956, 0.613, 0.443, 0.503 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   double expected[] = { -0.0956, 0.0613, 0.0443, 0.0503 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 143) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 143) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {1, 0};
Packit 67cb25
   double X[] = { -0.956, 0.613, 0.443, 0.503 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   double expected[] = { -0.956, 0.613, 0.443, 0.503 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 144) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 144) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {0, 0.1};
Packit 67cb25
   double X[] = { -0.956, 0.613, 0.443, 0.503 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   double expected[] = { -0.0613, -0.0956, -0.0503, 0.0443 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 145) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 145) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {0.1, 0.2};
Packit 67cb25
   double X[] = { -0.956, 0.613, 0.443, 0.503 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   double expected[] = { -0.2182, -0.1299, -0.0563, 0.1389 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 146) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 146) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {1, 0.3};
Packit 67cb25
   double X[] = { -0.956, 0.613, 0.443, 0.503 };
Packit 67cb25
   int incX = 1;
Packit 67cb25
   double expected[] = { -1.1399, 0.3262, 0.2921, 0.6359 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 147) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 147) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha = 0.0f;
Packit 67cb25
   float X[] = { 0.629f, -0.419f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.629f, -0.419f };
Packit 67cb25
   cblas_sscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], flteps, "sscal(case 148)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha = 0.1f;
Packit 67cb25
   float X[] = { 0.629f, -0.419f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.629f, -0.419f };
Packit 67cb25
   cblas_sscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], flteps, "sscal(case 149)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha = 1.0f;
Packit 67cb25
   float X[] = { 0.629f, -0.419f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.629f, -0.419f };
Packit 67cb25
   cblas_sscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], flteps, "sscal(case 150)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha = 0;
Packit 67cb25
   double X[] = { 0.398, -0.656 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.398, -0.656 };
Packit 67cb25
   cblas_dscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 151)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha = 0.1;
Packit 67cb25
   double X[] = { 0.398, -0.656 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.398, -0.656 };
Packit 67cb25
   cblas_dscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 152)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha = 1;
Packit 67cb25
   double X[] = { 0.398, -0.656 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.398, -0.656 };
Packit 67cb25
   cblas_dscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[i], expected[i], dbleps, "dscal(case 153)");
Packit 67cb25
     }
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {0.0f, 0.0f};
Packit 67cb25
   float X[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 154) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 154) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {0.1f, 0.0f};
Packit 67cb25
   float X[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 155) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 155) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {1.0f, 0.0f};
Packit 67cb25
   float X[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 156) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 156) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {0.0f, 0.1f};
Packit 67cb25
   float X[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 157) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 157) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {0.1f, 0.2f};
Packit 67cb25
   float X[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 158) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 158) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   float alpha[2] = {1.0f, 0.3f};
Packit 67cb25
   float X[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   float expected[] = { 0.736f, 0.331f, -0.318f, 0.622f };
Packit 67cb25
   cblas_cscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], flteps, "cscal(case 159) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], flteps, "cscal(case 159) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {0, 0};
Packit 67cb25
   double X[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 160) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 160) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {0.1, 0};
Packit 67cb25
   double X[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 161) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 161) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {1, 0};
Packit 67cb25
   double X[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 162) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 162) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {0, 0.1};
Packit 67cb25
   double X[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 163) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 163) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {0.1, 0.2};
Packit 67cb25
   double X[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 164) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 164) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
  {
Packit 67cb25
   int N = 2;
Packit 67cb25
   double alpha[2] = {1, 0.3};
Packit 67cb25
   double X[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   int incX = -1;
Packit 67cb25
   double expected[] = { 0.521, -0.811, 0.556, -0.147 };
Packit 67cb25
   cblas_zscal(N, alpha, X, incX);
Packit 67cb25
   {
Packit 67cb25
     int i;
Packit 67cb25
     for (i = 0; i < 2; i++) {
Packit 67cb25
       gsl_test_rel(X[2*i], expected[2*i], dbleps, "zscal(case 165) real");
Packit 67cb25
       gsl_test_rel(X[2*i+1], expected[2*i+1], dbleps, "zscal(case 165) imag");
Packit 67cb25
     };
Packit 67cb25
   };
Packit 67cb25
  };
Packit 67cb25
Packit 67cb25
Packit 67cb25
}