VALUE rb_blas_ixamax(int argc, VALUE *argv, VALUE self)
{
Matrix *dx;
int incx;
int incy;
int n;
//char error_msg[64];
VALUE n_value, incx_value;
rb_scan_args(argc, argv, "02", &incx_value, &n_value);
Data_Get_Struct(self, Matrix, dx);
if(incx_value == Qnil)
incx = 1;
else
incx = NUM2INT(incx_value);
if(n_value == Qnil)
n = dx->nrows;
else
n = NUM2INT(n_value);
if(dx == NULL || dx->ncols != 1)
{ //sprintf(error_msg, );
rb_raise(rb_eRuntimeError, "Self is not a Vector");
}
switch(dx->data_type)
{
case Single_t: //s
return INT2FIX(cblas_isamax(n , (float *)dx->data, incx));
case Double_t: //d
return INT2FIX(cblas_idamax(n , (double *)dx->data, incx));
case Complex_t: //c
return INT2FIX(cblas_icamax(n , dx->data, incx));
case Double_Complex_t: //z
return INT2FIX(cblas_izamax(n , dx->data, incx));
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);
return Qnil; //Never reaches here.
}
}
void
test_amax (void) {
{
int N = 1;
float X[] = { -0.388f };
int incX = -1;
int expected = 0;
int k;
k = cblas_isamax(N, X, incX);
gsl_test_int(k, expected, "samax(case 52)");
};
{
int N = 1;
double X[] = { 0.247 };
int incX = -1;
int expected = 0;
int k;
k = cblas_idamax(N, X, incX);
gsl_test_int(k, expected, "damax(case 53)");
};
{
int N = 1;
float X[] = { 0.704f, 0.665f };
int incX = -1;
int expected = 0;
int k;
k = cblas_icamax(N, X, incX);
gsl_test_int(k, expected, "camax(case 54)");
};
{
int N = 1;
double X[] = { -0.599, -0.758 };
int incX = -1;
int expected = 0;
int k;
k = cblas_izamax(N, X, incX);
gsl_test_int(k, expected, "zamax(case 55)");
};
{
int N = 2;
float X[] = { 0.909f, 0.037f };
int incX = 1;
int expected = 0;
int k;
k = cblas_isamax(N, X, incX);
gsl_test_int(k, expected, "samax(case 56)");
};
{
int N = 2;
double X[] = { 0.271, -0.426 };
int incX = 1;
int expected = 1;
int k;
k = cblas_idamax(N, X, incX);
gsl_test_int(k, expected, "damax(case 57)");
};
{
int N = 2;
float X[] = { -0.648f, 0.317f, 0.62f, 0.392f };
int incX = 1;
int expected = 1;
int k;
k = cblas_icamax(N, X, incX);
gsl_test_int(k, expected, "camax(case 58)");
};
{
int N = 2;
double X[] = { -0.789, 0.352, 0.562, 0.697 };
int incX = 1;
int expected = 1;
int k;
k = cblas_izamax(N, X, incX);
gsl_test_int(k, expected, "zamax(case 59)");
};
{
int N = 2;
float X[] = { 0.487f, 0.918f };
int incX = -1;
int expected = 0;
int k;
k = cblas_isamax(N, X, incX);
gsl_test_int(k, expected, "samax(case 60)");
};
{
int N = 2;
double X[] = { 0.537, 0.826 };
int incX = -1;
int expected = 0;
int k;
k = cblas_idamax(N, X, incX);
gsl_test_int(k, expected, "damax(case 61)");
};
{
int N = 2;
float X[] = { 0.993f, 0.172f, -0.825f, 0.873f };
int incX = -1;
int expected = 0;
int k;
k = cblas_icamax(N, X, incX);
gsl_test_int(k, expected, "camax(case 62)");
};
{
int N = 2;
double X[] = { 0.913, -0.436, -0.134, 0.129 };
int incX = -1;
int expected = 0;
int k;
k = cblas_izamax(N, X, incX);
gsl_test_int(k, expected, "zamax(case 63)");
};
}
//
// Overloaded function for dispatching to
// * CBLAS backend, and
// * complex<double> value-type.
//
inline std::ptrdiff_t iamax( const int n, const std::complex<double>* x,
const int incx ) {
return cblas_izamax( n, x, incx );
}
int F77_izamax(const int *N, const void *X, const int *incX)
{
if (*N < 1 || *incX < 1) return(0);
return(cblas_izamax(*N, X, *incX)+1);
}
