Mercurial > repos > public > sbplib

comparison +time/SBPInTimeSecondOrderFormImplicit.m @ 890:c70131daaa6e feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with feature/poroelastic.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 21 Nov 2018 18:29:29 -0800
parents 66eb4a2bbb72
children 8894e9c49e40
comparison
equal deleted inserted replaced
885:18e10217dca9 890:c70131daaa6e
17 function obj = SBPInTimeSecondOrderFormImplicit(A, B, C, f, k, t0, v0, v0t, do_scaling, TYPE, order, blockSize) 17 function obj = SBPInTimeSecondOrderFormImplicit(A, B, C, f, k, t0, v0, v0t, do_scaling, TYPE, order, blockSize)
18 default_arg('f', []); 18 default_arg('f', []);
19 default_arg('TYPE', []); 19 default_arg('TYPE', []);
20 default_arg('order', []); 20 default_arg('order', []);
21 default_arg('blockSize',[]); 21 default_arg('blockSize',[]);
22 default_arg('do_scaling', true); 22 default_arg('do_scaling', false);
23 23
24 m = length(v0); 24 m = length(v0);
25 25
26 default_arg('A', sparse(m, m)); 26 default_arg('A', speye(m, m));
27 default_arg('B', sparse(m, m)); 27 default_arg('B', sparse(m, m));
28 default_arg('C', sparse(m, m)); 28 default_arg('C', sparse(m, m));
29 29
30 I = speye(m); 30 I = speye(m);
31 O = sparse(m,m); 31 O = sparse(m,m);