Mercurial > repos > public > sbplib
comparison +time/CdiffImplicit.m @ 383:151b08366d63 feature/beams
Fix bug in CdiffImplicit.
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Mon, 02 Jan 2017 09:22:22 +0100 |
parents | 280ae4dbf93b |
children | b361a04894e3 |
comparison
equal
deleted
inserted
replaced
382:8b8672134be8 | 383:151b08366d63 |
---|---|
101 % % Backslash | 101 % % Backslash |
102 % obj.v = obj.AA\b; | 102 % obj.v = obj.AA\b; |
103 | 103 |
104 % LU with column pivot | 104 % LU with column pivot |
105 y = obj.L\b(obj.p); | 105 y = obj.L\b(obj.p); |
106 Qx = obj.U\y; | 106 z = obj.U\y; |
107 obj.v = Qx(obj.q); | 107 obj.v(obj.q) = z; |
108 | 108 |
109 % Update time | 109 % Update time |
110 obj.t = obj.t + obj.k; | 110 obj.t = obj.t + obj.k; |
111 obj.n = obj.n + 1; | 111 obj.n = obj.n + 1; |
112 end | 112 end |