Mercurial > forge
changeset 11935:95a0358ce980 octave-forge
fixed bug in (and cleaned) secs1d_mobility_model_noscale.m
author | dvd7587 |
---|---|
date | Mon, 08 Jul 2013 14:12:17 +0000 |
parents | 7b91ea5611f4 |
children | c6a399c9a1b6 |
files | extra/secs1d/inst/secs1d_impact_ionization_noscale.m extra/secs1d/inst/secs1d_mobility_model_noscale.m |
diffstat | 2 files changed, 0 insertions(+), 3 deletions(-) [+] |
line wrap: on
line diff
--- a/extra/secs1d/inst/secs1d_impact_ionization_noscale.m Mon Jul 08 11:02:19 2013 +0000 +++ b/extra/secs1d/inst/secs1d_impact_ionization_noscale.m Mon Jul 08 14:12:17 2013 +0000 @@ -29,7 +29,6 @@ [Ex,Ey] = bim2c_pde_gradient (device.msh, Fp); Fava_p = sqrt(Ex .^ 2 + Ey .^ 2)(:); elseif (rows (device.msh.t) == 5) - %% error("no 3d gradient computation routine available") [Ex,Ey,Ez] = bim3c_pde_gradient (device.msh, Fn); Fava_n = sqrt(Ex .^ 2 + Ey .^ 2 + Ez .^ 2)(:); [Ex,Ey,Ez] = bim3c_pde_gradient (device.msh, Fp);
--- a/extra/secs1d/inst/secs1d_mobility_model_noscale.m Mon Jul 08 11:02:19 2013 +0000 +++ b/extra/secs1d/inst/secs1d_mobility_model_noscale.m Mon Jul 08 14:12:17 2013 +0000 @@ -31,7 +31,6 @@ end t = 300 * ones (size (n)); %[K] t_300 = 300; %[K] - %t_elem = .5 * t(2:end) + .5 * t(1:end-1); t_elem = sum(t(meshelems(1:end-1, :)),1)(:) / (rows (meshelems) - 1); if (carrier =='n') @@ -64,7 +63,6 @@ error ("Mobility models only defined for electons (carrier=\'n\') or holes (carrier=\'p\')") endif - %%muph = .5 * mu_ph_nodes(2:end) + .5 * mu_ph_nodes(1:end-1); muph = sum(mu_ph_nodes(meshelems(1:end-1, :)),1)(:) / (rows (meshelems) - 1); u = (muph*(alpha +1)) ./ ... (alpha + (1 + (((alpha + 1) * muph .* abs(E)) ./ vsat_m) .^ beta_m) .^