diff +time/Cdiff.m @ 978:1a30dbe99c7c

Refactor CdiffImplicit to take input arguments in the right order
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 07 Jan 2019 16:15:49 +0100
parents c9009d5a3101
children 7a5e770974ed
line wrap: on
line diff
--- a/+time/Cdiff.m	Mon Jan 07 15:56:35 2019 +0100
+++ b/+time/Cdiff.m	Mon Jan 07 16:15:49 2019 +0100
@@ -5,8 +5,7 @@
         G
         k
         t
-        v
-        v_prev
+        v, v_prev
         n
     end
 
@@ -32,9 +31,9 @@
             obj.G = G;
 
             % Rewrite as AA*v_(n+1) + BB*v_n + CC*v_(n-1) = G(t_n)
-            obj.AA = A/k^2 + B/2/k;
+            obj.AA =    A/k^2 + B/(2*k);
             obj.BB = -2*A/k^2 + C;
-            obj.CC = A/k^2 - B/2/k;
+            obj.CC =    A/k^2 - B/(2*k);
 
             obj.k = k;
             obj.v_prev = v0;
@@ -70,4 +69,4 @@
             obj.n = obj.n + 1;
         end
     end
-end
\ No newline at end of file
+end