lu.c File Reference

#include "spConfig.h"
#include "spmatrix.h"
#include "spDefs.h"
#include "machine.h"

Include dependency graph for lu.c:

Go to the source code of this file.

Defines
#define	spINSIDE_SPARSE
Functions
static void	spFixThresold (char *eMatrix, double eps, double releps)
static void	spGetNumRank (char *eMatrix, int *n)
void	cerro ()
void C2F()	lufact1 (double *val, int *lln, int *col, int *n, int *nel, long *fmat, double *eps, double *releps, int *nrank, int *ierr)
void	Cout (char *str)
void C2F()	lusolve1 (long *fmat, double *b, double *x)
void C2F()	ludel1 (long *fmat)
static void	spSize (char *eMatrix, int *n)
void C2F()	lusize (long *fmat, int *n)
void	GetSigD (MatrixPtr Matrix, int indsigd[], double sigd[])
void	GetSigG (MatrixPtr Matrix, int indsigg[], double sigg[])
static void	spLuget (char *eMatrix, int *indP, double *P, int *indl, double *l, int *indu, double *u, int *indQ, double *Q)
void C2F()	luget1 (long *fmat, int *indP, double *P, int *indl, double *l, int *indu, double *u, int *indQ, double *Q)
static void	spLusiz (char *eMatrix, int *lsize, int *usize)
void C2F()	lusiz1 (long *fmat, int *lsize, int *usize)

---

Define Documentation

#define spINSIDE_SPARSE

Definition at line 31 of file lu.c.

---

Function Documentation

void cerro

(

)

void Cout

(

char *

str

)

Definition at line 7 of file cout.c.

References C2F, int, l, and out().

{
  int l;
  l = strlen(str) + 1;
  C2F(out)(str,l);
}

Here is the call graph for this function:

void GetSigD	(	MatrixPtr	Matrix,
		int	indsigd[],
		double	sigd[]
	)

Definition at line 190 of file lu.c.

References I, and last.

Referenced by spLuget().

{
  int I,J,mc=0,last=0;
  int Size=Matrix->Size;
for (I = 1; I <= Size; I++)
  {
    indsigd[I-1]=1;
    indsigd[Size+I-1]=  Matrix->IntToExtColMap[I];
    sigd[I-1]=1.0;
  }
/* counting missing colums */
for (I = 1; I <= Size; I++) 
  if (Matrix->ExtToIntColMap[I]== -1) mc++;
/* filling missing colums */
if (mc != 0)
  {
    for (I = Size -(mc)+1  ; I <= Size; I++)
      {
        for ( J=last+1; J <=Size; J++)
          {
            if (Matrix->ExtToIntColMap[J]==-1) 
              { 
                last=J;break;
              }
          }
        indsigd[I-1]=1;
        indsigd[Size+I-1]= last;
      }
  }
}

Here is the caller graph for this function:

void GetSigG	(	MatrixPtr	Matrix,
		int	indsigg[],
		double	sigg[]
	)

Definition at line 224 of file lu.c.

References I, and last.

Referenced by spLuget().

{
  int Size=Matrix->Size;
  int I,J,mc=0,last=0;
  /* counting missing Rows*/
  for (I = 1; I <= Size; I++) 
    if (Matrix->ExtToIntRowMap[I]== -1) mc++;

  for (I = 1; I <= Size-mc ; I++)
  {
        indsigg[I-1]= 1;
        indsigg[Size+Matrix->IntToExtRowMap[I]-1]=I;
        sigg[I-1]=1;
  }
/* filling missing Rows */
if (mc != 0)
  {
    for (I = Size -(mc)+1  ; I <= Size; I++)
      {
        for ( J=last+1; J <=Size; J++)
          {
            if (Matrix->ExtToIntRowMap[J]==-1) 
              { 

                last=J;break;
              }
          }
        indsigg[I-1]= 1;
        indsigg[Size+last-1]=I;
        sigg[I-1]=1;
      }
  }
}

Here is the caller graph for this function:

void C2F() ludel1

(

long *

fmat

)

Definition at line 153 of file lu.c.

References spDestroy().

{
  spDestroy((char*) *fmat);
}

Here is the call graph for this function:

void C2F() lufact1	(	double *	val,
		int*	lln,
		int*	col,
		int *	n,
		int*	nel,
		long *	fmat,
		double*	eps,
		double*	releps,
		int*	nrank,
		int*	ierr
	)

Definition at line 78 of file lu.c.

References cerro(), eps, error, i, i1, ierr, j, long, n, nel, spADD_REAL_ELEMENT, spCreate(), spFactor(), spFixThresold(), spGetElement(), spGetNumRank(), spNO_MEMORY, spOKAY, spREAL, spSINGULAR, spSMALL_PIVOT, spZERO_DIAG, and val.

