Read invariant torus on SSM

Contents

function sol = SSM_tor_read_solution(runid, varargin)
        

SSM_TOR_READ_SOLUTION

This function reads contination solutions stored in disk. The solutions here are obtained from continuation of 2D tori (quasi orbits) in run with runid. Each 2D torus in the slow time scale corresponds to a 3D torus or quasi-periodic orbit in the regular time dynamics. By specifying the label of solution in varargin, only this solution will be included in sol. Otherwise, all solutions will be included in sol.

FRC = SSM_TOR_READ_SOLUTION(RUN,VARARGIN)

See also: SSM_EP_READ_SOLUTION

sfname = 'SSMtor.mat';
runid  = [runid,'.tor'];
sfname = coco_fname(runid, sfname);
FRC = load(sfname);
FRC = FRC.FRC;
if isempty(varargin)
    sol = struct();
    sol.om = FRC.om;
    sol.ep = FRC.ep;
    sol.lab = FRC.lab;
else
    label = varargin{1};
    idx = find(FRC.lab==label);
    sol = struct();
    sol.om  = FRC.om(idx);
    sol.ep  = FRC.ep(idx);
    sol.tTr = FRC.tTr{idx};
    sol.xTr = FRC.zTr{idx};
    sol.z0  = FRC.Z0_frc(:,idx);
    torsol  = tor_read_solution('', runid,label);
    assert(abs(torsol.p(1)-sol.om)<=1e-5*abs(sol.om),...
        'the solution read from po is inconsistent with the one from SSM');
    sol.ttor = torsol.tbp;
    sol.xtor = torsol.xbp;
    sol.ptor = torsol.p;
end

info = struct();
info.SSMorder   = FRC.SSMorder;
info.SSMnonAuto = FRC.SSMnonAuto;
info.SSMispolar = FRC.SSMispolar;
sol.info = info;
        

end
        

Published with MATLAB® R2023a