JNIEXPORT void JNICALL Java_uncomplicate_neanderthal_CBLAS_srotm
(JNIEnv *env, jclass clazz, jint N,
jobject X, jint offsetX, jint incX,
jobject Y, jint offsetY, jint incY,
jobject param) {
float *cX = (float *) (*env)->GetDirectBufferAddress(env, X);
float *cY = (float *) (*env)->GetDirectBufferAddress(env, Y);
float *c_param = (float *) (*env)->GetDirectBufferAddress(env, param);
cblas_srotm(N, cX + offsetX, incX, cY + offsetY, incY, c_param);
};
JNIEXPORT void JNICALL Java_MKL_affineTransform (JNIEnv *env, jclass cl, jint n, jfloatArray jxs, jfloatArray jys, jfloat a, jfloat b, jfloat c, jfloat d, jfloat e, jfloat f) {
jfloat *xs = (*env)->GetFloatArrayElements(env, jxs, NULL);
jfloat *ys = (*env)->GetFloatArrayElements(env, jys, NULL);
float p[] = {-1, a, b, c, d};
cblas_srotm(n, xs, 1, ys, 1, p);
cblas_saxpy(n, 1, &e, 0, xs, 1);
cblas_saxpy(n, 1, &f, 0, ys, 1);
(*env)->ReleaseFloatArrayElements(env, jxs, xs, 0);
(*env)->ReleaseFloatArrayElements(env, jys, ys, 0);
}
VALUE rb_blas_xrotm_mod(int argc, VALUE *argv, VALUE self)
{
Matrix *dx, *dy;
int incx;
int incy;
int n;
sRotmg *srotmg;
dRotmg *drotmg;
//char error_msg[64];
VALUE vector_y, rotmg, n_value, incx_value, incy_value;
rb_scan_args(argc, argv, "23", &vector_y, &rotmg, &incx_value, &incy_value, &n_value);
Data_Get_Struct(self, Matrix, dx);
Data_Get_Struct(vector_y, Matrix, dy);
if(incx_value == Qnil)
incx = 1;
else
incx = NUM2INT(incx_value);
if(incy_value == Qnil)
incy = 1;
else
incy = NUM2INT(incy_value);
if(n_value == Qnil)
n = dx->nrows;
else
n = NUM2INT(n_value);
if(dx == NULL || dx->ncols != 1)
{ //sprintf(error_msg, "Self is not a Vector");
rb_raise(rb_eRuntimeError, "Self is not a Vector");
}
if(dy == NULL || dy->ncols != 1)
{ //sprintf(error_msg, "Argument dy is not a Vector");
rb_raise(rb_eRuntimeError, "Argument dy is not a Vector");
}
if(dy->data_type != dx->data_type)
{ //sprintf(error_msg, "Vectors are different data_types");
rb_raise(rb_eRuntimeError, "Vectors are different data_types");
}
if(rotmg == Qnil)
{ //sprintf(error_msg, "[SD]rotg argument is nil?");
rb_raise(rb_eRuntimeError, "[SD]rotg argument is nil?");
}
switch(dx->data_type)
{
case Single_t: //s
Data_Get_Struct(rotmg, sRotmg, srotmg);
cblas_srotm(n , (float *)dx->data, incx, (float *)dy->data, incy, srotmg->param );
break;
case Double_t: //d
Data_Get_Struct(rotmg, dRotmg, drotmg);
cblas_drotm(n , (double *)dx->data, incx, (double *)dy->data, incy, drotmg->param );
break;
default:
//sprintf(error_msg, "Invalid data_type (%d) in Matrix", dx->data_type);
rb_raise(rb_eRuntimeError, "Invalid data_type (%d) in Matrix", dx->data_type);
break; //Never reaches here.
}
return self;
}
JNIEXPORT jobject JNICALL Java_MKL_apply(JNIEnv *env, jclass cl, jint n, jobjectArray jts) {
int m = (*env)->GetArrayLength(env, jts);
float **p = (float **)malloc(m * sizeof (float *));
jclass classTransform = (*env)->FindClass(env, "Transform");
jfieldID idA = (*env)->GetFieldID(env, classTransform, "a", "F");
jfieldID idB = (*env)->GetFieldID(env, classTransform, "b", "F");
jfieldID idC = (*env)->GetFieldID(env, classTransform, "c", "F");
jfieldID idD = (*env)->GetFieldID(env, classTransform, "d", "F");
jfieldID idE = (*env)->GetFieldID(env, classTransform, "e", "F");
jfieldID idF = (*env)->GetFieldID(env, classTransform, "f", "F");
int i;
jobject t;
for (i = 0; i < m; ++i) {
t = (*env)->GetObjectArrayElement(env, jts, i);
p[i] = (float *)malloc(7 * sizeof (float));
p[i][0] = -1;
p[i][1] = (*env)->GetFloatField(env, t, idA);
p[i][2] = (*env)->GetFloatField(env, t, idB);
p[i][3] = (*env)->GetFloatField(env, t, idC);
p[i][4] = (*env)->GetFloatField(env, t, idD);
p[i][5] = (*env)->GetFloatField(env, t, idE);
p[i][6] = (*env)->GetFloatField(env, t, idF);
}
int k = ipow(m, n);
jfloatArray jxs = (*env)->NewFloatArray(env, k);
jfloatArray jys = (*env)->NewFloatArray(env, k);
jfloat *xs = (*env)->GetFloatArrayElements(env, jxs, NULL);
jfloat *ys = (*env)->GetFloatArrayElements(env, jys, NULL);
xs[0] = 0;
ys[0] = 0;
int j;
int w;
for (i = 0; i < n; ++i) {
w = ipow(m, i);
for (j = 1; j < m; ++j) {
cblas_scopy(w, xs, 1, xs + w * j, 1);
cblas_scopy(w, ys, 1, ys + w * j, 1);
}
for (j = 0; j < m; ++j) {
cblas_srotm(w, xs + w * j, 1, ys + w * j, 1, p[j]);
cblas_saxpy(w, 1, &(p[j][5]), 0, xs + w * j, 1);
cblas_saxpy(w, 1, &(p[j][6]), 0, ys + w * j, 1);
}
}
for (i = 0; i < m; ++i) {
free(p[i]);
}
free(p);
(*env)->ReleaseFloatArrayElements(env, jxs, xs, 0);
(*env)->ReleaseFloatArrayElements(env, jys, ys, 0);
jclass classPoints = (*env)->FindClass(env, "Points");
jmethodID idConstructorPoints = (*env)->GetMethodID(env, classPoints, "<init>", "([F[F)V");
jobject ps = (*env)->NewObject(env, classPoints, idConstructorPoints, jxs, jys);
return ps;
}
void
test_rotm (void) {
const double flteps = 1e-4, dbleps = 1e-6;
{
int N = 1;
float h[] = { -1.0f, -4.44982e+03f, -15.5826f, 7.091334e+04f, 2.95912e+04f };
float X[] = { -0.034f };
int incX = 1;
float Y[] = { -0.56f };
int incY = -1;
float x_expected[] = { -3.956017e+04f };
float y_expected[] = { -1.657054e+04f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 654)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 655)");
}
};
};
{
int N = 1;
float h[] = { 0.0f, 15.9728f, 6.400638e+03f, 1.733082e-05f, 1.524511e-04f };
float X[] = { -0.034f };
int incX = 1;
float Y[] = { -0.56f };
int incY = -1;
float x_expected[] = { -0.0340097f };
float y_expected[] = { -218.182f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 656)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 657)");
}
};
};
{
int N = 1;
float h[] = { 1.0f, 5.688411e+04f, 5.914789e+03f, 0.00210473f, 0.0231019f };
float X[] = { -0.034f };
int incX = 1;
float Y[] = { -0.56f };
int incY = -1;
float x_expected[] = { -1.93462e+03f };
float y_expected[] = { 0.0210629f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 658)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 659)");
}
};
};
{
int N = 1;
float h[] = { -2.0f, -0.582083f, 0.00103161f, -3.429851e-05f, 7.411469e-05f };
float X[] = { -0.034f };
int incX = 1;
float Y[] = { -0.56f };
int incY = -1;
float x_expected[] = { -0.034f };
float y_expected[] = { -0.56f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 660)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 661)");
}
};
};
{
int N = 1;
float h[] = { -1.0f, 115.163f, -6.715448e+04f, -258.695f, -16.2552f };
float X[] = { -0.034f };
int incX = 1;
float Y[] = { -0.56f };
int incY = -1;
float x_expected[] = { 140.954f };
float y_expected[] = { 2.292355e+03f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 662)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 663)");
}
};
};
{
int N = 1;
float h[] = { 0.0f, -3.314862e+03f, -442.976f, -214.586f, -25.9716f };
float X[] = { -0.034f };
int incX = 1;
float Y[] = { -0.56f };
int incY = -1;
float x_expected[] = { 120.134f };
float y_expected[] = { 14.5012f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 664)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 665)");
}
};
};
{
int N = 1;
float h[] = { 1.0f, -1.177304e+03f, -1.236662e-04f, -0.186585f, 1.15841f };
float X[] = { -0.034f };
int incX = 1;
float Y[] = { -0.56f };
int incY = -1;
float x_expected[] = { 39.4683f };
float y_expected[] = { -0.614711f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 666)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 667)");
}
};
};
{
int N = 1;
float h[] = { -2.0f, -88.9796f, 0.808226f, 1.106582e-05f, -0.00862288f };
float X[] = { -0.034f };
int incX = 1;
float Y[] = { -0.56f };
int incY = -1;
float x_expected[] = { -0.034f };
float y_expected[] = { -0.56f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 668)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 669)");
}
};
};
{
int N = 1;
float h[] = { -1.0f, -0.00225865f, 8.338551e+04f, -1.98282f, -2.409905e-05f };
float X[] = { -0.034f };
int incX = 1;
float Y[] = { -0.56f };
int incY = -1;
float x_expected[] = { 1.11046f };
float y_expected[] = { -2.835107e+03f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 670)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 671)");
}
};
};
{
int N = 1;
float h[] = { 0.0f, 0.258779f, 74.2802f, 0.923299f, 4.847128e+03f };
float X[] = { -0.034f };
int incX = 1;
float Y[] = { -0.56f };
int incY = -1;
float x_expected[] = { -0.551048f };
float y_expected[] = { -3.08553f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 672)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 673)");
}
};
};
{
int N = 1;
double h[] = { -1.0, -8.00850735044, 0.0204647351647, 1.898461360078e-04, -4.32701487194 };
double X[] = { 0.84 };
int incX = 1;
double Y[] = { -0.711 };
int incY = -1;
double x_expected[] = { -6.72728115497 };
double y_expected[] = { 3.09369795149 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 674)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 675)");
}
};
};
{
int N = 1;
double h[] = { 0.0, 1.230610998905e+04, 210.056650134, 9.20757074452, 2.072879691524e+03 };
double X[] = { 0.84 };
int incX = 1;
double Y[] = { -0.711 };
int incY = -1;
double x_expected[] = { -5.70658279935 };
double y_expected[] = { 175.736586112 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 676)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 677)");
}
};
};
{
int N = 1;
double h[] = { 1.0, -1.244580625511e+03, 1.11154682624, 2.269384716089e-05, -0.0143785338883 };
double X[] = { 0.84 };
int incX = 1;
double Y[] = { -0.711 };
int incY = -1;
double x_expected[] = { -1.046158725429e+03 };
double y_expected[] = { -0.829776862405 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 678)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 679)");
}
};
};
{
int N = 1;
double h[] = { -2.0, 293.927527276, -2.614737743134e+03, 10.3164975867, -7.947030813329e+03 };
double X[] = { 0.84 };
int incX = 1;
double Y[] = { -0.711 };
int incY = -1;
double x_expected[] = { 0.84 };
double y_expected[] = { -0.711 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 680)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 681)");
}
};
};
{
int N = 1;
double h[] = { -1.0, -0.0178609251786, 0.00983044958941, 105.944529127, 1.687350579234e-05 };
double X[] = { 0.84 };
int incX = 1;
double Y[] = { -0.711 };
int incY = -1;
double x_expected[] = { -75.3415633866 };
double y_expected[] = { 0.00824558059248 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 682)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 683)");
}
};
};
{
int N = 1;
double h[] = { 0.0, 6.241999071283e-05, 2.495425882445e+03, 304.604891146, 1.604644714854e+04 };
double X[] = { 0.84 };
int incX = 1;
double Y[] = { -0.711 };
int incY = -1;
double x_expected[] = { -215.734077605 };
double y_expected[] = { 2.095446741254e+03 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 684)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 685)");
}
};
};
{
int N = 1;
double h[] = { 1.0, -0.058097639487, 8.386083625428e+03, -10.5233229994, 184.653245391 };
double X[] = { 0.84 };
int incX = 1;
double Y[] = { -0.711 };
int incY = -1;
double x_expected[] = { -0.759802017169 };
double y_expected[] = { -132.128457473 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 686)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 687)");
}
};
};
{
int N = 1;
double h[] = { -2.0, -92.8754629217, 1.467547244529e-04, -3.197881072301e-04, -1.89874629713 };
double X[] = { 0.84 };
int incX = 1;
double Y[] = { -0.711 };
int incY = -1;
double x_expected[] = { 0.84 };
double y_expected[] = { -0.711 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 688)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 689)");
}
};
};
{
int N = 1;
double h[] = { -1.0, -0.0961996230646, -2.248344186185e-05, -316.856396787, 1.663969157848e+03 };
double X[] = { 0.84 };
int incX = 1;
double Y[] = { -0.711 };
int incY = -1;
double x_expected[] = { 225.204090432 };
double y_expected[] = { -1.183082090116e+03 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 690)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 691)");
}
};
};
{
int N = 1;
double h[] = { 0.0, -201.862043128, 4.999906166451e-04, -0.0653365534487, 586.454083328 };
double X[] = { 0.84 };
int incX = 1;
double Y[] = { -0.711 };
int incY = -1;
double x_expected[] = { 0.886454289502 };
double y_expected[] = { -0.710580007882 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 692)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 693)");
}
};
};
{
int N = 1;
float h[] = { -1.0f, 162.86f, 1.379231e-04f, 9.67285f, 0.929218f };
float X[] = { 0.629f };
int incX = -1;
float Y[] = { 0.386f };
int incY = 1;
float x_expected[] = { 106.173f };
float y_expected[] = { 0.358765f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 694)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 695)");
}
};
};
{
int N = 1;
float h[] = { 0.0f, 537.387f, -21.6404f, -1.017074e+03f, -1.730546e-05f };
float X[] = { 0.629f };
int incX = -1;
float Y[] = { 0.386f };
int incY = 1;
float x_expected[] = { -391.961f };
float y_expected[] = { -13.2258f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 696)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 697)");
}
};
};
{
int N = 1;
float h[] = { 1.0f, -1.339977e-05f, 0.00522784f, 2.020352e-05f, -0.0654088f };
float X[] = { 0.629f };
int incX = -1;
float Y[] = { 0.386f };
int incY = 1;
float x_expected[] = { 0.385992f };
float y_expected[] = { -0.654248f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 698)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 699)");
}
};
};
{
int N = 1;
float h[] = { -2.0f, -50.922f, 31.5261f, -0.194913f, 0.206417f };
float X[] = { 0.629f };
int incX = -1;
float Y[] = { 0.386f };
int incY = 1;
float x_expected[] = { 0.629f };
float y_expected[] = { 0.386f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 700)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 701)");
}
};
};
{
int N = 1;
float h[] = { -1.0f, 1.15659f, 2.599832e+04f, 435.891f, 1.546671e+03f };
float X[] = { 0.629f };
int incX = -1;
float Y[] = { 0.386f };
int incY = 1;
float x_expected[] = { 168.981f };
float y_expected[] = { 1.694996e+04f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 702)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 703)");
}
};
};
{
int N = 1;
float h[] = { 0.0f, 3.359889e-04f, -0.00134822f, -12.9136f, -5.655622e+04f };
float X[] = { 0.629f };
int incX = -1;
float Y[] = { 0.386f };
int incY = 1;
float x_expected[] = { -4.35566f };
float y_expected[] = { 0.385152f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 704)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 705)");
}
};
};
{
int N = 1;
float h[] = { 1.0f, 2.75119e-05f, 1.70314f, 18.4063f, 185.731f };
float X[] = { 0.629f };
int incX = -1;
float Y[] = { 0.386f };
int incY = 1;
float x_expected[] = { 0.386017f };
float y_expected[] = { 71.063f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 706)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 707)");
}
};
};
{
int N = 1;
float h[] = { -2.0f, -1.031009e-04f, -3.378602e+04f, 7.869358e-05f, 157.303f };
float X[] = { 0.629f };
int incX = -1;
float Y[] = { 0.386f };
int incY = 1;
float x_expected[] = { 0.629f };
float y_expected[] = { 0.386f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 708)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 709)");
}
};
};
{
int N = 1;
float h[] = { -1.0f, 0.00207419f, -89.9374f, -1.40414f, -25.1433f };
float X[] = { 0.629f };
int incX = -1;
float Y[] = { 0.386f };
int incY = 1;
float x_expected[] = { -0.540694f };
float y_expected[] = { -66.276f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 710)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 711)");
}
};
};
{
int N = 1;
float h[] = { 0.0f, -4.972562e+04f, 3.65698e-05f, 632.116f, 0.195207f };
float X[] = { 0.629f };
int incX = -1;
float Y[] = { 0.386f };
int incY = 1;
float x_expected[] = { 244.626f };
float y_expected[] = { 0.386023f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 712)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 713)");
}
};
};
{
int N = 1;
double h[] = { -1.0, 8.64768339859, -105.906731008, -347.053994991, -1.28802789909 };
double X[] = { -0.674 };
int incX = -1;
double Y[] = { -0.645 };
int incY = 1;
double x_expected[] = { 218.021288159 };
double y_expected[] = { 72.2119146942 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 714)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 715)");
}
};
};
{
int N = 1;
double h[] = { 0.0, 0.926057152065, 3.315158944851e-04, -1.203638835886e+03, 0.00197484344868 };
double X[] = { -0.674 };
int incX = -1;
double Y[] = { -0.645 };
int incY = 1;
double x_expected[] = { 775.673049147 };
double y_expected[] = { -0.645223441713 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 716)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 717)");
}
};
};
{
int N = 1;
double h[] = { 1.0, -9.404298701289e-05, -0.00380843381223, -0.0767212569647, -3.66628238398 };
double X[] = { -0.674 };
int incX = -1;
double Y[] = { -0.645 };
int incY = 1;
double x_expected[] = { -0.644936615027 };
double y_expected[] = { 3.03875213767 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 718)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 719)");
}
};
};
{
int N = 1;
double h[] = { -2.0, 0.0900662226146, 0.00250500071094, 6.46624826995, -2.159443948633e-05 };
double X[] = { -0.674 };
int incX = -1;
double Y[] = { -0.645 };
int incY = 1;
double x_expected[] = { -0.674 };
double y_expected[] = { -0.645 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 720)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 721)");
}
};
};
{
int N = 1;
double h[] = { -1.0, 8.011686652935e+03, -23.8989526115, -1.104879849207e+04, 0.108740065261 };
double X[] = { -0.674 };
int incX = -1;
double Y[] = { -0.645 };
int incY = 1;
double x_expected[] = { 1.726598223305e+03 };
double y_expected[] = { 16.0377567181 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 722)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 723)");
}
};
};
{
int N = 1;
double h[] = { 0.0, 5.162681717012e-05, 48.059409562, -4.701209666609e+04, -6.80333644488e+04 };
double X[] = { -0.674 };
int incX = -1;
double Y[] = { -0.645 };
int incY = 1;
double x_expected[] = { 3.032212834963e+04 };
double y_expected[] = { -33.0370420448 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 724)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 725)");
}
};
};
{
int N = 1;
double h[] = { 1.0, -5.554806445579e-04, 5.101973060197e+04, -5.932040237374e+03, 3.91045757161 };
double X[] = { -0.674 };
int incX = -1;
double Y[] = { -0.645 };
int incY = 1;
double x_expected[] = { -0.644625606046 };
double y_expected[] = { -1.84824513369 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 726)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 727)");
}
};
};
{
int N = 1;
double h[] = { -2.0, -1.697234884626e-05, 101.466514367, 5.772202675851e+03, -6.884724590773e-04 };
double X[] = { -0.674 };
int incX = -1;
double Y[] = { -0.645 };
int incY = 1;
double x_expected[] = { -0.674 };
double y_expected[] = { -0.645 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 728)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 729)");
}
};
};
{
int N = 1;
double h[] = { -1.0, -0.0199779342753, 13.013123509, -17.8393347684, 0.129333249919 };
double X[] = { -0.674 };
int incX = -1;
double Y[] = { -0.645 };
int incY = 1;
double x_expected[] = { 11.5198360534 };
double y_expected[] = { -8.85426519126 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 730)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 731)");
}
};
};
{
int N = 1;
double h[] = { 0.0, -6.673799053773e+04, 587.759435538, 3.493966594965e+04, 2.098374142331e-05 };
double X[] = { -0.674 };
int incX = -1;
double Y[] = { -0.645 };
int incY = 1;
double x_expected[] = { -2.253675853752e+04 };
double y_expected[] = { -396.794859553 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 732)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 733)");
}
};
};
{
int N = 1;
float h[] = { -1.0f, 0.070033f, 0.034824f, -0.00740144f, -0.153474f };
float X[] = { -0.111f };
int incX = -1;
float Y[] = { -0.103f };
int incY = -1;
float x_expected[] = { -0.00701131f };
float y_expected[] = { 0.0119423f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 734)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 735)");
}
};
};
{
int N = 1;
float h[] = { 0.0f, 7.618016e-04f, -0.00396806f, -92.8408f, -0.0018571f };
float X[] = { -0.111f };
int incX = -1;
float Y[] = { -0.103f };
int incY = -1;
float x_expected[] = { 9.4516f };
float y_expected[] = { -0.10256f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 736)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 737)");
}
};
};
{
int N = 1;
float h[] = { 1.0f, -5.833806e+03f, 0.00265668f, -587.573f, 0.0972416f };
float X[] = { -0.111f };
int incX = -1;
float Y[] = { -0.103f };
int incY = -1;
float x_expected[] = { 647.449f };
float y_expected[] = { 0.100984f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 738)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 739)");
}
};
};
{
int N = 1;
float h[] = { -2.0f, -8.93339e+04f, -5.16022e-05f, 2.589784e-05f, -7.52586f };
float X[] = { -0.111f };
int incX = -1;
float Y[] = { -0.103f };
int incY = -1;
float x_expected[] = { -0.111f };
float y_expected[] = { -0.103f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 740)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 741)");
}
};
};
{
int N = 1;
float h[] = { -1.0f, 0.125135f, 0.00586453f, 1.100694e-05f, -0.0137436f };
float X[] = { -0.111f };
int incX = -1;
float Y[] = { -0.103f };
int incY = -1;
float x_expected[] = { -0.0138912f };
float y_expected[] = { 7.64631e-04f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 742)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 743)");
}
};
};
{
int N = 1;
float h[] = { 0.0f, -0.0996414f, 0.00505806f, 1.321441e-05f, 1.151406e-04f };
float X[] = { -0.111f };
int incX = -1;
float Y[] = { -0.103f };
int incY = -1;
float x_expected[] = { -0.111001f };
float y_expected[] = { -0.103561f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 744)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 745)");
}
};
};
{
int N = 1;
float h[] = { 1.0f, 8.18165f, 169.902f, -1.453316e-05f, 1.539957e+03f };
float X[] = { -0.111f };
int incX = -1;
float Y[] = { -0.103f };
int incY = -1;
float x_expected[] = { -1.01116f };
float y_expected[] = { -158.505f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 746)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 747)");
}
};
};
{
int N = 1;
float h[] = { -2.0f, 1.827623e-04f, -0.0528808f, 24.7305f, 328.39f };
float X[] = { -0.111f };
int incX = -1;
float Y[] = { -0.103f };
int incY = -1;
float x_expected[] = { -0.111f };
float y_expected[] = { -0.103f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 748)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 749)");
}
};
};
{
int N = 1;
float h[] = { -1.0f, -0.0876053f, 7.858704e+04f, -4.758389e+03f, -0.0114841f };
float X[] = { -0.111f };
int incX = -1;
float Y[] = { -0.103f };
int incY = -1;
float x_expected[] = { 490.124f };
float y_expected[] = { -8.72316e+03f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 750)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 751)");
}
};
};
{
int N = 1;
float h[] = { 0.0f, 0.00192188f, -1.031412e-05f, -0.00123957f, 0.312197f };
float X[] = { -0.111f };
int incX = -1;
float Y[] = { -0.103f };
int incY = -1;
float x_expected[] = { -0.110872f };
float y_expected[] = { -0.102999f };
cblas_srotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], flteps, "srotm(case 752)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], flteps, "srotm(case 753)");
}
};
};
{
int N = 1;
double h[] = { -1.0, -0.0253351881542, -0.105247702585, -7.18405641016, -5.409804811228e+04 };
double X[] = { 0.203 };
int incX = -1;
double Y[] = { -0.03 };
int incY = -1;
double x_expected[] = { 0.21037864911 };
double y_expected[] = { 1.622920078085e+03 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 754)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 755)");
}
};
};
{
int N = 1;
double h[] = { 0.0, 8.503080247483e+03, -6.186691885896e-05, -0.201279925805, -5.810746179529e-05 };
double X[] = { 0.203 };
int incX = -1;
double Y[] = { -0.03 };
int incY = -1;
double x_expected[] = { 0.209038397774 };
double y_expected[] = { -0.0300125589845 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 756)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 757)");
}
};
};
{
int N = 1;
double h[] = { 1.0, 0.351101212426, 64.9574703355, 3.015315809025e-05, -5.291308403203e-04 };
double X[] = { 0.203 };
int incX = -1;
double Y[] = { -0.03 };
int incY = -1;
double x_expected[] = { 0.0412735461225 };
double y_expected[] = { -0.202984126075 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 758)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 759)");
}
};
};
{
int N = 1;
double h[] = { -2.0, 0.0220262018719, -0.00311338149392, -70.6413298654, 31.8952671416 };
double X[] = { 0.203 };
int incX = -1;
double Y[] = { -0.03 };
int incY = -1;
double x_expected[] = { 0.203 };
double y_expected[] = { -0.03 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 760)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 761)");
}
};
};
{
int N = 1;
double h[] = { -1.0, 1.549812806922e+04, -4.868519165134e+04, -5.230242596804e+04, 1.58043443456e+04 };
double X[] = { 0.203 };
int incX = -1;
double Y[] = { -0.03 };
int incY = -1;
double x_expected[] = { 4.715192777093e+03 };
double y_expected[] = { -1.035722423559e+04 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 762)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 763)");
}
};
};
{
int N = 1;
double h[] = { 0.0, -3.30917942895, -0.0100316602276, -0.0222191220411, -0.0881815578726 };
double X[] = { 0.203 };
int incX = -1;
double Y[] = { -0.03 };
int incY = -1;
double x_expected[] = { 0.203666573661 };
double y_expected[] = { -0.0320364270262 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 764)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 765)");
}
};
};
{
int N = 1;
double h[] = { 1.0, 5.68327898035, 1.646867755046e-04, -0.106527931872, -28.2458905362 };
double X[] = { 0.203 };
int incX = -1;
double Y[] = { -0.03 };
int incY = -1;
double x_expected[] = { 1.12370563301 };
double y_expected[] = { 0.644376716086 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 766)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 767)");
}
};
};
{
int N = 1;
double h[] = { -2.0, 2.20585352008, 1.117638462348e+03, -0.116329468158, 0.00362096329059 };
double X[] = { 0.203 };
int incX = -1;
double Y[] = { -0.03 };
int incY = -1;
double x_expected[] = { 0.203 };
double y_expected[] = { -0.03 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 768)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 769)");
}
};
};
{
int N = 1;
double h[] = { -1.0, -0.00182683798892, -2.288460066516e-05, -37.55844708, -9.54075659826e-05 };
double X[] = { 0.203 };
int incX = -1;
double Y[] = { -0.03 };
int incY = -1;
double x_expected[] = { 1.12638256429 };
double y_expected[] = { -1.783346955549e-06 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 770)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 771)");
}
};
};
{
int N = 1;
double h[] = { 0.0, 1.02690456955e-04, -20.1292302013, -1.703870486677e-04, 5.17477399477 };
double X[] = { 0.203 };
int incX = -1;
double Y[] = { -0.03 };
int incY = -1;
double x_expected[] = { 0.203005111611 };
double y_expected[] = { -4.11623373087 };
cblas_drotm(N, X, incX, Y, incY, h);
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(X[i], x_expected[i], dbleps, "drotm(case 772)");
}
};
{
int i;
for (i = 0; i < 1; i++) {
gsl_test_rel(Y[i], y_expected[i], dbleps, "drotm(case 773)");
}
};
};
}