|
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_sbmv (void) {
|
|
Packit |
67cb25 |
const double flteps = 1e-4, dbleps = 1e-6;
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int uplo = 121;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = 0.0f;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
float A[] = { 0.627f, -0.312f, 0.031f, 0.308f, 0.323f, -0.578f, 0.797f, 0.545f, -0.476f };
|
|
Packit |
67cb25 |
float X[] = { -0.542f, 0.606f, 0.727f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.755f, 0.268f, -0.99f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { -0.236236f, -0.215242f, 0.266757f };
|
|
Packit |
67cb25 |
cblas_ssbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "ssbmv(case 1102)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int uplo = 121;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = 0.0f;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
float A[] = { 0.627f, -0.312f, 0.031f, 0.308f, 0.323f, -0.578f, 0.797f, 0.545f, -0.476f };
|
|
Packit |
67cb25 |
float X[] = { -0.542f, 0.606f, 0.727f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.755f, 0.268f, -0.99f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { -0.236236f, -0.215242f, 0.266757f };
|
|
Packit |
67cb25 |
cblas_ssbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "ssbmv(case 1103)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int uplo = 122;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = 0.0f;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
float A[] = { 0.627f, -0.312f, 0.031f, 0.308f, 0.323f, -0.578f, 0.797f, 0.545f, -0.476f };
|
|
Packit |
67cb25 |
float X[] = { -0.542f, 0.606f, 0.727f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.755f, 0.268f, -0.99f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.187592f, -0.01232f, -0.040176f };
|
|
Packit |
67cb25 |
cblas_ssbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "ssbmv(case 1104)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int uplo = 122;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = 0.0f;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
float A[] = { 0.627f, -0.312f, 0.031f, 0.308f, 0.323f, -0.578f, 0.797f, 0.545f, -0.476f };
|
|
Packit |
67cb25 |
float X[] = { -0.542f, 0.606f, 0.727f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.755f, 0.268f, -0.99f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.187592f, -0.01232f, -0.040176f };
|
|
Packit |
67cb25 |
cblas_ssbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "ssbmv(case 1105)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int uplo = 121;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = 0.0f;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
float A[] = { 0.627f, -0.312f, 0.031f, 0.308f, 0.323f, -0.578f, 0.797f, 0.545f, -0.476f };
|
|
Packit |
67cb25 |
float X[] = { -0.542f, 0.606f, 0.727f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.755f, 0.268f, -0.99f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.187592f, -0.01232f, -0.040176f };
|
|
Packit |
67cb25 |
cblas_ssbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "ssbmv(case 1106)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int uplo = 121;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = 0.0f;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
float A[] = { 0.627f, -0.312f, 0.031f, 0.308f, 0.323f, -0.578f, 0.797f, 0.545f, -0.476f };
|
|
Packit |
67cb25 |
float X[] = { -0.542f, 0.606f, 0.727f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.755f, 0.268f, -0.99f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { 0.187592f, -0.01232f, -0.040176f };
|
|
Packit |
67cb25 |
cblas_ssbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "ssbmv(case 1107)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int uplo = 122;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = 0.0f;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
float A[] = { 0.627f, -0.312f, 0.031f, 0.308f, 0.323f, -0.578f, 0.797f, 0.545f, -0.476f };
|
|
Packit |
67cb25 |
float X[] = { -0.542f, 0.606f, 0.727f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.755f, 0.268f, -0.99f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { -0.236236f, -0.215242f, 0.266757f };
|
|
Packit |
67cb25 |
cblas_ssbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "ssbmv(case 1108)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int uplo = 122;
|
|
Packit |
67cb25 |
float alpha = 1.0f;
|
|
Packit |
67cb25 |
float beta = 0.0f;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
float A[] = { 0.627f, -0.312f, 0.031f, 0.308f, 0.323f, -0.578f, 0.797f, 0.545f, -0.476f };
|
|
Packit |
67cb25 |
float X[] = { -0.542f, 0.606f, 0.727f };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
float Y[] = { 0.755f, 0.268f, -0.99f };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
float y_expected[] = { -0.236236f, -0.215242f, 0.266757f };
|
|
Packit |
67cb25 |
cblas_ssbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], flteps, "ssbmv(case 1109)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int uplo = 121;
|
|
Packit |
67cb25 |
double alpha = 0;
|
|
Packit |
67cb25 |
double beta = 1;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
double A[] = { 0.83, -0.568, -0.888, 0.281, -0.779, -0.148, 0.138, 0.053, -0.757 };
|
|
Packit |
67cb25 |
double X[] = { 0.166, 0.808, 0.723 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
cblas_dsbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dsbmv(case 1110)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int uplo = 121;
|
|
Packit |
67cb25 |
double alpha = 0;
|
|
Packit |
67cb25 |
double beta = 1;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
double A[] = { 0.83, -0.568, -0.888, 0.281, -0.779, -0.148, 0.138, 0.053, -0.757 };
|
|
Packit |
67cb25 |
double X[] = { 0.166, 0.808, 0.723 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
cblas_dsbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dsbmv(case 1111)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int uplo = 122;
|
|
Packit |
67cb25 |
double alpha = 0;
|
|
Packit |
67cb25 |
double beta = 1;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
double A[] = { 0.83, -0.568, -0.888, 0.281, -0.779, -0.148, 0.138, 0.053, -0.757 };
|
|
Packit |
67cb25 |
double X[] = { 0.166, 0.808, 0.723 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
cblas_dsbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dsbmv(case 1112)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 101;
|
|
Packit |
67cb25 |
int uplo = 122;
|
|
Packit |
67cb25 |
double alpha = 0;
|
|
Packit |
67cb25 |
double beta = 1;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
double A[] = { 0.83, -0.568, -0.888, 0.281, -0.779, -0.148, 0.138, 0.053, -0.757 };
|
|
Packit |
67cb25 |
double X[] = { 0.166, 0.808, 0.723 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
cblas_dsbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dsbmv(case 1113)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int uplo = 121;
|
|
Packit |
67cb25 |
double alpha = 0;
|
|
Packit |
67cb25 |
double beta = 1;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
double A[] = { 0.83, -0.568, -0.888, 0.281, -0.779, -0.148, 0.138, 0.053, -0.757 };
|
|
Packit |
67cb25 |
double X[] = { 0.166, 0.808, 0.723 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
cblas_dsbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dsbmv(case 1114)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int uplo = 121;
|
|
Packit |
67cb25 |
double alpha = 0;
|
|
Packit |
67cb25 |
double beta = 1;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
double A[] = { 0.83, -0.568, -0.888, 0.281, -0.779, -0.148, 0.138, 0.053, -0.757 };
|
|
Packit |
67cb25 |
double X[] = { 0.166, 0.808, 0.723 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
cblas_dsbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dsbmv(case 1115)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int uplo = 122;
|
|
Packit |
67cb25 |
double alpha = 0;
|
|
Packit |
67cb25 |
double beta = 1;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
double A[] = { 0.83, -0.568, -0.888, 0.281, -0.779, -0.148, 0.138, 0.053, -0.757 };
|
|
Packit |
67cb25 |
double X[] = { 0.166, 0.808, 0.723 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
cblas_dsbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dsbmv(case 1116)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int order = 102;
|
|
Packit |
67cb25 |
int uplo = 122;
|
|
Packit |
67cb25 |
double alpha = 0;
|
|
Packit |
67cb25 |
double beta = 1;
|
|
Packit |
67cb25 |
int N = 3;
|
|
Packit |
67cb25 |
int k = 1;
|
|
Packit |
67cb25 |
int lda = 3;
|
|
Packit |
67cb25 |
double A[] = { 0.83, -0.568, -0.888, 0.281, -0.779, -0.148, 0.138, 0.053, -0.757 };
|
|
Packit |
67cb25 |
double X[] = { 0.166, 0.808, 0.723 };
|
|
Packit |
67cb25 |
int incX = -1;
|
|
Packit |
67cb25 |
double Y[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
int incY = -1;
|
|
Packit |
67cb25 |
double y_expected[] = { 0.9, 0.99, -0.578 };
|
|
Packit |
67cb25 |
cblas_dsbmv(order, uplo, N, k, alpha, A, lda, X, incX, beta, Y, incY);
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int i;
|
|
Packit |
67cb25 |
for (i = 0; i < 3; i++) {
|
|
Packit |
67cb25 |
gsl_test_rel(Y[i], y_expected[i], dbleps, "dsbmv(case 1117)");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
}
|