Mathc complexes/a72

Application

Le code de la fonction :

/* ------------------------------------ */
/* ------------------------------------ */
double **c_A_b_Ab_mZ(
double **A,
double **b,
double **Ab
)
{
int r;
int c;
int c_A;

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

            Ab[r][c] = A[r][c];

             c_A = --c;

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

            Ab[r][c+c_A] = b[r][c];
            
return(Ab);
}
/* ------------------------------------ */
/* ------------------------------------ */

Copier A et b dans Ab. Les deux premières boucles "for" introduisent la matrice A. Les deux dernières boucle "for", introduisent b. 

 //La matrice Ab doit être initialiser avec la fonction :
 
 double **Ab = i_Abr_Ac_bc_mZ( RAb, CA, Cb); 
 
 // Cette fonction installe la bonne taille pour la matrices Ab.
 
 double **i_Abr_Ac_bc_mZ(
 int     Ab_r,                  // nombre de lignes   de Ab
 int      A_c,                  // nombre de colonnes de A
 int      b_c                   // nombre de colonnes de b
 )

Installer et compiler ce fichier dans votre répertoire de travail.

	c00a.c

/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A  = r_mZ(i_mZ(r,r ),999);
double **b  = r_mZ(i_mZ(r,C1),999);
double **Ab = i_Abr_Ac_bc_mZ(r,    r,   C1);
/*            i_Abr_Ac_bc_mZ(R_Ab, C_A, C_b));  */


  clrscrn();
  printf(" A :");
  p_mZ(A,S5,P0,S5,P0,C7);
  
  printf(" b :");
  p_mZ(b,S5,P0,S5,P0,C7);

  printf(" Ab :  c_A_b_Ab_mZ(A,b, Ab)");
  p_mZ(c_A_b_Ab_mZ(A,b, Ab),S5,P0,S5,P0,C7);

  f_mZ(Ab);
  f_mZ(b);
  f_mZ(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
do
{
  fun(rp_I(RC2)+R1); /* +R1 -> row >= R2 */


} while(stop_w());

  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */

Exemple de sortie écran :

 A :
 +998 -499i   +81 -495i  -331 -424i 
 -651 +864i  -797 +749i  -404 +464i 
 +292 +494i  +593 +517i  -196 +764i 

 b :
 +866 -384i 
 -451  +61i 
 -261 -216i 

 Ab :  c_A_b_Ab_mZ(A,b, Ab)
 +998 -499i   +81 -495i  -331 -424i  +866 -384i 
 -651 +864i  -797 +749i  -404 +464i  -451  +61i 
 +292 +494i  +593 +517i  -196 +764i  -261 -216i 


 Press return to continue
 Press X      to stop