{
  int error,i,i0,i1,k,j;
  spREAL *pelement;
  *ierr = 0;
  *fmat = (long)spCreate(*n,0,&error);
    if (error != spOKAY) {
        *ierr = 1;
        return;
    }


  i0=0;
  i1=i0;
  i=1;
  for (k = 0 ;k < *nel; k++) {
      i0=i0+1;
      while (i0-i1 > lln[i-1]) {
          i1=i0;
          i=i+1;
          i0=i0+1;
      }
      j=col[k];
      pelement = spGetElement((char*) *fmat,i,j);
      if (pelement == 0) {
          *ierr=2;
          return;
      }
      spADD_REAL_ELEMENT(pelement,(spREAL)(val[k]));

  }
  /* Fix the AbsThresold with scilex %eps */
  spFixThresold((char*) *fmat,*eps,*releps);
  /* spPrint((char *) *fmat,1,1,1); */
  error = spFactor((char*) *fmat);
  spGetNumRank((char *) *fmat,nrank);
  switch (error) {
  case spZERO_DIAG:
    cerro("zero_diag: A zero was encountered on the diagonal the matrix ");
    break;
  case spNO_MEMORY:
    *ierr=3;
    break;
  case spSINGULAR:
    *ierr=-1; /*Singular matrix" */
    break;
  case spSMALL_PIVOT:
    *ierr=-2; /* matrix is singular at precision level */
    break;
  }
}

Here is the call graph for this function:

void C2F() luget1	(	long *	fmat,
		int *	indP,
		double *	P,
		int*	indl,
		double*	l,
		int*	indu,
		double*	u,
		int*	indQ,
		double*	Q
	)

Definition at line 315 of file lu.c.

References l, P, and spLuget().

{
  spLuget((char *) *fmat,indP,P,indl,l,indu,u,indQ,Q);
}

Here is the call graph for this function:

void C2F() lusiz1	(	long *	fmat,
		int *	lsize,
		int*	usize
	)

Definition at line 355 of file lu.c.

References spLusiz().

{
  spLusiz((char *) *fmat,lsize,usize);
}

Here is the call graph for this function:

void C2F() lusize	(	long *	fmat,
		int *	n
	)

Definition at line 173 of file lu.c.

References n, and spSize().

{
  spSize((char *) *fmat,n);
}

Here is the call graph for this function:

void C2F() lusolve1	(	long *	fmat,
		double *	b,
		double *	x
	)

Definition at line 141 of file lu.c.

References b, spREAL, spSolve(), and x.

{
    spSolve((char*) *fmat,(spREAL*)b,(spREAL*)x);
}

Here is the call graph for this function:

static void spFixThresold	(	char *	eMatrix,
		double	eps,
		double	releps
	)			[static]

Definition at line 56 of file lu.c.

Referenced by lufact1().

{
  MatrixPtr  Matrix = (MatrixPtr)eMatrix;
  Matrix->AbsThreshold = eps;
  Matrix->RelThreshold = releps;
}

Here is the caller graph for this function:

static void spGetNumRank	(	char *	eMatrix,
		int *	n
	)			[static]

Definition at line 67 of file lu.c.

Referenced by lufact1().

{
  MatrixPtr  Matrix = (MatrixPtr)eMatrix;
  *n = Matrix->NumRank;
}

Here is the caller graph for this function:

static void spLuget	(	char *	eMatrix,
		int *	indP,
		double *	P,
		int*	indl,
		double*	l,
		int*	indu,
		double*	u,
		int*	indQ,
		double*	Q
	)			[static]

Definition at line 260 of file lu.c.

References MatrixElement::Col, GetSigD(), GetSigG(), I, MatrixElement::NextInRow, NULL, and MatrixElement::Real.

Referenced by luget1().

{
int I,J;
int lsize,usize;

MatrixPtr Matrix = (MatrixPtr) eMatrix;
ElementPtr  pElement;
int Size;
Size = Matrix->Size;
GetSigD(Matrix,indQ,Q);
GetSigG(Matrix,indP,P);
for (J = 1; J <= Size ; J++)
  {
    indl[J-1] = 0;
    indu[J-1] = 0;
}
lsize=0;
usize=0;
for (I = 1; I <= Size ; I++)
  {
    indu[I-1]=indu[I-1]+1;
    indu[Size+usize]=I;
    u[usize]=1.0;
    usize=usize+1;

    pElement = Matrix->FirstInRow[I];
    while ( pElement != NULL )
      {
          J = pElement->Col;
          if (I >= J) {
              indl[I-1] = indl[I-1]+1;
              indl[Size+lsize]=J;
              l[lsize]=(double) pElement->Real ;
              lsize=lsize+1;

          }
          else {
              indu[I-1] = indu[I-1]+1;
              indu[Size+usize]=J;
              u[usize]=(double) pElement->Real ;
              usize=usize+1;

          }
        pElement = pElement->NextInRow;
      };
  };
}

Here is the call graph for this function:

Here is the caller graph for this function:

static void spLusiz	(	char *	eMatrix,
		int *	lsize,
		int*	usize
	)			[static]

Definition at line 331 of file lu.c.

References MatrixElement::NextInCol, NULL, and MatrixElement::Row.

Referenced by lusiz1().

{
int J;
MatrixPtr Matrix = (MatrixPtr) eMatrix;
ElementPtr  pElement;
int Size;
Size = Matrix->Size;
*lsize=0;
*usize=Size;
for (J = 1; J <= Size ; J++)
  {
    pElement = Matrix->FirstInCol[J];
    while ( pElement != NULL )
      {
        if (pElement->Row >= J) 
            *lsize=*lsize+1;
        else
            *usize=*usize+1;
        pElement = pElement->NextInCol;
    };
  };
}  

Here is the caller graph for this function:

static void spSize	(	char *	eMatrix,
		int *	n
	)			[static]

Definition at line 165 of file lu.c.

Referenced by lusize().

{
  MatrixPtr  Matrix = (MatrixPtr)eMatrix;
  *n=Matrix->Size;
}

Here is the caller graph for this function:

---