vector.f

Go to the documentation of this file.

      SUBROUTINE vector(Q,P,W)
*
*
*       KS matrix of Q times P.
*       -----------------------
*
      IMPLICIT REAL*8 (a-h,o-z)
      REAL*8 q(4),p(4),w(4)
*
*
*       Evaluate W = L(Q)P; the KS matrix of Q times P.
      w(1)=(+q(1)*p(1)+q(2)*p(2)+q(3)*p(3)+q(4)*p(4))
      w(2)=(-q(2)*p(1)+q(1)*p(2)+q(4)*p(3)-q(3)*p(4))
      w(3)=(-q(3)*p(1)-q(4)*p(2)+q(1)*p(3)+q(2)*p(4))
      w(4)=(+q(4)*p(1)-q(3)*p(2)+q(2)*p(3)-q(1)*p(4))
*
      RETURN
      END