rD.h

Go to the documentation of this file.

#ifndef R_D_H
#define R_D_H  "rD.h"


//-----------------------------------------------------------------------------
// Define second order difference approximations:
// The include file cgux2a.h defines second order difference approximations
//
// Notes:
//  For efficiency we define UR2, US2 and UT2 to be the arrays that we have precomputed
//  the parametric derivatives in. (Otherwise the derivatives would get
//  recomputed for all derivatives). UR2, US2 and UT2 are used in cgux2a.h
//  to define UX2, UY2, UXX2, ...
//-----------------------------------------------------------------------------
#define U(I1,I2,I3,N)   u(I1,I2,I3,N)
#define UR2(I1,I2,I3,KD) ( (U(I1+1,I2,I3,KD)-U(I1-1,I2,I3,KD))*d12(axis1) )
#define US2(I1,I2,I3,KD) ( (U(I1,I2+1,I3,KD)-U(I1,I2-1,I3,KD))*d12(axis2) )
#define UT2(I1,I2,I3,KD) ( (U(I1,I2,I3+1,KD)-U(I1,I2,I3-1,KD))*d12(axis3) )
#include "cgux2a.h"    // define 2nd order difference approximations

//-----------------------------------------------------------------------------
// Define fourth order difference approximations:
// The include file cgux4a.h defines fourth order difference approximations
//
// Notes : see above
//----------------------------------------------------------------------------
#define UR4(I1,I2,I3,KD) ( (8.*(U(I1+1,I2,I3,KD)-U(I1-1,I2,I3,KD))  \
                              -(U(I1+2,I2,I3,KD)-U(I1-2,I2,I3,KD)))*d14(axis1) )
#define US4(I1,I2,I3,KD) ( (8.*(U(I1,I2+1,I3,KD)-U(I1,I2-1,I3,KD))  \
                              -(U(I1,I2+2,I3,KD)-U(I1,I2-2,I3,KD)))*d14(axis2) )
#define UT4(I1,I2,I3,KD) ( (8.*(U(I1,I2,I3+1,KD)-U(I1,I2,I3-1,KD)) \
                              -(U(I1,I2,I3+2,KD)-U(I1,I2,I3-2,KD)))*d14(axis3) )
#include "cgux4a.h"    // define 4th order difference approximations

//------------------------------------------------------------------------------------
// This function defines derivatives for all derivative types
//   derivativeType : which derivative to compute
//   ux  : array to put the results in
//   4: order of accuracy 2 or 4
//   .nd   : number of dimensions 2 or 3
// Notes:
//  This macro calls one of the macros UX22, UX23, UX42, UX43 etc,
//  (defined in the include file cgux2a.h or cgux4a.h) to compute a 
//  derivative to first or second order accuracy
//------------------------------------------------------------------------------------


#endif