From 52928065276bfc0a1eb1870dccf21e51c8596100 Mon Sep 17 00:00:00 2001
From: pnenzi <pnenzi>
Date: Sat, 20 Sep 2003 20:21:08 +0000
Subject: [PATCH] Added a new function to take into account dtemp.

---
 src/spicelib/analysis/nevalsrc.c | 40 ++++++++++++++++++++++++++++++++
 1 file changed, 40 insertions(+)

diff --git a/src/spicelib/analysis/nevalsrc.c b/src/spicelib/analysis/nevalsrc.c
index 14fb6d952..e17f9974e 100644
--- a/src/spicelib/analysis/nevalsrc.c
+++ b/src/spicelib/analysis/nevalsrc.c
@@ -49,3 +49,43 @@ NevalSrc (double *noise, double *lnNoise, CKTcircuit *ckt, int type, int node1,
 
     }
 }
+
+
+/*
+PN 2003:
+The following function includes instance dtemp in 
+thermal noise calculation. 
+It will replace NevalSrc as soon as all devices 
+will implement dtemp feature.
+*/
+
+void
+NevalSrc2 (double *noise, double *lnNoise, CKTcircuit *ckt, int type, 
+           int node1, int node2, double param, double param2)
+{
+    double realVal;
+    double imagVal;
+    double gain;
+
+    realVal = *((ckt->CKTrhs) + node1) - *((ckt->CKTrhs) + node2);
+    imagVal = *((ckt->CKTirhs) + node1) - *((ckt->CKTirhs) + node2);
+    gain = (realVal*realVal) + (imagVal*imagVal);
+    switch (type) {
+
+    case SHOTNOISE:
+	*noise = gain * 2 * CHARGE * fabs(param);          /* param is the dc current in a semiconductor */
+  	*lnNoise = log( MAX(*noise,N_MINLOG) ); 
+	break;
+
+    case THERMNOISE:
+	*noise = gain * 4 * CONSTboltz * (ckt->CKTtemp + param2)  /* param2 is the instance temperature difference */ 
+	         * param;                                         /* param is the conductance of a resistor */
+        *lnNoise = log( MAX(*noise,N_MINLOG) );
+        break;
+
+    case N_GAIN:
+	*noise = gain;
+	break;
+
+    }
+}