runAnalyzer.m

function runtime = runAnalyzer(self, options)
% RUNTIME = RUN()
% Run the solver
 
tic;
if nargin<2
    options = self.getOptions;
end
line_debug(options, 'LQNS: starting (method=%s, multiserver=%s)', options.method, options.config.multiserver);
 
% Resolve the execution backend before choosing the temp dir. options.config.container:
%   ''         -> native binary; fall back to Docker only if lqns is missing
%   'auto'/true-> prefer Docker when a local imperialqore/lqns image exists
%   '<image>'  -> use the named image if present locally
% Docker sidesteps the MATLAB libstdc++/GLIBCXX clash noted below and pins a
% reproducible lqns version, so tests opt into 'auto'.
dockerImage = '';
useRemote = isfield(options.config, 'remote') && options.config.remote;
if isunix && ~useRemote
    reqContainer = '';
    if isfield(options.config, 'container')
        reqContainer = options.config.container;
    end
    if isempty(reqContainer) && ~SolverLQNS.hasLocalBinary()
        % Native lqns absent: transparently fall back to a Docker image.
        reqContainer = 'auto';
    end
    if ~isempty(reqContainer)
        dockerImage = SolverLQNS.getDockerImage(reqContainer);
    end
end
 
% Snap-confined Docker cannot bind-mount the system temp dir; stage the model
% under HOME when dispatching through Docker so the volume mount is visible.
dirpath = lineTempName('lqns', ~isempty(dockerImage));
filename = [dirpath,filesep,'model.lqnx'];
self.model.writeXML(filename);
 
%self.runAnalyzerChecks(options);
Solver.resetRandomGeneratorSeed(options.seed);
 
if options.verbose
    %verbose = '-v';
    verbose = '';
else
    verbose = '-a -w';
end
 
multiserver_praqma = '';
switch options.method
    case 'lqsim'
        %no-op
    otherwise
        switch options.config.multiserver
            case 'conway'
                multiserver_praqma='-Pmultiserver=conway';
            case 'rolia'
                multiserver_praqma='-Pmultiserver=rolia';
            case 'zhou'
                multiserver_praqma='-Pmultiserver=zhou';
            case 'suri'
                multiserver_praqma='-Pmultiserver=suri';
            case 'reiser'
                multiserver_praqma='-Pmultiserver=reiser';
            case 'schmidt'
                multiserver_praqma='-Pmultiserver=schmidt';
            case 'default'
                multiserver_praqma='-Pmultiserver=rolia';
        end
end
 
if isunix
    %                 switch options.method
    %                     case {'default','lqns'}
    %                         cmd=['lqns ',verbose,' ',multiserver_praqma,' --iteration-limit=',num2str(options.iter_max),' -Pstop-on-message-loss=false -x ',filename]);
    %                     case {'srvn'}
    %                         cmd=['lqns ',verbose,' ',multiserver_praqma,' --iteration-limit=',num2str(options.iter_max),' -Playering=srvn -Pstop-on-message-loss=false -x ',filename]);
    %                     case {'exactmva'}
    %                         cmd=['lqns ',verbose,' ',multiserver_praqma,' --iteration-limit=',num2str(options.iter_max),' -Pmva=exact -Pstop-on-message-loss=false -x ',filename]);
    %                     case {'srvnexact'}
    %                         cmd=['lqns ',verbose,' ',multiserver_praqma,' --iteration-limit=',num2str(options.iter_max),' -Playering=srvn -Pmva=exact -Pstop-on-message-loss=false -x ',filename]);
    %                     case {'sim','lqsim'}
    %                         cmd=['lqsim ',verbose,' ',multiserver_praqma,' -A ',num2str(options.samples),',3  -Pstop-on-message-loss=false -x ',filename]);
    %                     case {'lqnsdefault'}
    %                         cmd=['lqns ',verbose,' ',multiserver_praqma,' -x ',filename]);
    %                     otherwise
    %                         cmd=['lqns ',verbose,' ',multiserver_praqma,' --iteration-limit=',num2str(options.iter_max),' -Pstop-on-message-loss=false -x ',filename]);
    %                 end
 
    % --iteration-limit seems faulty as of 6.2.27
    if options.verbose
        switch options.method
            case {'default','lqns'}
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -Pstop-on-message-loss=false -x ',filename];
            case {'srvn'}
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -Playering=srvn -Pstop-on-message-loss=false -x ',filename];
            case {'exactmva'}
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -Pmva=exact -Pstop-on-message-loss=false -x ',filename];
            case {'srvn.exactmva'}
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -Playering=srvn -Pmva=exact -Pstop-on-message-loss=false -x ',filename];
            case {'sim','lqsim'}
                cmd=['lqsim ',verbose,' ',multiserver_praqma,' -A ',num2str(options.samples),',3  -Pstop-on-message-loss=false -x ',filename];
            case {'lqnsdefault'}
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -x ',filename];
            otherwise
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -Pstop-on-message-loss=false -x ',filename];
        end
    else
        switch options.method
            case {'default','lqns'}
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -Pstop-on-message-loss=false -x ',filename,' 2>&1'];
            case {'srvn'}
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -Playering=srvn -Pstop-on-message-loss=false -x ',filename,' 2>&1'];
            case {'exactmva'}
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -Pmva=exact -Pstop-on-message-loss=false -x ',filename,' 2>&1'];
            case {'srvn.exactmva'}
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -Playering=srvn -Pmva=exact -Pstop-on-message-loss=false -x ',filename,' 2>&1'];
            case {'sim','lqsim'}
                cmd=['lqsim ',verbose,' ',multiserver_praqma,' -A ',num2str(options.samples),' -Pstop-on-message-loss=false -x ',filename,' 2>&1'];
            case {'lqnsdefault'}
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -x ',filename,' 2>&1'];
            otherwise
                cmd=['lqns ',verbose,' ',multiserver_praqma,' -Pstop-on-message-loss=false -x ',filename,' 2>&1'];
        end
    end
