Mathc complexes/h10m

Bibliothèque

Installer ce fichier dans votre répertoire de travail.

	wrm_sing.h

/* ------------------------------------ */
/*  Save as :   wrm_sing.h              */
/* ------------------------------------ */
double **r_singular_Ab_mZ(
double **A,
int      n
)
{

  r_mZ(A,n);

  if(A[R_SIZE][C0]==3) c_r1Ar2_mZ(1,A, 2);
  if(A[R_SIZE][C0]==4){c_r1Ar2_mZ(  rp_I(A[R_SIZE][C0]-1),
                                  A, 
                                    rp_I(A[R_SIZE][C0]-1));
                       c_r1Ar2_mZ(  rp_I(A[R_SIZE][C0]-1),
                                  A,
                                    rp_I(A[R_SIZE][C0]-1));
                      }
  if(A[R_SIZE][C0]>4){  c_r1Ar2_mZ( rp_I(A[R_SIZE][C0]-1),
                                   A,
                                     rp_I(A[R_SIZE][C0]-1));
                        c_r1Ar2_mZ(  rp_I(A[R_SIZE][C0]-1),
                                   A,
                                     rp_I(A[R_SIZE][C0]-1));
                        c_r1Ar2_mZ(  rp_I(A[R_SIZE][C0]-1),
                                   A,
                                     rp_I(A[R_SIZE][C0]-1));
                        c_r1Ar2_mZ(  rp_I(A[R_SIZE][C0]-1),
                                   A,
                                     rp_I(A[R_SIZE][C0]-1));
                      }
        
return(A);
}
/* ------------------------------------ */
/* ------------------------------------ */

Création de matrices singulières.