discfunc/ldg/ldg_evaluate.m

function res = ldg_evaluate(df, eindices, lcoord)
%function res = ldg_evaluate(df, eindices, lcoord)
%
% method evaluating a ldg function in local coordinates in the point
% `\hat x` = lcoord in several elements simultaneously given by 
% indices eindices.
% res is a length(eindices) x dimrange vector 

% Bernard Haasdonk 28.1.2009

% evaluate all reference basis functions in the lcoords, 
% is a ndofs_per_element x dimrange matrix
basis_values = ldg_evaluate_basis(lcoord,df);

%res = zeros(params.dimrange,length(eindices));
%res = zeros(length(eindices),params.dimrange);
% linear combination with DOFS to get function values

d = reshape(df.dofs,df.ndofs_per_element,df.nelements);

% small loop over basis functions
%for i = 1:params.ndofs_per_element
%%  res = res + basis_values(:,i) * d(i,eindices);
%  res = res + d(eindices,i)* basis_values(i,:);
%end;

%res = basis_values' * d(:,eindices); 
res = d(:,eindices)' * basis_values; 

%ldg_evaluate_basis([0,0],params)'*dofs(1,:)'
%| \docupdate