aa_calc_system_mat

 AA_CALC_SYSTEM_MAT: SS= aa_calc_system_mat( fwd_model, img)
 Calc system matrix for Andy Adler's EIT code
 fwd_model = forward model
 img       = image background for system matrix calc
 s_mat.E = CC' * SS * conductivites * CC;
 where:
   SS  = Unconnected system Matrix
   CC  = Connectivity Matrix

0001 function s_mat= aa_calc_system_mat( fwd_model, img)
0002 % AA_CALC_SYSTEM_MAT: SS= aa_calc_system_mat( fwd_model, img)
0003 % Calc system matrix for Andy Adler's EIT code
0004 % fwd_model = forward model
0005 % img       = image background for system matrix calc
0006 % s_mat.E = CC' * SS * conductivites * CC;
0007 % where:
0008 %   SS  = Unconnected system Matrix
0009 %   CC  = Connectivity Matrix
0010 
0011 % (C) 2005 Andy Adler. License: GPL version 2 or version 3
0012 % $Id: aa_calc_system_mat.m 3087 2012-06-07 18:11:42Z aadler $
0013 
0014 warning('EIDORS:deprecated','AA_CALC_SYSTEM_MAT is deprecated as of 07-Jun-2012. Use SYSTEM_MAT_1ST_ORDER instead.');
0015 
0016 FC= system_mat_fields( fwd_model);
0017 lFC= size(FC,1);
0018 
0019 elem_sigma = kron( img.elem_data(:), ones(elem_dim(fwd_model),1) );
0020 
0021 ES= ones(lFC,1);
0022 ES(1:length(elem_sigma))= elem_sigma;
0023 ES= spdiags(ES,0,lFC,lFC);
0024 
0025 s_mat.E= FC' * ES * FC;
0026