/*@@
  @header   DYDK_declare.h
  @date     Jul 98
  @author   Gabrielle Allen
  @desc
  Declarations for macro @seefile DYDG_guts.h to compute first 
  derivatives of the extrinsic curvature with respect to y
  @enddesc
@@*/

#ifndef DYDK_DECLARE

#ifdef FCODE

/* Input variables */
#undef  DYDK_KXX_JP 
#define DYDK_KXX_JP hxx(i,j+1,k)
#undef  DYDK_KXX_JM
#define DYDK_KXX_JM hxx(i,j-1,k)
#undef  DYDK_KXY_JP 
#define DYDK_KXY_JP hxy(i,j+1,k)
#undef  DYDK_KXY_JM 
#define DYDK_KXY_JM hxy(i,j-1,k)
#undef  DYDK_KXZ_JP 
#define DYDK_KXZ_JP hxz(i,j+1,k)
#undef  DYDK_KXZ_JM
#define DYDK_KXZ_JM hxz(i,j-1,k)
#undef  DYDK_KYY_JP 
#define DYDK_KYY_JP hyy(i,j+1,k)
#undef  DYDK_KYY_JM 
#define DYDK_KYY_JM hyy(i,j-1,k)
#undef  DYDK_KYZ_JP 
#define DYDK_KYZ_JP hyz(i,j+1,k)
#undef  DYDK_KYZ_JM 
#define DYDK_KYZ_JM hyz(i,j-1,k)
#undef  DYDK_KZZ_JP 
#define DYDK_KZZ_JP hzz(i,j+1,k)
#undef  DYDK_KZZ_JM
#define DYDK_KZZ_JM hzz(i,j-1,k)

/* Output variables */ 
#undef  DYDK_DYDKXX
#define DYDK_DYDKXX dydk_dydkxx
#undef  DYDK_DYDKXY
#define DYDK_DYDKXY dydk_dydkxy
#undef  DYDK_DYDKXZ
#define DYDK_DYDKXZ dydk_dydkxz
#undef  DYDK_DYDKYY
#define DYDK_DYDKYY dydk_dydkyy
#undef  DYDK_DYDKYZ
#define DYDK_DYDKYZ dydk_dydkyz
#undef  DYDK_DYDKZZ
#define DYDK_DYDKZZ dydk_dydkzz 

/* Internal variables */
#undef  DYDK_DY
#define DYDK_DY    dy
#undef  DYDK_OO2DY   
#define DYDK_OO2DY dydk_oo2dy

/* Declare internal variables */
      REAL DYDK_OO2DY

/* Declare output variables */
      REAL DYDK_DYDKXX
      REAL DYDK_DYDKXY
      REAL DYDK_DYDKXZ
      REAL DYDK_DYDKYY
      REAL DYDK_DYDKYZ
      REAL DYDK_DYDKZZ

#endif


#ifdef CCODE

/* Input variables */
#undef  DYDK_KXX_JP 
#define DYDK_KXX_JP hxx[ dj+ijk]
#undef  DYDK_KXX_JM
#define DYDK_KXX_JM hxx[-dj+ijk]
#undef  DYDK_KXY_JP 
#define DYDK_KXY_JP hxy[ dj+ijk]
#undef  DYDK_KXY_JM 
#define DYDK_KXY_JM hxy[-dj+ijk]
#undef  DYDK_KXZ_JP 
#define DYDK_KXZ_JP hxz[ dj+ijk]
#undef  DYDK_KXZ_JM
#define DYDK_KXZ_JM hxz[-dj+ijk]
#undef  DYDK_KYY_JP 
#define DYDK_KYY_JP hyy[ dj+ijk]
#undef  DYDK_KYY_JM 
#define DYDK_KYY_JM hyy[-dj+ijk]
#undef  DYDK_KYZ_JP 
#define DYDK_KYZ_JP hyz[ dj+ijk]
#undef  DYDK_KYZ_JM 
#define DYDK_KYZ_JM hyz[-dj+ijk]
#undef  DYDK_KZZ_JP 
#define DYDK_KZZ_JP hzz[ dj+ijk]
#undef  DYDK_KZZ_JM
#define DYDK_KZZ_JM hzz[-dj+ijk]

/* Output variables */ 
#undef  DYDK_DYDKXX
#define DYDK_DYDKXX dydk_dydkxx
#undef  DYDK_DYDKXY
#define DYDK_DYDKXY dydk_dydkxy
#undef  DYDK_DYDKXZ
#define DYDK_DYDKXZ dydk_dydkxz
#undef  DYDK_DYDKYY
#define DYDK_DYDKYY dydk_dydkyy
#undef  DYDK_DYDKYZ
#define DYDK_DYDKYZ dydk_dydkyz
#undef  DYDK_DYDKZZ
#define DYDK_DYDKZZ dydk_dydkzz 

/* Internal variables */
#undef  DYDK_OO2DY   
#define DYDK_OO2DY dydk_oo2dy

/* Declare internal variables */
double DYDK_OO2DY;

/* Declare output variables */
double DYDK_DYDKXX;
double DYDK_DYDKXY;
double DYDK_DYDKXZ;
double DYDK_DYDKYY;
double DYDK_DYDKYZ;
double DYDK_DYDKZZ;

#endif

#define DYDK_DECLARE

#endif