Mathc complexes/Fichiers h : vsm

Bibliothèque

Installer ce fichier dans votre répertoire de travail.

	wsm.h

/* ------------------------------------ */
/*  Save as :   wi_sm.h                 */
/* ------------------------------------ */
double **m0_mZ(
double **Zer
)
{
int r;
int c;

  for  (r=R1; r<Zer[R_SIZE][C0]; r++)
    for(c=R1; c<Zer[C_SIZE][C0]; c+=C2)
          {
           Zer[r][c   ] = 0.;
           Zer[r][c+C1] = 0.;
          }
          
return(Zer);
}
/* ------------------------------------ */
double **m1_mZ(
double **One
)
{
int r;
int c;

  for  (r=R1; r<One[R_SIZE][C0]; r++)
    for(c=R1; c<One[C_SIZE][C0]; c+=C2)
          {
           One[r][c   ] = 1.;
           One[r][c+C1] = 0.;
          }
          
return(One);
}
/* ------------------------------------ */
/* ------------------------------------ */

/* ------------------------------------ */
/* ------------------------------------ */
double **erase_i_mZ(
double **Zer
)
{
int r;
int c;

  for  (r=R1; r<Zer[R_SIZE][C0];   r++)
    for(c=R1; c<Zer[C_SIZE][C0]; c+=C2)
           
           Zer[r][c+C1] = 0.;
                
return(Zer);
}
/* ------------------------------------ */
/* ------------------------------------ */

/* ------------------------------------ */
/* ------------------------------------ */
double **eye_mZ(
double **ID
)
{
int r;
int c;

  isquare_mZ(ID,"eye_mZ();","(ID)");

  for (   r=R1; r<ID[R_SIZE][C0]; r++)
    for ( c=R1; c<ID[C_SIZE][C0]; c+=C2)

        if(((r-R1)*C2)==c-C1)
            {
                 ID[r][c   ] = 1.;
                 ID[r][c+C1] = 0.;
            }
        else{
                 ID[r][c   ] = 0.;
                 ID[r][c+C1] = 0.;
            }
            
return(ID);
}
/* ------------------------------------ */
/* ------------------------------------ */
double **clean_eyes_mZ(
double **ID
)
{
int r;
int c;

  isquare_mZ(ID,"clean_eyes_mZ;","(ID)");

  for (   r=R1; r<ID[R_SIZE][C0]; r++)
    for ( c=R1; c<ID[C_SIZE][C0]; c+=C2)

        if(((r-R1)*C2)!=c-C1)
            {
                 ID[r][c   ] = 0.;
                 ID[r][c+C1] = 0.;
            }
            
return(ID);
}
/* ------------------------------------ */
/* ------------------------------------ */

Dans ce fichier h on introduit la construction de quelques matrices spécifiques, la matrice zéro, la matrice identité.