add free_charge() to free memory used to interpolate the charge distributions

1 files changed, 13 insertions, 0 deletions

diff --git a/src/sno_charge.c b/src/sno_charge.c
index f96655b..1f39706 100644
--- a/src/sno_charge.c
+++ b/src/sno_charge.c
@@ -243,3 +243,16 @@ void init_charge(void)
     free(x);
     free(y);
 }
+
+void free_charge(void)
+{
+    size_t i;
+
+    for (i = 0; i < MAX_PE; i++) {
+        gsl_spline_free(splines[i]);
+    }
+
+    gsl_interp_accel_free(acc);
+
+    initialized = 0;
+}