tests for utils folder added

Author: holgerteichgraeber

src/utils/utils.jl

@@ -1,16 +1,3 @@
-"""
-    get_EUR_to_USD(region::String)
-convert Euro to US dollars
-introduced because the clusters generated by the python script are in EUR for GER
-"""
-function get_EUR_to_USD(region::String)
-   if region =="GER"
-     ret = 1.109729
-   else
-     ret =1
-   end
-   return ret
-end
 
 """
     sort_centers(centers::Array,weights::Array)
@@ -229,20 +216,21 @@ function calc_medoids(data::Array,assignments::Array)
   return medoids
 end
 
-"""
-     resize_medoids(data::Array,centers::Array,weights::Array,assignments::Array)
-Takes in centers (typically medoids) and normalizes them such that for all clusters the average of the cluster is the same as the average of the respective original data that belongs to that cluster.
-In order to use this method of the resize function, add assignments to the function call (e.g. clustids[5,1]).
-"""
-function resize_medoids(data::Array,centers::Array,weights::Array,assignments::Array)
-    new_centers = zeros(centers)
-    for k=1:size(centers)[2] # number of clusters
-       is_in_k = assignments.==k
-       n = sum(is_in_k)
-       new_centers[:,k]=resize_medoids(reshape(data[:,is_in_k],:,n),reshape(centers[:,k] , : ,1),[1.0])# reshape is used for the side case with only one vector, so that resulting vector is 24x1 instead of 24-element
-    end
-    return new_centers
-end
+#"""
+# Not used in literature. Only uncomment if test added.
+#     resize_medoids(data::Array,centers::Array,weights::Array,assignments::Array)
+#Takes in centers (typically medoids) and normalizes them such that for all clusters the average of the cluster is the same as the average of the respective original data that belongs to that cluster.
+#In order to use this method of the resize function, add assignments to the function call (e.g. clustids[5,1]).
+#"""
+#function resize_medoids(data::Array,centers::Array,weights::Array,assignments::Array)#
+#    new_centers = zeros(centers)
+#    for k=1:size(centers)[2] # number of clusters
+#       is_in_k = assignments.==k
+#       n = sum(is_in_k)
+#       new_centers[:,k]=resize_medoids(reshape(data[:,is_in_k],:,n),reshape(centers[:,k] , : ,1),[1.0])# reshape is used for the side case with only one vector, so that resulting vector is 24x1 instead of 24-element
+#    end
+#    return new_centers
+#end
 
 
 """

test/load_data.jl

@@ -0,0 +1,21 @@
+
+@testset "load_timeseries_data" begin
+    # load one and check if all the fields are correct
+    ts_input_data = load_timeseries_data(:CEP_GER1)
+    @test ts_input_data.K == 366
+    @test ts_input_data.T == 24
+    @test all(ts_input_data.weights .== 1.)
+
+    # load all four by name - just let them run and see if they run without error
+    ts_input_data = load_timeseries_data(:DAM_GER)
+    ts_input_data = load_timeseries_data(:DAM_CA)
+    ts_input_data = load_timeseries_data(:CEP_GER1)
+    ts_input_data = load_timeseries_data(:CEP_GER18)
+
+    #load a folder by path
+    ts_input_data = load_timeseries_data(normpath(joinpath(dirname(@__FILE__),"..","data","TS_GER_1")))
+
+    # load single file by path
+    ts_input_data = load_timeseries_data(normpath(joinpath(dirname(@__FILE__),"..","data","TS_GER_1","solar.csv")))
+    @test all(["solar-germany"] .== keys(ts_input_data.data) )
+end

test/runtests.jl

@@ -8,4 +8,6 @@ include("test_utils.jl")
 
 include("run_clust.jl")
 include("extreme_vals.jl")
-include("z_normalize.jl")
+include("datastructs.jl")
+include("load_data.jl")
+include("utils.jl")

test/test_workflow_introduction.jl

@@ -1,6 +1,37 @@
-using Test
-using StatsBase
-using ClustForOpt
+@testset "sort centers" begin
+  c = [1 2 3;4 5 6]
+  w = [0.2, 0.5, 0.3]
+  c_s, w_s = sort_centers(c,w)
+  @test all(c_s .== c[:,[2,3,1]])
+  @test all(w_s .== w[[2,3,1]])
+end
+
+@testset "sakoe_chiba_band" begin
+  r=0
+  l=3
+  i2min,i2max = sakoe_chiba_band(r,l)
+  @test all(i2min .== [1,2,3])
+  @test all(i2max .== [1,2,3])
+  r=1
+  l=3
+  i2min,i2max = sakoe_chiba_band(r,l)
+  @test all(i2min .== [1,1,2])
+  @test all(i2max .== [2,3,3])
+  r=3
+  l=3
+  i2min,i2max = sakoe_chiba_band(r,l)
+  @test all(i2min .== [1,1,1])
+  @test all(i2max .== [3,3,3])
+end
+
+@testset "calc SSE" begin
+  c = [1 2 3;4 5 6]
+  a = [1,1,1]
+  sse = calc_SSE(c,a)
+  @test sse ≈ 2 + 2
+end
+
+# resize medoids
 
 @testset "z_normalize" begin
   #srand(1991)