invert_4x4 Subroutine

private subroutine invert_4x4(a, ainv)

Inverse of a 4x4 matrix via LAPACK LU (dgetrf) + inversion (dgetri).

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp),	intent(in)			::	a(4,4)
real(kind=wp),	intent(out)			::	ainv(4,4)