bpmnr/gsl_block.c

Go to the documentation of this file.

#include <bpm/bpm_messages.h>
#include <bpm/bpm_nr.h>

gsl_block *gsl_block_alloc(const size_t n)
{
  gsl_block * b;

  if (n == 0)
    {
      bpm_error("block length n must be positive integer in gsl_block_alloc(...)",
                __FILE__, __LINE__);
      return NULL;
    }

  b = (gsl_block *) malloc (sizeof (gsl_block));

  if (b == 0)
    {
      bpm_error("failed to allocate space for block struct in gsl_block_alloc(...)",
                __FILE__, __LINE__);
      return NULL;
    }

  b->data = (double *) malloc (n * sizeof (double));

  if (b->data == 0)
    {
      free (b);         /* exception in constructor, avoid memory leak */
      
      bpm_error("failed to allocate space for block data in gsl_block_alloc(...)",
                __FILE__, __LINE__);
      return NULL;
    }

  b->size = n;

  return b;
}

void gsl_block_free (gsl_block * b)
{

  free (b->data);
  free (b);
}

---