ideget.m

function o = ideget(options,name,default,flag)
%DDEGET  Get DDE OPTIONS parameters.
%   VAL = IDEGET(OPTIONS,'NAME') extracts the value of the named property
%   from integrator options structure OPTIONS, returning an empty matrix if
%   the property value is not specified in OPTIONS. It is sufficient to type
%   only the leading characters that uniquely identify the property. Case is
%   ignored for property names. [] is a valid OPTIONS argument.
%   
%   VAL = IDEGET(OPTIONS,'NAME',DEFAULT) extracts the named property as
%   above, but returns VAL = DEFAULT if the named property is not specified
%   in OPTIONS. For example
%   
%       val = ideget(opts,'RelTol',1e-4);
%   
%   returns val = 1e-4 if the RelTol property is not specified in opts.
%   
%   See also FCRKSET, DDE23, DDESD, DDENSD.

% undocumented usage for fast access with no error checking
if (nargin == 4) && isequal(char(flag),'fast')
  o = getknownfield(options,name,default);
  return
end

if nargin < 2
  error(message('MATLAB:fcrkget:NotEnoughInputs'));
end
if nargin < 3
  default = [];
end
if isstring(name) && isscalar(name)
  name = char(name);
end
if ~isempty(options) && ~isa(options,'struct')
  error(message('MATLAB:fcrkget:Arg1NotStruct'));
end

if isempty(options)
  o = default;
  return;
end

Names = { 'AbsTol', 'BreakPoints', 'BPOrders', 'Events', 'InitialStep',... 
    'InitialY', 'IntEqs', 'Jumps', 'MaxStep', 'NormControl','NumFlags', 'OutputFcn', ...
    'OutputSel', 'Refine', 'RelTol', 'Stats' };

j = strncmpi(name, Names, length(name));
if ~any(j)               % if no matches
  error(message('MATLAB:fcrkget:InvalidPropName', name));
elseif nnz(j) > 1            % if more than one match
  % No names are subsets of others, so there will be no exact match
  msg = strjoin(Names(j), ', ');
  error(message('MATLAB:fcrkget:AmbiguousPropName', name, msg));
end

if any(strcmp(fieldnames(options),Names{j}))
  o = options.(Names{j});
  if isempty(o)
    o = default;
  end
else
  o = default;
end

% --------------------------------------------------------------------------
function v = getknownfield(s, f, d)
%GETKNOWNFIELD  Get field f from struct s, or else yield default d.

if isfield(s,f)   % s could be empty.
  v = s.(f);
  if isempty(v)
    v = d;
  end
else
  v = d;
end