From 56c51c7476beeab5f1b701da8339f9bc1422cd52 Mon Sep 17 00:00:00 2001
From: Dominic Steinitz <dominic@steinitz.org>
Date: Sun, 8 Apr 2018 10:38:59 +0100
Subject: Tidy and add picture

---
 packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs | 18 ++----------------
 1 file changed, 2 insertions(+), 16 deletions(-)

(limited to 'packages/sundials/src/Numeric/Sundials/ARKode')

diff --git a/packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs b/packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs
index 8f83fe7..2577b8e 100644
--- a/packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs
+++ b/packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs
@@ -35,6 +35,8 @@
 -- main = mplot (ts : toColumns sol)
 -- @
 --
+-- <<diagrams/brusselator.png#diagram=brusselator&height=400&width=500>>
+--
 -- KVAERNO_4_2_3
 --
 -- \[
@@ -70,7 +72,6 @@
 module Numeric.Sundials.ARKode.ODE ( odeSolve
                                    , odeSolveV
                                    , odeSolveVWith
-                                   , odeSolve'
                                    , getButcherTable
                                    , getBT
                                    , btGet
@@ -244,21 +245,6 @@ odeSolve f y0 ts =
     nC = length y0
     g t x0 = V.fromList $ f t (V.toList x0)
 
-odeSolve' :: ODEMethod
-         -> (Double -> [Double] -> [Double]) -- ^ The RHS of the system \(\dot{y} = f(t,y)\)
-         -> [Double]                         -- ^ initial conditions
-         -> Vector Double                    -- ^ desired solution times
-         -> Matrix Double                    -- ^ solution
-odeSolve' method f y0 ts =
-  case odeSolveVWith method (XX' 1.0e-6 1.0e-10 1 1) g (V.fromList y0) (V.fromList $ toList ts) of
-    Left c -> error $ show c -- FIXME
-    Right (v, d) -> trace (show d) $ (nR >< nC) (V.toList v)
-  where
-    us = toList ts
-    nR = length us
-    nC = length y0
-    g t x0 = V.fromList $ f t (V.toList x0)
-
 odeSolveVWith ::
   ODEMethod
   -> StepControl
-- 
cgit v1.2.3