diff --git a/src/spicelib/devices/nbjt/Makefile.am b/src/spicelib/devices/nbjt/Makefile.am
index 4a5c3ee8a..f6d1287d9 100644
--- a/src/spicelib/devices/nbjt/Makefile.am
+++ b/src/spicelib/devices/nbjt/Makefile.am
@@ -24,6 +24,17 @@ libnbjt_la_SOURCES =	\
 	nbjttrun.c
 
 
+if KLU_WANTED
+libnbjt_la_SOURCES += nbjtbindCSC.c
+endif
+
+if SuperLU_WANTED
+libnbjt_la_SOURCES += nbjtbindCSC.c
+endif
+
+if UMFPACK_WANTED
+libnbjt_la_SOURCES += nbjtbindCSC.c
+endif
 
 AM_CPPFLAGS = @AM_CPPFLAGS@ -I$(top_srcdir)/src/include
 AM_CFLAGS = $(STATIC)
diff --git a/src/spicelib/devices/nbjt/nbjtbindCSC.c b/src/spicelib/devices/nbjt/nbjtbindCSC.c
new file mode 100644
index 000000000..c910971ab
--- /dev/null
+++ b/src/spicelib/devices/nbjt/nbjtbindCSC.c
@@ -0,0 +1,211 @@
+/**********
+Author: 2013 Francesco Lannutti
+**********/
+
+#include "ngspice/ngspice.h"
+#include "ngspice/cktdefs.h"
+#include "nbjtdefs.h"
+#include "ngspice/sperror.h"
+
+#include <stdlib.h>
+
+static
+int
+BindCompare (const void *a, const void *b)
+{
+    BindElement *A, *B ;
+    A = (BindElement *)a ;
+    B = (BindElement *)b ;
+
+    return ((int)(A->Sparse - B->Sparse)) ;
+}
+
+int
+NBJTbindCSC (GENmodel *inModel, CKTcircuit *ckt)
+{
+    NBJTmodel *model = (NBJTmodel *)inModel ;
+    NBJTinstance *here ;
+    double *i ;
+    BindElement *matched, *BindStruct ;
+    size_t nz ;
+
+    BindStruct = ckt->CKTmatrix->CKTbindStruct ;
+    nz = (size_t)ckt->CKTmatrix->CKTklunz ;
+
+    /* loop through all the NBJT models */
+    for ( ; model != NULL ; model = model->NBJTnextModel)
+    {
+        /* loop through all the instances of the model */
+        for (here = model->NBJTinstances ; here != NULL ; here = here->NBJTnextInstance)
+        {
+            if ((here-> NBJTcolNode != 0) && (here-> NBJTcolNode != 0))
+            {
+                i = here->NBJTcolColPtr ;
+                matched = (BindElement *) bsearch (&i, BindStruct, nz, sizeof(BindElement), BindCompare) ;
+                here->NBJTcolColStructPtr = matched ;
+                here->NBJTcolColPtr = matched->CSC ;
+            }
+
+            if ((here-> NBJTbaseNode != 0) && (here-> NBJTbaseNode != 0))
+            {
+                i = here->NBJTbaseBasePtr ;
+                matched = (BindElement *) bsearch (&i, BindStruct, nz, sizeof(BindElement), BindCompare) ;
+                here->NBJTbaseBaseStructPtr = matched ;
+                here->NBJTbaseBasePtr = matched->CSC ;
+            }
+
+            if ((here-> NBJTemitNode != 0) && (here-> NBJTemitNode != 0))
+            {
+                i = here->NBJTemitEmitPtr ;
+                matched = (BindElement *) bsearch (&i, BindStruct, nz, sizeof(BindElement), BindCompare) ;
+                here->NBJTemitEmitStructPtr = matched ;
+                here->NBJTemitEmitPtr = matched->CSC ;
+            }
+
+            if ((here-> NBJTcolNode != 0) && (here-> NBJTbaseNode != 0))
+            {
+                i = here->NBJTcolBasePtr ;
+                matched = (BindElement *) bsearch (&i, BindStruct, nz, sizeof(BindElement), BindCompare) ;
+                here->NBJTcolBaseStructPtr = matched ;
+                here->NBJTcolBasePtr = matched->CSC ;
+            }
+
+            if ((here-> NBJTcolNode != 0) && (here-> NBJTemitNode != 0))
+            {
+                i = here->NBJTcolEmitPtr ;
+                matched = (BindElement *) bsearch (&i, BindStruct, nz, sizeof(BindElement), BindCompare) ;
+                here->NBJTcolEmitStructPtr = matched ;
+                here->NBJTcolEmitPtr = matched->CSC ;
+            }
+
+            if ((here-> NBJTbaseNode != 0) && (here-> NBJTcolNode != 0))
+            {
+                i = here->NBJTbaseColPtr ;
+                matched = (BindElement *) bsearch (&i, BindStruct, nz, sizeof(BindElement), BindCompare) ;
+                here->NBJTbaseColStructPtr = matched ;
+                here->NBJTbaseColPtr = matched->CSC ;
+            }
+
+            if ((here-> NBJTbaseNode != 0) && (here-> NBJTemitNode != 0))
+            {
+                i = here->NBJTbaseEmitPtr ;
+                matched = (BindElement *) bsearch (&i, BindStruct, nz, sizeof(BindElement), BindCompare) ;
+                here->NBJTbaseEmitStructPtr = matched ;
+                here->NBJTbaseEmitPtr = matched->CSC ;
+            }
+
+            if ((here-> NBJTemitNode != 0) && (here-> NBJTcolNode != 0))
+            {
+                i = here->NBJTemitColPtr ;
+                matched = (BindElement *) bsearch (&i, BindStruct, nz, sizeof(BindElement), BindCompare) ;
+                here->NBJTemitColStructPtr = matched ;
+                here->NBJTemitColPtr = matched->CSC ;
+            }
+
+            if ((here-> NBJTemitNode != 0) && (here-> NBJTbaseNode != 0))
+            {
+                i = here->NBJTemitBasePtr ;
+                matched = (BindElement *) bsearch (&i, BindStruct, nz, sizeof(BindElement), BindCompare) ;
+                here->NBJTemitBaseStructPtr = matched ;
+                here->NBJTemitBasePtr = matched->CSC ;
+            }
+
+        }
+    }
+
+    return (OK) ;
+}
+
+int
+NBJTbindCSCComplex (GENmodel *inModel, CKTcircuit *ckt)
+{
+    NBJTmodel *model = (NBJTmodel *)inModel ;
+    NBJTinstance *here ;
+
+    NG_IGNORE (ckt) ;
+
+    /* loop through all the NBJT models */
+    for ( ; model != NULL ; model = model->NBJTnextModel)
+    {
+        /* loop through all the instances of the model */
+        for (here = model->NBJTinstances ; here != NULL ; here = here->NBJTnextInstance)
+        {
+            if ((here-> NBJTcolNode != 0) && (here-> NBJTcolNode != 0))
+                here->NBJTcolColPtr = here->NBJTcolColStructPtr->CSC_Complex ;
+
+            if ((here-> NBJTbaseNode != 0) && (here-> NBJTbaseNode != 0))
+                here->NBJTbaseBasePtr = here->NBJTbaseBaseStructPtr->CSC_Complex ;
+
+            if ((here-> NBJTemitNode != 0) && (here-> NBJTemitNode != 0))
+                here->NBJTemitEmitPtr = here->NBJTemitEmitStructPtr->CSC_Complex ;
+
+            if ((here-> NBJTcolNode != 0) && (here-> NBJTbaseNode != 0))
+                here->NBJTcolBasePtr = here->NBJTcolBaseStructPtr->CSC_Complex ;
+
+            if ((here-> NBJTcolNode != 0) && (here-> NBJTemitNode != 0))
+                here->NBJTcolEmitPtr = here->NBJTcolEmitStructPtr->CSC_Complex ;
+
+            if ((here-> NBJTbaseNode != 0) && (here-> NBJTcolNode != 0))
+                here->NBJTbaseColPtr = here->NBJTbaseColStructPtr->CSC_Complex ;
+
+            if ((here-> NBJTbaseNode != 0) && (here-> NBJTemitNode != 0))
+                here->NBJTbaseEmitPtr = here->NBJTbaseEmitStructPtr->CSC_Complex ;
+
+            if ((here-> NBJTemitNode != 0) && (here-> NBJTcolNode != 0))
+                here->NBJTemitColPtr = here->NBJTemitColStructPtr->CSC_Complex ;
+
+            if ((here-> NBJTemitNode != 0) && (here-> NBJTbaseNode != 0))
+                here->NBJTemitBasePtr = here->NBJTemitBaseStructPtr->CSC_Complex ;
+
+        }
+    }
+
+    return (OK) ;
+}
+
+int
+NBJTbindCSCComplexToReal (GENmodel *inModel, CKTcircuit *ckt)
+{
+    NBJTmodel *model = (NBJTmodel *)inModel ;
+    NBJTinstance *here ;
+
+    NG_IGNORE (ckt) ;
+
+    /* loop through all the NBJT models */
+    for ( ; model != NULL ; model = model->NBJTnextModel)
+    {
+        /* loop through all the instances of the model */
+        for (here = model->NBJTinstances ; here != NULL ; here = here->NBJTnextInstance)
+        {
+            if ((here-> NBJTcolNode != 0) && (here-> NBJTcolNode != 0))
+                here->NBJTcolColPtr = here->NBJTcolColStructPtr->CSC ;
+
+            if ((here-> NBJTbaseNode != 0) && (here-> NBJTbaseNode != 0))
+                here->NBJTbaseBasePtr = here->NBJTbaseBaseStructPtr->CSC ;
+
+            if ((here-> NBJTemitNode != 0) && (here-> NBJTemitNode != 0))
+                here->NBJTemitEmitPtr = here->NBJTemitEmitStructPtr->CSC ;
+
+            if ((here-> NBJTcolNode != 0) && (here-> NBJTbaseNode != 0))
+                here->NBJTcolBasePtr = here->NBJTcolBaseStructPtr->CSC ;
+
+            if ((here-> NBJTcolNode != 0) && (here-> NBJTemitNode != 0))
+                here->NBJTcolEmitPtr = here->NBJTcolEmitStructPtr->CSC ;
+
+            if ((here-> NBJTbaseNode != 0) && (here-> NBJTcolNode != 0))
+                here->NBJTbaseColPtr = here->NBJTbaseColStructPtr->CSC ;
+
+            if ((here-> NBJTbaseNode != 0) && (here-> NBJTemitNode != 0))
+                here->NBJTbaseEmitPtr = here->NBJTbaseEmitStructPtr->CSC ;
+
+            if ((here-> NBJTemitNode != 0) && (here-> NBJTcolNode != 0))
+                here->NBJTemitColPtr = here->NBJTemitColStructPtr->CSC ;
+
+            if ((here-> NBJTemitNode != 0) && (here-> NBJTbaseNode != 0))
+                here->NBJTemitBasePtr = here->NBJTemitBaseStructPtr->CSC ;
+
+        }
+    }
+
+    return (OK) ;
+}
diff --git a/src/spicelib/devices/nbjt/nbjtdefs.h b/src/spicelib/devices/nbjt/nbjtdefs.h
index 886841130..0db007074 100644
--- a/src/spicelib/devices/nbjt/nbjtdefs.h
+++ b/src/spicelib/devices/nbjt/nbjtdefs.h
@@ -83,6 +83,19 @@ typedef struct sNBJTinstance {
   unsigned NBJTicFileGiven:1;	/* flag to indicate init. cond. file given */
   unsigned NBJTprintGiven:1;	/* flag to indicate if print was given */
   unsigned NBJTtempGiven:1;	/* flag to indicate if temp was given */
+
+#ifdef KLU
+    BindElement *NBJTcolColStructPtr ;
+    BindElement *NBJTbaseBaseStructPtr ;
+    BindElement *NBJTemitEmitStructPtr ;
+    BindElement *NBJTcolBaseStructPtr ;
+    BindElement *NBJTcolEmitStructPtr ;
+    BindElement *NBJTbaseColStructPtr ;
+    BindElement *NBJTbaseEmitStructPtr ;
+    BindElement *NBJTemitColStructPtr ;
+    BindElement *NBJTemitBaseStructPtr ;
+#endif
+
 } NBJTinstance;
 
 /* per model data */