Mercurial > repos > public > sbplib

comparison +time/Cdiff.m @ 978:1a30dbe99c7c

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Refactor CdiffImplicit to take input arguments in the right order
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 07 Jan 2019 16:15:49 +0100
parents c9009d5a3101
children 7a5e770974ed
comparison
equal deleted inserted replaced
977:c9009d5a3101 978:1a30dbe99c7c
3 A, B, C 3 A, B, C
4 AA, BB, CC 4 AA, BB, CC
5 G 5 G
6 k 6 k
7 t 7 t
8 v 8 v, v_prev
9 v_prev
10 n 9 n
11 end 10 end
12 11
13 12
14 methods 13 methods
30 obj.B = B; 29 obj.B = B;
31 obj.C = C; 30 obj.C = C;
32 obj.G = G; 31 obj.G = G;
33 32
34 % Rewrite as AA*v_(n+1) + BB*v_n + CC*v_(n-1) = G(t_n) 33 % Rewrite as AA*v_(n+1) + BB*v_n + CC*v_(n-1) = G(t_n)
35 obj.AA = A/k^2 + B/2/k; 34 obj.AA = A/k^2 + B/(2*k);
36 obj.BB = -2*A/k^2 + C; 35 obj.BB = -2*A/k^2 + C;
37 obj.CC = A/k^2 - B/2/k; 36 obj.CC = A/k^2 - B/(2*k);
38 37
39 obj.k = k; 38 obj.k = k;
40 obj.v_prev = v0; 39 obj.v_prev = v0;
41 obj.v = v0 + k*v0t; 40 obj.v = v0 + k*v0t;
42 obj.t = t0+k; 41 obj.t = t0+k;