/usr/share/octave/packages/interval-3.1.0/@infsup/triu.m is in octave-interval 3.1.0-5.
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 | ## Copyright 2016 Oliver Heimlich
##
## This program 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.
##
## This program 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 this program; if not, see <http://www.gnu.org/licenses/>.
## -*- texinfo -*-
## @documentencoding UTF-8
## @defmethod {@@infsup} triu (@var{A})
## @defmethodx {@@infsup} triu (@var{A}, @var{k})
## @defmethodx {@@infsup} triu (@var{A}, @var{k}, "pack")
## Return a new matrix formed by extracting the upper triangular part of the
## matrix @var{A}, and setting all other elements to zero.
##
## The second argument is optional, and specifies how many diagonals above or
## below the main diagonal should also be set to zero.
##
## If the option @option{pack} is given as third argument, the extracted
## elements are not inserted into a matrix, but rather stacked column-wise one
## above other.
##
## @seealso{@@infsup/tril, @@infsup/diag}
## @end defmethod
## Author: Oliver Heimlich
## Keywords: interval
## Created: 2016-08-04
function A = triu (A, varargin)
if (nargin >= 2 && isa (varargin{1}, 'infsup'))
error ('triu: invalid second argument; it must not be an interval');
endif
if (nargin >= 3 && isa (varargin{2}, 'infsup'))
error ('triu: invalid third argument; it must not be an interval');
endif
A.inf = triu (A.inf, varargin{:});
A.sup = triu (A.sup, varargin{:});
A.inf(A.inf == 0) = -0;
endfunction
%!assert (triu (infsup (magic (10))) == triu (magic (10)));
%!assert (triu (infsup (magic (10)), 1) == triu (magic (10), 1));
%!assert (triu (infsup (magic (10)), -1) == triu (magic (10), -1));
%!assert (triu (infsup (magic (10)), 0, "pack") == triu (magic (10), 0, "pack"));
|