diff options
Diffstat (limited to 'packages/sundials/src/Main.hs')
-rw-r--r-- | packages/sundials/src/Main.hs | 12 |
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)) |