testLoadTemplate :: Assertion
testLoadTemplate = do
    t <- getModelTemplate (ModelByFile "exampleData/fourPop.template.txt")
    assertEqual "testLoadTemplate" t template

-- testGetNewParams :: Assertion
-- testGetNewParams = do
--     let p = V.fromList [1,1,1,1,1,0.02,0.04,0.041,0.05,1,2,3,4]
--         newP = getNewParams template p 6 0.002
--     assertEqual "switch time with overlapping join" (V.fromList [1, 1, 1, 1, 1, 0.02, 0.042, 0.041, 0.05, 1, 3, 3, 4]) newP
--     let newP' = getNewParams template p 7 0.002
--         diffVec = V.zipWith (\a b -> abs (a - b)) (V.fromList [1, 1, 1, 1, 1, 0.02, 0.04, 0.043, 0.05, 1, 2, 3, 4]) newP'
--     assertBool "no overlapping join" $ V.all (<1.0e-8) diffVec
testLoadTemplate :: Assertion
testLoadTemplate = do
    t <- getModelTemplate (ModelByFile "exampleData/fourPop.template.txt")
    assertEqual "testLoadTemplate" t template

-- testGetNewParams :: Assertion
-- testGetNewParams = do
--     let p = V.fromList [1,1,1,1,1,0.02,0.04,0.041,0.05,1,2,3,4]
--         newP = getNewParams template p 6 0.002
--     assertEqual "switch time with overlapping join" (V.fromList [1, 1, 1, 1, 1, 0.02, 0.042, 0.041, 0.05, 1, 3, 3, 4]) newP
--     let newP' = getNewParams template p 7 0.002
--         diffVec = V.zipWith (\a b -> abs (a - b)) (V.fromList [1, 1, 1, 1, 1, 0.02, 0.04, 0.043, 0.05, 1, 2, 3, 4]) newP'
--     assertBool "no overlapping join" $ V.all (<1.0e-8) diffVec
testLoadTemplate = do
    t <- getModelTemplate (ModelByFile "exampleData/fourPop.template.txt")
    assertEqual "testLoadTemplate" t template

-- testGetNewParams :: Assertion
-- testGetNewParams = do
--     let p = V.fromList [1,1,1,1,1,0.02,0.04,0.041,0.05,1,2,3,4]
--         newP = getNewParams template p 6 0.002
--     assertEqual "switch time with overlapping join" (V.fromList [1, 1, 1, 1, 1, 0.02, 0.042, 0.041, 0.05, 1, 3, 3, 4]) newP
--     let newP' = getNewParams template p 7 0.002
--         diffVec = V.zipWith (\a b -> abs (a - b)) (V.fromList [1, 1, 1, 1, 1, 0.02, 0.04, 0.043, 0.05, 1, 2, 3, 4]) newP'
--     assertBool "no overlapping join" $ V.all (<1.0e-8) diffVec