func_phi_product.m

function res = func_phi_product(lcoord,func,params)
%function res = func_phi_product(lcoord,func,params)
%
% auxiliary function computing the product, which is used in l2
% projecting a given analytical function. 
% lcoord is a 2-vector of local coordinates, func the analytical
% function with calling syntax func(einds,lcoord,grid,params)
% res is a vector of projected values on all elements. This
% function is used for integration and results in the dof values 
% of the l2-projected function.
% params.evaluate_basis is supposed to be a function doing
% point-evaluation of basis functions

% Bernard Haasdonk 2.2.2009

% Phi will be a ndofs_per_element x dimrange matrix
Phi = params.evaluate_basis(lcoord,params);
%Phitrans = Phi';
nel = params.grid.nelements;
% F will be a nelements x dimrange matrix
%F = func(1:nel,ones(nel,1)*lcoord(:)',grid,params);
F = func(1:nel,lcoord(:)',params.grid,params);
% small loop over dofs

res = Phi * F';


%| \docupdate