BalanceRules.hs

-----------------------------------------------------------------------------
-- Copyright 2015, Open Universiteit Nederland. This file is distributed
-- under the terms of the GNU General Public License. For more information,
-- see the file "LICENSE.txt", which is included in the distribution.
-----------------------------------------------------------------------------
-- |
-- Maintainer  :  bastiaan.heeren@ou.nl
-- Stability   :  provisional
-- Portability :  portable (depends on ghc)
--
-----------------------------------------------------------------------------
--  $Id: BalanceRules.hs 7527 2015-04-08 07:58:06Z bastiaan $

module Domain.Math.Equation.BalanceRules
   ( plusRule, minusRule, timesRule, divisionRule
   ) where

import Control.Monad
import Domain.Math.Data.Relation
import Domain.Math.Expr
import Domain.Math.Numeric.Views
import Ideas.Common.Library

plusRule :: Functor f => ParamTrans Expr (f Expr)
plusRule = parameter1 "term" $ \a ->
   makeTrans $ Just . fmap (:+: a)

minusRule :: Functor f => ParamTrans Expr (f Expr)
minusRule = parameter1 "term" $ \a ->
   makeTrans $ Just . fmap (:-: a)

timesRule :: Functor f => ParamTrans Expr (f Expr)
timesRule = parameter1 "factor" $ \a ->
   makeTrans $ unlessZero a . fmap (a :*:)

divisionRule :: ParamTrans Expr (Equation Expr)
divisionRule = parameter1 "factor" $ \a ->
   makeTrans $ unlessZero a . fmap (:/: a)

unlessZero :: Expr -> a -> Maybe a
unlessZero e a = do
   r <- matchM rationalView e
   guard (r /= 0)
   return a