1 function res = fem_evaluate_scalar_basis_derivative(df,lcoord)
2 %
function res = fem_evaluate_scalar_basis_derivative(df,lcoord)
4 % evaluation of scalar fem reference basis
5 % derivatives `\nabla(\hat \phi_i), i=1,\ldots,m` in point
6 % `\hat x` =
'lcoord'. The argument
'df' can either be a ldg discfunc or
7 % a structure with fields pdeg and dimrange.
8 % Res is a nbasefunctions x dimworld array,
'Res(i,:)' = `\nabla \hat \phi_i`
9 %
'lcoord' is a 2-vector with the barycentric coordinates in the
12 % with V being the fem_weight_matrix we have
14 % ``(\hat \phi_1, \ldots, \hat \phi_{\mbox{nbasefct}}) =
16 % (c_1(x),c_2(x),c_k(x)) ``
18 % and `c_i(x)= v_i p(x) `
20 % with `v_i` = i-th row of V and p(x) the powervector
22 % with k = number of scalar base functions = dim(powervector2)
24 % Hence the derivatives also have similar structure.
26 % `\nabla \hat \phi_1 = [w_1 * D p(x)]' `
29 % perhaps the general vectorial case (tensor-product basis) can be
30 % generalized from this version later...
32 % Bernard Haasdonk 12.1.2011
35 warning('caution: scalar basis evaluated, but df is vectorial!');
40 V = fem_basis_weight_matrix(df.pdeg);
41 % nrep=(df.pdeg+1)*(df.pdeg+2)/2; %[3,6,10,15];
43 %res = zeros(2,nbasefct);
44 Dp = power_vector2_derivative(lcoord,df.pdeg);