Mathc initiation/a476

Sommaire

Installer ce fichier dans votre répertoire de travail.

	x_juvw.h

/* --------------------------------- */
/* save as x_juvw.h                  */
/* --------------------------------- */
double J_x_uvw(

double (*P_f)(double x, double y, double z),

double (*P_X)(double x, double y, double z),
double (*P_Y)(double x, double y, double z),
double (*P_Z)(double x, double y, double z),

double (*Pux)(double y, double z),
double (*Pvx)(double y, double z),
   int nx,
   
double y,
   
double z
)
{
pt3d   p;

double x = 0;

   int i = 0;
double m = 0.;
double M = 0.;

 for(i = 0; i <= nx; i++)
 {
       if(i ==0 || i== nx){m = 1.;}
  else if(fmod(i,2) == 0){m = 2.;}
  else                   {m = 4.;}

  x = ((*Pux)(y,z)) + i*(((*Pvx)(y,z))-((*Pux)(y,z)))/nx;
 
  p.x = x;
  p.y = y;   
  p.z = z;    
     
  M += m * (*P_f)( ((*P_X)(x,y,z)), 
                   ((*P_Y)(x,y,z)), 
                   ((*P_Z)(x,y,z)))
                          *
                   /* Jacobian */   
                      
                    fabs(
                    
              fxyz_x((*P_X),H,p) 
                       * 
   ( fxyz_y((*P_Y),H,p)*fxyz_z((*P_Z),H,p)
                       -
     fxyz_z((*P_Y),H,p)*fxyz_y((*P_Z),H,p) ) 
     
                       - 
                                       
              fxyz_y((*P_X),H,p) 
                       * 
   ( fxyz_x((*P_Y),H,p)*fxyz_z((*P_Z),H,p)
                       -
     fxyz_z((*P_Y),H,p)*fxyz_x((*P_Z),H,p) )
     
                       + 
                                       
              fxyz_z((*P_X),H,p) 
                       * 
   ( fxyz_x((*P_Y),H,p)*fxyz_y((*P_Z),H,p)
                       -
     fxyz_y((*P_Y),H,p)*fxyz_x((*P_Z),H,p) )
                        );  
 }

  return( ((((*Pvx)(y,z)) -((*Pux)(y,z)))*M) / (3*nx) );
}
/* --------------------------------- */
double J_y_uvw(
double (*P_f)(double x, double y, double z),

double (*P_X)(double x, double y, double z),
double (*P_Y)(double x, double y, double z),
double (*P_Z)(double x, double y, double z),
   
double (*Pux)(double y, double z),
double (*Pvx)(double y, double z),
   int   nx,

double (*Psy)(double x),
double (*Pty)(double x),
   int   ny,
   
double z
)
{
   int i = 0;
double m = 0.;
double M = 0.;

 for(i = 0; i <= ny; i++)
 {
       if(i ==0 || i== ny){m = 1.;}
  else if(fmod(i,2) == 0){m = 2.;}
  else                   {m = 4.;}

  M += m * J_x_uvw( (*P_f),
                    (*P_X),(*P_Y),(*P_Z),                             
                    (*Pux),(*Pvx),nx,
                    (((*Psy)(z))+i*(((*Pty)(z))-((*Psy)(z)))/ny),
                             z);
 }

 return( ((((*Pty)(z)) -((*Psy)(z)))*M) / (3*ny) );
}
/* --------------------------------- */
double J_uvw(
double (*P_f)(double x, double y, double z),

double (*P_X)(double x, double y, double z),
double (*P_Y)(double x, double y, double z),
double (*P_Z)(double x, double y, double z),
   
double (*Pux)(double y, double z),
double (*Pvx)(double y, double z),
   int   nx,
   
double (*Psy)(double x),
double (*Pty)(double x),
   int   ny,

double az,
double bz,
   int nz
)
{
   int i = 0;
double m = 0.;
double M = 0.;

 for(i = 0; i <= nz; i++)
 {
       if(i ==0 || i== nz){m = 1.;}
  else if(fmod(i,2) == 0){m = 2.;}
  else                   {m = 4.;}

  M += m * J_y_uvw( (*P_f),
                    (*P_X),(*P_Y),(*P_Z),                                
                    (*Pux),(*Pvx),nx,
                    (*Psy),(*Pty),ny,
                    (az+i*(bz-az)/nz));
 }

  return( ((bz -az)*M) / (3*nz) );
}
/* --------------------------------- */
/* --------------------------------- */

----