PDIPM_ABS inverse solver for absolute data using Primal/Dual interior point method img= pdipm_abs( inv_model, data); img => output image (or vector of images) inv_model => inverse model struct data => vector of eit data inv_model.pdipm_abs.norm_data 1 or 2 (DEFAULT 2) inv_model.pdipm_abs.norm_prior 1 or 2 (DEFAULT 2) inv_model.pdipm_abs.beta (default 1e-6) Parameters: max_iter = inv_model.parameters.max_iteration (default 10) Max number of iterations before stopping min change = inv_model.parameters.min_change (default 0) Min Change in objective fcn (norm(y-Jx)^2 + hp*TV(x)) before stopping beta is the parameter that smooths the TV functional
0001 function img=pdipm_abs( varargin ) 0002 % PDIPM_ABS inverse solver for absolute data using Primal/Dual interior point method 0003 % img= pdipm_abs( inv_model, data); 0004 % img => output image (or vector of images) 0005 % inv_model => inverse model struct 0006 % data => vector of eit data 0007 % 0008 % inv_model.pdipm_abs.norm_data 1 or 2 (DEFAULT 2) 0009 % inv_model.pdipm_abs.norm_prior 1 or 2 (DEFAULT 2) 0010 % inv_model.pdipm_abs.beta (default 1e-6) 0011 % 0012 % Parameters: 0013 % max_iter = inv_model.parameters.max_iteration (default 10) 0014 % Max number of iterations before stopping 0015 % min change = inv_model.parameters.min_change (default 0) 0016 % Min Change in objective fcn (norm(y-Jx)^2 + hp*TV(x)) before stopping 0017 % beta is the parameter that smooths the TV functional 0018 0019 % (C) 2010 Andrea Borsic + Andy Adler. License: GPL v2 or v3 0020 % $Id: pdipm_abs.m 3289 2012-07-01 10:40:31Z aadler $ 0021 0022 0023 warning('EIDORS:deprecated','PDIPM_ABS is deprecated as of 08-Jun-2012. Use INV_SOLVE_ABS_PDIPM instead.'); 0024 0025 if isfield(inv_model,'pdipm_abs'); 0026 inv_model.inv_solve_abs_pdipm = inv_model.pdipm_abs; 0027 end 0028 0029 img = inv_solve_abs_pdipm(varargin{:});