rbmatlab  1.16.09
 All Classes Namespaces Files Functions Variables Modules 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 
function r = verbose(level, message, messageId)
This function displays messages depending on a message-id and/or a level. Aditionally you can set/res...
Definition: verbose.m:17
function INC = fv_conv_explicit_space(U, NU_ind,gridbase grid, params)
fv_conv_explicit_space(U,NU_ind,grid,params) function applying an FV-space-discretization operator st...
function [ INC , fluxes ] = fv_space_operator(model, model_data, U, NU_ind, weights)
fv_space_operator(model,model_data,U,NU_ind,grid,weights)