1 function [scalar_dofs, scalar_df_info] = ldg_scalar_component(df,ncomp)
2 %
function [scalar_dofs, scalar_df_info] = ldg_scalar_component(df,ncomp)
4 % extract single component ncomp of vectorial discrete
function and
5 % generate
new scalar ldg
function dof vector of same degree.
7 % params must provide params.pdeg, params.nelements, params.dimrange
10 % Bernard Haasdonk 27.1.2009
12 indices = (0:(df.ndofs/df.dimrange-1))*df.dimrange + ncomp;
13 scalar_dofs = df.dofs(indices)
';
15 scalar_params.dimrange = 1;
16 scalar_params.pdeg = df.pdeg;
18 scalar_df_info = ldginfo(scalar_params,df.grid);
19 %scalar_params.nelements = params.nelements;
20 %scalar_params.ndofs_per_element = ldg_ndofs_per_element(scalar_params);
21 %scalar_params.ndofs = ldg_ndofs(scalar_params);