From c0c3470dff407cfe82030720103d21fbdd94486e Mon Sep 17 00:00:00 2001
From: Holger Vogt <holger.vogt@uni-due.de>
Date: Mon, 2 Apr 2018 15:48:51 +0200
Subject: [PATCH] add the simple vdmos capacitance model instead of Meyer's
 model

---
 src/include/ngspice/devdefs.h          |  2 ++
 src/spicelib/devices/devsup.c          | 17 +++++++++++++
 src/spicelib/devices/vdmos/vdmosload.c | 35 +++++++++++---------------
 3 files changed, 34 insertions(+), 20 deletions(-)

diff --git a/src/include/ngspice/devdefs.h b/src/include/ngspice/devdefs.h
index 0e338e348..b4be7a2b6 100644
--- a/src/include/ngspice/devdefs.h
+++ b/src/include/ngspice/devdefs.h
@@ -20,6 +20,8 @@ void DEVcmeyer(double,double,double,double,double,double,double,double,double,
         double,double,double*,double*,double*,double,double,double,double);
 void DEVqmeyer(double,double,double,double,double,double*,double*,double*,
         double,double);
+void DevCapVDMOS(double, double, double, double, double,
+                 double*, double*, double*);
 double DEVpred(CKTcircuit*,int);
 
 /* Cider integration */
diff --git a/src/spicelib/devices/devsup.c b/src/spicelib/devices/devsup.c
index eb57ca868..ac2737707 100644
--- a/src/spicelib/devices/devsup.c
+++ b/src/spicelib/devices/devsup.c
@@ -614,6 +614,23 @@ DEVcmeyer(double vgs0,		/* initial voltage gate-source */
     *cgb = *cgb *.5 + covlgb;
 }
 
+/* model according to
+http://ltwiki.org/index.php5?title=Undocumented_LTspice#VDMOS:_Breakdown_and_Sub-threshold_Enhancements
+*/
+void
+DevCapVDMOS(double vgd, double cgdmin,
+            double cgdmax, double a, double cgs,
+            double *capgs, double *capgd, double *capgb)
+{
+    double s = (cgdmax - cgdmin) / (1 + M_PI / 2);
+    double y = cgdmax - s;
+    if (vgd > 0)
+        *capgd = s * tanh(a * vgd) + y;
+    else
+        *capgd = s * atan(a * vgd) + y;
+    *capgs = cgs;
+    *capgb = 0;
+}
 
 /* Compute the MOS overlap capacitances as functions of the device
  * terminal voltages 
diff --git a/src/spicelib/devices/vdmos/vdmosload.c b/src/spicelib/devices/vdmos/vdmosload.c
index f56a62b27..779a9a48f 100644
--- a/src/spicelib/devices/vdmos/vdmosload.c
+++ b/src/spicelib/devices/vdmos/vdmosload.c
@@ -28,7 +28,6 @@ VDMOSload(GENmodel *inModel, CKTcircuit *ckt)
     double GateBulkOverlapCap;
     double GateDrainOverlapCap;
     double GateSourceOverlapCap;
-    double OxideCap;
     double SourceSatCur;
     double arg;
     double cbhat;
@@ -83,6 +82,11 @@ VDMOSload(GENmodel *inModel, CKTcircuit *ckt)
 
     /*  loop through all the VDMOS device models */
     for( ; model != NULL; model = VDMOSnextModel(model)) {
+        /* VDMOS capacitance parameters */
+        const double cgdmin = model->VDMOScgdmin;
+        const double cgdmax = model->VDMOScgdmax;
+        const double a = model->VDMOSa;
+        const double cgs = model->VDMOScgs;
 
         /* loop through all the instances of the model */
         for (here = VDMOSinstances(model); here != NULL ;
@@ -120,8 +124,6 @@ VDMOSload(GENmodel *inModel, CKTcircuit *ckt)
                     here->VDMOSm * EffectiveLength;
             Beta = here->VDMOStTransconductance * here->VDMOSm *
                     here->VDMOSw/EffectiveLength;
-            OxideCap = model->VDMOSoxideCapFactor * EffectiveLength * 
-                    here->VDMOSm * here->VDMOSw;
            
             /* 
              * ok - now to do the start-up operations
@@ -652,32 +654,25 @@ VDMOSload(GENmodel *inModel, CKTcircuit *ckt)
 */
 
             /*
-             *     meyer's capacitor model
+             * vdmos capacitor model
              */
             if ( ckt->CKTmode & (MODETRAN | MODETRANOP | MODEINITSMSIG) ) {
                 /*
-                 *     calculate meyer's capacitors
+                 * calculate gate - drain, gate - source capacitors
+                 * drain-source capacitor is evaluated with the bulk diode below
                  */
-                /* 
-                 * new cmeyer - this just evaluates at the current time,
+                /*
+                 * this just evaluates at the current time,
                  * expects you to remember values from previous time
                  * returns 1/2 of non-constant portion of capacitance
                  * you must add in the other half from previous time
                  * and the constant part
                  */
-                if (here->VDMOSmode > 0){
-                    DEVqmeyer (vgs,vgd,vgb,von,vdsat,
-			       (ckt->CKTstate0 + here->VDMOScapgs),
-			       (ckt->CKTstate0 + here->VDMOScapgd),
-			       (ckt->CKTstate0 + here->VDMOScapgb),
-			       here->VDMOStPhi,OxideCap);
-                } else {
-                    DEVqmeyer (vgd,vgs,vgb,von,vdsat,
-			       (ckt->CKTstate0 + here->VDMOScapgd),
-			       (ckt->CKTstate0 + here->VDMOScapgs),
-			       (ckt->CKTstate0 + here->VDMOScapgb),
-			       here->VDMOStPhi,OxideCap);
-                }
+                DevCapVDMOS(vgd, cgdmin, cgdmax, a, cgs,
+                            (ckt->CKTstate0 + here->VDMOScapgs),
+                            (ckt->CKTstate0 + here->VDMOScapgd),
+                            (ckt->CKTstate0 + here->VDMOScapgb));
+
                 vgs1 = *(ckt->CKTstate1 + here->VDMOSvgs);
                 vgd1 = vgs1 - *(ckt->CKTstate1 + here->VDMOSvds);
                 vgb1 = vgs1 - *(ckt->CKTstate1 + here->VDMOSvbs);