sampleSysAggr.m

function tranSysState = sampleSysAggr(self, numEvents)
% TRANSYSSTATE = SAMPLESYSAGGR(NUMSAMPLES)
options = self.getOptions;
options.force = true;
if isempty(self.result) || ~isfield(self.result,'infGen')
    runAnalyzer(self);
end
[infGen, eventFilt] = getGenerator(self);
[stateSpace, localStateSpace] = getStateSpace(self);
stateSpaceAggr = getStateSpaceAggr(self);
 
sn = self.getStruct;
initState = sn.state;
% Pad each stateful node's initial state (compact representation) with
% leading zeros up to the enumerated per-node state-space width, so that
% s0 lines up column-wise with the joint stateSpace. Phase-type stations
% (e.g. Erlang/Coxian service) enumerate phase-augmented rows that are much
% wider than the compact initial state, so an unpadded concatenation would
% not match any stateSpace row. Mirrors solver_ctmc_transient_analyzer.
spaceWidths = zeros(1,length(initState));
s0 = [];
for isf=1:length(initState)
    w = size(localStateSpace{isf},2);
    spaceWidths(isf) = w;
    s0 = [s0, zeros(1,w-length(initState{isf})), initState{isf}];
end
nst = cumsum([1,spaceWidths]);
 
% set initial state
pi0 = zeros(1,size(stateSpace,1));
pi0(matchrow(stateSpace,s0))=1;
 
% filter all CTMC events as a marked Markovian arrival process
D1 = cellsum(eventFilt);
D0 = infGen-D1;
MMAP = mmap_normalize([{D0},{D1},eventFilt(:)']);
 
% now sampel the MMAP
[sjt,event,~,~,sts] = mmap_sample(MMAP,numEvents, pi0);
 
sn = self.getStruct;
tranSysState = struct();
tranSysState.handle = self.model.getStatefulNodes';
tranSysState.t = cumsum([0,sjt(1:end-1)']');
for isf=1:length(initState)
    tranSysState.state{isf} = stateSpace(sts,(nst(isf):nst(isf+1)-1));
    [~,tranSysState.state{isf}] = State.toMarginal(sn,sn.statefulToNode(isf),tranSysState.state{isf});
end
 
tranSysState.event = {};
for e = 1:length(event)
    for a=1:length(sn.sync{event(e)}.active)
        tranSysState.event{end+1} = sn.sync{event(e)}.active{a};
        tranSysState.event{end}.t = tranSysState.t(e);
    end
    for p=1:length(sn.sync{event(e)}.passive)
        tranSysState.event{end+1} = sn.sync{event(e)}.passive{p};
        tranSysState.event{end}.t = tranSysState.t(e);
    end
end
tranSysState.isaggregate = true;
end