#include <config.h>

#include <gsl/gsl_math.h>

#include <gsl/gsl_errno.h>

#include <gsl/gsl_sum.h>

gsl_sum_levin_u_workspace * 

gsl_sum_levin_u_alloc (size_t n)

{

  gsl_sum_levin_u_workspace * w;

  if (n == 0)

    {

      GSL_ERROR_VAL ("length n must be positive integer", GSL_EDOM, 0);

    }

  w = (gsl_sum_levin_u_workspace *) malloc(sizeof(gsl_sum_levin_u_workspace));

  if (w == NULL)

    {

      GSL_ERROR_VAL ("failed to allocate struct", GSL_ENOMEM, 0);

    }

  w->q_num = (double *) malloc (n * sizeof (double));

  if (w->q_num == NULL)

    {

      free(w) ; /* error in constructor, prevent memory leak */

      GSL_ERROR_VAL ("failed to allocate space for q_num", GSL_ENOMEM, 0);

    }

  w->q_den = (double *) malloc (n * sizeof (double));

  if (w->q_den == NULL)

    {

      free (w->q_num);

      free (w) ; /* error in constructor, prevent memory leak */

      GSL_ERROR_VAL ("failed to allocate space for q_den", GSL_ENOMEM, 0);

    }

  w->dq_num = (double *) malloc (n * n * sizeof (double));

  if (w->dq_num == NULL)

    {

      free (w->q_den);

      free (w->q_num);

      free(w) ; /* error in constructor, prevent memory leak */

      GSL_ERROR_VAL ("failed to allocate space for dq_num", GSL_ENOMEM, 0);

    }

  w->dq_den = (double *) malloc (n * n * sizeof (double));

  if (w->dq_den == NULL)

    {

      free (w->dq_num);

      free (w->q_den);

      free (w->q_num);

      free (w) ; /* error in constructor, prevent memory leak */

      GSL_ERROR_VAL ("failed to allocate space for dq_den", GSL_ENOMEM, 0);

    }

  w->dsum = (double *) malloc (n * sizeof (double));

  if (w->dsum == NULL)

    {

      free (w->dq_den);

      free (w->dq_num);

      free (w->q_den);

      free (w->q_num);

      free (w) ; /* error in constructor, prevent memory leak */

      GSL_ERROR_VAL ("failed to allocate space for dsum", GSL_ENOMEM, 0);

    }

  w->size = n;

  w->terms_used = 0;

  w->sum_plain = 0;

  return w;

}

void

gsl_sum_levin_u_free (gsl_sum_levin_u_workspace * w)

{

  RETURN_IF_NULL (w);

  free (w->dsum);

  free (w->dq_den);

  free (w->dq_num);

  free (w->q_den);

  free (w->q_num);

  free (w);

}