From 38a45d34f4c96e92a2a34ccd97383449fd5207ee Mon Sep 17 00:00:00 2001
From: tlatorre <tlatorre@uchicago.edu>
Date: Wed, 4 Jul 2018 17:13:00 -0400
Subject: update

---
 quantum_efficiency.c | 11 +++++++++++
 1 file changed, 11 insertions(+)

diff --git a/quantum_efficiency.c b/quantum_efficiency.c
index 02aea06..f8ae972 100644
--- a/quantum_efficiency.c
+++ b/quantum_efficiency.c
@@ -2,12 +2,17 @@
 #include <errno.h>
 #include <string.h>
 #include <stdlib.h>
+#include <gsl/gsl_errno.h>
+#include <gsl/gsl_spline.h>
 
 static int initialized = 0;
 
 static double *x, *y;
 static size_t size;
 
+gsl_interp_accel *acc;
+gsl_spline *spline;
+
 static int init()
 {
     int i;
@@ -94,6 +99,10 @@ static int init()
 
     fclose(f);
 
+    acc = gsl_interp_accel_alloc();
+    spline = gsl_spline_alloc(gsl_interp_linear, size);
+    gsl_spline_init(spline, x, y, size);
+
     initialized = 1;
 
     return 0;
@@ -111,4 +120,6 @@ double get_quantum_efficiency(double wavelength)
             exit(1);
         }
     }
+
+    return gsl_spline_eval(spline, wavelength, acc);
 }
-- 
cgit