Mathc complexes/Fichiers h : wel cop

Installer ce fichier dans votre répertoire de travail.

	`wel_cop.h`

/* ------------------------------------ */
/*  Save as :  wel_cop.h                */
/* ------------------------------------ */
/* ------------------------------------ */
int TheCorrect_C(int C)
{
int c = C;

     if(c!=C1) c= (C-C1)*C2+C1;

return(c);
}
/* ------------------------------------ */
/* ------------------------------------ */
double **swapC_mZ(
double **A,
int c1,
int c2)
{
double T;
int r;

    for ( r=R0; r<A[R_SIZE][C0]; r++)
        {
         T        = A[r][c1];
         A[r][c1] = A[r][c2];
         A[r][c2] = T;

         T           = A[r][c1+C1];
         A[r][c1+C1] = A[r][c2+C1];
         A[r][c2+C1] = T;
        }
        
return(A);
}
/* ------------------------------------ */
double **mulC_mZ(
double **A,
nb_Z z,
int  c
)
{
nb_Z T;
int r;

    for ( r=R1; r<A[R_SIZE][C0]; r++)
        {
          T = mul_Z( i_Z(A[r][c],
                         A[r][c+C1]),
                     z);

          A[r][c] =    T.r;
          A[r][c+C1] = T.i;
         }
        
return(A);
}
/* ------------------------------------ */
double **addC_mZ(
double **A,
nb_Z z,
int c1,
int c2
)
{
int      r;
nb_Z T;

    for ( r=R1; r<A[R_SIZE][C0]; r++)
         {
          T = add_Z(
                     mul_Z(i_Z(A[r][c1],A[r][c1+C1]), z),
                           i_Z(A[r][c2],A[r][c2+C1]));

          A[r][c2]    = T.r;
          A[r][c2+C1] = T.i;
         }
        
return(A);
}
/* ------------------------------------ */
/* ------------------------------------ */

Il y a les trois opérations élémentaires sur les colonnes. Il faut utiliser la fonction TheCorrect_C(); pour selectionner les bonnes colonnes avec les fonctions de ce fichier.

Si on veut travailler sur la colonne trois, il faut écrire c = TheCorrect_C(C3); et utiliser c au lieu de C3.