Return in b int inverse of a. Uses utility_inv3
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=dp), | intent(in) | :: | a(3,3) | |||
| real(kind=dp), | intent(out) | :: | b(3,3) |
subroutine utility_inverse_mat(a, b) !================================================! ! !! Return in b int inverse of a. Uses utility_inv3 ! !================================================ implicit none real(kind=dp), intent(in) :: a(3, 3) real(kind=dp), intent(out) :: b(3, 3) real(kind=dp) :: det call utility_inv3(a, b, det) b = b/det return end subroutine utility_inverse_mat