summaryrefslogtreecommitdiff
path: root/packages/sundials/src/Main.hs
diff options
context:
space:
mode:
Diffstat (limited to 'packages/sundials/src/Main.hs')
-rw-r--r--packages/sundials/src/Main.hs12
1 files changed, 6 insertions, 6 deletions
diff --git a/packages/sundials/src/Main.hs b/packages/sundials/src/Main.hs
index 6d6a397..ac19e7f 100644
--- a/packages/sundials/src/Main.hs
+++ b/packages/sundials/src/Main.hs
@@ -133,17 +133,17 @@ main = do
133 -- putStrLn $ show res 133 -- putStrLn $ show res
134 -- putStrLn $ butcherTableauTex res 134 -- putStrLn $ butcherTableauTex res
135 135
136 let res1 = odeSolve' (SDIRK_5_3_4 brussJac) brusselator [1.2, 3.1, 3.0] (fromList [0.0, 0.1 .. 10.0]) 136 let res1 = odeSolve brusselator [1.2, 3.1, 3.0] (fromList [0.0, 0.1 .. 10.0])
137 renderRasterific "diagrams/brusselator.png" 137 renderRasterific "diagrams/brusselator.png"
138 (D.dims2D 500.0 500.0) 138 (D.dims2D 500.0 500.0)
139 (renderAxis $ lSaxis $ [0.0, 0.1 .. 10.0]:(toLists $ tr res1)) 139 (renderAxis $ lSaxis $ [0.0, 0.1 .. 10.0]:(toLists $ tr res1))
140 140
141 let res1a = odeSolve' (SDIRK_5_3_4') brusselator [1.2, 3.1, 3.0] (fromList [0.0, 0.1 .. 10.0]) 141 let res1a = odeSolve brusselator [1.2, 3.1, 3.0] (fromList [0.0, 0.1 .. 10.0])
142 renderRasterific "diagrams/brusselatorA.png" 142 renderRasterific "diagrams/brusselatorA.png"
143 (D.dims2D 500.0 500.0) 143 (D.dims2D 500.0 500.0)
144 (renderAxis $ lSaxis $ [0.0, 0.1 .. 10.0]:(toLists $ tr res1a)) 144 (renderAxis $ lSaxis $ [0.0, 0.1 .. 10.0]:(toLists $ tr res1a))
145 145
146 let res2 = odeSolve' (SDIRK_5_3_4 stiffJac) stiffish [0.0] (fromList [0.0, 0.1 .. 10.0]) 146 let res2 = odeSolve stiffish [0.0] (fromList [0.0, 0.1 .. 10.0])
147 putStrLn $ show res2 147 putStrLn $ show res2
148 renderRasterific "diagrams/stiffish.png" 148 renderRasterific "diagrams/stiffish.png"
149 (D.dims2D 500.0 500.0) 149 (D.dims2D 500.0 500.0)
@@ -153,17 +153,17 @@ main = do
153 putStrLn "Lower tolerances" 153 putStrLn "Lower tolerances"
154 putStrLn $ show res2a 154 putStrLn $ show res2a
155 155
156 let res3 = odeSolve' (SDIRK_5_3_4 lorenzJac) lorenz [-5.0, -5.0, 1.0] (fromList [0.0, 0.01 .. 10.0]) 156 let res3 = odeSolve lorenz [-5.0, -5.0, 1.0] (fromList [0.0, 0.01 .. 10.0])
157 putStrLn $ show $ last ((toLists $ tr res3)!!0) 157 putStrLn $ show $ last ((toLists $ tr res3)!!0)
158 158
159 let res3 = odeSolve' (SDIRK_5_3_4') lorenz [-5.0, -5.0, 1.0] (fromList [0.0, 0.01 .. 10.0]) 159 let res3 = odeSolve lorenz [-5.0, -5.0, 1.0] (fromList [0.0, 0.01 .. 10.0])
160 putStrLn $ show $ last ((toLists $ tr res3)!!0) 160 putStrLn $ show $ last ((toLists $ tr res3)!!0)
161 161
162 renderRasterific "diagrams/lorenz.png" 162 renderRasterific "diagrams/lorenz.png"
163 (D.dims2D 500.0 500.0) 163 (D.dims2D 500.0 500.0)
164 (renderAxis $ kSaxis $ zip ((toLists $ tr res3)!!0) ((toLists $ tr res3)!!1)) 164 (renderAxis $ kSaxis $ zip ((toLists $ tr res3)!!0) ((toLists $ tr res3)!!1))
165 165
166 let res3a = odeSolve' (SDIRK_5_3_4') lorenz [-5.0, -5.0, 1.0] (fromList [0.0, 0.01 .. 10.0]) 166 let res3a = odeSolve lorenz [-5.0, -5.0, 1.0] (fromList [0.0, 0.01 .. 10.0])
167 renderRasterific "diagrams/lorenzA.png" 167 renderRasterific "diagrams/lorenzA.png"
168 (D.dims2D 500.0 500.0) 168 (D.dims2D 500.0 500.0)
169 (renderAxis $ kSaxis $ zip ((toLists $ tr res3)!!0) ((toLists $ tr res3a)!!1)) 169 (renderAxis $ kSaxis $ zip ((toLists $ tr res3)!!0) ((toLists $ tr res3a)!!1))