htswanalysis/MACS/lib/gsl/gsl-1.11/doc/examples/matrixw.c

   1 #include <stdio.h>
   2 #include <gsl/gsl_matrix.h>
   3
   4 int
   5 main (void)
   6 {
   7   int i, j, k = 0;
   8   gsl_matrix * m = gsl_matrix_alloc (100, 100);
   9   gsl_matrix * a = gsl_matrix_alloc (100, 100);
  10
  11   for (i = 0; i < 100; i++)
  12     for (j = 0; j < 100; j++)
  13       gsl_matrix_set (m, i, j, 0.23 + i + j);
  14
  15   {
  16      FILE * f = fopen ("test.dat", "wb");
  17      gsl_matrix_fwrite (f, m);
  18      fclose (f);
  19   }
  20
  21   {
  22      FILE * f = fopen ("test.dat", "rb");
  23      gsl_matrix_fread (f, a);
  24      fclose (f);
  25   }
  26
  27   for (i = 0; i < 100; i++)
  28     for (j = 0; j < 100; j++)
  29       {
  30         double mij = gsl_matrix_get (m, i, j);
  31         double aij = gsl_matrix_get (a, i, j);
  32         if (mij != aij) k++;
  33       }
  34
  35   gsl_matrix_free (m);
  36   gsl_matrix_free (a);
  37
  38   printf ("differences = %d (should be zero)\n", k);
  39   return (k > 0);
  40 }