better use transpose
This commit is contained in:
parent
e25b4d6f17
commit
dc57f6aba0
16
src/math.f90
16
src/math.f90
|
@ -97,7 +97,6 @@ module math
|
|||
math_mul3333xx33, &
|
||||
math_mul3333xx3333, &
|
||||
math_exp33 , &
|
||||
math_transpose33, &
|
||||
math_inv33, &
|
||||
math_invert33, &
|
||||
math_invSym3333, &
|
||||
|
@ -631,21 +630,6 @@ pure function math_exp33(A,n)
|
|||
end function math_exp33
|
||||
|
||||
|
||||
!--------------------------------------------------------------------------------------------------
|
||||
!> @brief transposition of a 33 matrix
|
||||
!--------------------------------------------------------------------------------------------------
|
||||
pure function math_transpose33(A)
|
||||
|
||||
implicit none
|
||||
real(pReal),dimension(3,3) :: math_transpose33
|
||||
real(pReal),dimension(3,3),intent(in) :: A
|
||||
integer :: i,j
|
||||
|
||||
forall(i=1:3, j=1:3) math_transpose33(i,j) = A(j,i)
|
||||
|
||||
end function math_transpose33
|
||||
|
||||
|
||||
!--------------------------------------------------------------------------------------------------
|
||||
!> @brief Cramer inversion of 33 matrix (function)
|
||||
!> @details Direct Cramer inversion of matrix A. Returns all zeroes if not possible, i.e.
|
||||
|
|
Loading…
Reference in New Issue