38 const SquareMatrix<Type>& matrix,
42 Type diagProduct = pTraits<Type>::one;
44 for (
label i = 0; i < matrix.m(); ++i)
46 diagProduct *= matrix(i, i);
49 return sign*diagProduct;
54 Foam::scalar
Foam::det(
const SquareMatrix<Type>& matrix)
56 SquareMatrix<Type> matrixTmp = matrix;
67 Foam::scalar
Foam::det(SquareMatrix<Type>& matrix)
80 template<
class AnyType>
85 for (
label i=0; i < this->
n(); ++i)