intdsyev.c File Reference

#include "stack-c.h"

Include dependency graph for intdsyev.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	intdsyev (char *fname)

---

Define Documentation

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

Definition at line 5 of file intdsyev.c.

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

Definition at line 6 of file intdsyev.c.

---

Function Documentation

int intdsyev

(

char*

fname

)

Definition at line 8 of file intdsyev.c.

References A, C2F, CheckLhs, CheckRhs, CreateVar, Error, errorinfo(), GetRhsVar, L, Lhs, LhsVar, lw, M, MAX, maxvol(), N, sciprint(), stk, and WORK.

Referenced by C2F().

{ 
  int M,N,lA;
  int NLHS;
  int mU, nU, lU, U;
  int un=1; int WORK, lWORK, LWORKMIN, LWORK, INFO;
  int LDA, lw;
  int A=1, W=2;

  static int minlhs=1, minrhs=1, maxlhs=2, maxrhs=1;
  CheckRhs(minrhs,maxrhs) ;  CheckLhs(minlhs,maxlhs) ;

  /*--------------------A---------------------------*/
  GetRhsVar(A, "d", &M, &N, &lA);
  /*  To be done: Check A symmetric and Real*/
  NLHS=Lhs;
  switch ( NLHS ) {
  case 2:        
    /*  [w,U]=dsyev(A)     */
    CreateVar(W, "d", &N, &un, &lw);
    WORK=3;
    LWORKMIN = MAX(1, 3*N-1);
    LWORK = LWORKMIN;
    /* LWORK=C2F(maxvol)(&WORK, "d", 1L);   max memory currently available */
    if (LWORK < LWORKMIN) { 
      sciprint("not enough memory (use stacksize)");
      Error(9999); };
    CreateVar(3, "d", &LWORK, &un, &lWORK);
    LDA = MAX(1,N);
    C2F(dsyev)("V", "L", &N, stk(lA), &LDA, stk(lw), stk(lWORK), &LWORK, &INFO);
    if (INFO != 0) C2F(errorinfo)("dsyev ", &INFO, 5L);
    LhsVar(1)=A; LhsVar(2)=W;
    return 0;
    break;
  case 1:
    /*   w=dsyev(A)    */
    CreateVar(W, "d", &N, &un, &lw);
    WORK=3;
    LWORKMIN = MAX(1, 3*N-1);
    LWORK=C2F(maxvol)(&WORK, "d", 1L);  /* max memory currently available */
    if (LWORK < LWORKMIN) { 
      sciprint("not enough memory (use stacksize)");
      Error(9999); };
    CreateVar(3, "d", &LWORK, &un, &lWORK);
    LDA = MAX(1,N);
    C2F(dsyev)("N", "L", &N, stk(lA), &LDA, stk(lw), stk(lWORK), &LWORK, &INFO);
    if (INFO != 0) C2F(errorinfo)("dsyev ", &INFO, 5L);
    LhsVar(1)=W;
    return 0;
    break;
  default:
    return 0;
  }
  return(0);
}

Here is the call graph for this function:

Here is the caller graph for this function:

---