ksphys.f

Go to the documentation of this file.

      SUBROUTINE ksphys(Q,P,XR,PR)
*
*
*       Transformation from KS to physical variables.
*       ---------------------------------------------
*
      IMPLICIT REAL*8 (a-h,o-z)
      REAL*8 q(4),p(4),xr(3),pr(3)
*
*
*       Obtain physical coordinates & momenta from KS.
      xr(1)=q(1)**2-q(2)**2-q(3)**2+q(4)**2
      xr(2)=2.d0*(q(1)*q(2)-q(3)*q(4))
      xr(3)=2.d0*(q(1)*q(3)+q(2)*q(4))
*
      r=q(1)**2+q(2)**2+q(3)**2+q(4)**2
      a=0.5d0/r
      pr(1)=(q(1)*p(1)-q(2)*p(2)-q(3)*p(3)+q(4)*p(4))*a
      pr(2)=(q(2)*p(1)+q(1)*p(2)-q(4)*p(3)-q(3)*p(4))*a
      pr(3)=(q(3)*p(1)+q(4)*p(2)+q(1)*p(3)+q(2)*p(4))*a
*
      RETURN
      END