|
Packit |
67cb25 |
/* matrix/test_complex_source.c
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* Copyright (C) 1996, 1997, 1998, 1999, 2000, 2007 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) (const size_t M, const size_t N);
|
|
Packit |
67cb25 |
void FUNCTION (test, ops) (const size_t P, const size_t Q);
|
|
Packit |
67cb25 |
void FUNCTION (test, trap) (const size_t M, const size_t N);
|
|
Packit |
67cb25 |
void FUNCTION (test, text) (const size_t M, const size_t N);
|
|
Packit |
67cb25 |
void FUNCTION (test, binary) (const size_t M, const size_t N);
|
|
Packit |
67cb25 |
void FUNCTION (test, binary_noncontiguous) (const size_t M, const size_t N);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#define TEST(expr,desc) gsl_test((expr), NAME(gsl_matrix) desc " M=%d, N=%d", M, N)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, func) (const size_t M, const size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
size_t i, j;
|
|
Packit |
67cb25 |
int k = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * m = FUNCTION (gsl_matrix, alloc) (M, N);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test (m->data == 0, NAME (gsl_matrix) "_alloc returns valid pointer");
|
|
Packit |
67cb25 |
gsl_test (m->size1 != M, NAME (gsl_matrix) "_alloc returns valid size1");
|
|
Packit |
67cb25 |
gsl_test (m->size2 != N, NAME (gsl_matrix) "_alloc returns valid size2");
|
|
Packit |
67cb25 |
gsl_test (m->tda != N, NAME (gsl_matrix) "_alloc returns valid tda");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = ZERO;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
GSL_REAL (z) = (ATOMIC) k;
|
|
Packit |
67cb25 |
GSL_IMAG (z) = (ATOMIC) (k + 1000);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (m, i, j, z);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
if (m->data[2 * (i * N + j)] != k ||
|
|
Packit |
67cb25 |
m->data[2 * (i * N + j) + 1] != k + 1000)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_set writes into array");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = FUNCTION (gsl_matrix, get) (m, i, j);
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
if (GSL_REAL (z) != k || GSL_IMAG (z) != k + 1000)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_get reads from array");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (m); /* free whatever is in m */
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
m = FUNCTION (gsl_matrix, calloc) (M, N);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = (FUNCTION(gsl_matrix,isnull)(m) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on calloc matrix");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,ispos)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on calloc matrix");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isneg)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on calloc matrix");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = ZERO;
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (m, i, j, z);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isnull)(m) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on null matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,ispos)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on null matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isneg)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on null matrix") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = ZERO;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
GSL_REAL (z) = (ATOMIC) (k % 10);
|
|
Packit |
67cb25 |
GSL_IMAG (z) = (ATOMIC) ((k + 5) % 10);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (m, i, j, z);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isnull)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on non-negative matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,ispos)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on non-negative matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isneg)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on non-negative matrix") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = ZERO;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
GSL_REAL (z) = (ATOMIC) ((k % 10) - 5);
|
|
Packit |
67cb25 |
GSL_IMAG (z) = (ATOMIC) (((k + 5) % 10) - 5);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (m, i, j, z);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isnull)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on mixed matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,ispos)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on mixed matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isneg)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on mixed matrix") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = ZERO;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
GSL_REAL (z) = -(ATOMIC) (k % 10);
|
|
Packit |
67cb25 |
GSL_IMAG (z) = -(ATOMIC) ((k + 5) % 10);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (m, i, j, z);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isnull)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on non-positive matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,ispos)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on non-positive matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isneg)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on non-positive matrix") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = ZERO;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
GSL_REAL (z) = (ATOMIC) (k % 10 + 1);
|
|
Packit |
67cb25 |
GSL_IMAG (z) = (ATOMIC) ((k + 5) % 10 + 1);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (m, i, j, z);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isnull)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on positive matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,ispos)(m) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on positive matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isneg)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on positive matrix") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = ZERO;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
GSL_REAL (z) = -(ATOMIC) (k % 10 + 1);
|
|
Packit |
67cb25 |
GSL_IMAG (z) = -(ATOMIC) ((k + 5) % 10 + 1);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (m, i, j, z);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isnull)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_isnull" DESC " on negative matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,ispos)(m) != 0);
|
|
Packit |
67cb25 |
TEST (status, "_ispos" DESC " on negative matrix") ;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = (FUNCTION(gsl_matrix,isneg)(m) != 1);
|
|
Packit |
67cb25 |
TEST (status, "_isneg" DESC " on negative matrix") ;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (m); /* free whatever is in m */
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#if !(USES_LONGDOUBLE && !HAVE_PRINTF_LONGDOUBLE)
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, text) (const size_t M, const size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * m = FUNCTION (gsl_matrix, alloc) (M, N);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
size_t i, j, k;
|
|
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 |
/* write file */
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FILE *f = fopen(filename, "w");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
GSL_REAL (z) = (ATOMIC) k;
|
|
Packit |
67cb25 |
GSL_IMAG (z) = (ATOMIC) (k + 1000);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (m, i, j, z);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, fprintf) (f, m, OUT_FORMAT);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
fclose(f);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* read file */
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FILE *f = fopen(filename, "r");
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * mm = FUNCTION (gsl_matrix, alloc) (M, N);
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, fscanf) (f, mm);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
if (mm->data[2 * (i * N + j)] != k
|
|
Packit |
67cb25 |
|| mm->data[2 * (i * N + j) + 1] != k + 1000)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_fprintf and fscanf");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (mm);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
fclose (f);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (m);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
#endif
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, binary) (const size_t M, const size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * m = FUNCTION (gsl_matrix, alloc) (M, N);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
size_t i, j, k;
|
|
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 |
/* write file */
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FILE *f = fopen(filename, "wb");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = ZERO;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
GSL_REAL (z) = (ATOMIC) k;
|
|
Packit |
67cb25 |
GSL_IMAG (z) = (ATOMIC) (k + 1000);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (m, i, j, z);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, fwrite) (f, m);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
fclose(f);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* read file */
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FILE *f = fopen(filename, "rb");
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * mm = FUNCTION (gsl_matrix, alloc) (M, N);
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, fread) (f, mm);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
if (mm->data[2 * (i * N + j)] != k
|
|
Packit |
67cb25 |
|| mm->data[2 * (i * N + j) + 1] != k + 1000)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_write and read");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (mm);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
fclose (f);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (m);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, binary_noncontiguous) (const size_t M, const size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * l = FUNCTION (gsl_matrix, calloc) (M+1, N+1);
|
|
Packit |
67cb25 |
VIEW (gsl_matrix, view) m = FUNCTION (gsl_matrix, submatrix) (l, 0, 0, M, N);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
size_t i, j, k;
|
|
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 |
/* write file */
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FILE *f = fopen(filename, "wb");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = ZERO;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
GSL_REAL (z) = (ATOMIC) k;
|
|
Packit |
67cb25 |
GSL_IMAG (z) = (ATOMIC) (k + 1000);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (&m.matrix, i, j, z);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, fwrite) (f, &m.matrix);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
fclose(f);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* read file */
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FILE *f = fopen(filename, "rb");
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * ll = FUNCTION (gsl_matrix, alloc) (M+1, N+1);
|
|
Packit |
67cb25 |
VIEW (gsl_matrix, view) mm = FUNCTION (gsl_matrix, submatrix) (ll, 0, 0, M, N);
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, fread) (f, &mm.matrix);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
for (i = 0; i < M; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < N; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = FUNCTION (gsl_matrix, get) (&mm.matrix, i, j);
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
if (GSL_REAL (z) != k || GSL_IMAG (z) != k + 1000)
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_write and read (noncontiguous)");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (ll);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
fclose (f);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (l);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, trap) (const size_t M, const size_t N)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * mc = FUNCTION (gsl_matrix, alloc) (M, N);
|
|
Packit |
67cb25 |
size_t i = 0, j = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
BASE z = { {(ATOMIC) 1.2, (ATOMIC) 3.4} };
|
|
Packit |
67cb25 |
BASE z1;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (mc, i - 1, j, z);
|
|
Packit |
67cb25 |
gsl_test (!status,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_set traps 1st index below lower bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (mc, i, j - 1, z);
|
|
Packit |
67cb25 |
gsl_test (!status,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_set traps 2nd index below lower bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (mc, M + 1, 0, z);
|
|
Packit |
67cb25 |
gsl_test (!status,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_set traps 1st index above upper bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (mc, 0, N + 1, z);
|
|
Packit |
67cb25 |
gsl_test (!status,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_set traps 2nd index above upper bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (mc, M, 0, z);
|
|
Packit |
67cb25 |
gsl_test (!status, NAME (gsl_matrix) "_set traps 1st index at upper bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (mc, 0, N, z);
|
|
Packit |
67cb25 |
gsl_test (!status, NAME (gsl_matrix) "_set traps 2nd index at upper bound");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
z1 = FUNCTION (gsl_matrix, get) (mc, i - 1, 0);
|
|
Packit |
67cb25 |
gsl_test (!status,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get traps 1st index below lower bound");
|
|
Packit |
67cb25 |
gsl_test (GSL_REAL (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero real for 1st index below l.b.");
|
|
Packit |
67cb25 |
gsl_test (GSL_IMAG (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero imag for 1st index below l.b.");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
z1 = FUNCTION (gsl_matrix, get) (mc, 0, j - 1);
|
|
Packit |
67cb25 |
gsl_test (!status,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get traps 2nd index below lower bound");
|
|
Packit |
67cb25 |
gsl_test (GSL_REAL (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero real for 2nd index below l.b.");
|
|
Packit |
67cb25 |
gsl_test (GSL_IMAG (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero imag for 2nd index below l.b.");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
z1 = FUNCTION (gsl_matrix, get) (mc, M + 1, 0);
|
|
Packit |
67cb25 |
gsl_test (!status,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get traps 1st index above upper bound");
|
|
Packit |
67cb25 |
gsl_test (GSL_REAL (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero real for 1st index above u.b.");
|
|
Packit |
67cb25 |
gsl_test (GSL_IMAG (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero imag for 1st index above u.b.");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
z1 = FUNCTION (gsl_matrix, get) (mc, 0, N + 1);
|
|
Packit |
67cb25 |
gsl_test (!status,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get traps 2nd index above upper bound");
|
|
Packit |
67cb25 |
gsl_test (GSL_REAL (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero real for 2nd index above u.b.");
|
|
Packit |
67cb25 |
gsl_test (GSL_IMAG (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero imag for 2nd index above u.b.");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
z1 = FUNCTION (gsl_matrix, get) (mc, M, 0);
|
|
Packit |
67cb25 |
gsl_test (!status, NAME (gsl_matrix) "_get traps 1st index at upper bound");
|
|
Packit |
67cb25 |
gsl_test (GSL_REAL (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero real for 1st index at u.b.");
|
|
Packit |
67cb25 |
gsl_test (GSL_IMAG (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero imag for 1st index at u.b.");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
z1 = FUNCTION (gsl_matrix, get) (mc, 0, N);
|
|
Packit |
67cb25 |
gsl_test (!status, NAME (gsl_matrix) "_get traps 2nd index at upper bound");
|
|
Packit |
67cb25 |
gsl_test (GSL_REAL (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero real for 2nd index at u.b.");
|
|
Packit |
67cb25 |
gsl_test (GSL_IMAG (z1) != 0,
|
|
Packit |
67cb25 |
NAME (gsl_matrix) "_get, zero imag for 2nd index at u.b.");
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (mc);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
FUNCTION (test, ops) (const size_t P, const size_t Q)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * a = FUNCTION (gsl_matrix, alloc) (P, Q);
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * b = FUNCTION (gsl_matrix, alloc) (P, Q);
|
|
Packit |
67cb25 |
TYPE (gsl_matrix) * m = FUNCTION (gsl_matrix, alloc) (P, Q);
|
|
Packit |
67cb25 |
size_t i, j;
|
|
Packit |
67cb25 |
size_t k = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
size_t status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < Q; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z, z1;
|
|
Packit |
67cb25 |
GSL_REAL (z) = (ATOMIC) k;
|
|
Packit |
67cb25 |
GSL_IMAG (z) = (ATOMIC) (k + 10);
|
|
Packit |
67cb25 |
GSL_REAL (z1) = (ATOMIC) (k + 5);
|
|
Packit |
67cb25 |
GSL_IMAG (z1) = (ATOMIC) (k + 20);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (a, i, j, z);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, set) (b, i, j, z1);
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = (FUNCTION(gsl_matrix,equal) (a,b) != 0);
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_equal matrix unequal");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, memcpy) (m, a);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = (FUNCTION(gsl_matrix,equal) (a,m) != 1);
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_equal matrix equal");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, add) (m, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < Q; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = FUNCTION (gsl_matrix, get) (m, i, j);
|
|
Packit |
67cb25 |
if (GSL_REAL (z) != (ATOMIC) (2 * k + 5) ||
|
|
Packit |
67cb25 |
GSL_IMAG (z) != (ATOMIC) (2 * k + 30))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_add matrix addition");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, memcpy) (m, a);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, sub) (m, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < Q; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE z = FUNCTION (gsl_matrix, get) (m, i, j);
|
|
Packit |
67cb25 |
if (GSL_REAL (z) != (ATOMIC) (-5)
|
|
Packit |
67cb25 |
|| GSL_IMAG (z) != (ATOMIC) (-10))
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_sub matrix subtraction");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, memcpy) (m, a);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, mul_elements) (m, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < Q; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
ATOMIC real = -(ATOMIC) (25 * k + 200);
|
|
Packit |
67cb25 |
ATOMIC imag = (ATOMIC) (2 * k * k + 35 * k + 50);
|
|
Packit |
67cb25 |
BASE z = FUNCTION (gsl_matrix, get) (m, i, j);
|
|
Packit |
67cb25 |
if (fabs (GSL_REAL (z) - real) > 100 * BASE_EPSILON ||
|
|
Packit |
67cb25 |
fabs (GSL_IMAG (z) - imag) > 100 * BASE_EPSILON)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
#ifdef DEBUG
|
|
Packit |
67cb25 |
printf ("%d\t%d\n", i, j);
|
|
Packit |
67cb25 |
printf (OUT_FORMAT "\n",
|
|
Packit |
67cb25 |
GSL_REAL (z) + (ATOMIC) (25 * (ATOMIC) k + 200));
|
|
Packit |
67cb25 |
printf (OUT_FORMAT "\n",
|
|
Packit |
67cb25 |
GSL_IMAG (z) - (ATOMIC) (2 * k * k + 35 * k + 50));
|
|
Packit |
67cb25 |
printf ("\n");
|
|
Packit |
67cb25 |
#endif
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_mul_elements multiplication");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, memcpy) (m, a);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, div_elements) (m, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < Q; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
ATOMIC denom = (2 * k * k + 50 * k + 425);
|
|
Packit |
67cb25 |
ATOMIC real = (ATOMIC) (2 * k * k + 35 * k + 200) / denom;
|
|
Packit |
67cb25 |
ATOMIC imag = ((ATOMIC) (50) - (ATOMIC) (5 * k)) / denom;
|
|
Packit |
67cb25 |
BASE z = FUNCTION (gsl_matrix, get) (m, i, j);
|
|
Packit |
67cb25 |
if (fabs (GSL_REAL (z) - real) > 100 * BASE_EPSILON ||
|
|
Packit |
67cb25 |
fabs (GSL_IMAG (z) - imag) > 100 * BASE_EPSILON)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
#ifdef DEBUG
|
|
Packit |
67cb25 |
printf (OUT_FORMAT "\t",
|
|
Packit |
67cb25 |
GSL_REAL (z) - (ATOMIC) (2 * k * k + 35 * k +
|
|
Packit |
67cb25 |
200) / denom);
|
|
Packit |
67cb25 |
printf (OUT_FORMAT "\n",
|
|
Packit |
67cb25 |
GSL_IMAG (z) - ((ATOMIC) (50) -
|
|
Packit |
67cb25 |
(ATOMIC) (5 * k)) / denom);
|
|
Packit |
67cb25 |
#endif
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_div_elements division");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE s;
|
|
Packit |
67cb25 |
GSL_SET_COMPLEX(&s, 2.0, 3.0);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, memcpy) (m, a);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, scale) (m, s);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < Q; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
ATOMIC real = (ATOMIC) (-(ATOMIC)k - 30);
|
|
Packit |
67cb25 |
ATOMIC imag = (ATOMIC) (5 * (ATOMIC)k + 20);
|
|
Packit |
67cb25 |
BASE z = FUNCTION (gsl_matrix, get) (m, i, j);
|
|
Packit |
67cb25 |
if (GSL_REAL (z) != real || GSL_IMAG (z) != imag)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_scale");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE s;
|
|
Packit |
67cb25 |
GSL_SET_COMPLEX(&s, 2.0, 3.0);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, memcpy) (m, a);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, add_constant) (m, s);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < Q; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
ATOMIC real = (ATOMIC) ((ATOMIC)k + 2);
|
|
Packit |
67cb25 |
ATOMIC imag = (ATOMIC) ((ATOMIC)k + 10 + 3);
|
|
Packit |
67cb25 |
BASE z = FUNCTION (gsl_matrix, get) (m, i, j);
|
|
Packit |
67cb25 |
if (GSL_REAL (z) != real || GSL_IMAG (z) != imag)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_add_constant");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE s;
|
|
Packit |
67cb25 |
GSL_SET_COMPLEX(&s, 2.0, 3.0);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, memcpy) (m, a);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, add_diagonal) (m, s);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < Q; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
ATOMIC real = (ATOMIC) ((ATOMIC)k + ((i==j) ? 2 : 0));
|
|
Packit |
67cb25 |
ATOMIC imag = (ATOMIC) ((ATOMIC)k + 10 +((i==j) ? 3 : 0));
|
|
Packit |
67cb25 |
BASE z = FUNCTION (gsl_matrix, get) (m, i, j);
|
|
Packit |
67cb25 |
if (GSL_REAL (z) != real || GSL_IMAG (z) != imag)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_add_diagonal");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, swap) (a, b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
k = 0;
|
|
Packit |
67cb25 |
status = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < Q; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
BASE x = FUNCTION (gsl_matrix, get) (a, i, j);
|
|
Packit |
67cb25 |
BASE y = FUNCTION (gsl_matrix, get) (b, i, j);
|
|
Packit |
67cb25 |
if (GSL_REAL (x) != (ATOMIC) (k + 5) || GSL_IMAG (x) != (ATOMIC) (k + 20) ||
|
|
Packit |
67cb25 |
GSL_REAL (y) != (ATOMIC) (k) || GSL_IMAG (y) != (ATOMIC) (k + 10))
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
status = 1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
k++;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
gsl_test (status, NAME (gsl_matrix) "_swap");
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (a);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (b);
|
|
Packit |
67cb25 |
FUNCTION (gsl_matrix, free) (m);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
}
|