Numeric.HiGHS.LP.Monad

The monadic interface to CLP allows to optimize with respect to multiple objectives, successively.

simple solver with warm restart

data T sh a

Instances

Applicative (T sh)

Defined in Numeric.HiGHS.LP.Monad

Methods

pure :: a -> T sh a

(<*>) :: T sh (a -> b) -> T sh a -> T sh b

liftA2 :: (a -> b -> c) -> T sh a -> T sh b -> T sh c

(*>) :: T sh a -> T sh b -> T sh b

(<*) :: T sh a -> T sh b -> T sh a

Functor (T sh)

Defined in Numeric.HiGHS.LP.Monad

Methods

fmap :: (a -> b) -> T sh a -> T sh b

(<$) :: a -> T sh b -> T sh a

Monad (T sh)

Defined in Numeric.HiGHS.LP.Monad

Methods

(>>=) :: T sh a -> (a -> T sh b) -> T sh b

(>>) :: T sh a -> T sh b -> T sh b

return :: a -> T sh a

run :: (Indexed sh, Index sh ~ ix) => sh -> Bounds ix -> T sh a -> a

solve :: (Eq sh, Indexed sh, Index sh ~ ix) => Method -> Constraints Double ix -> (Direction, Objective sh) -> T sh (Result sh)

Add new constraints to an existing problem and run with a new direction and objective.

>>> case Shape.indexTupleFromShape tripletShape of
      (x,y,z) ->
         fmap (mapSnd Array.toTuple) $ snd $
         LP.run tripletShape []
            (LP.solve LP.simplex
               [[2.*x, 1.*y] <=. 10, [1.*y, (5::Double).*z] <=. 20]
               (LP.Maximize, Array.fromTuple (4,-3,2)
                                 :: Array.Array TripletShape Double))
Just (28.0,(5.0,0.0,4.0))

forAllMethod $ \method ->
   TestLP.forAllOrigin $ \origin ->
   TestLP.forAllProblem origin $ \bounds constrs ->
   QC.forAll (TestLP.genObjective origin) $ \(dir,obj) ->
   case (CLP.solve method bounds constrs (dir,obj),
         LP.run (Array.shape origin) bounds $
            LP.solve method constrs (dir,obj)) of
      ((_, Just (optA,_)), (_, Just (optB,_))) ->
         TestLP.approxReal 0.1 optA optB; _ -> False

forAllMethod $ \method ->
   TestLP.forAllOrigin $ \origin ->
   TestLP.forAllProblem origin $ \bounds constrs ->
   TestLP.forAllObjectives origin $ \objs_ ->
   case TestLP.successiveObjectives origin 0.01 objs_ of
      (dirObj, objs) ->
         either
            (\msg -> QC.counterexample (show msg) False)
            (const $ QC.property True) $
         runSuccessive method (Array.shape origin) bounds (constrs,dirObj) objs

forAllMethod $ \method ->
   TestLP.forAllOrigin $ \origin ->
   TestLP.forAllProblem origin $ \bounds constrs ->
   TestLP.forAllObjectives origin $ \objs_ ->
   let shape = Array.shape origin in
   case TestLP.successiveObjectives origin 0.01 objs_ of
      (dirObj, objs) ->
         approxSuccession 0.01
            (solveSuccessiveWarm method shape bounds (constrs,dirObj) objs)
            (solveSuccessiveGen method shape bounds (constrs,dirObj) objs)

data Direction

Constructors

Minimize

Maximize

Instances

Bounded Direction

Defined in Numeric.LinearProgramming.Common

Methods

minBound :: Direction

maxBound :: Direction

Enum Direction

Defined in Numeric.LinearProgramming.Common

Methods

succ :: Direction -> Direction

pred :: Direction -> Direction

toEnum :: Int -> Direction

fromEnum :: Direction -> Int

enumFrom :: Direction -> [Direction]

enumFromThen :: Direction -> Direction -> [Direction]

enumFromTo :: Direction -> Direction -> [Direction]

enumFromThenTo :: Direction -> Direction -> Direction -> [Direction]

Show Direction

Defined in Numeric.LinearProgramming.Common

Methods

showsPrec :: Int -> Direction -> ShowS

show :: Direction -> String

showList :: [Direction] -> ShowS

Eq Direction

Defined in Numeric.LinearProgramming.Common

Methods

(==) :: Direction -> Direction -> Bool

(/=) :: Direction -> Direction -> Bool

data Method

Instances

Show Method

Defined in Numeric.HiGHS.LP.Private

Methods

showsPrec :: Int -> Method -> ShowS

show :: Method -> String

showList :: [Method] -> ShowS

simplex :: Method

choose :: Method

ipm :: Method

data ModelStatus

Instances

Show ModelStatus

Defined in Numeric.HiGHS.LP.Enumeration

Methods

showsPrec :: Int -> ModelStatus -> ShowS

show :: ModelStatus -> String

showList :: [ModelStatus] -> ShowS

Eq ModelStatus

Defined in Numeric.HiGHS.LP.Enumeration

Methods

(==) :: ModelStatus -> ModelStatus -> Bool

(/=) :: ModelStatus -> ModelStatus -> Bool

type Result sh = (ModelStatus, Maybe (Double, Array sh Double))

solve with extra queries on the result

solveWith :: (Eq sh, Indexed sh, Index sh ~ ix) => Query sh result -> Method -> Constraints Double ix -> (Direction, Objective sh) -> T sh (ModelStatus, Maybe result)

data Query sh a

Instances

Applicative (Query sh)

Defined in Numeric.HiGHS.LP.Private

Methods

pure :: a -> Query sh a

(<*>) :: Query sh (a -> b) -> Query sh a -> Query sh b

liftA2 :: (a -> b -> c) -> Query sh a -> Query sh b -> Query sh c

(*>) :: Query sh a -> Query sh b -> Query sh b

(<*) :: Query sh a -> Query sh b -> Query sh a

Functor (Query sh)

Defined in Numeric.HiGHS.LP.Private

Methods

fmap :: (a -> b) -> Query sh a -> Query sh b

(<$) :: a -> Query sh b -> Query sh a

Monad (Query sh)

Defined in Numeric.HiGHS.LP.Private

Methods

(>>=) :: Query sh a -> (a -> Query sh b) -> Query sh b

(>>) :: Query sh a -> Query sh b -> Query sh b

return :: a -> Query sh a

getObjectiveValue :: Query sh Double

getOptimalVector :: C sh => Query sh (Array sh Double)

getSolutionVectors :: C sh => Query sh ((Array sh Double, Array sh Double), (Array ShapeInt Double, Array ShapeInt Double))

getBasisStatus :: C sh => Query sh (Array sh BasisStatus, Array ShapeInt BasisStatus)

data BasisStatus

Instances

Storable BasisStatus

Defined in Numeric.HiGHS.LP.FFI

Methods

sizeOf :: BasisStatus -> Int

alignment :: BasisStatus -> Int

peekElemOff :: Ptr BasisStatus -> Int -> IO BasisStatus

pokeElemOff :: Ptr BasisStatus -> Int -> BasisStatus -> IO ()

peekByteOff :: Ptr b -> Int -> IO BasisStatus

pokeByteOff :: Ptr b -> Int -> BasisStatus -> IO ()

peek :: Ptr BasisStatus -> IO BasisStatus

poke :: Ptr BasisStatus -> BasisStatus -> IO ()

Show BasisStatus

Defined in Numeric.HiGHS.LP.FFI

Methods

showsPrec :: Int -> BasisStatus -> ShowS

show :: BasisStatus -> String

showList :: [BasisStatus] -> ShowS