mrespA

 MRESPA   Get measurement response from an averaging kernel matrix

 FORMAT   [mresp,As] = mrespA( A, ji )
        
 OUT   mresp   Measurement response.
       As      Array of matrices. One matrix for each retrieval quantity. 
               As returned by *splitA*.
 IN    A       Complete averaging kernel matrix.
       ji      Matrix of two columns, where the columns hold start and end 
               index for each retrieval quantity.

0001 % MRESPA   Get measurement response from an averaging kernel matrix
0002 %
0003 % FORMAT   [mresp,As] = mrespA( A, ji )
0004 %
0005 % OUT   mresp   Measurement response.
0006 %       As      Array of matrices. One matrix for each retrieval quantity.
0007 %               As returned by *splitA*.
0008 % IN    A       Complete averaging kernel matrix.
0009 %       ji      Matrix of two columns, where the columns hold start and end
0010 %               index for each retrieval quantity.
0011 
0012 % 2006-09-29   Created by Patrick Eriksson.
0013 
0014 
0015 function [mresp,As] = mrespA(A,ji)
0016 
0017 % Checks are done in *splitA*
0018 
0019 As = splitA( A, ji );
0020 
0021 mresp = zeros( size(A,1), 1 );
0022 
0023 i0 = 0;
0024 
0025 for i = 1:length(As)
0026   ind        = i0+1:i0+size(As{i},1);
0027   mresp(ind) = sum( A(ind,ind)' );
0028   i0         = ind(end);
0029 end
0030 
0031 
0032