Contents

Initialize continuation of equilibrium bifurcation (fold or Hopf)

function [bifbranch,suc]=SetupStstBifurcation(funcs,branch,ind,type,varargin)

Inputs

Important name-value pair inputs

All other named arguments are passed on to fields of bifbranch

Outputs

Parameter limits for bifbranch etc are inherited from branch, unless overridden by optional input arguments.

$Id$

process options

default={'contpar',[],'correc',true,'dir',[],'step',1e-3,'excludefreqs',[]};
[options,pass_on]=dde_set_options(default,varargin,'pass_on');
% initialize branch of bifurcations (bifbranch)
bifbranch=branch;
bifbranch=replace_branch_pars(bifbranch,options.contpar,pass_on);
point=branch.point(ind);
if ~isfield(point,'stability') || isempty(point.stability)
    point.stability=p_stabil(funcs,point,branch.method.stability);
end

create initial guess for correction

switch type
    case 'hopf'
        pini0=p_tohopf(funcs,point,options.excludefreqs);
    case 'fold'
        pini0=p_tofold(funcs,point);
    otherwise
        warning('SetupStstBifurcation:type',...
            'SetupStstBifurcation: type %s not supported',type);
        suc=0;
        return
end

correct and add 2nd point if desired

[bifbranch,suc]=correct_ini(funcs,bifbranch,pini0,...
    options.dir,options.step,options.correc);
end