qforce.f

Go to the documentation of this file.

      SUBROUTINE qforce(Q,F,QF)
*
*
*       KS transpose of Q times 2F.
*       ---------------------------
*
      IMPLICIT REAL*8 (a-h,o-z)
      REAL*8 q(4),f(3),qf(4)
*
*
      qf(1)=2.d0*(+q(1)*f(1)+q(2)*f(2)+q(3)*f(3))
      qf(2)=2.d0*(-q(2)*f(1)+q(1)*f(2)+q(4)*f(3))
      qf(3)=2.d0*(-q(3)*f(1)-q(4)*f(2)+q(1)*f(3))
      qf(4)=2.d0*(+q(4)*f(1)-q(3)*f(2)+q(2)*f(3))
*
      RETURN
      END