1 files changed, 94 insertions, 0 deletions
diff --git a/gsl-1.9/sum/work_u.c b/gsl-1.9/sum/work_u.c
new file mode 100644
index 0000000..3f7052b
--- /dev/null
+++ b/gsl-1.9/sum/work_u.c
@@ -0,0 +1,94 @@
+#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)
+{
+  free (w->dsum);
+  free (w->dq_den);
+  free (w->dq_num);
+  free (w->q_den);
+  free (w->q_num);
+  free (w);
+}