|
Packit |
67cb25 |
/* vector/test_source.c
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* Copyright (C) 1996, 1997, 1998, 1999, 2000, 2007, 2010 Gerard Jungman, Brian Gough
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* This program is free software; you can redistribute it and/or modify
|
|
Packit |
67cb25 |
* it under the terms of the GNU General Public License as published by
|
|
Packit |
67cb25 |
* the Free Software Foundation; either version 3 of the License, or (at
|
|
Packit |
67cb25 |
* your option) any later version.
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* This program is distributed in the hope that it will be useful, but
|
|
Packit |
67cb25 |
* WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
Packit |
67cb25 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
Packit |
67cb25 |
* General Public License for more details.
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* You should have received a copy of the GNU General Public License
|
|
Packit |
67cb25 |
* along with this program; if not, write to the Free Software
|
|
Packit |
67cb25 |
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void FUNCTION (test, func) (size_t stride, size_t N);
|
|
Packit |
67cb25 |
void FUNCTION (test, ops) (size_t stride1, size_t stride2, size_t N);
|
|
Packit |
67cb25 |
void FUNCTION (test, file) (size_t stride, size_t N);
|
|
Packit |
67cb25 |
void FUNCTION (test, text) (size_t stride, size_t N);
|
|
Packit |
67cb25 |
void FUNCTION (test, trap) (size_t stride, size_t N);
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * FUNCTION(create, vector) (size_t stride, size_t N);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#define TEST(expr,desc) gsl_test((expr), NAME(gsl_vector) desc " stride=%d, N=%d", stride, N)
|
|
Packit |
67cb25 |
#define TEST2(expr,desc) gsl_test((expr), NAME(gsl_vector) desc " stride1=%d, stride2=%d, N=%d", stride1, stride2, N)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TYPE (gsl_vector) *
|
|
Packit |
67cb25 |
FUNCTION(create, vector) (size_t stride, size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * v = FUNCTION (gsl_vector, calloc) (N*stride);
|
|
Packit |
67cb25 |
v->stride = stride;
|
|
Packit |
67cb25 |
v->size = N;
|
|
Packit |
67cb25 |
return v;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, func) (size_t stride, size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * v0;
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * v;
|
|
Packit |
67cb25 |
QUALIFIED_VIEW(gsl_vector,view) view;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
size_t i, j;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
if (stride == 1)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
v = FUNCTION (gsl_vector, calloc) (N);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST(v->data == 0, "_calloc pointer");
|
|
Packit |
67cb25 |
TEST(v->size != N, "_calloc size");
|
|
Packit |
67cb25 |
TEST(v->stride != 1, "_calloc stride");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = (FUNCTION(gsl_vector,isnull)(v) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on calloc vector");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,ispos)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on calloc vector");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isneg)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on calloc vector");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnonneg)(v) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_isnonneg" DESC " on calloc vector");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, free) (v); /* free whatever is in v */
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
if (stride == 1)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
v = FUNCTION (gsl_vector, alloc) (N);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST(v->data == 0, "_alloc pointer");
|
|
Packit |
67cb25 |
TEST(v->size != N, "_alloc size");
|
|
Packit |
67cb25 |
TEST(v->stride != 1, "_alloc stride");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, free) (v); /* free whatever is in v */
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
if (stride == 1)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
v0 = FUNCTION (gsl_vector, alloc) (N);
|
|
Packit |
67cb25 |
view = FUNCTION (gsl_vector, subvector) (v0, 0, N);
|
|
Packit |
67cb25 |
v = &view.vector;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
else
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
v0 = FUNCTION (gsl_vector, alloc) (N * stride);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N*stride; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
v0->data[i] = i;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
view = FUNCTION (gsl_vector, subvector_with_stride) (v0, 0, stride, N);
|
|
Packit |
67cb25 |
v = &view.vector;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, (ATOMIC) i);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (v->data[i*stride] != (ATOMIC) (i))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST(status,"_set" DESC " writes into array");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (v, i) != (ATOMIC) (i))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_get" DESC " reads from array");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, ptr) (v, i) != v->data + i*stride)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_ptr" DESC " access to array");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, const_ptr) (v, i) != v->data + i*stride)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_const_ptr" DESC " access to array");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, (ATOMIC) 0);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnull)(v) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on null vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,ispos)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on null vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isneg)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on null vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnonneg)(v) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_isnonneg" DESC " on null vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, (ATOMIC) (i % 10));
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnull)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on non-negative vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,ispos)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on non-negative vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isneg)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on non-negative vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnonneg)(v) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_isnonneg" DESC " on non-negative vector") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#ifndef UNSIGNED
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
ATOMIC vi = (i % 10) - (ATOMIC) 5;
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, vi);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnull)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on mixed vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,ispos)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on mixed vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isneg)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on mixed vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnonneg)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnonneg" DESC " on mixed vector") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, -(ATOMIC) (i % 10));
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnull)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on non-positive vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,ispos)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on non-positive vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isneg)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on non-positive non-null vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnonneg)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnonneg" DESC " on non-positive non-null vector") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
#endif
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, (ATOMIC) (i % 10 + 1));
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnull)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on positive vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,ispos)(v) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on positive vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isneg)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on positive vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnonneg)(v) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_isnonneg" DESC " on positive vector") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#if (!defined(UNSIGNED) && !defined(BASE_CHAR))
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, -(ATOMIC) (i % 10 + 1));
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnull)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on negative vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,ispos)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on negative vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isneg)(v) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on negative vector") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,isnonneg)(v) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnonneg" DESC " on negative vector") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
#endif
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set_zero) (v);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (v, i) != (ATOMIC)0)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_setzero" DESC " on non-null vector") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set_all) (v, (ATOMIC)27);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (v, i) != (ATOMIC) (27))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_setall" DESC " to non-zero value") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set_basis) (v, i);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (i == j)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (v, j) != (ATOMIC)1)
|
|
Packit |
67cb25 |
status = 1 ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
else
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (v, j) != (ATOMIC)(0))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_setbasis" DESC " over range") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, (ATOMIC) i);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, scale) (v, 2.0);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (v, i) != (ATOMIC) ((ATOMIC)i*(ATOMIC)2.0))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_scale" DESC " by 2") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, add_constant) (v, (ATOMIC)7);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (v, i) != (ATOMIC) ((ATOMIC)i*(ATOMIC)2.0 + (ATOMIC)7))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_add_constant" DESC) ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, (ATOMIC) i);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector,swap_elements) (v, 2, 5) ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_vector,get)(v,2) != 5) ;
|
|
Packit |
67cb25 |
status |= (FUNCTION(gsl_vector,get)(v,5) != 2) ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector,swap_elements) (v, 2, 5) ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status |= (FUNCTION(gsl_vector,get)(v,2) != 2) ;
|
|
Packit |
67cb25 |
status |= (FUNCTION(gsl_vector,get)(v,5) != 5) ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_swap_elements" DESC " (2,5)") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector,reverse) (v) ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status |= (FUNCTION (gsl_vector, get) (v, i) != (ATOMIC) (N - i - 1));
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_reverse" DESC " reverses elements") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
QUALIFIED_VIEW(gsl_vector,view) v1 = FUNCTION(gsl_vector, view_array) (v->data, N*stride);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (&v1.vector, i*stride) != FUNCTION (gsl_vector, get) (v, i))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_view_array" DESC);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
QUALIFIED_VIEW(gsl_vector,view) v1 = FUNCTION(gsl_vector, view_array_with_stride) (v->data, stride, N*stride);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (&v1.vector, i) != FUNCTION (gsl_vector, get) (v, i))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_view_array_with_stride" DESC);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
QUALIFIED_VIEW(gsl_vector,view) v1 = FUNCTION(gsl_vector, subvector) (v, N/3, N/2);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N/2; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (&v1.vector, i) != FUNCTION (gsl_vector, get) (v, (N/3) + i))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_view_subvector" DESC);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
QUALIFIED_VIEW(gsl_vector,view) v1 = FUNCTION(gsl_vector, subvector_with_stride) (v, N/5, 3, N/4);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N/4; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (FUNCTION (gsl_vector, get) (&v1.vector, i) != FUNCTION (gsl_vector, get) (v, (N/5) + 3*i))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_view_subvector_with_stride" DESC);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE exp_max = FUNCTION(gsl_vector,get)(v, 0);
|
|
Packit |
67cb25 |
BASE exp_min = FUNCTION(gsl_vector,get)(v, 0);
|
|
Packit |
67cb25 |
size_t exp_imax = 0, exp_imin = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE k = FUNCTION(gsl_vector, get) (v, i) ;
|
|
Packit |
67cb25 |
if (k < exp_min) {
|
|
Packit |
67cb25 |
exp_min = FUNCTION(gsl_vector, get) (v, i);
|
|
Packit |
67cb25 |
exp_imin = i;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE k = FUNCTION(gsl_vector, get) (v, i) ;
|
|
Packit |
67cb25 |
if (k > exp_max) {
|
|
Packit |
67cb25 |
exp_max = FUNCTION(gsl_vector, get) (v, i) ;
|
|
Packit |
67cb25 |
exp_imax = i;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE max = FUNCTION(gsl_vector, max) (v) ;
|
|
Packit |
67cb25 |
TEST (max != exp_max, "_max returns correct maximum value");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE min = FUNCTION(gsl_vector, min) (v) ;
|
|
Packit |
67cb25 |
TEST (min != exp_min, "_min returns correct minimum value");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE min, max;
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, minmax) (v, &min, &max;;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (max != exp_max, "_minmax returns correct maximum value");
|
|
Packit |
67cb25 |
TEST (min != exp_min, "_minmax returns correct minimum value");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
size_t imax = FUNCTION(gsl_vector, max_index) (v) ;
|
|
Packit |
67cb25 |
TEST (imax != exp_imax, "_max_index returns correct maximum i");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
size_t imin = FUNCTION(gsl_vector, min_index) (v) ;
|
|
Packit |
67cb25 |
TEST (imin != exp_imin, "_min_index returns correct minimum i");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
size_t imin, imax;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, minmax_index) (v, &imin, &imax;;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (imax != exp_imax, "_minmax_index returns correct maximum i");
|
|
Packit |
67cb25 |
TEST (imin != exp_imin, "_minmax_index returns correct minimum i");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#if FP
|
|
Packit |
67cb25 |
i = N/2;
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, set) (v, i, GSL_NAN);
|
|
Packit |
67cb25 |
exp_max = GSL_NAN; exp_min = GSL_NAN;
|
|
Packit |
67cb25 |
exp_imax = i; exp_imin = i;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE max = FUNCTION(gsl_vector, max) (v) ;
|
|
Packit |
67cb25 |
gsl_test_abs (max, exp_max, 0, "_max returns correct maximum value for NaN");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE min = FUNCTION(gsl_vector, min) (v) ;
|
|
Packit |
67cb25 |
gsl_test_abs (min, exp_min, 0, "_min returns correct minimum value for NaN");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE min, max;
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, minmax) (v, &min, &max;;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test_abs (max, exp_max, 0, "_minmax returns correct maximum value for NaN");
|
|
Packit |
67cb25 |
gsl_test_abs (min, exp_min, 0, "_minmax returns correct minimum value for NaN");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
size_t imax = FUNCTION(gsl_vector, max_index) (v) ;
|
|
Packit |
67cb25 |
TEST (imax != exp_imax, "_max_index returns correct maximum i for NaN");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
size_t imin = FUNCTION(gsl_vector, min_index) (v) ;
|
|
Packit |
67cb25 |
TEST (imin != exp_imin, "_min_index returns correct minimum i for NaN");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
size_t imin, imax;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, minmax_index) (v, &imin, &imax;;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (imax != exp_imax, "_minmax_index returns correct maximum i for NaN");
|
|
Packit |
67cb25 |
TEST (imin != exp_imin, "_minmax_index returns correct minimum i for NaN");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
#endif
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, free) (v0); /* free whatever is in v */
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, ops) (size_t stride1, size_t stride2, size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
size_t i;
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * a = FUNCTION (create, vector) (stride1, N);
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * b = FUNCTION (create, vector) (stride2, N);
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * v = FUNCTION (create, vector) (stride1, N);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (a, i, (BASE)(3 + i));
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (b, i, (BASE)(3 + 2 * i));
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = (FUNCTION(gsl_vector,equal) (a,b) != 0);
|
|
Packit |
67cb25 |
TEST2 (status, "_equal vectors unequal");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, memcpy) (v, a);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = (FUNCTION(gsl_vector,equal) (a,v) != 1);
|
|
Packit |
67cb25 |
TEST2 (status, "_equal vectors equal");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, add) (v, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE r = FUNCTION(gsl_vector,get) (v,i);
|
|
Packit |
67cb25 |
BASE x = FUNCTION(gsl_vector,get) (a,i);
|
|
Packit |
67cb25 |
BASE y = FUNCTION(gsl_vector,get) (b,i);
|
|
Packit |
67cb25 |
BASE z = x + y;
|
|
Packit |
67cb25 |
if (r != z)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
TEST2 (status, "_add vector addition");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, swap) (a, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status |= (FUNCTION (gsl_vector, get) (a, i) != (BASE)(3 + 2 * i));
|
|
Packit |
67cb25 |
status |= (FUNCTION (gsl_vector, get) (b, i) != (BASE)(3 + i));
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, swap) (a, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status |= (FUNCTION (gsl_vector, get) (a, i) != (BASE)(3 + i));
|
|
Packit |
67cb25 |
status |= (FUNCTION (gsl_vector, get) (b, i) != (BASE)(3 + 2 * i));
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST2 (status, "_swap exchange vectors");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, memcpy) (v, a);
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, sub) (v, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE r = FUNCTION(gsl_vector,get) (v,i);
|
|
Packit |
67cb25 |
BASE x = FUNCTION(gsl_vector,get) (a,i);
|
|
Packit |
67cb25 |
BASE y = FUNCTION(gsl_vector,get) (b,i);
|
|
Packit |
67cb25 |
BASE z = x - y;
|
|
Packit |
67cb25 |
if (r != z)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST2 (status, "_sub vector subtraction");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, memcpy) (v, a);
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, mul) (v, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE r = FUNCTION(gsl_vector,get) (v,i);
|
|
Packit |
67cb25 |
BASE x = FUNCTION(gsl_vector,get) (a,i);
|
|
Packit |
67cb25 |
BASE y = FUNCTION(gsl_vector,get) (b,i);
|
|
Packit |
67cb25 |
BASE z = x * y;
|
|
Packit |
67cb25 |
if (r != z)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST2 (status, "_mul multiplication");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, memcpy) (v, a);
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, div) (v, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE r = FUNCTION(gsl_vector,get) (v,i);
|
|
Packit |
67cb25 |
BASE x = FUNCTION(gsl_vector,get) (a,i);
|
|
Packit |
67cb25 |
BASE y = FUNCTION(gsl_vector,get) (b,i);
|
|
Packit |
67cb25 |
BASE z = x / y;
|
|
Packit |
67cb25 |
if (fabs(r - z) > 2 * GSL_FLT_EPSILON * fabs(z))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
TEST2 (status, "_div division");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, free) (a);
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, free) (b);
|
|
Packit |
67cb25 |
FUNCTION(gsl_vector, free) (v);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, file) (size_t stride, size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * v = FUNCTION (create, vector) (stride, N);
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * w = FUNCTION (create, vector) (stride, N);
|
|
Packit |
67cb25 |
size_t i;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#ifdef NO_INLINE
|
|
Packit |
67cb25 |
char filename[] = "test_static.dat";
|
|
Packit |
67cb25 |
#else
|
|
Packit |
67cb25 |
char filename[] = "test.dat";
|
|
Packit |
67cb25 |
#endif
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
/* write file */
|
|
Packit |
67cb25 |
FILE *f = fopen(filename, "wb");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, (ATOMIC) (N - i));
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, fwrite) (f, v);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
fclose(f);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
/* read file */
|
|
Packit |
67cb25 |
FILE *f = fopen(filename, "rb");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, fread) (f, w);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (w->data[i*stride] != (ATOMIC) (N - i))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TEST (status, "_write and read");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
fclose(f);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, free) (v); /* free whatever is in v */
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, free) (w); /* free whatever is in w */
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#if USES_LONGDOUBLE && ! HAVE_PRINTF_LONGDOUBLE
|
|
Packit |
67cb25 |
/* skip this test */
|
|
Packit |
67cb25 |
#else
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, text) (size_t stride, size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * v = FUNCTION (create, vector) (stride, N);
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * w = FUNCTION (create, vector) (stride, N);
|
|
Packit |
67cb25 |
size_t i;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#ifdef NO_INLINE
|
|
Packit |
67cb25 |
char filename[] = "test_static.dat";
|
|
Packit |
67cb25 |
#else
|
|
Packit |
67cb25 |
char filename[] = "test.dat";
|
|
Packit |
67cb25 |
#endif
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
/* write file */
|
|
Packit |
67cb25 |
FILE *f = fopen(filename, "w");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, i, (ATOMIC) i);
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, fprintf) (f, v, OUT_FORMAT);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
fclose(f);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
/* read file */
|
|
Packit |
67cb25 |
FILE *f = fopen(filename, "r");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, fscanf) (f, w);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (w->data[i*stride] != (ATOMIC) i)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_vector) "_fprintf and fscanf");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
fclose (f);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, free) (v);
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, free) (w);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
#endif
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, trap) (size_t stride, size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double x;
|
|
Packit |
67cb25 |
size_t j = 0;
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * v = FUNCTION (create, vector) (stride, N);
|
|
Packit |
67cb25 |
v->size = N;
|
|
Packit |
67cb25 |
v->stride = stride;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, j - 1, (ATOMIC)0);
|
|
Packit |
67cb25 |
TEST (!status, "_set traps index below lower bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, N + 1, (ATOMIC)0);
|
|
Packit |
67cb25 |
TEST (!status, "_set traps index above upper bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, set) (v, N, (ATOMIC)0);
|
|
Packit |
67cb25 |
TEST (!status, "_set traps index at upper bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
x = FUNCTION (gsl_vector, get) (v, j - 1);
|
|
Packit |
67cb25 |
TEST (!status, "_get traps index below lower bound");
|
|
Packit |
67cb25 |
TEST (x != 0, "_get returns zero for index below lower bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
x = FUNCTION (gsl_vector, get) (v, N + 1);
|
|
Packit |
67cb25 |
TEST (!status, "_get traps index above upper bound");
|
|
Packit |
67cb25 |
TEST (x != 0, "_get returns zero for index above upper bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
x = FUNCTION (gsl_vector, get) (v, N);
|
|
Packit |
67cb25 |
TEST (!status, "_get traps index at upper bound");
|
|
Packit |
67cb25 |
TEST (x != 0, "_get returns zero for index at upper bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, free) (v); /* free whatever is in v */
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, alloc_zero_length) (void)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * b = FUNCTION (gsl_vector, alloc) (0);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test (b == 0, NAME (gsl_vector) "_alloc permits zero length");
|
|
Packit |
67cb25 |
gsl_test (b->size != 0, NAME (gsl_vector) "_alloc reflects zero length");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, free) (b);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, calloc_zero_length) (void)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_vector) * b = FUNCTION (gsl_vector, calloc) (0);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test (b == 0, NAME (gsl_vector) "_calloc permits zero length");
|
|
Packit |
67cb25 |
gsl_test (b->size != 0, NAME (gsl_vector) "_calloc reflects zero length");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_vector, free) (b);
|
|
Packit |
67cb25 |
}
|