Numeric.LAPACK.Example.EconomicAllocation

Do not import this module. It is only for demonstration purposes.

Documentation

>>> import Numeric.LAPACK.Example.EconomicAllocation
>>> import Test.Utility (approxVector)
>>> 
>>> import qualified Numeric.LAPACK.Vector as Vector
>>> import Numeric.LAPACK.Vector ((+++))
>>> 
>>> import qualified Data.Array.Comfort.Storable as Array

type ZeroInt2 = ShapeInt ::+ ShapeInt

type Vector sh = Vector sh Double

type Matrix height width = General height width Double

type SquareMatrix size = Square size Double

balances0 :: Vector ZeroInt2

expenses0 :: Matrix ShapeInt ZeroInt2

normalize :: (Eq height, C height, C width) => Matrix height width -> Matrix height width

normalizeSplit :: (C sh0, C sh1, Eq sh1) => Matrix sh1 (sh0 ::+ sh1) -> (Matrix sh0 sh1, SquareMatrix sh1)

completeIdSquare :: (C sh0, Eq sh0, C sh1, Eq sh1) => Matrix sh1 (sh0 ::+ sh1) -> SquareMatrix (sh0 ::+ sh1)

iterated :: (C sh0, Eq sh0, C sh1, Eq sh1) => Matrix sh1 (sh0 ::+ sh1) -> Vector (sh0 ::+ sh1) -> Vector (sh0 ::+ sh1)

compensated :: (C sh0, Eq sh0, C sh1, Eq sh1) => Matrix sh1 (sh0 ::+ sh1) -> Vector (sh0 ::+ sh1) -> Vector sh0

main :: IO ()

let result = iterated expenses0 balances0 in approxVector result $ compensated expenses0 balances0 +++ Vector.zero (Array.shape $ Vector.takeRight result)