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 [vr,vc] = size(vtx);
0023 [sr,sc] = size(simp);
0024 
0025 el_no = size(elec,1);
0026 
0027 if sum(df)~= size(v_f,2);
0028    error('Mismatched data input');
0029 end
0030 
0031 [E,D,Ela,pp] = fem_master_full(vtx,simp,mat_ref,gnd_ind,elec,zc,perm_sym);
0032 
0033 [V] = forward_solver(E,I,tol,pp);
0034 
0035 [J] = jacobian_3d_fields(V,Ela,D,elec,vtx,simp,mat_ref,v_f,df);
0036 
0037 
0038 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0039 % This is part of the EIDORS suite.
0040 % Copyright (c) N. Polydorides and W.R.B. Lionheart 2003
0041 % Copying permitted under terms of GNU GPL
0042 % See enclosed file gpl.html for details.
0043 % EIDORS 3D version 2.0
0044 % MATLAB version 5.3 R11
0045 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%