annotate +util/semidiscrete2function.m @ 1012:1e437c9e5132
feature/advectionRV
Create residual viscosity package +rv and generalize the ResidualViscosity class
- Generalize residual viscosity, by passing user-defined flux and calculating the time derivative outside of the update.
- Create separate RungekuttaRV specifically using interior RV updates
- Separate the artifical dissipation operator from the scheme AdvectionRV1D so that the same scheme can be reused for creating the diff op used by the ResidualViscosity class
| author |
Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
| date |
Wed, 05 Dec 2018 13:44:10 +0100 |
| parents |
48b6fb693025 |
| children |
|
| rev |
line source |
|
0
|
1 function F = semidiscrete2function(D,S)
|
|
|
2 % TODO: Find threshhold when sparse makes sense.
|
|
|
3 % TODO: Why does having the spase block slow down convergence.m (or does it?)
|
|
|
4 if size(D,1) > 500
|
|
|
5 D = sparse(D);
|
|
|
6 end
|
|
|
7
|
|
|
8 if ~isa(S,'function_handle')
|
|
|
9 F = @(w,t)(D*w + S);
|
|
|
10 else
|
|
|
11 F = @(w,t)(D*w + S(t));
|
|
|
12 end
|
|
|
13 end
|
