get_AdtQA_V_index

 Sparse matrix index for each vertex & each delay

0001 function [sp_dl_indx] = get_AdtQA_V_index(Nvact, Ndelta, NL_indx, AdtQA_V)
0002 
0003 % Sparse matrix index for each vertex & each delay
0004 sp_dl_indx = cell(Nvact,Ndelta,Ndelta);
0005 
0006 % number of elements
0007 %Nindx_nl = zeros(Nvact,Ndelta);
0008 %Nindx_sp = zeros(Nvact,Ndelta);
0009 
0010 for nv = 1:Nvact
0011     for nl = 1:Ndelta
0012          if ~isempty(NL_indx{nv,nl})
0013         for nld = 1:Ndelta
0014          if ~isempty(NL_indx{nv,nld})
0015           if nl ~= nld
0016            [J_indx, I_indx] = meshgrid(NL_indx{nv,nl},NL_indx{nv,nld});
0017           sp_dl_indx{nv,nld,nl}  = ...
0018               sp_calc_index(I_indx(:), J_indx(:), AdtQA_V{nld,nl});
0019           end
0020         end
0021       end
0022     end
0023   end
0024 end
0025 
0026 
0027 %NLN_indx = sum(sum(Nindx_nl))
0028 %NL2_indx = sum(sum(Nindx_nl.^2))
0029 %Nsp_dl_indx = sum(sum(Nindx_sp))
0030 
0031 return