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