Mercurial > repos > public > sbplib
diff +time/CdiffImplicit.m @ 704:111fcbcff2e9 feature/optim
merg with featuew grids
| author | Ylva Rydin <ylva.rydin@telia.com> |
|---|---|
| date | Fri, 03 Nov 2017 10:53:15 +0100 |
| parents | d5bce13ece23 |
| children | d6ede7f5cbf9 |
line wrap: on
line diff
--- a/+time/CdiffImplicit.m Fri Nov 03 10:43:27 2017 +0100 +++ b/+time/CdiffImplicit.m Fri Nov 03 10:53:15 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.')