1 files changed, 121 insertions, 0 deletions
diff --git a/gsl-1.9/block/test_source.c b/gsl-1.9/block/test_source.c
new file mode 100644
index 0000000..8eab871
--- /dev/null
+++ b/gsl-1.9/block/test_source.c
@@ -0,0 +1,121 @@
+/* block/test_source.c
+ * 
+ * Copyright (C) 1996, 1997, 1998, 1999, 2000 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 2 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, trap) (void)
+{
+  TYPE (gsl_block) * b = FUNCTION (gsl_block, alloc) (0);
+
+  gsl_test (b != 0, NAME (gsl_block) "_alloc traps zero length");
+}