Mathc complexes/h09n

Bibliothèque

Installer ce fichier dans votre répertoire de travail.

	wgjinv.h

/* ------------------------------------ */
/*  Save as :   wgjinv.h                */
/* ------------------------------------ */
void invgj_T_mZ(
double **Ab
)
{
int r;
int c;
double pivot;

      r = c = C1;

    while((r<Ab[R_SIZE][C0]) &&
          (c<Ab[C_SIZE_A][C0])
         )
       {
        pivot = pivotbestAId_mZ(Ab,r,c);
        if(fabs(pivot)>ERROR_E)  zero_under_pivot_mZ(Ab,r,c);
        r++;
        c+=C2;
       }

    r--;
    c-=C2;

    while( (r>R1) )
       {
        zero_below_pivot_gj_mZ(Ab,r,c);
        r--;
        c-=C2;
       }
}
/* ------------------------------------ */
double **invgj_mZ(
double **A,
double **invA
)
{
double **ID = i_RC_mZ(A[R_SIZE][C0],A[C_SIZE][C0]);
double **AID;
int r = A[R_SIZE][C0];
   
 --r;

  AID = i_Abr_Ac_bc_mZ(r,r,r);

         eye_mZ(ID);
  c_A_b_Ab_mZ(A,ID,AID);
  invgj_T_mZ(AID);
   sort_c_mZ(AID);
   sort_r_mZ(AID);
  c_Inv_A_mZ(AID,invA);

  f_mZ(AID);
  f_mZ(ID);
  
return(invA);
}
/* ------------------------------------ */
/* ------------------------------------ */

Déclaration des fichiers h.