This file is indexed.

/usr/share/octave/packages/nurbs-1.3.10/nrbkntins.m is in octave-nurbs 1.3.10-1.

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
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
function inurbs = nrbkntins(nurbs,iknots)
% 
% NRBKNTINS: Insert a single or multiple knots into a NURBS curve,
%            surface or volume.
% 
% Calling Sequence:
% 
%   icrv = nrbkntins(crv,iuknots);
%   isrf = nrbkntins(srf,{iuknots ivknots});
%   ivol = nrbkntins(vol,{iuknots ivknots iwknots});
% 
% INPUT:
% 
%   crv		: NURBS curve, see nrbmak.
% 
%   srf		: NURBS surface, see nrbmak.
% 
%   srf		: NURBS volume, see nrbmak.
% 
%   iuknots	: Knots to be inserted along U direction.
% 
%   ivknots	: Knots to be inserted along V direction.
% 
%   iwknots	: Knots to be inserted along W direction.
% 
% OUTPUT:
% 
%   icrv	: new NURBS structure for a curve with knots inserted.
% 
%   isrf	: new NURBS structure for a surface with knots inserted.
% 
%   ivol	: new NURBS structure for a volume with knots inserted.
% 
% Description:
% 
%   Inserts knots into the NURBS data structure, these can be knots at
%   new positions or at the location of existing knots to increase the
%   multiplicity. Note that the knot multiplicity cannot be increased
%   beyond the order of the spline. Knots along the V direction can only
%   inserted into NURBS surfaces, not curve that are always defined along
%   the U direction. This function use the B-Spline function bspkntins,
%   which interfaces to an internal 'C' routine.
% 
% Examples:
% 
%   Insert two knots into a curve, one at 0.3 and another
%   twice at 0.4
%
%   icrv = nrbkntins(crv, [0.3 0.4 0.4])
% 
%   Insert into a surface two knots as (1) into the U knot
%   sequence and one knot into the V knot sequence at 0.5.
%
%   isrf = nrbkntins(srf, {[0.3 0.4 0.4] [0.5]})
% 
% See also:
% 
%   bspkntins
%
% Note:
%
%   No knot multiplicity will be increased beyond the order of the spline.
%
%    Copyright (C) 2000 Mark Spink, 2010 Rafael Vazquez
%
%    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/>.

if nargin < 2
  error('Input argument must include the NURBS and knots to be inserted');
end

if ~isstruct(nurbs)
  error('NURBS representation is not structure!');
end

if ~strcmp(nurbs.form,'B-NURBS')
  error('Not a recognised NURBS representation');
end

degree = nurbs.order-1;

if iscell(nurbs.knots)
 fmax = @(x,y) any (y > max(x)); fmin = @(x,y) any (y < min(x));
 if (any(cellfun(fmax, nurbs.knots, iknots)) || any(cellfun(fmin, nurbs.knots, iknots)))
   error ('Trying to insert a knot outside the interval of definition')
 end
    
 if size(nurbs.knots,2)==3
  % NURBS represents a volume
  num1 = nurbs.number(1);
  num2 = nurbs.number(2);
  num3 = nurbs.number(3);

  % Insert knots along the w direction
  if isempty(iknots{3})
    coefs = nurbs.coefs;
    knots{3} = nurbs.knots{3};
  else
    coefs = reshape(nurbs.coefs,4*num1*num2,num3);
    [coefs,knots{3}] = bspkntins(degree(3),coefs,nurbs.knots{3},iknots{3});
    num3 = size(coefs,2);
    coefs = reshape(coefs,[4 num1 num2 num3]);
  end

  % Insert knots along the v direction
  if isempty(iknots{2})
    knots{2} = nurbs.knots{2};
  else
    coefs = permute(coefs,[1 2 4 3]);
    coefs = reshape(coefs,4*num1*num3,num2);
    [coefs,knots{2}] = bspkntins(degree(2),coefs,nurbs.knots{2},iknots{2});
    num2 = size(coefs,2);
    coefs = reshape(coefs,[4 num1 num3 num2]);
    coefs = permute(coefs,[1 2 4 3]);
  end

  % Insert knots along the u direction
  if isempty(iknots{1})
    knots{1} = nurbs.knots{1};
  else   
    coefs = permute(coefs,[1 3 4 2]);
    coefs = reshape(coefs,4*num2*num3,num1);
    [coefs,knots{1}] = bspkntins(degree(1),coefs,nurbs.knots{1},iknots{1});
    coefs = reshape(coefs,[4 num2 num3 size(coefs,2)]);
    coefs = permute(coefs,[1 4 2 3]);
  end
     
 elseif size(nurbs.knots,2)==2
  % NURBS represents a surface
  num1 = nurbs.number(1);
  num2 = nurbs.number(2);

  % Insert knots along the v direction
  if isempty(iknots{2})
    coefs = nurbs.coefs;
    knots{2} = nurbs.knots{2};
  else
    coefs = reshape(nurbs.coefs,4*num1,num2);
    [coefs,knots{2}] = bspkntins(degree(2),coefs,nurbs.knots{2},iknots{2});
    num2 = size(coefs,2);
    coefs = reshape(coefs,[4 num1 num2]);
  end

  % Insert knots along the u direction
  if isempty(iknots{1})
    knots{1} = nurbs.knots{1};
  else   
    coefs = permute(coefs,[1 3 2]);
    coefs = reshape(coefs,4*num2,num1);
    [coefs,knots{1}] = bspkntins(degree(1),coefs,nurbs.knots{1},iknots{1});
    coefs = reshape(coefs,[4 num2 size(coefs,2)]);
    coefs = permute(coefs,[1 3 2]);
  end
 end
else

  if (any(iknots > max(nurbs.knots)) || any(iknots < min(nurbs.knots)))
    error ('Trying to insert a knot outside the interval of definition')
  end
  % NURBS represents a curve
  if isempty(iknots)
    coefs = nurbs.coefs;
    knots = nurbs.knots;
  else
    [coefs,knots] = bspkntins(degree,nurbs.coefs,nurbs.knots,iknots);  
  end

end

% construct new NURBS
inurbs = nrbmak(coefs,knots); 

end

%!demo
%! crv = nrbtestcrv;
%! plot(crv.coefs(1,:),crv.coefs(2,:),'bo')
%! title('Knot insertion along test curve: curve and control polygons.');
%! hold on;
%! plot(crv.coefs(1,:),crv.coefs(2,:),'b--');
%!
%! nrbplot(crv,48);
%!
%! icrv = nrbkntins(crv,[0.125 0.375 0.625 0.875] );
%! plot(icrv.coefs(1,:),icrv.coefs(2,:),'ro')
%! plot(icrv.coefs(1,:),icrv.coefs(2,:),'r--');
%! hold off