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

   1 #include <gsl/gsl_rng.h>
   2 #include <gsl/gsl_sort_double.h>
   3
   4 int
   5 main (void)
   6 {
   7   const gsl_rng_type * T;
   8   gsl_rng * r;
   9
  10   size_t i, k = 5, N = 100000;
  11
  12   double * x = malloc (N * sizeof(double));
  13   double * small = malloc (k * sizeof(double));
  14
  15   gsl_rng_env_setup();
  16
  17   T = gsl_rng_default;
  18   r = gsl_rng_alloc (T);
  19
  20   for (i = 0; i < N; i++)
  21     {
  22       x[i] = gsl_rng_uniform(r);
  23     }
  24
  25   gsl_sort_smallest (small, k, x, 1, N);
  26
  27   printf ("%d smallest values from %d\n", k, N);
  28
  29   for (i = 0; i < k; i++)
  30     {
  31       printf ("%d: %.18f\n", i, small[i]);
  32     }
  33
  34   free (x);
  35   free (small);
  36   gsl_rng_free (r);
  37   return 0;
  38 }