jacobian_3d

function [J] = jacobian_3d(I,elec,vtx,simp,gnd_ind,mat_ref,zc,v_f,df,tol,perm_sym);

This function calculates the Jacobian (sensitivity) matrix of the system at
mat_ref. 



I        = The currents used
elec     = the electrodes matrix
vtx      = The vertices matrix
simp     = The simplices matrix
gnd_ind  = The ground index (node)
mat_ref  = The reference conductivity vector
zc       = The electrode contact impedance vector
IntGrad  = The integrals of the gradients
v_f      = The measurement fields
df       = Measurements per current pattern as used in v_f 
tol      = Tolerance 
J        = The Jacobian (sensitivity) matrix with respect to conductivity

0001 function [J] = jacobian_3d(I,elec,vtx,simp,gnd_ind,mat_ref,zc,v_f,df,tol,perm_sym);
0002 %function [J] = jacobian_3d(I,elec,vtx,simp,gnd_ind,mat_ref,zc,v_f,df,tol,perm_sym);
0003 %
0004 %This function calculates the Jacobian (sensitivity) matrix of the system at
0005 %mat_ref.
0006 %
0007 %
0008 %
0009 %I        = The currents used
0010 %elec     = the electrodes matrix
0011 %vtx      = The vertices matrix
0012 %simp     = The simplices matrix
0013 %gnd_ind  = The ground index (node)
0014 %mat_ref  = The reference conductivity vector
0015 %zc       = The electrode contact impedance vector
0016 %IntGrad  = The integrals of the gradients
0017 %v_f      = The measurement fields
0018 %df       = Measurements per current pattern as used in v_f
0019 %tol      = Tolerance
0020 %J        = The Jacobian (sensitivity) matrix with respect to conductivity
0021 
0022 warning('EIDORS:deprecated','JACOBIAN_3D is deprecated as of 07-Jun-2012. ');
0023 
0024 [vr,vc] = size(vtx);
0025 [sr,sc] = size(simp);
0026 
0027 el_no = size(elec,1);
0028 
0029 if sum(df)~= size(v_f,2);
0030    error('Mismatched data input');
0031 end
0032 
0033 [E,D,Ela,pp] = fem_master_full(vtx,simp,mat_ref,gnd_ind,elec,zc,perm_sym);
0034 
0035 [V] = forward_solver(E,I,tol,pp);
0036 
0037 [J] = jacobian_3d_fields(V,Ela,D,elec,vtx,simp,mat_ref,v_f,df);
0038 
0039 
0040 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0041 % This is part of the EIDORS suite.
0042 % Copyright (c) N. Polydorides and W.R.B. Lionheart 2003
0043 % Copying permitted under terms of GNU GPL
0044 % See enclosed file gpl.html for details.
0045 % EIDORS 3D version 2.0
0046 % MATLAB version 5.3 R11
0047 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%