fem_rhs_neumann_integral_kernel.m

function res = fem_rhs_neumann_integral_kernel(x,model,df_info,i,neumann_elind,neumann_edgeind);
%function res = fem_rhs_neumann_integral_kernel(x,model,df_info,i,neumann_elind,neumann_edgeind);
%
% auxiliary function for integral kernel for r_neumann
% integral kernel on all neumann-edges simultaneously
% here x is a scalar value on unit interval
% f(x) = g_N * hatphi_i  
% multiplication with |edge| is realized in caller after quadrature.
% function can handle cell-array valued data

% B. Haasdonk 22.2.2011
% I. Maier 24.03.2011

res = zeros(0,1);
for local_edge_id = 1:3
  lcoord = llocal2local(df_info.grid,local_edge_id,x); % local coord for x on edge 1    
  hat_phi_i = fem_evaluate_basis_function(df_info,lcoord,i);
  inds = find(neumann_edgeind==local_edge_id);
  g_N = model.neumann_values(df_info.grid,neumann_elind(inds), ...
                             local_edge_id,x,model); % local model!
  if ~iscell(g_N)
    res = [res;...
           g_N * hat_phi_i; ...
          ];
  else
    if ~iscell(res)
      res = cell(1,length(g_N));
      for q = 1:length(g_N)
        res{q} = zeros(0,1);
      end;
    end;
    for q = 1:length(g_N)
      res{q} = [res{q};...
                g_N{q} * hat_phi_i; ...
               ];
    end;
  end;
end;