/* block/test_source.c

 * 

 * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2007 Gerard Jungman, Brian Gough

 * 

 * This program is free software; you can redistribute it and/or modify

 * it under the terms of the GNU General Public License as published by

 * the Free Software Foundation; either version 3 of the License, or (at

 * your option) any later version.

 * 

 * This program is distributed in the hope that it will be useful, but

 * WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU

 * General Public License for more details.

 * 

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.

 */

void FUNCTION (test, func) (void);

void FUNCTION (test, binary) (void);

void FUNCTION (test, trap) (void);

void

FUNCTION (test, func) (void)

{

  TYPE (gsl_block) * v;

  ATOMIC * data;

  size_t i, size;

  v = FUNCTION (gsl_block, alloc) (N);

  gsl_test (v->data == 0, NAME (gsl_block) "_alloc returns valid pointer");

  gsl_test (v->size != N, NAME (gsl_block) "_alloc returns valid size");

  data = FUNCTION(gsl_block, data) (v);

  size = FUNCTION(gsl_block, size) (v);

  gsl_test (data == 0, NAME (gsl_block) "_data returns valid pointer");

  gsl_test (size != N, NAME (gsl_block) "_size returns valid size");

  FUNCTION (gsl_block, free) (v);       /* free whatever is in v */

  v = FUNCTION (gsl_block, calloc) (N);

  gsl_test (v->data == 0, NAME (gsl_block) "_calloc returns valid pointer");

  gsl_test (v->size != N, NAME (gsl_block) "_calloc returns valid size");

  data = FUNCTION(gsl_block, data) (v);

  size = FUNCTION(gsl_block, size) (v);

  gsl_test (data == 0, NAME (gsl_block) "_data returns valid pointer from calloc");

  gsl_test (size != N, NAME (gsl_block) "_size returns valid size from calloc");

  status = 0;

  for (i = 0; i < N; i++)

    {

      if (v->data[i] != 0.0)

        status = 1;

    };

  gsl_test (status, NAME (gsl_block) "_calloc initializes array to zero");

  FUNCTION (gsl_block, free) (v);       /* free whatever is in v */

}

void

FUNCTION (test, binary) (void)

{

  size_t i;

  {

    TYPE (gsl_block) * v = FUNCTION (gsl_block, calloc) (N);

    FILE *f = fopen ("test.dat", "wb");

    for (i = 0; i < N; i++)

      {

        v->data[i] = (ATOMIC)(N - i);

      };

    FUNCTION (gsl_block, fwrite) (f, v);

    fclose (f);

    FUNCTION (gsl_block, free) (v);

  }

  {

    TYPE (gsl_block) * w = FUNCTION (gsl_block, calloc) (N);

    FILE *f = fopen ("test.dat", "rb");

    FUNCTION (gsl_block, fread) (f, w);

    status = 0;

    for (i = 0; i < N; i++)

      {

        if (w->data[i] != (ATOMIC) (N - i))

          status = 1;

      };

    fclose (f);

    FUNCTION (gsl_block, free) (w);

  }

  gsl_test (status, NAME (gsl_block) "_write and read");

}

void

FUNCTION (test, alloc_zero_length) (void)

{

  TYPE (gsl_block) * b = FUNCTION (gsl_block, alloc) (0);

  gsl_test (b == 0, NAME (gsl_block) "_alloc permits zero length");

  gsl_test (b->size != 0, NAME (gsl_block) "_alloc reflects zero length");

  FUNCTION (gsl_block, free) (b);

}

void

FUNCTION (test, calloc_zero_length) (void)

{

  TYPE (gsl_block) * b = FUNCTION (gsl_block, calloc) (0);

  gsl_test (b == 0, NAME (gsl_block) "_calloc permits zero length");

  gsl_test (b->size != 0, NAME (gsl_block) "_calloc reflects zero length");

  FUNCTION (gsl_block, free) (b);

}