/*@@
  @header   DKDT_declare.h
  @date     Jul 98
  @author   Gabrielle Allen
  @desc
  Declarations for macro @seefile DKDT_guts.h to calculate the source term 
  in the evolution equation for the extrinsic curvature. That is

  d K_ij/dt = alpha*(R_ij-2 K_ik K^k_j + K_ij trK] - nabla_i nabla_j alpha

       + Lie_beta K_ij
 
  @enddesc
@@*/

#ifndef DKDT_DECLARE
#define DKDT_DECLARE

#include "../../packages/CactusEinstein/Einstein/src/macro/RICCI_declare.h"
#include "../../packages/CactusEinstein/Einstein/src/macro/KK_declare.h"
#include "../../packages/CactusEinstein/Einstein/src/macro/TRK_declare.h"
#include "../../packages/CactusEinstein/Einstein/src/macro/CDCDA_declare.h"
#include "../../packages/CactusEinstein/Einstein/src/macro/LIEK_declare.h"


#ifdef FCODE

/* Input variables */
#undef  DKDT_A
#define DKDT_A alp(i,j,k)
#undef  DKDT_KXX 
#define DKDT_KXX hxx(i,j,k)
#undef  DKDT_KXY 
#define DKDT_KXY hxy(i,j,k)
#undef  DKDT_KXZ 
#define DKDT_KXZ hxz(i,j,k)
#undef  DKDT_KYY 
#define DKDT_KYY hyy(i,j,k)
#undef  DKDT_KYZ 
#define DKDT_KYZ hyz(i,j,k)
#undef  DKDT_KZZ 
#define DKDT_KZZ hzz(i,j,k)

/* Output variables */ 
#undef  DKDT_DKXXDT
#define DKDT_DKXXDT dkdt_dkxxdt
#undef  DKDT_DKXYDT
#define DKDT_DKXYDT dkdt_dkxydt
#undef  DKDT_DKXZDT
#define DKDT_DKXZDT dkdt_dkxzdt
#undef  DKDT_DKYYDT
#define DKDT_DKYYDT dkdt_dkyydt
#undef  DKDT_DKYZDT
#define DKDT_DKYZDT dkdt_dkyzdt
#undef  DKDT_DKZZDT
#define DKDT_DKZZDT dkdt_dkzzdt

/* Declare output variables */
      REAL DKDT_DKXXDT
      REAL DKDT_DKXYDT
      REAL DKDT_DKXZDT
      REAL DKDT_DKYYDT
      REAL DKDT_DKYZDT
      REAL DKDT_DKZZDT

#endif


#ifdef CCODE

/* Input variables */
#undef  DKDT_A
#define DKDT_A alp[ijk]
#undef  DKDT_KXX 
#define DKDT_KXX hxx[ijk]
#undef  DKDT_KXY 
#define DKDT_KXY hxy[ijk]
#undef  DKDT_KXZ 
#define DKDT_KXZ hxz[ijk]
#undef  DKDT_KYY 
#define DKDT_KYY hyy[ijk]
#undef  DKDT_KYZ 
#define DKDT_KYZ hyz[ijk]
#undef  DKDT_KZZ 
#define DKDT_KZZ hzz[ijk]

/* Output variables */ 
#undef  DKDT_DKXXDT
#define DKDT_DKXXDT dkdt_dkxxdt
#undef  DKDT_DKXYDT
#define DKDT_DKXYDT dkdt_dkxydt
#undef  DKDT_DKXZDT
#define DKDT_DKXZDT dkdt_dkxzdt
#undef  DKDT_DKYYDT
#define DKDT_DKYYDT dkdt_dkyydt
#undef  DKDT_DKYZDT
#define DKDT_DKYZDT dkdt_dkyzdt
#undef  DKDT_DKZZDT
#define DKDT_DKZZDT dkdt_dkzzdt

/* Declare output variables */
      double DKDT_DKXXDT;
      double DKDT_DKXYDT;
      double DKDT_DKXZDT;
      double DKDT_DKYYDT;
      double DKDT_DKYZDT;
      double DKDT_DKZZDT;

#endif

#endif