|
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_gemv (void) {
|
|
Packit |
67cb25 |
const double flteps = 1e-4, dbleps = 1e-6;
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int trans = 111;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = -0.3f;
|
|
Packit |
67cb25 |
float A[] = { -0.805f };
|
|
Packit |
67cb25 |
float X[] = { -0.965f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.537f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.615725f };
|
|
Packit |
67cb25 |
cblas_sgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "sgemv(case 774)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int trans = 111;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = -0.3f;
|
|
Packit |
67cb25 |
float A[] = { -0.805f };
|
|
Packit |
67cb25 |
float X[] = { -0.965f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.537f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.615725f };
|
|
Packit |
67cb25 |
cblas_sgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "sgemv(case 775)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int trans = 112;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = 0.0f;
|
|
Packit |
67cb25 |
float A[] = { -0.805f };
|
|
Packit |
67cb25 |
float X[] = { -0.965f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.537f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.776825f };
|
|
Packit |
67cb25 |
cblas_sgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "sgemv(case 776)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int trans = 112;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = 0.0f;
|
|
Packit |
67cb25 |
float A[] = { -0.805f };
|
|
Packit |
67cb25 |
float X[] = { -0.965f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.537f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.776825f };
|
|
Packit |
67cb25 |
cblas_sgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "sgemv(case 777)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int trans = 111;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
double alpha = -0.3;
|
|
Packit |
67cb25 |
double beta = -1;
|
|
Packit |
67cb25 |
double A[] = { -0.047 };
|
|
Packit |
67cb25 |
double X[] = { 0.672 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.554 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { -0.5445248 };
|
|
Packit |
67cb25 |
cblas_dgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dgemv(case 778)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int trans = 111;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
double alpha = -0.3;
|
|
Packit |
67cb25 |
double beta = -1;
|
|
Packit |
67cb25 |
double A[] = { -0.047 };
|
|
Packit |
67cb25 |
double X[] = { 0.672 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.554 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { -0.5445248 };
|
|
Packit |
67cb25 |
cblas_dgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dgemv(case 779)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int trans = 112;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
double alpha = -1;
|
|
Packit |
67cb25 |
double beta = 1;
|
|
Packit |
67cb25 |
double A[] = { -0.047 };
|
|
Packit |
67cb25 |
double X[] = { 0.672 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.554 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.585584 };
|
|
Packit |
67cb25 |
cblas_dgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dgemv(case 780)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int trans = 112;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
double alpha = -1;
|
|
Packit |
67cb25 |
double beta = 1;
|
|
Packit |
67cb25 |
double A[] = { -0.047 };
|
|
Packit |
67cb25 |
double X[] = { 0.672 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.554 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.585584 };
|
|
Packit |
67cb25 |
cblas_dgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dgemv(case 781)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int trans = 111;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
float alpha[2] = {0.0f, 0.1f};
|
|
Packit |
67cb25 |
float beta[2] = {0.0f, 1.0f};
|
|
Packit |
67cb25 |
float A[] = { 0.629f, 0.801f };
|
|
Packit |
67cb25 |
float X[] = { 0.778f, -0.073f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { -0.976f, -0.682f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.624274f, -0.921216f };
|
|
Packit |
67cb25 |
cblas_cgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "cgemv(case 782) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "cgemv(case 782) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int trans = 111;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
float alpha[2] = {0.0f, 0.1f};
|
|
Packit |
67cb25 |
float beta[2] = {0.0f, 1.0f};
|
|
Packit |
67cb25 |
float A[] = { 0.629f, 0.801f };
|
|
Packit |
67cb25 |
float X[] = { 0.778f, -0.073f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { -0.976f, -0.682f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.624274f, -0.921216f };
|
|
Packit |
67cb25 |
cblas_cgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "cgemv(case 783) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "cgemv(case 783) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int trans = 112;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
float alpha[2] = {0.0f, 1.0f};
|
|
Packit |
67cb25 |
float beta[2] = {-0.3f, 0.1f};
|
|
Packit |
67cb25 |
float A[] = { 0.629f, 0.801f };
|
|
Packit |
67cb25 |
float X[] = { 0.778f, -0.073f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { -0.976f, -0.682f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { -0.216261f, 0.654835f };
|
|
Packit |
67cb25 |
cblas_cgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "cgemv(case 784) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "cgemv(case 784) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int trans = 112;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
float alpha[2] = {0.0f, 1.0f};
|
|
Packit |
67cb25 |
float beta[2] = {-0.3f, 0.1f};
|
|
Packit |
67cb25 |
float A[] = { 0.629f, 0.801f };
|
|
Packit |
67cb25 |
float X[] = { 0.778f, -0.073f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { -0.976f, -0.682f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { -0.216261f, 0.654835f };
|
|
Packit |
67cb25 |
cblas_cgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "cgemv(case 785) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "cgemv(case 785) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int trans = 113;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
float alpha[2] = {0.0f, 0.1f};
|
|
Packit |
67cb25 |
float beta[2] = {-0.3f, 0.1f};
|
|
Packit |
67cb25 |
float A[] = { 0.629f, 0.801f };
|
|
Packit |
67cb25 |
float X[] = { 0.778f, -0.073f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { -0.976f, -0.682f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.427909f, 0.150089f };
|
|
Packit |
67cb25 |
cblas_cgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "cgemv(case 786) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "cgemv(case 786) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int trans = 113;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
float alpha[2] = {0.0f, 0.1f};
|
|
Packit |
67cb25 |
float beta[2] = {-0.3f, 0.1f};
|
|
Packit |
67cb25 |
float A[] = { 0.629f, 0.801f };
|
|
Packit |
67cb25 |
float X[] = { 0.778f, -0.073f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { -0.976f, -0.682f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.427909f, 0.150089f };
|
|
Packit |
67cb25 |
cblas_cgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "cgemv(case 787) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "cgemv(case 787) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int trans = 111;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
double alpha[2] = {0, 0.1};
|
|
Packit |
67cb25 |
double beta[2] = {1, 0};
|
|
Packit |
67cb25 |
double A[] = { 0.932, -0.724 };
|
|
Packit |
67cb25 |
double X[] = { 0.334, -0.317 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.348, 0.07 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.401726, 0.078178 };
|
|
Packit |
67cb25 |
cblas_zgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zgemv(case 788) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zgemv(case 788) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int trans = 111;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
double alpha[2] = {0, 0.1};
|
|
Packit |
67cb25 |
double beta[2] = {1, 0};
|
|
Packit |
67cb25 |
double A[] = { 0.932, -0.724 };
|
|
Packit |
67cb25 |
double X[] = { 0.334, -0.317 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.348, 0.07 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.401726, 0.078178 };
|
|
Packit |
67cb25 |
cblas_zgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zgemv(case 789) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zgemv(case 789) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int trans = 112;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
double alpha[2] = {-0.3, 0.1};
|
|
Packit |
67cb25 |
double beta[2] = {0, 1};
|
|
Packit |
67cb25 |
double A[] = { 0.932, -0.724 };
|
|
Packit |
67cb25 |
double X[] = { 0.334, -0.317 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.348, 0.07 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { -0.040808, 0.517356 };
|
|
Packit |
67cb25 |
cblas_zgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zgemv(case 790) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zgemv(case 790) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int trans = 112;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
double alpha[2] = {-0.3, 0.1};
|
|
Packit |
67cb25 |
double beta[2] = {0, 1};
|
|
Packit |
67cb25 |
double A[] = { 0.932, -0.724 };
|
|
Packit |
67cb25 |
double X[] = { 0.334, -0.317 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.348, 0.07 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { -0.040808, 0.517356 };
|
|
Packit |
67cb25 |
cblas_zgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zgemv(case 791) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zgemv(case 791) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int trans = 113;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
double alpha[2] = {1, 0};
|
|
Packit |
67cb25 |
double beta[2] = {0, 0};
|
|
Packit |
67cb25 |
double A[] = { 0.932, -0.724 };
|
|
Packit |
67cb25 |
double X[] = { 0.334, -0.317 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.348, 0.07 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.540796, -0.053628 };
|
|
Packit |
67cb25 |
cblas_zgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zgemv(case 792) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zgemv(case 792) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int trans = 113;
|
|
Packit |
67cb25 |
int M = 1;
|
|
Packit |
67cb25 |
int N = 1;
|
|
Packit |
67cb25 |
int lda = 1;
|
|
Packit |
67cb25 |
double alpha[2] = {1, 0};
|
|
Packit |
67cb25 |
double beta[2] = {0, 0};
|
|
Packit |
67cb25 |
double A[] = { 0.932, -0.724 };
|
|
Packit |
67cb25 |
double X[] = { 0.334, -0.317 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.348, 0.07 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.540796, -0.053628 };
|
|
Packit |
67cb25 |
cblas_zgemv(order, trans, M, N, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 1; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zgemv(case 793) real");
|
|
Packit |
67cb25 |
gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zgemv(case 793) imag");
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
}
|