else
    switch options.method
        %         case {'default','lqns'}
        %             cmd=['lqns ',verbose,' ',multiserver_praqma,' -i ',num2str(options.iter_max),' -Pstop-on-message-loss=false -x ',filename]);
        %         case {'srvn'}
        %             cmd=['lqns ',verbose,' ',multiserver_praqma,' -i ',num2str(options.iter_max),' -Playering=srvn -Pstop-on-message-loss=false -x ',filename]);
        %         case {'exactmva'}
        %             cmd=['lqns ',verbose,' ',multiserver_praqma,' -i ',num2str(options.iter_max),' -Pmva=exact -Pstop-on-message-loss=false -x ',filename]);
        %         case {'srvnexact'}
        %             cmd=['lqns ',verbose,' ',multiserver_praqma,' -i ',num2str(options.iter_max),' -Playering=srvn -Pmva=exact -Pstop-on-message-loss=false -x ',filename]);
        %         case {'sim','lqsim'}
        %             cmd=['lqsim ',verbose,' ',multiserver_praqma,' -A ',num2str(options.samples),',3  -Pstop-on-message-loss=false -x ',filename]);
        %         case {'lqnsdefault'}
        %             cmd=['lqns ',verbose,' ',multiserver_praqma,' -x ',filename]);
        %         otherwise
        %             cmd=['lqns ',verbose,' ',multiserver_praqma,' -i ',num2str(options.iter_max),' -Pstop-on-message-loss=false -x ',filename]);
        case {'default','lqns'}
            cmd=['lqns ',verbose,' ',multiserver_praqma,' -Pstop-on-message-loss=false -x ',filename];
        case {'srvn'}
            cmd=['lqns ',verbose,' ',multiserver_praqma,' -Playering=srvn -Pstop-on-message-loss=false -x ',filename];
        case {'exactmva'}
            cmd=['lqns ',verbose,' ',multiserver_praqma,' -Pmva=exact -Pstop-on-message-loss=false -x ',filename];
        case {'srvn.exactmva'}
            cmd=['lqns ',verbose,' ',multiserver_praqma,' -Playering=srvn -Pmva=exact -Pstop-on-message-loss=false -x ',filename];
        case {'sim','lqsim'}
            %cmd=['lqsim ',verbose,' ',multiserver_praqma,' -Pstop-on-message-loss=false -x ',filename];
            cmd=['lqsim ',verbose,' ',multiserver_praqma,' -A ',num2str(options.samples),' -Pstop-on-message-loss=false -x ',filename];
        case {'lqnsdefault'}
            cmd=['lqns ',verbose,' ',multiserver_praqma,' -x ',filename];
        otherwise
            cmd=['lqns ',verbose,' ',multiserver_praqma,' -Pstop-on-message-loss=false -x ',filename];
    end
end
if ~isempty(dockerImage)
    % Rewrite the local invocation to run inside the container. The temp dir
    % is bind-mounted at the same absolute path so the model/result file paths
    % embedded in cmd resolve unchanged; --user matches the host uid/gid so the
    % .lqxo output stays host-owned and rmdir cleanup succeeds. The image's
    % 'lqns' is a REST wrapper, so call the real 'lqns.original' binary.
    [~, hostUid] = unix('id -u'); hostUid = strtrim(hostUid);
    [~, hostGid] = unix('id -g'); hostGid = strtrim(hostGid);
    dockerPrefix = ['docker run --rm --user ', hostUid, ':', hostGid, ...
        ' -v ', dirpath, ':', dirpath, ' -w ', dirpath, ' ', dockerImage, ' '];
    if strncmp(cmd, 'lqns ', 5)
        cmd = [dockerPrefix, 'lqns.original ', cmd(6:end)];
    elseif strncmp(cmd, 'lqsim ', 6)
        cmd = [dockerPrefix, 'lqsim ', cmd(7:end)];
    end
end
 
if options.verbose
%    line_printf('\nLQNS model: %s',filename);
    line_printf('\nLQNS command: %s\n',cmd);
end
 
% Check for remote execution
if useRemote
    line_debug(options, 'LQNS: using remote execution at %s', options.config.remote_url);
    if options.verbose
        line_printf('\nUsing remote LQNS at: %s\n', options.config.remote_url);
    end
    self.runRemoteLQNS(filename, options);
else
    line_debug(options, 'LQNS: using local execution, command: %s', cmd);
    % Local execution.
    %
    % MATLAB prepends its bundled runtime libraries to LD_LIBRARY_PATH before
    % spawning child processes. Its libstdc++ (e.g. GLIBCXX_3.4.28 in R2024a)
    % is older than the one lqns/lqsim require (GLIBCXX_3.4.29+), so the
    % binaries abort with a "version `GLIBCXX_...' not found" error before
    % writing their .lqxo output. That later surfaces as a FileNotFoundException
    % in parseXMLResults. Strip LD_LIBRARY_PATH for the child so it resolves the
    % system libstdc++ instead of MATLAB's.
    if isunix && isempty(dockerImage)
        cmd = ['env -u LD_LIBRARY_PATH ', cmd];
    end
    system(cmd);
end
self.parseXMLResults(filename);
 
if ~options.keep
    rmdir(dirpath,'s');
end
runtime = toc;
end