Mercurial > repos > public > sbplib

annotate +time/+rk4/rungekutta_4RV.m @ 844:9e4e0576ca0f feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Only compute dissipation operator when requested
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 18 Sep 2018 10:32:00 +0200
parents f63b99f0729d
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
835
008496ca38f3 Compute the residual in between each runge-kutta stage.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
1 % Takes one time step of size k using the rungekutta method
008496ca38f3 Compute the residual in between each runge-kutta stage.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
2 % starting from v_0 and where the function F(v,t) gives the
008496ca38f3 Compute the residual in between each runge-kutta stage.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
3 % time derivatives.
008496ca38f3 Compute the residual in between each runge-kutta stage.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
4 function v = rungekutta_4RV(v, t , k, F, RV)
842
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
5
843
f63b99f0729d Minor rewrite of RK4RV
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 842
diff changeset
6 v1 = v + k/2*F(v,t, RV.getViscosity());
842
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
7
843
f63b99f0729d Minor rewrite of RK4RV
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 842
diff changeset
8 RV.update(v1, v, k/2);
f63b99f0729d Minor rewrite of RK4RV
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 842
diff changeset
9 v2 = v + k/2*F(v1,t+k/2, RV.getViscosity());
842
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
10
843
f63b99f0729d Minor rewrite of RK4RV
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 842
diff changeset
11 RV.update(v2, v, k/2);
f63b99f0729d Minor rewrite of RK4RV
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 842
diff changeset
12 v3 = v + k*F(v2,t+k/2, RV.getViscosity());
835
008496ca38f3 Compute the residual in between each runge-kutta stage.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
13
843
f63b99f0729d Minor rewrite of RK4RV
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 842
diff changeset
14 RV.update(v3,v,k);
f63b99f0729d Minor rewrite of RK4RV
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 842
diff changeset
15 v4 = v + k*F(v3,t+k, RV.getViscosity());
f63b99f0729d Minor rewrite of RK4RV
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 842
diff changeset
16
f63b99f0729d Minor rewrite of RK4RV
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 842
diff changeset
17 v_next = 1/6*(-3*v + 2*v1 + 4*v2 + 2*v3 + v4);
842
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
18 RV.update(v_next,v,k);
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
19 v = v_next;
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
20
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
21 % k1 = F(v, t, RV.getViscosity());
835
008496ca38f3 Compute the residual in between each runge-kutta stage.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
22
842
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
23 % RV.update(v+0.5*k*k1, v, 0.5*k);
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
24 % k2 = F(v+0.5*k*k1, t+0.5*k, RV.getViscosity());
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
25
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
26 % RV.update(v+0.5*k*k2, v, 0.5*k);
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
27 % k3 = F(v+0.5*k*k2, t+0.5*k, RV.getViscosity());
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
28
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
29 % RV.update(v+k*k3, v, k);
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
30 % k4 = F(v+k*k3,t+k, RV.getViscosity());
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
31
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
32 % RV.update(v + (1/6)*(k1+2*(k2+k3)+k4)*k, v, k);
619561e9ec0e Reformulate the RK4 time stages in order to easier see how the residual should be updated.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 835
diff changeset
33 % v = v + (1/6)*(k1+2*(k2+k3)+k4)*k;
835
008496ca38f3 Compute the residual in between each runge-kutta stage.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
34 end