pdipm_diff

 PDIPM_DIFF inverse solver for difference data using Primal/Dual interior point method
 img= ab_pdipm( inv_model, data1, data2)
 img        => output image (or vector of images)
 inv_model  => inverse model struct
 data1      => differential data at earlier time
 data2      => differential data at later time

  inv_model.pdipm_diff.norm_data  1 or 2 (DEFAULT 2)
  inv_model.pdipm_diff.norm_image 1 or 2 (DEFAULT 2)
  inv_model.pdipm_diff.beta     (default 1e-6)

 Parameters:
  max_iters =  inv_model.parameters.max_iterations (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_diff( varargin )
0002 % PDIPM_DIFF inverse solver for difference data using Primal/Dual interior point method
0003 % img= ab_pdipm( inv_model, data1, data2)
0004 % img        => output image (or vector of images)
0005 % inv_model  => inverse model struct
0006 % data1      => differential data at earlier time
0007 % data2      => differential data at later time
0008 %
0009 %  inv_model.pdipm_diff.norm_data  1 or 2 (DEFAULT 2)
0010 %  inv_model.pdipm_diff.norm_image 1 or 2 (DEFAULT 2)
0011 %  inv_model.pdipm_diff.beta     (default 1e-6)
0012 %
0013 % Parameters:
0014 %  max_iters =  inv_model.parameters.max_iterations (default 10)
0015 %      Max number of iterations before stopping
0016 %  min change = inv_model.parameters.min_change   (default 0)
0017 %      Min Change in objective fcn (norm(y-Jx)^2 + hp*TV(x)) before stopping
0018 % beta is the parameter that smooths the TV functional
0019 
0020 % (C) 2008 Andrea Borsic. License: GPL version 2 or version 3
0021 % $Id: pdipm_diff.m 3544 2012-07-25 12:51:47Z aadler $
0022 
0023 
0024 warning('EIDORS:deprecated','PDIPM_DIFF is deprecated as of 08-Jun-2012. Use INV_SOLVE_DIFF_PDIPM instead.');
0025 
0026 if isfield(varargin{1},'pdipm_diff');
0027   varargin{1}.inv_solve_diff_pdipm = varargin{1}.pdipm_diff;
0028 end
0029 
0030 img = inv_solve_diff_pdipm(varargin{:});