From cc4bffc3f70eb73fc40f5771b6a8ce973eddcb54 Mon Sep 17 00:00:00 2001
From: Martin Craig <martin.craig@eng.ox.ac.uk>
Date: Wed, 12 Feb 2020 15:23:11 +0000
Subject: [PATCH] Remove commented out old code

---
 fwdmodel_cest.cc | 27 ---------------------------
 1 file changed, 27 deletions(-)

diff --git a/fwdmodel_cest.cc b/fwdmodel_cest.cc
index 7f166dd..b488c15 100644
--- a/fwdmodel_cest.cc
+++ b/fwdmodel_cest.cc
@@ -374,33 +374,6 @@ void CESTFwdModel::EvaluateCestRstar(const ColumnVector &params, ColumnVector &r
     //LOG << "pool " << pool << endl;
     //LOG << "frac " << pool/water << endl;
 
-    // Need to use a spline interpolation
-    /* vector<double> vec_wvec;
-    vector<double> vec_water_only;
-    vector<double> vec_with_pool;
-    
-    vec_wvec.reserve(wvec.Nrows());
-    vec_water_only.reserve(wvec.Nrows());
-    vec_with_pool.reserve(wvec.Nrows());
-    
-    for (int ii = 1; ii <= wvec.Nrows(); ++ii)
-    {
-        vec_wvec.push_back(wvec(ii));
-        vec_water_only.push_back(water_only(ii));
-        vec_with_pool.push_back(with_pool(ii));
-    }
-    
-    NaturalSplineInterpolator interp_water_only(vec_wvec, vec_water_only);
-    NaturalSplineInterpolator interp_with_pool(vec_wvec, vec_with_pool);
-    
-    water = interp_water_only(ppm_eval * wlam / 1e6);
-    pool = interp_with_pool(ppm_eval * wlam / 1e6);
-    LOG << "Spline interp" << endl;
-    LOG << "water " << water << endl;
-    LOG << "pool " << pool << endl;
-    LOG << "frac " << pool/water << endl;
-    */
-    
     result.ReSize(1);
     result(1) = 100 * (water - pool) / params(1);
     // LOG << "res " << result.t() << endl;