c--- File written by FORM program dkdrwz.frm on Thu Jul 26 20:15:58 CDT 2012
      subroutine dkdrwz(p1,p2,p3,p4,p5,p6,p7,Fa,Fb,Fc,Fd)
      implicit none
      include 'constants.f'
      include 'sprods_com.f'
      include 'zprods_com.f'
      integer p1,p2,p3,p4,p5,p6,p7
      double complex zab,iza,izb,Fa(2,2),Fb(2,2),Fc(2,2),Fd(2,2)
      double precision s12,s34,s134,s234,s567
      zab(p1,p2,p3,p4)=za(p1,p2)*zb(p2,p4)+za(p1,p3)*zb(p3,p4)
      iza(p1,p2)=cone/za(p1,p2)
      izb(p1,p2)=cone/zb(p1,p2)
      s134=s(p1,p3)+s(p1,p4)+s(p3,p4)
      s234=s(p2,p3)+s(p2,p4)+s(p3,p4)
      s567=s(p5,p6)+s(p5,p7)+s(p6,p7)
      s34=s(p3,p4)
      s12=s(p1,p2)
c--- First label is the quark helicity
c--- Second label is the gluon helicity
      Fa(1,1)= + izb(p5,p7)*s567**(-1)*s34**(-1) * ( za(p2,p3)*za(p2,p7
     &    )*zb(p1,p6)*zb(p2,p4)*s234**(-1) + za(p2,p3)*za(p3,p7)*zb(p1,
     &    p6)*zb(p3,p4)*s234**(-1) )
      Fa(1,1) = Fa(1,1) + izb(p5,p7)*izb(p6,p7)*s567**(-1)*s34**(-1)
     &  * (  - za(p2,p3)*za(p2,p5)*zb(p1,p6)*zb(p2,p4)*zb(p5,p6)*
     &    s234**(-1) - za(p2,p3)*za(p3,p5)*zb(p1,p6)*zb(p3,p4)*zb(p5,p6
     &    )*s234**(-1) )

      Fa(1,2)= + iza(p5,p7)*iza(p6,p7)*s567**(-1)*s34**(-1) * ( za(p2,
     &    p3)*za(p2,p5)*za(p5,p6)*zb(p1,p6)*zb(p2,p4)*s234**(-1) + za(
     &    p2,p3)*za(p3,p5)*za(p5,p6)*zb(p1,p6)*zb(p3,p4)*s234**(-1) )
      Fa(1,2) = Fa(1,2) + iza(p6,p7)*s567**(-1)*s34**(-1) * ( za(p2,p3)
     &    *za(p2,p5)*zb(p1,p7)*zb(p2,p4)*s234**(-1) + za(p2,p3)*za(p3,
     &    p5)*zb(p1,p7)*zb(p3,p4)*s234**(-1) )

      Fa(2,1)= + izb(p5,p7)*izb(p6,p7)*s567**(-1)*s34**(-1) * (  - za(
     &    p2,p3)*za(p2,p6)*zb(p1,p5)*zb(p2,p4)*zb(p5,p6)*s234**(-1) - 
     &    za(p2,p3)*za(p3,p6)*zb(p1,p5)*zb(p3,p4)*zb(p5,p6)*s234**(-1)
     &     )
      Fa(2,1) = Fa(2,1) + izb(p6,p7)*s567**(-1)*s34**(-1) * (  - za(p2,
     &    p3)*za(p2,p7)*zb(p1,p5)*zb(p2,p4)*s234**(-1) - za(p2,p3)*za(
     &    p3,p7)*zb(p1,p5)*zb(p3,p4)*s234**(-1) )

      Fa(2,2)= + iza(p5,p7)*s567**(-1)*s34**(-1) * (  - za(p2,p3)*za(p2
     &    ,p6)*zb(p1,p7)*zb(p2,p4)*s234**(-1) - za(p2,p3)*za(p3,p6)*zb(
     &    p1,p7)*zb(p3,p4)*s234**(-1) )
      Fa(2,2) = Fa(2,2) + iza(p5,p7)*iza(p6,p7)*s567**(-1)*s34**(-1)
     &  * ( za(p2,p3)*za(p2,p6)*za(p5,p6)*zb(p1,p5)*zb(p2,p4)*
     &    s234**(-1) + za(p2,p3)*za(p3,p6)*za(p5,p6)*zb(p1,p5)*zb(p3,p4
     &    )*s234**(-1) )

      Fb(1,1)= + izb(p5,p7)*s567**(-1)*s34**(-1) * (  - za(p1,p3)*za(p2
     &    ,p7)*zb(p1,p4)*zb(p1,p6)*s134**(-1) + za(p2,p7)*za(p3,p4)*zb(
     &    p1,p4)*zb(p4,p6)*s134**(-1) )
      Fb(1,1) = Fb(1,1) + izb(p5,p7)*izb(p6,p7)*s567**(-1)*s34**(-1)
     &  * ( za(p1,p3)*za(p2,p5)*zb(p1,p4)*zb(p1,p6)*zb(p5,p6)*
     &    s134**(-1) - za(p2,p5)*za(p3,p4)*zb(p1,p4)*zb(p4,p6)*zb(p5,p6
     &    )*s134**(-1) )

      Fb(1,2)= + iza(p5,p7)*iza(p6,p7)*s567**(-1)*s34**(-1) * (  - za(
     &    p1,p3)*za(p2,p5)*za(p5,p6)*zb(p1,p4)*zb(p1,p6)*s134**(-1) + 
     &    za(p2,p5)*za(p3,p4)*za(p5,p6)*zb(p1,p4)*zb(p4,p6)*s134**(-1)
     &     )
      Fb(1,2) = Fb(1,2) + iza(p6,p7)*s567**(-1)*s34**(-1) * (  - za(p1,
     &    p3)*za(p2,p5)*zb(p1,p4)*zb(p1,p7)*s134**(-1) + za(p2,p5)*za(
     &    p3,p4)*zb(p1,p4)*zb(p4,p7)*s134**(-1) )

      Fb(2,1)= + izb(p5,p7)*izb(p6,p7)*s567**(-1)*s34**(-1) * ( za(p1,
     &    p3)*za(p2,p6)*zb(p1,p4)*zb(p1,p5)*zb(p5,p6)*s134**(-1) - za(
     &    p2,p6)*za(p3,p4)*zb(p1,p4)*zb(p4,p5)*zb(p5,p6)*s134**(-1) )
      Fb(2,1) = Fb(2,1) + izb(p6,p7)*s567**(-1)*s34**(-1) * ( za(p1,p3)
     &    *za(p2,p7)*zb(p1,p4)*zb(p1,p5)*s134**(-1) - za(p2,p7)*za(p3,
     &    p4)*zb(p1,p4)*zb(p4,p5)*s134**(-1) )

      Fb(2,2)= + iza(p5,p7)*s567**(-1)*s34**(-1) * ( za(p1,p3)*za(p2,p6
     &    )*zb(p1,p4)*zb(p1,p7)*s134**(-1) - za(p2,p6)*za(p3,p4)*zb(p1,
     &    p4)*zb(p4,p7)*s134**(-1) )
      Fb(2,2) = Fb(2,2) + iza(p5,p7)*iza(p6,p7)*s567**(-1)*s34**(-1)
     &  * (  - za(p1,p3)*za(p2,p6)*za(p5,p6)*zb(p1,p4)*zb(p1,p5)*
     &    s134**(-1) + za(p2,p6)*za(p3,p4)*za(p5,p6)*zb(p1,p4)*zb(p4,p5
     &    )*s134**(-1) )

      Fc(1,1)= + izb(p5,p7)*s567**(-1)*s34**(-1)*s12**(-1) * ( za(p2,p3
     &    )*zb(p1,p4)*zab(p7,p1,p2,p6) - za(p2,p7)*zb(p1,p6)*zab(p3,p1,
     &    p2,p4) + za(p3,p7)*zb(p4,p6)*zab(p2,p3,p4,p1) )
      Fc(1,1) = Fc(1,1) + izb(p5,p7)*izb(p6,p7)*s567**(-1)*s34**(-1)*
     & s12**(-1) * (  - za(p2,p3)*zb(p1,p4)*zb(p5,p6)*zab(p5,p1,p2,p6)
     &     + za(p2,p5)*zb(p1,p6)*zb(p5,p6)*zab(p3,p1,p2,p4) - za(p3,p5)
     &    *zb(p4,p6)*zb(p5,p6)*zab(p2,p3,p4,p1) )

      Fc(1,2)= + iza(p5,p7)*iza(p6,p7)*s567**(-1)*s34**(-1)*s12**(-1)
     &  * ( za(p2,p3)*za(p5,p6)*zb(p1,p4)*zab(p5,p1,p2,p6) - za(p2,p5)*
     &    za(p5,p6)*zb(p1,p6)*zab(p3,p1,p2,p4) + za(p3,p5)*za(p5,p6)*
     &    zb(p4,p6)*zab(p2,p3,p4,p1) )
      Fc(1,2) = Fc(1,2) + iza(p6,p7)*s567**(-1)*s34**(-1)*s12**(-1)
     &  * ( za(p2,p3)*zb(p1,p4)*zab(p5,p1,p2,p7) - za(p2,p5)*zb(p1,p7)*
     &    zab(p3,p1,p2,p4) + za(p3,p5)*zb(p4,p7)*zab(p2,p3,p4,p1) )

      Fc(2,1)= + izb(p5,p7)*izb(p6,p7)*s567**(-1)*s34**(-1)*s12**(-1)
     &  * (  - za(p2,p3)*zb(p1,p4)*zb(p5,p6)*zab(p6,p1,p2,p5) + za(p2,
     &    p6)*zb(p1,p5)*zb(p5,p6)*zab(p3,p1,p2,p4) - za(p3,p6)*zb(p4,p5
     &    )*zb(p5,p6)*zab(p2,p3,p4,p1) )
      Fc(2,1) = Fc(2,1) + izb(p6,p7)*s567**(-1)*s34**(-1)*s12**(-1)
     &  * (  - za(p2,p3)*zb(p1,p4)*zab(p7,p1,p2,p5) + za(p2,p7)*zb(p1,
     &    p5)*zab(p3,p1,p2,p4) - za(p3,p7)*zb(p4,p5)*zab(p2,p3,p4,p1) )

      Fc(2,2)= + iza(p5,p7)*s567**(-1)*s34**(-1)*s12**(-1) * (  - za(p2
     &    ,p3)*zb(p1,p4)*zab(p6,p1,p2,p7) + za(p2,p6)*zb(p1,p7)*zab(p3,
     &    p1,p2,p4) - za(p3,p6)*zb(p4,p7)*zab(p2,p3,p4,p1) )
      Fc(2,2) = Fc(2,2) + iza(p5,p7)*iza(p6,p7)*s567**(-1)*s34**(-1)*
     & s12**(-1) * ( za(p2,p3)*za(p5,p6)*zb(p1,p4)*zab(p6,p1,p2,p5) - 
     &    za(p2,p6)*za(p5,p6)*zb(p1,p5)*zab(p3,p1,p2,p4) + za(p3,p6)*
     &    za(p5,p6)*zb(p4,p5)*zab(p2,p3,p4,p1) )

      Fd(:,:)=-Fc(:,:)
      return
      end