1 # Created by Octave 3.6.1, Thu Mar 22 16:18:37 2012 UTC <root@t61>
13 # name: <cell-element>
20 bp = B(x)=x/(exp(x)-1)
25 # name: <cell-element>
34 # name: <cell-element>
41 # name: <cell-element>
49 b = B(x)=x/(exp(x)-1) if sg==1
50 b = B(-x)=x+B(x) if sg==0
51 also works if x is a vector
55 # name: <cell-element>
63 b = B(x)=x/(exp(x)-1) if sg=
67 # name: <cell-element>
74 # name: <cell-element>
78 R = compconst (nodes,Nnodes,elements,Nelements,D,C);
82 # name: <cell-element>
86 R = compconst (nodes,Nnodes,elements,Nelements,D,C);
91 # name: <cell-element>
98 # name: <cell-element>
102 L = Ucomplap (nodes,Nnode,elements,Nelements,coeff)
103 Computes the P1 finite element approximation of the
104 differential operator - d ( coeff d (.)\dx)\dx
108 # name: <cell-element>
112 L = Ucomplap (nodes,Nnode,elements,Nelements,coeff)
113 Computes the P1 finite
117 # name: <cell-element>
124 # name: <cell-element>
128 Bmat = Ucompmass (nodes,Nnodes,elements,Nelements,Bvect,Cvect);
132 # name: <cell-element>
136 Bmat = Ucompmass (nodes,Nnodes,elements,Nelements,Bvect,Cvect);
141 # name: <cell-element>
148 # name: <cell-element>
152 A=Udriftdiffusion(x,psi,coeff)
154 Builds the Scharfetter-Gummel approximation
155 of the differential operator - (coeff (n' - n psi'))'
160 # name: <cell-element>
164 A=Udriftdiffusion(x,psi,coeff)
169 # name: <cell-element>
176 # name: <cell-element>
180 m = mediaarmonica(w,x);
181 returns the harmonic mean value of w in each of the intervals x_i , x_i+1
186 # name: <cell-element>
190 m = mediaarmonica(w,x);
191 returns the harmonic mean value of w in each of the in
195 # name: <cell-element>
202 # name: <cell-element>
206 A=Uscharfettergummel(nodes,Nnodes,elements,Nelements,acoeff,bcoeff,v)
208 Builds the Scharfetter-Gummel matrix for the
209 the discretization of the LHS
210 of the Drift-Diffusion equation:
212 $ -(a(x) (u' - b v'(x) u))'= f $
214 where a(x) is piecewise constant
215 and v(x) is piecewise linear, so that
216 v'(x) is still piecewise constant
217 b is a constant independent of x
223 # name: <cell-element>
227 A=Uscharfettergummel(nodes,Nnodes,elements,Nelements,acoeff,bcoeff,v)
232 # name: <cell-element>
239 # name: <cell-element>
245 SECS1D - A 1-D Drift--Diffusion Semiconductor Device Simulator
246 -------------------------------------------------------------------
247 Copyright (C) 2004-2007 Carlo de Falco
251 SECS1D is free software; you can redistribute it and/or modify
252 it under the terms of the GNU General Public License as published by
253 the Free Software Foundation; either version 2 of the License, or
254 (at your option) any later version.
256 SECS1D is distributed in the hope that it will be useful,
257 but WITHOUT ANY WARRANTY; without even the implied warranty of
258 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
259 GNU General Public License for more details.
261 You should have received a copy of the GNU General Public License
262 along with SECS1D; If not, see <http://www.gnu.org/licenses/>.
266 # name: <cell-element>