/usr/share/octave/packages/secs2d-0.0.8/DDGOX/DDGOXddcurrent.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 42 43 44 45 46 47 48 49 50 | function [current,divrg]=DDGOXddcurrent(mesh,Sinodes,data,contacts);
% [current,divrg]=DDGOXddcurrent(mesh,Sinodes,data,contacts);
% This file is part of
%
% SECS2D - A 2-D Drift--Diffusion Semiconductor Device Simulator
% -------------------------------------------------------------------
% Copyright (C) 2004-2006 Carlo de Falco
%
%
%
% SECS2D 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 2 of the License, or
% (at your option) any later version.
%
% SECS2D 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 SECS2D; If not, see <http://www.gnu.org/licenses/>.
load (file_in_path(path,'constants.mat'))
Nelements = size(mesh.t,2);
mob = Ufielddepmob(mesh,data.un,data.Fn,data.vsatn,data.mubn);
An = Uscharfettergummel(mesh,data.V(Sinodes),mob);
mob = Ufielddepmob(mesh,data.up,data.Fp,data.vsatp,data.mubp);
Ap = Uscharfettergummel(mesh,-data.V(Sinodes),mob);
divrg = An * data.n + Ap * data.p;
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;
|