Mercurial > repos > public > sbplib
diff +time/CdiffImplicit.m @ 707:0de70ec8bf60 feature/quantumTriangles
merge with feature/optim
| author | Ylva Rydin <ylva.rydin@telia.com> |
|---|---|
| date | Fri, 10 Nov 2017 14:22:56 +0100 |
| parents | d5bce13ece23 |
| children | d6ede7f5cbf9 |
line wrap: on
line diff
--- a/+time/CdiffImplicit.m Thu Oct 05 18:04:23 2017 +0200 +++ b/+time/CdiffImplicit.m Fri Nov 10 14:22:56 2017 +0100 @@ -63,7 +63,9 @@ v_prev = f1; I = speye(m); - v = (1/k^2*A)\((1/k^2*A - 1/2*B)*f1 + (1/k*I - 1/2*C)*f2 + 1/2*G(0)); + % v = (1/k^2*A)\((1/k^2*A - 1/2*B)*f1 + (1/k*I - 1/2*C)*f2 + 1/2*G(0)); + v = f1 + k*f2; + if ~issparse(A) || ~issparse(B) || ~issparse(C) error('LU factorization with full pivoting only works for sparse matrices.')