Multiply a p-vector by the transpose of an r-matrix.
Status: vector/matrix support routine.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in), | dimension(3,3) | :: | r | r-matrix |
|
real(kind=wp), | intent(in), | dimension(3) | :: | p | p-vector |
|
real(kind=wp), | intent(out), | dimension(3) | :: | trp | R * P |
subroutine TRXP ( r, p, trp )
implicit none
real(wp),dimension(3,3),intent(in) :: r !! r-matrix
real(wp),dimension(3),intent(in) :: p !! p-vector
real(wp),dimension(3),intent(out) :: trp !! R * P
real(wp) :: ri(3,3)
! Transpose of matrix R.
call TR ( r, ri )
! Matrix RI * vector P -> vector TRP.
call RXP ( ri, p, trp )
end subroutine TRXP