hgeometry-0.14: src/Algorithms/Geometry/SoS/Determinant.hs

module Algorithms.Geometry.SoS.Determinant where

import           Algorithms.Geometry.SoS.Sign
import           Algorithms.Geometry.SoS.Symbolic
import           Data.Geometry.Matrix


-- | pre: computes the sign of the determinant
signDet   :: (HasDeterminant d, Ord i, Num r, Ord r) => Matrix d d (Symbolic i r) -> Sign
signDet m = case det m `compare` 0 of
              LT -> Negative
              GT -> Positive
              EQ -> error "signDet: determinant is zero! this should not happen!"