/usr/share/octave/packages/secs2d-0.0.8/METLINES/METLINEScapcomp.m is in octave-secs2d 0.0.8-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 | function [C,L,Lii,LiI,LII,Dnodes,varnodes]=...
METLINEScapcomp(imesh,epsilon,contacts)
##
##
## [C,L,Lii,LiI,LII,Dnodes,varnodes]=METLINEScapcomp(imesh,epsilon,contacts)
##
##
Ncontacts = length(contacts);
Nnodes = columns(imesh.p);
varnodes = [1:Nnodes];
for ii=1:Ncontacts
Dnodes{ii}=Unodesonside(imesh,contacts{ii});
varnodes = setdiff(varnodes,Dnodes{ii});
end
L = Ucomplap (imesh,epsilon);
for ii=1:Ncontacts
Lii{ii} = L(Dnodes{ii},Dnodes{ii});
LiI{ii} = L(Dnodes{ii},varnodes);
end
LII = L(varnodes,varnodes);
for ii=1:Ncontacts
for jj=ii:Ncontacts
if ii==jj
C(ii,jj)=sum(sum(Lii{ii}-LiI{ii}*(LII\(LiI{ii})')));
else
C(ii,jj)=sum(sum(-LiI{ii}*(LII\(LiI{jj})')));
end
end
end
for ii=1:4
for jj=1:ii-1
C(ii,jj)=C(jj,ii);
end
end
|