0001 function [AdtQA_V] = set_AdtQA_V ...
0002 (Nvact, Ndelta, ax_z, mar, Vmar, nl_indx, NL_indx, sp_dl_indx, AdtQA_V0)
0003
0004
0005 AdtQA_V = cell(Ndelta,Ndelta);
0006 for nld = 1:Ndelta
0007 for nl = 1:Ndelta
0008 if nl ~= nld
0009 AdtQA_V{nld,nl} = sp_init_value(AdtQA_V0{nld,nl},0);
0010 end
0011 end
0012 end
0013
0014 for nv = 1:Nvact
0015 for nld = 1:Ndelta
0016 if ~isempty(nl_indx{nv,nld})
0017 for nl = 1:Ndelta
0018 if ~isempty(nl_indx{nv,nl})
0019 if nl ~= nld
0020 Vmar_fz = Vmar{nv}(nl_indx{nv,nld},nl_indx{nv,nl})/ax_z(nv);
0021
0022 sp_add_value(sp_dl_indx{nv,nld,nl}, Vmar_fz(:), AdtQA_V{nld,nl});
0023 end
0024 end
0025 end
0026 end
0027 end
0028 end