Continuation Core and Toolboxes (COCO) Code

function [varargout] = coco(varargin)
%COCO  Top-level entry point to continuation toolboxes.
% 
%   BD = COCO(TOOLBOX, NAME, RUN, IST, ST, ...)
%   BD = COCO(OPTS, TOOLBOX, NAME, RUN, IST, ST, ...) is the main wrapper
%   for the continuation package. This call initiates the continuation of
%   solutions of a type as defined by TOOLBOX. NAME is the name of the
%   problem, typically the name of an m-file or a function handle for a
%   system of equations or an ODE. RUN is a unique name assigned to the
%   results produced during the continuation. IST specifies the type of the
%   initial solution and ST the type of solutions to be continued from the
%   starting point of type IST. COCO initialises an options structure used
%   by all continuation toolboses and some variables and then calls the
%   entry function to TOOLBOX. The name of this function is defined as
%   SPRINTF('%s_%s2%s',TOOLBOX,IST,ST). All remaining arguments after ST
%   are passed to this entry function. The optional first argument OPTS is
%   an options structure created with coco_set. This argument allows to
%   adjust the behaviour of the algorithms.
%
%   See also: coco_set
%

current_spparms = spparms();
spparms('umfpack', 0);

opts = coco_set_defaults();
argidx = 1;

if isstruct(varargin{argidx})
	opts = coco_set(opts, varargin{argidx});
	argidx = argidx + 1;
end

Toolbox = varargin{argidx  };
FName   = varargin{argidx+1};
run     = varargin{argidx+2};
from_st = varargin{argidx+3};
to_st   = varargin{argidx+4};

if (exist(FName, 'file') == 2) || (exist(FName, 'file') == 3)
	data_dir = sprintf('data/%s/%s', FName, run);
	if ~ exist(data_dir, 'dir')
		[s,msg,msgid]=mkdir(data_dir);
	end
	if ~ exist(data_dir, 'dir')
		disp(msg);
		error('%s: could not create directory ''%s''', ...
			mfilename, data_dir);
	end
	delete([data_dir '/*']);
else
	error('%s: file ''%s.m'' not found', mfilename, FName);
end

opts = coco_set(opts, 'toolbox', 'name', Toolbox);

opts = coco_set(opts, 'run', 'run',       run             );
opts = coco_set(opts, 'run', 'mfname',    FName           );
opts = coco_set(opts, 'run', 'dir',       data_dir        );
opts = coco_set(opts, 'run', 'bdfname',   [data_dir '/bd']);
opts = coco_set(opts, 'run', 'isol_type', from_st         );
opts = coco_set(opts, 'run', 'sol_type',  to_st           );

cfname = sprintf('%s_%s2%s', Toolbox, from_st, to_st);
cfhan  = str2func(cfname);

[varargout{1:nargout}] = cfhan(opts, varargin{argidx+5:end});

spparms(current_spparms);