Mercurial > repos > public > sbplib

changeset 380:280ae4dbf93b feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
CdiffImplicit: Fixed a bug. Switched to LU.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 12 Dec 2016 11:35:23 +0100
parents ca73ee0623e5
children d32c749408cb
files +time/CdiffImplicit.m
diffstat 1 files changed, 9 insertions(+), 10 deletions(-) [+]
line wrap: on
line diff
diff -r ca73ee0623e5 -r 280ae4dbf93b +time/CdiffImplicit.m
--- a/+time/CdiffImplicit.m	Fri Dec 09 16:03:30 2016 +0100
+++ b/+time/CdiffImplicit.m	Mon Dec 12 11:35:23 2016 +0100
@@ -36,7 +36,7 @@
             end
 
             if isempty(G)
-                G = @(t) sparse(m,m);
+                G = @(t) sparse(m,1);
             end
 
             if isempty(f1)
@@ -96,18 +96,17 @@
 
         function obj = step(obj)
             b = obj.G(obj.t) - obj.BB*obj.v - obj.CC*obj.v_prev;
+            obj.v_prev = obj.v;
 
-            % Backslash
-            obj.v_prev = obj.v;
-            obj.v = obj.AA\b;
+            % % Backslash
+            % obj.v = obj.AA\b;
 
-            % % LU with column pivot
-            % y = obj.L\b(obj.p);
-            % Qx = U\y;
-            % v = Qx(q);
+            % LU with column pivot
+            y = obj.L\b(obj.p);
+            Qx = obj.U\y;
+            obj.v = Qx(obj.q);
 
-
-
+            % Update time
             obj.t = obj.t + obj.k;
             obj.n = obj.n + 1;
         end