htswanalysis/MACS/lib/gsl/gsl-1.11/block/test_source.c

   1 /* block/test_source.c
   2  *
   3  * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2007 Gerard Jungman, Brian Gough
   4  *
   5  * This program is free software; you can redistribute it and/or modify
   6  * it under the terms of the GNU General Public License as published by
   7  * the Free Software Foundation; either version 3 of the License, or (at
   8  * your option) any later version.
   9  *
  10  * This program is distributed in the hope that it will be useful, but
  11  * WITHOUT ANY WARRANTY; without even the implied warranty of
  12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13  * General Public License for more details.
  14  *
  15  * You should have received a copy of the GNU General Public License
  16  * along with this program; if not, write to the Free Software
  17  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
  18  */
  19
  20 void FUNCTION (test, func) (void);
  21 void FUNCTION (test, binary) (void);
  22 void FUNCTION (test, trap) (void);
  23
  24 void
  25 FUNCTION (test, func) (void)
  26 {
  27   TYPE (gsl_block) * v;
  28   ATOMIC * data;
  29   size_t i, size;
  30
  31   v = FUNCTION (gsl_block, alloc) (N);
  32
  33   gsl_test (v->data == 0, NAME (gsl_block) "_alloc returns valid pointer");
  34   gsl_test (v->size != N, NAME (gsl_block) "_alloc returns valid size");
  35
  36   data = FUNCTION(gsl_block, data) (v);
  37   size = FUNCTION(gsl_block, size) (v);
  38
  39   gsl_test (data == 0, NAME (gsl_block) "_data returns valid pointer");
  40   gsl_test (size != N, NAME (gsl_block) "_size returns valid size");
  41
  42
  43   FUNCTION (gsl_block, free) (v);       /* free whatever is in v */
  44
  45   v = FUNCTION (gsl_block, calloc) (N);
  46
  47   gsl_test (v->data == 0, NAME (gsl_block) "_calloc returns valid pointer");
  48   gsl_test (v->size != N, NAME (gsl_block) "_calloc returns valid size");
  49
  50   data = FUNCTION(gsl_block, data) (v);
  51   size = FUNCTION(gsl_block, size) (v);
  52
  53   gsl_test (data == 0, NAME (gsl_block) "_data returns valid pointer from calloc");
  54   gsl_test (size != N, NAME (gsl_block) "_size returns valid size from calloc");
  55
  56   status = 0;
  57
  58   for (i = 0; i < N; i++)
  59     {
  60       if (v->data[i] != 0.0)
  61         status = 1;
  62     };
  63
  64   gsl_test (status, NAME (gsl_block) "_calloc initializes array to zero");
  65
  66   FUNCTION (gsl_block, free) (v);       /* free whatever is in v */
  67 }
  68
  69
  70 void
  71 FUNCTION (test, binary) (void)
  72 {
  73   size_t i;
  74
  75   {
  76     TYPE (gsl_block) * v = FUNCTION (gsl_block, calloc) (N);
  77
  78     FILE *f = fopen ("test.dat", "wb");
  79
  80     for (i = 0; i < N; i++)
  81       {
  82         v->data[i] = (ATOMIC)(N - i);
  83       };
  84
  85     FUNCTION (gsl_block, fwrite) (f, v);
  86
  87     fclose (f);
  88
  89     FUNCTION (gsl_block, free) (v);
  90   }
  91
  92   {
  93     TYPE (gsl_block) * w = FUNCTION (gsl_block, calloc) (N);
  94
  95     FILE *f = fopen ("test.dat", "rb");
  96
  97     FUNCTION (gsl_block, fread) (f, w);
  98
  99     status = 0;
 100
 101     for (i = 0; i < N; i++)
 102       {
 103         if (w->data[i] != (ATOMIC) (N - i))
 104           status = 1;
 105       };
 106
 107     fclose (f);
 108
 109     FUNCTION (gsl_block, free) (w);
 110   }
 111
 112   gsl_test (status, NAME (gsl_block) "_write and read");
 113 }
 114
 115 void
 116 FUNCTION (test, trap) (void)
 117 {
 118   TYPE (gsl_block) * b = FUNCTION (gsl_block, alloc) (0);
 119
 120   gsl_test (b != 0, NAME (gsl_block) "_alloc traps zero length");
 121 }