correct init state vector for qth integration

6 months ago · 931e223493
1 changed files with 4 additions and 0 deletions
--- a/src/spicelib/devices/dio/dioload.c
+++ b/src/spicelib/devices/dio/dioload.c
@ -492,6 +492,10 @@ next1:
                    }
                    if (selfheat)
                    {
+                        if (ckt->CKTmode & MODEINITTRAN) {
+                            *(ckt->CKTstate1 + here->DIOqth) =
+                                    *(ckt->CKTstate0 + here->DIOqth);
+                        }
                        error = NIintegrate(ckt, &gcTt, &ceqqth, model->DIOcth0, here->DIOqth);
                        if (error) return(error);
                        if (ckt->CKTmode & MODEINITTRAN) {