luke
/
ngspice
mirror of https://git.code.sf.net/p/ngspice/ngspice
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

New Binding Method extended to NBJT device

This commit is contained in:
Francesco Lannutti 2013-04-01 22:50:18 +02:00
parent 6c2704e400
commit 18e0a7ac5d
3 changed files with 235 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
src/spicelib/devices/nbjt/Makefile.am
View File

@ -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)

211
src/spicelib/devices/nbjt/nbjtbindCSC.c Normal file
View File

@ -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) ;
}

13
src/spicelib/devices/nbjt/nbjtdefs.h
View File

@ -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 */