test_ldgfunc.m

function OK = test_ldgfunc
% function testing ldgfunc abilities
%
% currently testing:
%   l2projection of vectorial function
%   extraction of scalar component
%   local evaluation of vectorial function
%   ldg_l2error

% Bernard Haasdonk 18.1.2010

OK = 1;

% initialize grid
grid = triagrid();

%params.nelements = grid.nelements; 
params.pdeg = 1;
params.dimrange = 2; % vectorial function
%params.ndofs = ldg_ndofs(params);
%params.ndofs_per_element = ldg_ndofs_per_element(params);
df_info = ldginfo(params,grid);

% function constant [1.0, 1.0]
f = @(einds,loc,grid,params) ones(length(einds),1)*[1.0, 1.0];

%f = @(einds,loc,grid,params) f_global(...
%    local2global(grid,einds,loc,params),params); 

qdeg = 4;
df = ldgdiscfunc(df_info);
df.dofs = ldg_l2project(f,qdeg,df_info);

df_scalar = scalar_component(df,1);
%[dofs1, params1 ] = ldg_scalar_component(dofs,1,params);
%disp('local evaluation should be constant [1 1]');

res = df(1:10,[0,0]);
if max(max(abs(res-ones(10,df.dimrange))))>1e-6
  disp('test_ldgfunc: local evaluation result of projected dofs not [1,1]!!');
  OK = 0;
end;

res = df_scalar(1:10,[0,0]);
if max(max(abs(res-ones(10,1))))>1e-6
  disp('test_ldgfunc: local evaluation of scalar component not 1!!');
  OK = 0;
end;

res = ldg_l2error(df,f,qdeg);
if abs(res)>1e-6
  disp('test_ldgfunc: ldg_l2error not 0!!');
  OK = 0;
end;

%disp('plot should be constant 1');
%plot(df_scalar,[]);