/*@@ @header TRKK_guts.h @date Aug 98 @author Gabrielle Allen @desc Macro to calculate the trace of the K_lj K^j_i @enddesc @@*/ #ifndef TRKK_GUTS #define TRKK_GUTS #include "UPPERMET_guts.h" #include "KK_guts.h" #ifdef FCODE TRKK_TRKK = UPPERMET_UXX*KK_KKXX + UPPERMET_UYY*KK_KKYY \ + UPPERMET_UZZ*KK_KKZZ + 2D0*(UPPERMET_UXY*KK_KKXY \ + UPPERMET_UXZ*KK_KKXZ + UPPERMET_UYZ*KK_KKYZ) #endif #ifdef CCODE TRKK_TRKK = UPPERMET_UXX*KK_KKXX + UPPERMET_UYY*KK_KKYY + UPPERMET_UZZ*KK_KKZZ + 2.0*(UPPERMET_UXY*KK_KKXY + UPPERMET_UXZ*KK_KKXZ + UPPERMET_UYZ*KK_KKYZ); #endif #endif