rbmatlab  1.13.10
 All Classes Namespaces Files Functions Variables Groups Pages
fv_explicit_space.m
1 function [INC, b_I] = fv_explicit_space(model,model_data,U,NU_ind)
2 % function INC = fv_explicit_space(model, model_data, U, [NU_ind])
3 %
4 % function applying an FV-space-discretization operator starting from old
5 % values U corresponding to the geometry given in grid producing a
6 % new vector of elementwise scalars NU but only on for the
7 % subelements with numbers given in NU_ind. If NU_ind is empty, all
8 % new values NU are determined, i.e. length(NU) = length(U) =
9 % grid.nelements
10 %
11 % By this, the operator evaluation can be performed in a localized
12 % way, i.e. used for empirical interpolation in rb_nonlin_evol_simulation
13 %
14 % usual timestepping can be performed afterwards by (NU = Id - deltat *
15 % INC).
16 %
17 % required fields of model:
18 % verbose : a verbosity level
19 % fv_expl_conv_weight : factor for the explicit part of the convective
20 % discretization operator of the scheme
21 % fv_expl_diff_weight : factor for the explicit part of the diffusion
22 % discretization operator of the scheme
23 % fv_expl_react_weight : factor for the explicit part of the reaction term
24 % of the scheme
25 
26 % Martin Drohmann 9.12.2007 based on fv_conv_explicit_space by Bernard
27 % Haasdonk
28 
29 weights.diff_weight = model.fv_expl_diff_weight;
30 weights.conv_weight = model.fv_expl_conv_weight;
31 weights.react_weight = model.fv_expl_react_weight;
32 
33 INC = fv_space_operator(model,model_data,U,NU_ind,weights);
34 
35 b_I = [];
36 
37