/usr/share/octave/packages/secs2d-0.0.8/QDDGOX/QDDGOXddcurrent.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 | function [current,divrg]=QDDGOXddcurrent(mesh,Sinodes,data,contacts,options);
% [current,divrg]=QDDGOXddcurrent(mesh,Sinodes,data,contacts);
constants
Nelements = size(mesh.t,2);
mob = Ufielddepmob(mesh,data.un,data.Fn,data.vsatn,data.mubn);
if (~isfield("options","FD")|(options.FD==0))
data.FDn=0;
data.FDp=0;
end
An = Uscharfettergummel(mesh,data.V(Sinodes)+data.G+data.FDn,mob);
if (options.holes==1)
mob = Ufielddepmob(mesh,data.up,data.Fp,data.vsatp,data.mubp);
Ap = Uscharfettergummel(mesh,-data.V(Sinodes)-data.Gp-data.FDp,mob);
divrg = An * data.n + Ap * data.p;
else
divrg = An * data.n;
end
for con = 1:length(contacts)
cedges = [];
cedges=[cedges,find(mesh.e(5,:)==contacts(con))];
cnodes = mesh.e(1:2,cedges);
cnodes = [cnodes(1,:) cnodes(2,:)];
cnodes = unique(cnodes);
current(con) = sum(divrg(cnodes));
end
Is = q*data.us*data.Vs*data.ns;
current = current * Is;
|