line source
% GETMATRIXSELECTIONS allows to select data in a block-defined matrix
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% DESCRIPTION: GETMATRIXSELECTIONS allows to select data in a
% block-defined matrix with user i/o selection, and i/o data selection.
% Permuations matrices are returned
%
% CALL: [ao_out pl_out] = getModelMatrixSelection(sys, plist_inputs)
%
% INPUTS:
%
% OUTPUTS:
%
% <a href="matlab:utils.helper.displayMethodInfo('ssm', 'getModelMatrixSelection')">Parameters Description</a>
%
% VERSION: $Id: getMatrixSelection.m,v 1.5 2011/04/08 08:56:23 hewitson Exp $
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function varargout = getMatrixSelection(blockMat, oldColumns, newColumns, oldRows, newRows)
error('This function is deprecated, it will be modified or removed. Use double instead.')
utils.helper.msg(utils.const.msg.MNAME, ['running ', mfilename]);
%% making indexes out of i/o data
indexInputsIn = makePortIndex( oldColumns, newColumns );
colSizesIn = getsizes(indexInputsIn);
colSizes = oldColumns.Nports;
indexOutputsOut = makePortIndex( oldRows, newRows );
rowSizesOut = getsizes(indexOutputsOut);
rowSizes = oldRows.Nports;
%% pruning ssmBlocks
inputNamesFull = getcellstr(oldColumns, colSizes);
inputsIn = blocksPrune(oldColumns, indexInputsIn);
inputNames = getcellstr(inputsIn, colSizesIn);
inputUnits = getUnit(inputsIn, colSizesIn);
outputNamesFull = getcellstr(oldRows, rowSizes);
outputsOut = blocksPrune(oldRows, indexOutputsOut);
outputNames = getcellstr(outputsOut, rowSizesOut);
outputUnits = getUnit(outputsOut, rowSizesOut);
%% pruning matrix and turning it to a double
A = ssm.blockMatFusion( blockMat, rowSizes, colSizes);
%% finding input and output variable positions of data in option plists
if ischar(newColumns) && strcmp(newColumns, 'ALL') %case 'ALL' 'NONE'
inputsConv = eye(sum(colSizes),sum(colSizes));
elseif ischar(newColumns) && strcmp(newColumns, 'NONE') %case 'NONE'
inputsConv = eye(sum(colSizes),0);
else
if ischar(newColumns)
newColumns = {newColumns};
end
NinputsIn = numel(newColumns);
Ninputs = sum(colSizes);
inputsConv = zeros(Ninputs,NinputsIn);
for i=1:NinputsIn
position = reshape( strcmpi(inputNamesFull, newColumns{i}) , [Ninputs, 1]);
inputsConv(:,i) = double(position).';
if sum(position)==0
error(['No variable found named ', newColumns{i}]);
end
end
end
%% finding input and output variable positions of data in option plists
if ischar(newRows) && strcmp(newRows, 'ALL') %case 'ALL' 'NONE'
outputsConv = eye(sum(rowSizes),sum(rowSizes));
elseif ischar(newRows) && strcmp(newRows, 'NONE') %case 'NONE'
outputsConv = eye(0,sum(rowSizes));
else
if ischar(newRows)
newRows = {newRows};
end
NoutputsOut = numel(newRows);
Noutputs = sum(rowSizes);
outputsConv = zeros(NoutputsOut,Noutputs);
for i=1:NoutputsOut
position = reshape( strcmpi(newRows{i}, outputNamesFull ) , [1, Noutputs]);
outputsConv(i,:) = double(position);
if sum(position)==0
error(['No variable found named ', newRows{i}]);
end
end
end
%% parsing output matrices
varargout = {A, inputsConv, outputsConv, inputNames, outputNames, inputUnits, outputUnits};
end
function sizes = getsizes(idx)
sizes = zeros(1, numel(idx));
for i=1:numel(idx)
sizes(i) = sum(idx{i});
end
end
function names = getcellstr(block, sizes)
names = cell(1, sum(sizes));
k = 1;
for i=1:numel(block)
for j=1:sizes(i)
names{k} = block(i).ports(j).name;
k = k+1;
end
end
end
function units = getUnit(block, sizes)
units = unit.initObjectWithSize(1, sum(sizes));
k = 1;
for i=1:numel(block)
for j=1:sizes(i)
units(k) = block(i).ports(j).units;
k = k+1;
end
end
end