src/solver/qddm1e/jflux1q.h

   1 /*****************************************************************************/
   2 /*              8888888         88888888         88888888                    */
   3 /*            8                8                8                            */
   4 /*           8                 8                8                            */
   5 /*           8                  88888888         88888888                    */
   6 /*           8      8888                8                8                   */
   7 /*            8       8                 8                8                   */
   8 /*              888888         888888888        888888888                    */
   9 /*                                                                           */
  10 /*       A Two-Dimensional General Purpose Semiconductor Simulator.          */
  11 /*                                                                           */
  12 /*  GSS 0.4x                                                                 */
  13 /*  Last update: July 18, 2007                                               */
  14 /*                                                                           */
  15 /*  Gong Ding                                                                */
  16 /*  gdiso@ustc.edu                                                           */
  17 /*  NINT, No.69 P.O.Box, Xi'an City, China                                   */
  18 /*                                                                           */
  19 /*****************************************************************************/
  20
  21 #ifndef _jflux1q_h_
  22 #define _jflux1q_h_
  23 #include "petsc.h"
  24 #include "mathfunc.h"
  25
  26
  27 inline PetscScalar nmid(PetscScalar Vt, PetscScalar Vc1, PetscScalar Vc2, PetscScalar n1, PetscScalar n2)
  28 {
  29   return n1*aux2((Vc2-Vc1)/(2*Vt)) + n2*aux2((Vc1-Vc2)/(2*Vt));
  30 }
  31
  32 inline PetscScalar pmid(PetscScalar Vt, PetscScalar Vv1, PetscScalar Vv2, PetscScalar p1, PetscScalar p2)
  33 {
  34   return p1*aux2((Vv1-Vv2)/(2*Vt)) + p2*aux2((Vv2-Vv1)/(2*Vt));
  35 }
  36
  37 inline AutoDScalar nmid(PetscScalar Vt, const AutoDScalar &Vc1, const AutoDScalar &Vc2, const AutoDScalar &n1, const AutoDScalar &n2)
  38 {
  39   return n1*aux2((Vc2-Vc1)/(2*Vt)) + n2*aux2((Vc1-Vc2)/(2*Vt));
  40 }
  41
  42 inline AutoDScalar pmid(PetscScalar Vt, const AutoDScalar &Vv1, const AutoDScalar &Vv2, const AutoDScalar &p1, const AutoDScalar &p2)
  43 {
  44   return p1*aux2((Vv1-Vv2)/(2*Vt)) + p2*aux2((Vv2-Vv1)/(2*Vt));
  45 }
  46
  47
  48 //-----------------------------------------------------------------------------
  49 inline PetscScalar In(PetscScalar Vt,PetscScalar Vc1,PetscScalar Vc2,PetscScalar n1,PetscScalar n2,PetscScalar h)
  50 {
  51   PetscScalar E    = (Vc2-Vc1)/h;
  52   PetscScalar n    = n1*aux2((Vc2-Vc1)/(2*Vt)) + n2*aux2((Vc1-Vc2)/(2*Vt));
  53   PetscScalar dndx = aux1((Vc2-Vc1)/(2*Vt))*(n2-n1)/h;
  54   return (E*n+Vt*dndx);
  55 }
  56
  57 inline PetscScalar Ip(PetscScalar Vt,PetscScalar Vv1,PetscScalar Vv2,PetscScalar p1,PetscScalar p2,PetscScalar h)
  58 {
  59   PetscScalar E    = (Vv2-Vv1)/h;
  60   PetscScalar p    = p1*aux2((Vv1-Vv2)/(2*Vt)) + p2*aux2((Vv2-Vv1)/(2*Vt));
  61   PetscScalar dpdx = aux1((Vv2-Vv1)/(2*Vt))*(p2-p1)/h;
  62   return (E*p-Vt*dpdx);
  63 }
  64
  65 inline AutoDScalar In(PetscScalar Vt, const AutoDScalar &Vc1, const AutoDScalar &Vc2, const AutoDScalar &n1, const AutoDScalar &n2, PetscScalar h)
  66 {
  67   AutoDScalar E    = (Vc2-Vc1)/h;
  68   AutoDScalar n    = n1*aux2((Vc2-Vc1)/(2*Vt)) + n2*aux2((Vc1-Vc2)/(2*Vt));
  69   AutoDScalar dndx = aux1((Vc2-Vc1)/(2*Vt))*(n2-n1)/h;
  70   return (E*n+Vt*dndx);
  71 }
  72
  73 inline AutoDScalar Ip(PetscScalar Vt, const AutoDScalar &Vv1, const AutoDScalar &Vv2, const AutoDScalar &p1, const AutoDScalar &p2, PetscScalar h)
  74 {
  75   AutoDScalar E    = (Vv2-Vv1)/h;
  76   AutoDScalar p    = p1*aux2((Vv1-Vv2)/(2*Vt)) + p2*aux2((Vv2-Vv1)/(2*Vt));
  77   AutoDScalar dpdx = aux1((Vv2-Vv1)/(2*Vt))*(p2-p1)/h;
  78   return (E*p-Vt*dpdx);
  79 }
  80
  81
  82
  83
  84
  85 //-----------------------------------------------------------------------------
  86 inline PetscScalar In(PetscScalar Vt,PetscScalar dVc,PetscScalar n1,PetscScalar n2,PetscScalar h)
  87 {
  88   return Vt/h*(n2*bern(-dVc/Vt)-n1*bern(dVc/Vt));
  89 }
  90
  91 inline PetscScalar Ip(PetscScalar Vt,PetscScalar dVv,PetscScalar p1,PetscScalar p2,PetscScalar h)
  92 {
  93   return Vt/h*(p1*bern(-dVv/Vt)-p2*bern(dVv/Vt));
  94 }
  95
  96 inline AutoDScalar In(PetscScalar Vt,const AutoDScalar &dVc,const AutoDScalar &n1,const AutoDScalar &n2, PetscScalar h)
  97 {
  98   return Vt/h*(n2*bern(-dVc/Vt)-n1*bern(dVc/Vt));
  99 }
 100
 101 inline AutoDScalar Ip(PetscScalar Vt,const AutoDScalar &dVv,const AutoDScalar &p1,const AutoDScalar &p2, PetscScalar h)
 102 {
 103   return Vt/h*(p1*bern(-dVv/Vt)-p2*bern(dVv/Vt));
 104 }
 105
 106
 107
 108
 109 //-----------------------------------------------------------------------------
 110 inline PetscScalar qV(PetscScalar q1,PetscScalar q2,PetscScalar h)
 111 {
 112   PetscScalar x = 0.5*(q1-q2);
 113   if(x<0)
 114    return (1-exp(x))/h;
 115   else
 116    return (-x-x*x/2)/h;
 117 }
 118
 119 inline AutoDScalar qV(const AutoDScalar &q1, const AutoDScalar &q2, PetscScalar h)
 120 {
 121   AutoDScalar x = 0.5*(q1-q2);
 122   if(x<0)
 123    return (1-exp(x))/h;
 124   else
 125    return (-x-x*x/2)/h;
 126 }
 127
 128
 129 #endif