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

   1 #include <gsl/gsl_ntuple.h>
   2 #include <gsl/gsl_rng.h>
   3 #include <gsl/gsl_randist.h>
   4
   5 struct data
   6 {
   7   double x;
   8   double y;
   9   double z;
  10 };
  11
  12 int
  13 main (void)
  14 {
  15   const gsl_rng_type * T;
  16   gsl_rng * r;
  17
  18   struct data ntuple_row;
  19   int i;
  20
  21   gsl_ntuple *ntuple
  22     = gsl_ntuple_create ("test.dat", &ntuple_row,
  23                          sizeof (ntuple_row));
  24
  25   gsl_rng_env_setup ();
  26
  27   T = gsl_rng_default;
  28   r = gsl_rng_alloc (T);
  29
  30   for (i = 0; i < 10000; i++)
  31     {
  32       ntuple_row.x = gsl_ran_ugaussian (r);
  33       ntuple_row.y = gsl_ran_ugaussian (r);
  34       ntuple_row.z = gsl_ran_ugaussian (r);
  35
  36       gsl_ntuple_write (ntuple);
  37     }
  38
  39   gsl_ntuple_close (ntuple);
  40   gsl_rng_free (r);
  41
  42   return 0;
  43 }