0001 function img= inv_solve_TV_pdipm( inv_model, data1, data2)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023 p= get_params(inv_model);
0024
0025 dva = calc_difference_data( data1, data2, inv_model.fwd_model);
0026
0027 back_val = get_good_background(inv_model, data1);
0028 inv_model.jacobian_bkgnd.value= back_val;
0029
0030 sol= [];
0031 for i=1:size(dva,2)
0032 [soln,dual_x]=primaldual_tvrecon_lsearch(inv_model, dva(:,i), ...
0033 p.maxiter,p.alpha1,p.alpha2, p.tol, p.beta, p.min_change);
0034
0035 if ~p.keepiters
0036 soln=soln(:,end);
0037 end
0038
0039 sol=[sol, soln];
0040 end
0041
0042 img.name= 'solved by inv_solve_TV_pdipm';
0043 img.elem_data = sol;
0044 img.fwd_model= inv_model.fwd_model;
0045 try if inv_model.inv_solve_TV_pdipm.want_dual_variable
0046 img.dual_data = dual_x;
0047 end; end
0048
0049 function p = get_params(inv_model);
0050
0051 try
0052 p.alpha1= inv_model.inv_solve_TV_pdipm.alpha1;
0053 catch
0054 p.alpha1= 1e-2;
0055 end
0056
0057 try
0058 p.beta= inv_model.inv_solve_TV_pdipm.beta;
0059 catch
0060 p.beta= 1e-4;
0061 end
0062
0063 p.alpha2= calc_hyperparameter( inv_model);
0064
0065 try
0066 p.min_change = inv_model.parameters.min_change;
0067 catch
0068 p.min_change = 0;
0069 end
0070
0071 try
0072 p.maxiter = inv_model.parameters.max_iterations;
0073 catch
0074 p.maxiter= 10;
0075 end
0076
0077 try
0078 p.keepiters = inv_model.parameters.keep_iterations;
0079 catch
0080 p.keepiters= 0;
0081 end
0082
0083 p.tol = 0;
0084
0085 function back_val = get_good_background(inv_mdl, data1);
0086
0087
0088 IM= eidors_obj('image','');
0089 IM.fwd_model= inv_mdl.fwd_model;
0090 s= ones(size(IM.fwd_model.elems,1),1);
0091 IM.elem_data= s;
0092
0093 vsim= fwd_solve( IM);
0094 back_val=abs( data1\vsim.meas );
0095 back_val=1;
0096