annotate +time/+rk4/rungekutta_4.m @ 958:72cd29107a9a
feature/poroelastic
Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug!
| author |
Martin Almquist <malmquist@stanford.edu> |
| date |
Wed, 05 Dec 2018 18:58:10 -0800 |
| parents |
48b6fb693025 |
| children |
|
| rev |
line source |
|
0
|
1 % Takes one time step of size k using the rungekutta method
|
|
|
2 % starting from v_0 and where the function F(v,t) gives the
|
|
|
3 % time derivatives.
|
|
|
4 function v = rungekutta_4(v, t , k, F)
|
|
|
5 k1 = F(v ,t );
|
|
|
6 k2 = F(v+0.5*k*k1,t+0.5*k);
|
|
|
7 k3 = F(v+0.5*k*k2,t+0.5*k);
|
|
|
8 k4 = F(v+ k*k3,t+ k);
|
|
|
9 v = v + (1/6)*(k1+2*(k2+k3)+k4)*k;
|
|
|
10 end |
