DET33 Function

public pure function DET33(Mat)

compute determinant of 3x3 matrix

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: Mat(3,3)

input matrix

Return Value real(kind=wp)

determinant of the input matrix


Source Code

PURE FUNCTION DET33(Mat)
  ! MODULES
  IMPLICIT NONE
  ! INPUT/OUTPUT VARIABLES ------------------------------------------------------------------------------------------------------!
  REAL(wp),INTENT(IN)  :: Mat(3,3) !! input matrix
  REAL(wp)             :: DET33 !! determinant of the input matrix
  ! LOCAL VARIABLES -------------------------------------------------------------------------------------------------------------!
  ! CODE ------------------------------------------------------------------------------------------------------------------------!
  DET33=   ( Mat(1,1) * Mat(2,2) - Mat(1,2) * Mat(2,1) ) * Mat(3,3) &
           + ( Mat(1,2) * Mat(2,3) - Mat(1,3) * Mat(2,2) ) * Mat(3,1) &
           + ( Mat(1,3) * Mat(2,1) - Mat(1,1) * Mat(2,3) ) * Mat(3,2)
END FUNCTION DET33