/* ------------------------------------ */
/*  Save as :   wgc_m.h                  */
/* ------------------------------------ */
double **cR0_mZ(
double **A,
double **B
)
{
int c;
int r;

 	for    (r=R0;    r<A[R_SIZE][C0]; r++)
 	   for (c=C1; c<A[C_SIZE][C0]; c++)

            B[r][c] = A[r][c];
            
return(B);
}
/* ------------------------------------ */
void c_n_r_mZ(
double **A,
int r_n,
double **B
)
{
int r;
int c;

 	for    (r=R1; r<r_n; r++)
 	   for (c=C1; c<A[C_SIZE][C0]; c++)

            B[r][c] = A[r][c];
}
/* ------------------------------------ */
void c_r1Ar2_mZ(
double **A,
int r1,
int r2)
{
int c;

    for ( c=C1; c<A[C_SIZE][C0]; c++)

             A[r2][c] = A[r1][c];
}
/* ------------------------------------ */
void c_r1zr2_mZ(
double **A,
int r1,
int r2,
nb_Z z)
{
nb_Z T;
int c;

    for ( c=C1; c<A[C_SIZE][C0]; c+=C2)
     {
      T = mul_Z(i_Z(A[r1][c],A[r1][c+C1]),z);

             A[r2][c]    = T.r;
             A[r2][c+C1] = T.i;
     }
}
/* ------------------------------------ */
void cR0_c_mZ(
double **A,
int cA,
double **B,
int cB
)
{
int r;

    for(r=R0; r<A[R_SIZE][C0]; r++)
       {
        B[r][cB]    = A[r][cA];
        B[r][cB+C1] = A[r][cA+C1];
       }
}
/* ------------------------------------ */
/* ------------------------------------ */