niditer.c File Reference

#include "spice.h"
#include <stdio.h>
#include "cktdefs.h"
#include "util.h"
#include "sperror.h"
#include "cktext.h"
#include "niext.h"

Include dependency graph for niditer.c:

Go to the source code of this file.

Functions
int	NIdIter (CKTcircuit *ckt)

Function Documentation

int NIdIter

(

CKTcircuit *

ckt

)

Definition at line 26 of file niditer.c.

{
    int error;
    int ignore;
    double *temp;

    for (;;) {
        ckt->CKTnoncon = 0;

        if (ckt->CKTniState & NIACSHOULDREORDER) {
            spSetComplex(ckt->CKTmatrix);
            error = spOrderAndFactor(ckt->CKTmatrix,NULL,
                ckt->CKTpivotRelTol,ckt->CKTpivotAbsTol,1);
            ckt->CKTniState &= ~NIACSHOULDREORDER;
            if (error != 0) {
                /* either singular equations or no memory, in either case,
                 * let caller handle problem
                 */
                return (error);
            }
        }
        else {
            spSetComplex(ckt->CKTmatrix);
            error = spFactor(ckt->CKTmatrix);
            if (error != 0) {
                if (error == E_SINGULAR) {
                    /* the problem is that the matrix can't be solved with
                     * the current LU factorization.  Maybe if we reload and
                     * try to reorder again it will help...
                     */
                    ckt->CKTniState |= NIACSHOULDREORDER;
                    temp = ckt->CKTrhs;
                    ckt->CKTrhs = ckt->CKTrhsSpare;
                    ckt->CKTrhsSpare = temp;

                    temp = ckt->CKTirhs;
                    ckt->CKTirhs = ckt->CKTirhsSpare;
                    ckt->CKTirhsSpare = temp;

                    error = CKTacLoad(ckt);
                    if (error) return(error);
    
                    temp = ckt->CKTrhs;
                    ckt->CKTrhs = ckt->CKTrhsSpare;
                    ckt->CKTrhsSpare = temp;

                    temp = ckt->CKTirhs;
                    ckt->CKTirhs = ckt->CKTirhsSpare;
                    ckt->CKTirhsSpare = temp;
                    continue;
                }
                return (error); /* can't handle E_BADMATRIX, so let caller */
            }
        }
        break;
    }
    spSolve(ckt->CKTmatrix,ckt->CKTrhs,ckt->CKTrhs, 
            ckt->CKTirhs,ckt->CKTirhs);

    *ckt->CKTrhs = 0;
    *ckt->CKTrhsSpare = 0;
    *ckt->CKTrhsOld = 0;
    *ckt->CKTirhs = 0;
    *ckt->CKTirhsSpare = 0;
    *ckt->CKTirhsOld = 0;

    temp = ckt->CKTirhsOld;
    ckt->CKTirhsOld = ckt->CKTirhs;
    ckt->CKTirhs = temp;

    temp = ckt->CKTrhsOld;
    ckt->CKTrhsOld = ckt->CKTrhs;
    ckt->CKTrhs = temp;
    return (OK);
}