Multiply a pv-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,2) | :: | pv | pv-vector |
|
real(kind=wp), | intent(out), | dimension(3,2) | :: | trpv | R * PV |
subroutine TRXPV ( r, pv, trpv )
implicit none
real(wp),dimension(3,3),intent(in) :: r !! r-matrix
real(wp),dimension(3,2),intent(in) :: pv !! pv-vector
real(wp),dimension(3,2),intent(out) :: trpv !! R * PV
real(wp) :: ri(3,3)
! Transpose of matrix R.
call TR ( r, ri )
! Matrix RI * vector PV -> vector TRPV.
call RXPV ( ri, pv, trpv )
end subroutine TRXPV