change interpolation to linear

lanl-ansi · Sep 29, 2020 · 52d5bd9 · 52d5bd9 · kaarthiksundar · Sep 29, 2020
1 parent 31e990f
commit 52d5bd9
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/CHANGELOG.md b/CHANGELOG.md
@@ -4,6 +4,7 @@
 
 ## v0.8.2 
 
+- changed default interpolation to linear
 - bug fix in reservoir constraints 
 
 ## v0.8.1 

diff --git a/src/io/transient.jl b/src/io/transient.jl
@@ -383,17 +383,19 @@ at least 4 time series data points are available (and result in an error otherwi
     for (type, id, param) in fields
         if (additional_time > 0.0)
             start_val = interpolators[type][id][param]["reduced_data_points"][1]
+            #= remove cubic spline interpolation 
             end_val = interpolators[type][id][param]["reduced_data_points"][end]
             middle_time = total_time + additional_time / 2
             middle_val = (end_val + start_val) / 2
             push!(interpolators[type][id][param]["times"], middle_time)
             push!(interpolators[type][id][param]["reduced_data_points"], middle_val)
+            =# 
             push!(interpolators[type][id][param]["times"], end_time)
             push!(interpolators[type][id][param]["reduced_data_points"], start_val)
         end
         x = interpolators[type][id][param]["times"]
         y = interpolators[type][id][param]["reduced_data_points"]
-        interpolators[type][id][param]["itp"] = Spline1D(x, y, periodic = periodic)
+        interpolators[type][id][param]["itp"] = Spline1D(x, y, k = 1, periodic = periodic)
 
         if !haskey(time_series_block, type)
             time_series_block[type] = Dict{String,Any}()