intdgeesx.c File Reference

#include "stack-c.h"

Include dependency graph for intdgeesx.c:

Go to the source code of this file.

Defines
#define	MAX(x, y)   (((x)>(y))?(x):(y))
#define	MIN(x, y)   (((x)<(y))?(x):(y))
Functions
int	intdgeesx (char *fname)

---

Define Documentation

#define MAX	(	x,
		y		)	(((x)>(y))?(x):(y))

Definition at line 10 of file intdgeesx.c.

#define MIN	(	x,
		y		)	(((x)<(y))?(x):(y))

Definition at line 11 of file intdgeesx.c.

---

Function Documentation

int intdgeesx

(

char*

fname

)

Definition at line 20 of file intdgeesx.c.

References A, C2F, CheckLhs, CheckRhs, CheckSquare, CreateCVar, CreateVar, Errorinfo, GetRhsVar, it, Lhs, LhsVar, M, MAX, Maxvol, N, stk, and WORK.

{
        int A = 1, W = 2, VS = 3, WORK;
        int M, N, lA, lVS;
        int NLHS, it=1;
        int un=1, lWORK, LWORKMIN, LWORK, IWORK, LIWORK;
        int LDA, LDVS, lWR, lWI;
        int INFO, SDIM, BWORK;
        double RCONDE,RCONDV;
        char *JOBVS, *SORT, *SENSE;
        
        #ifndef WIN32
                extern select();
        #endif
        
        CheckRhs(1,1) ;
        CheckLhs(1,3) ;
        
        /*--------------------A---------------------------*/
        GetRhsVar(A, "d", &M, &N, &lA);
        CheckSquare(A, M, N);
        
        /*--------------------W---------------------------*/
        CreateCVar(W, "d", &it, &N, &un, &lWR, &lWI);
        
        NLHS=Lhs;
        LDA = MAX(1,N);
        LDVS=MAX(1,N);
        SORT="N";
        SENSE="N";
        
        switch ( NLHS ){
                case 3:     /* [W,VS,A]=dgeesx(A) */
                case 2:     /* [W,VS]=dgeesx(A)   */
                        printf("== 2 ==\r\n");
                        JOBVS="V";
                        /*--------------------VS---------------------------*/
                        CreateVar(VS, "d", &N, &N, &lVS);
                        /*--------------------WORK---------------------------*/
                        WORK=4; LWORKMIN = MAX(1, 3*N);
                        LWORK=Maxvol(WORK,"d");  /* max memory currently available */
                        if (LWORK < LWORKMIN) Scierror(999,"%s:  not enough memory (use stacksize) \r\n",fname);
                        
                        LIWORK = 1;
                        
                        CreateVar(WORK, "d", &LWORK, &un, &lWORK);
                        C2F(dgeesx)(JOBVS, SORT, select, SENSE, &N, stk(lA), &LDA, &SDIM, 
                                        stk(lWR), stk(lWI),
                                        stk(lVS), &LDVS, &RCONDE, &RCONDV, stk(lWORK), 
                                        &LWORK, &IWORK, &LIWORK, BWORK, &INFO);
                        
                        if (INFO != 0) Errorinfo("DGEESX ", INFO);
                        LhsVar(1)=W; LhsVar(2)=VS; LhsVar(3)=A;
                        break;
                case 1:   /*  [W]=dgeesx(A)    */
                        JOBVS="N";
                        /*--------------------WORK---------------------------*/
                        WORK=3; LWORKMIN = MAX(1, 3*N);
                        LWORK=Maxvol(WORK, "d"); /* max memory currently available */
                        if (LWORK < LWORKMIN)  Scierror(999,"%s:  not enough memory (use stacksize) \r\n",fname); 
                        CreateVar(WORK, "d", &LWORK, &un, &lWORK);
                        
                        C2F(dgeesx)(JOBVS, SORT, select, SENSE, &N, stk(lA), &LDA, &SDIM, 
                                        stk(lWR), stk(lWI),
                                        VS, &LDVS, &RCONDE, &RCONDV, stk(lWORK), 
                                        &LWORK, IWORK, LIWORK, BWORK, &INFO);
                        if (INFO != 0) Errorinfo("DGEESX ", INFO);
                        LhsVar(1)=W;
                        break;
                default:
                        return 0;
        }
        
        return 0;
}

---