https://github.com/florentrenaud/nbody6tt
Tip revision: 8aea50c213fd132d500c415511ae1e27eeabab80 authored by florent on 14 February 2015, 16:38:53 UTC
corrected mb typo in ttgalaxy
corrected mb typo in ttgalaxy
Tip revision: 8aea50c
matrix.f
SUBROUTINE MATRIX(U,A)
*
*
* Levi-Civita matrix.
* -------------------
*
REAL*8 U(4),A(3,4)
*
*
* Set current transformation matrix.
A(1,1) = U(1)
A(1,2) = -U(2)
A(1,3) = -U(3)
A(1,4) = U(4)
A(2,1) = U(2)
A(2,2) = U(1)
A(2,3) = -U(4)
A(2,4) = -U(3)
A(3,1) = U(3)
A(3,2) = U(4)
A(3,3) = U(1)
A(3,4) = U(2)
*
RETURN
*
END