/usr/share/octave/packages/control-2.6.2/@lti/__lti_prune__.m is in octave-control 2.6.2-1build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 | ## Copyright (C) 2009-2014 Lukas F. Reichlin
##
## This file is part of LTI Syncope.
##
## LTI Syncope is free software: you can redistribute it and/or modify
## it under the terms of the GNU General Public License as published by
## the Free Software Foundation, either version 3 of the License, or
## (at your option) any later version.
##
## LTI Syncope is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
## GNU General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with LTI Syncope. If not, see <http://www.gnu.org/licenses/>.
## -*- texinfo -*-
## Submodel extraction and reordering for LTI objects.
## This file is part of the Model Abstraction Layer.
## For internal use only.
## Author: Lukas Reichlin <lukas.reichlin@gmail.com>
## Created: September 2009
## Version: 0.2
function [lti, out_idx, in_idx] = __lti_prune__ (lti, out_idx, in_idx)
if (ischar (out_idx) && ! strcmp (out_idx, ":")) # sys("grp", :)
out_idx = {out_idx};
endif
if (ischar (in_idx) && ! strcmp (in_idx, ":")) # sys(:, "grp")
in_idx = {in_idx};
endif
if (iscell (out_idx)) # sys({"grp1", "grp2"}, :)
tmp = cellfun (@(x) __str2idx__ (lti.outgroup, lti.outname, x, "out"), out_idx, "uniformoutput", false);
out_idx = vertcat (tmp{:});
endif
if (iscell (in_idx)) # sys(:, {"grp1", "grp2"})
tmp = cellfun (@(x) __str2idx__ (lti.ingroup, lti.inname, x, "in"), in_idx, "uniformoutput", false);
in_idx = vertcat (tmp{:});
endif
if (nfields (lti.outgroup))
p = numel (lti.outname); # get size before pruning outnames!
[lti.outgroup, empty] = structfun (@(x) __group_prune__ (x, out_idx, p), lti.outgroup, "uniformoutput", false);
empty = cell2mat (struct2cell (empty));
fields = fieldnames (lti.outgroup);
lti.outgroup = rmfield (lti.outgroup, fields(empty));
endif
if (nfields (lti.ingroup))
m = numel (lti.inname);
[lti.ingroup, empty] = structfun (@(x) __group_prune__ (x, in_idx, m), lti.ingroup, "uniformoutput", false);
empty = cell2mat (struct2cell (empty));
fields = fieldnames (lti.ingroup);
lti.ingroup = rmfield (lti.ingroup, fields(empty));
endif
lti.outname = lti.outname(out_idx);
lti.inname = lti.inname(in_idx);
endfunction
function [group, empty] = __group_prune__ (group, idx, n)
lg = length (group);
group = sparse (group, 1:lg, 1, n, lg);
group = group(idx, :);
[group, ~] = find (group);
empty = isempty (group);
endfunction
|