% SPLITA   Splits up an averaging kernel matrix
%
% FORMAT   As = splitA( A, ji )
%        
% OUT   As   Array of matrices. One matrix for each retrieval quantity. 
% IN    A    Complete averaging kernel matrix.
%       ji   Matrix of two columns, where the columns hold start and end index
%            for each retrieval quantity.

% 2006-09-29   Created by Patrick Eriksson.


function As = splitA(A,ji)


if ~ismatrix(A)  |  size(A,1)~=size(A,2)
  error( 'Input argument *A* must be a square matrix.' );
end
if ~ismatrix(ji)  &  size(ji,2)
  error( 'Input argument *ji* must be a matrix with 2 columns.' );
end
if ji(end,2) ~= size(A,1)
  error( 'Inconsistency between *A* and *ji*.' );
end

for i = 1:size(ji,1)
  ind   = ji(i,1) : ji(i,2);
  As{i} = A(ind,ind);
end