1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187
|
#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_hpr2 (void) {
const double flteps = 1e-4, dbleps = 1e-6;
{
int order = 101;
int uplo = 121;
int N = 1;
float alpha[2] = {-1.0f, 0.0f};
float Ap[] = { 0.159f, -0.13f };
float X[] = { 0.854f, 0.851f };
int incX = -1;
float Y[] = { 0.526f, -0.267f };
int incY = -1;
float Ap_expected[] = { -0.284974f, 0.0f };
cblas_chpr2(order, uplo, N, alpha, X, incX, Y, incY, Ap);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Ap[2*i], Ap_expected[2*i], flteps, "chpr2(case 1458) real");
gsl_test_rel(Ap[2*i+1], Ap_expected[2*i+1], flteps, "chpr2(case 1458) imag");
};
};
};
{
int order = 101;
int uplo = 122;
int N = 1;
float alpha[2] = {-1.0f, 0.0f};
float Ap[] = { 0.159f, -0.13f };
float X[] = { 0.854f, 0.851f };
int incX = -1;
float Y[] = { 0.526f, -0.267f };
int incY = -1;
float Ap_expected[] = { -0.284974f, 0.0f };
cblas_chpr2(order, uplo, N, alpha, X, incX, Y, incY, Ap);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Ap[2*i], Ap_expected[2*i], flteps, "chpr2(case 1459) real");
gsl_test_rel(Ap[2*i+1], Ap_expected[2*i+1], flteps, "chpr2(case 1459) imag");
};
};
};
{
int order = 102;
int uplo = 121;
int N = 1;
float alpha[2] = {-1.0f, 0.0f};
float Ap[] = { 0.159f, -0.13f };
float X[] = { 0.854f, 0.851f };
int incX = -1;
float Y[] = { 0.526f, -0.267f };
int incY = -1;
float Ap_expected[] = { -0.284974f, 0.0f };
cblas_chpr2(order, uplo, N, alpha, X, incX, Y, incY, Ap);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Ap[2*i], Ap_expected[2*i], flteps, "chpr2(case 1460) real");
gsl_test_rel(Ap[2*i+1], Ap_expected[2*i+1], flteps, "chpr2(case 1460) imag");
};
};
};
{
int order = 102;
int uplo = 122;
int N = 1;
float alpha[2] = {-1.0f, 0.0f};
float Ap[] = { 0.159f, -0.13f };
float X[] = { 0.854f, 0.851f };
int incX = -1;
float Y[] = { 0.526f, -0.267f };
int incY = -1;
float Ap_expected[] = { -0.284974f, 0.0f };
cblas_chpr2(order, uplo, N, alpha, X, incX, Y, incY, Ap);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Ap[2*i], Ap_expected[2*i], flteps, "chpr2(case 1461) real");
gsl_test_rel(Ap[2*i+1], Ap_expected[2*i+1], flteps, "chpr2(case 1461) imag");
};
};
};
{
int order = 101;
int uplo = 121;
int N = 1;
double alpha[2] = {-0.3, 0.1};
double Ap[] = { 0.772, 0.997 };
double X[] = { -0.173, -0.839 };
int incX = -1;
double Y[] = { 0.941, -0.422 };
int incY = -1;
double Ap_expected[] = { 0.829742, 0.0 };
cblas_zhpr2(order, uplo, N, alpha, X, incX, Y, incY, Ap);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Ap[2*i], Ap_expected[2*i], dbleps, "zhpr2(case 1462) real");
gsl_test_rel(Ap[2*i+1], Ap_expected[2*i+1], dbleps, "zhpr2(case 1462) imag");
};
};
};
{
int order = 101;
int uplo = 122;
int N = 1;
double alpha[2] = {-0.3, 0.1};
double Ap[] = { 0.772, 0.997 };
double X[] = { -0.173, -0.839 };
int incX = -1;
double Y[] = { 0.941, -0.422 };
int incY = -1;
double Ap_expected[] = { 0.829742, 0.0 };
cblas_zhpr2(order, uplo, N, alpha, X, incX, Y, incY, Ap);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Ap[2*i], Ap_expected[2*i], dbleps, "zhpr2(case 1463) real");
gsl_test_rel(Ap[2*i+1], Ap_expected[2*i+1], dbleps, "zhpr2(case 1463) imag");
};
};
};
{
int order = 102;
int uplo = 121;
int N = 1;
double alpha[2] = {-0.3, 0.1};
double Ap[] = { 0.772, 0.997 };
double X[] = { -0.173, -0.839 };
int incX = -1;
double Y[] = { 0.941, -0.422 };
int incY = -1;
double Ap_expected[] = { 0.829742, 0.0 };
cblas_zhpr2(order, uplo, N, alpha, X, incX, Y, incY, Ap);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Ap[2*i], Ap_expected[2*i], dbleps, "zhpr2(case 1464) real");
gsl_test_rel(Ap[2*i+1], Ap_expected[2*i+1], dbleps, "zhpr2(case 1464) imag");
};
};
};
{
int order = 102;
int uplo = 122;
int N = 1;
double alpha[2] = {-0.3, 0.1};
double Ap[] = { 0.772, 0.997 };
double X[] = { -0.173, -0.839 };
int incX = -1;
double Y[] = { 0.941, -0.422 };
int incY = -1;
double Ap_expected[] = { 0.829742, 0.0 };
cblas_zhpr2(order, uplo, N, alpha, X, incX, Y, incY, Ap);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Ap[2*i], Ap_expected[2*i], dbleps, "zhpr2(case 1465) real");
gsl_test_rel(Ap[2*i+1], Ap_expected[2*i+1], dbleps, "zhpr2(case 1465) imag");
};
};
};
}
|