@brief Compute A = A + A^T
of a square matrix
@param[inout] a square NxN matrix
@param[in] n matrix dimension
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=dp), | intent(inout) | :: | a(n,*) | |||
integer, | intent(in) | :: | n |
subroutine symmetrize_matrix(a,n) use precision, only: dp real(kind=dp), intent(inout) :: a(n,*) integer, intent(in) :: n integer :: i do i = 1, n a(i:n,i) = a(i:n,i) + a(i,i:n) a(1:i-1,i) = a(i,1:i-1) end do end subroutine symmetrize_matrix