prior_tikhonov

 PRIOR_TIKHONOV calculate image prior
 Reg= prior_tikhonov( inv_model )
 Reg        => output regularization term
 inv_model  => inverse model struct

0001 function Reg= prior_tikhonov( inv_model );
0002 % PRIOR_TIKHONOV calculate image prior
0003 % Reg= prior_tikhonov( inv_model )
0004 % Reg        => output regularization term
0005 % inv_model  => inverse model struct
0006 
0007 % (C) 2005 Andy Adler. License: GPL version 2 or version 3
0008 % $Id: prior_tikhonov.m 4227 2013-06-20 12:57:24Z bgrychtol $
0009 
0010 %pp= fwd_model_parameters( inv_model.fwd_model );
0011 switch inv_model.type
0012   case 'inv_model'; fwd_model = inv_model.fwd_model;
0013   case 'fwd_model'; fwd_model = inv_model;
0014   otherwise; error('PRIOR_TIKHONOV requires input type of inv_model or fwd_model');
0015 end
0016 
0017 if isfield( fwd_model, 'coarse2fine' )
0018     no_dof = size(fwd_model.coarse2fine,2);
0019 else
0020     no_dof = size(fwd_model.elems,1);
0021 end
0022 
0023 Reg = speye( no_dof );
0024