(options, 
% options.empty_weight);|)
% 
% Huber prior (|options = computeWeights(options);| and |options = huberWeights(options);|)
% 
% L-filter (|options = computeWeights(options);|, |options = computeOffsets(options);| 
% and |options.a_L = lfilter_weights(options.Ndx, options.Ndy, options.Ndz, options.FOVa_x/options.Nx, 
% options.FOVa_y/options.Ny, options.axial_fov/options.Nz, options.oneD_weights);|)
% 
% FMH (|options = computeWeights(options);|, |options = computeOffsets(options);| 
% and |options