test_AtQA

0001 clear all
0002 
0003 udir = 'D:\MatData\Connection\Test\';
0004 %fn = 'input_Nvact10.mat';
0005 fn = 'input.mat';
0006 
0007 load([udir fn])
0008 
0009 Nvact = 1000;
0010 Nloop = 1;
0011 profile_on = 1;
0012 
0013 % Reduce dimension for debug
0014 Delta = Delta(1:Nvact, 1:Nvact);
0015 
0016 profile_start(profile_on);
0017 
0018 [Delta_set, Ndelta, indx, deltax, Nindx, nl_indx, NL_indx] = ...
0019     init_index(Delta);
0020 
0021 init_matrix;
0022                 
0023 [AtQA0] = init_AtQA(Nvact, Ndelta, NL_indx);
0024 
0025 [sp_indx] = get_AtQA_index(Nvact, Ndelta, NL_indx, AtQA0);
0026 
0027 if profile_on==0
0028     profile_end(profile_on);
0029     profile_start(profile_on);
0030 end
0031 
0032 for n=1:Nloop
0033     [AtQA1] = cal_AtQA ...
0034       (Nvact, Ndelta, Delta_set, ax_z, mar, Vmar, indx, deltax);
0035 end
0036 
0037 if profile_on==0
0038     profile_end(profile_on);
0039     profile_start(profile_on);
0040 end
0041 
0042 for n=1:Nloop
0043 %    [AtQA2] = set_AtQA2(Nvact, Ndelta, ax_z, mar, Vmar, nl_indx, NL_indx);
0044     
0045     [AtQA2] = set_AtQA(Nvact, Ndelta, ax_z, mar, Vmar, ...
0046         nl_indx, NL_indx, sp_indx, AtQA0);
0047 end
0048 
0049 profile_end(profile_on);
0050 
0051 [err2 ] = vb_err_cell(AtQA1, AtQA2)
0052 %[err4 ] = vb_err_cell(AtQA1, AtQA4)
0053 return