From ae78dd06ef762de83e1389c106d72b00a408e736 Mon Sep 17 00:00:00 2001
From: dwarning <dwarning>
Date: Tue, 25 Jun 2013 22:03:32 +0200
Subject: [PATCH] devices/mos*: prevent NaN in sqrt for negative value of model
 parameter `PHI'

---
 src/spicelib/devices/mos1/mos1temp.c | 6 ++++++
 src/spicelib/devices/mos2/mos2temp.c | 5 +++++
 src/spicelib/devices/mos3/mos3temp.c | 6 ++++++
 src/spicelib/devices/mos6/mos6temp.c | 6 ++++++
 src/spicelib/devices/mos9/mos9temp.c | 5 +++++
 5 files changed, 28 insertions(+)

diff --git a/src/spicelib/devices/mos1/mos1temp.c b/src/spicelib/devices/mos1/mos1temp.c
index 3764bc96a..dcf5853d6 100644
--- a/src/spicelib/devices/mos1/mos1temp.c
+++ b/src/spicelib/devices/mos1/mos1temp.c
@@ -52,6 +52,12 @@ MOS1temp(GENmodel *inModel, CKTcircuit *ckt)
 
     /* now model parameter preprocessing */
 
+        if (model->MOS1phi <= 0.0) {
+            SPfrontEnd->IFerror (ERR_FATAL,
+               "%s: Phi is not positive.", &model->MOS1modName);
+            return(E_BADPARM);
+        }
+
         if(!model->MOS1oxideThicknessGiven || model->MOS1oxideThickness == 0) {
             model->MOS1oxideCapFactor = 0;
         } else {
diff --git a/src/spicelib/devices/mos2/mos2temp.c b/src/spicelib/devices/mos2/mos2temp.c
index 982c784dd..4921af791 100644
--- a/src/spicelib/devices/mos2/mos2temp.c
+++ b/src/spicelib/devices/mos2/mos2temp.c
@@ -53,6 +53,11 @@ MOS2temp(GENmodel *inModel, CKTcircuit *ckt)
         arg1 = -egfet1/(kt1+kt1)+1.1150877/(CONSTboltz*(REFTEMP+REFTEMP));
         pbfact1 = -2*vtnom *(1.5*log(fact1)+CHARGE*arg1);
 
+        if (model->MOS2phi <= 0.0) {
+            SPfrontEnd->IFerror (ERR_FATAL,
+               "%s: Phi is not positive.", &model->MOS2modName);
+            return(E_BADPARM);
+        }
 
         if(!model->MOS2oxideThicknessGiven) {
             model->MOS2oxideThickness = 1e-7;
diff --git a/src/spicelib/devices/mos3/mos3temp.c b/src/spicelib/devices/mos3/mos3temp.c
index 7357b47ae..19d86d651 100644
--- a/src/spicelib/devices/mos3/mos3temp.c
+++ b/src/spicelib/devices/mos3/mos3temp.c
@@ -53,6 +53,12 @@ MOS3temp(GENmodel *inModel, CKTcircuit *ckt)
                                                  CONSTKoverQ);
         ni_temp=1.45e16*nifact;
 
+        if (model->MOS3phi <= 0.0) {
+            SPfrontEnd->IFerror (ERR_FATAL,
+               "%s: Phi is not positive.", &model->MOS3modName);
+            return(E_BADPARM);
+        }
+
         model->MOS3oxideCapFactor = 3.9 * 8.854214871e-12/
                 model->MOS3oxideThickness;
         if(!model->MOS3surfaceMobilityGiven) model->MOS3surfaceMobility=600;
diff --git a/src/spicelib/devices/mos6/mos6temp.c b/src/spicelib/devices/mos6/mos6temp.c
index 046f62c79..26442e8a2 100644
--- a/src/spicelib/devices/mos6/mos6temp.c
+++ b/src/spicelib/devices/mos6/mos6temp.c
@@ -51,6 +51,12 @@ MOS6temp(GENmodel *inModel, CKTcircuit *ckt)
 
     /* now model parameter preprocessing */
 
+        if (model->MOS6phi <= 0.0) {
+            SPfrontEnd->IFerror (ERR_FATAL,
+               "%s: Phi is not positive.", &model->MOS6modName);
+            return(E_BADPARM);
+        }
+
         if(!model->MOS6oxideThicknessGiven || model->MOS6oxideThickness == 0) {
             model->MOS6oxideCapFactor = 0;
         } else {
diff --git a/src/spicelib/devices/mos9/mos9temp.c b/src/spicelib/devices/mos9/mos9temp.c
index 0ffcd11df..237f09c0b 100644
--- a/src/spicelib/devices/mos9/mos9temp.c
+++ b/src/spicelib/devices/mos9/mos9temp.c
@@ -54,6 +54,11 @@ MOS9temp(GENmodel *inModel, CKTcircuit *ckt)
                                                                   CONSTKoverQ);
         ni_temp=1.45e16*nifact;
 
+        if (model->MOS9phi <= 0.0) {
+            SPfrontEnd->IFerror (ERR_FATAL,
+               "%s: Phi is not positive.", &model->MOS9modName);
+            return(E_BADPARM);
+        }
 
         model->MOS9oxideCapFactor = 3.9 * 8.854214871e-12/
                 model->MOS9oxideThickness;