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

   1 /* block/init_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 TYPE (gsl_block) *
  21 FUNCTION (gsl_block, alloc) (const size_t n)
  22 {
  23   TYPE (gsl_block) * b;
  24
  25   if (n == 0)
  26     {
  27       GSL_ERROR_VAL ("block length n must be positive integer",
  28                         GSL_EINVAL, 0);
  29     }
  30
  31   b = (TYPE (gsl_block) *) malloc (sizeof (TYPE (gsl_block)));
  32
  33   if (b == 0)
  34     {
  35       GSL_ERROR_VAL ("failed to allocate space for block struct",
  36                         GSL_ENOMEM, 0);
  37     }
  38
  39   b->data = (ATOMIC *) malloc (MULTIPLICITY * n * sizeof (ATOMIC));
  40
  41   if (b->data == 0)
  42     {
  43       free (b);         /* exception in constructor, avoid memory leak */
  44
  45       GSL_ERROR_VAL ("failed to allocate space for block data",
  46                         GSL_ENOMEM, 0);
  47     }
  48
  49   b->size = n;
  50
  51   return b;
  52 }
  53
  54 TYPE (gsl_block) *
  55 FUNCTION (gsl_block, calloc) (const size_t n)
  56 {
  57   size_t i;
  58
  59   TYPE (gsl_block) * b = FUNCTION (gsl_block, alloc) (n);
  60
  61   if (b == 0)
  62     return 0;
  63
  64   /* initialize block to zero */
  65
  66   for (i = 0; i < MULTIPLICITY * n; i++)
  67     {
  68       b->data[i] = 0;
  69     }
  70
  71   return b;
  72 }
  73
  74 void
  75 FUNCTION (gsl_block, free) (TYPE (gsl_block) * b)
  76 {
  77   free (b->data);
  78   free (b);
  79 }