{-# LANGUAGE DerivingStrategies #-}
{-# LANGUAGE FlexibleInstances #-}
{-# LANGUAGE GeneralizedNewtypeDeriving #-}
{-# LANGUAGE MultiParamTypeClasses #-}

-- |
-- Slower than Control.Monad.Bayes.Density.Free, so not used by default,
-- but more elementary to understand. Just uses standard
-- monad transformer techniques.
module Control.Monad.Bayes.Density.State where

import Control.Monad.Bayes.Class (MonadDistribution (random))
import Control.Monad.State (MonadState (get, put), StateT, evalStateT)
import Control.Monad.Writer

newtype DensityT m a = DensityT {forall (m :: * -> *) a.
DensityT m a -> WriterT [Double] (StateT [Double] m) a
getDensityT :: WriterT [Double] (StateT [Double] m) a} deriving newtype (forall a b. a -> DensityT m b -> DensityT m a
forall a b. (a -> b) -> DensityT m a -> DensityT m b
forall (m :: * -> *) a b.
Functor m =>
a -> DensityT m b -> DensityT m a
forall (m :: * -> *) a b.
Functor m =>
(a -> b) -> DensityT m a -> DensityT m b
forall (f :: * -> *).
(forall a b. (a -> b) -> f a -> f b)
-> (forall a b. a -> f b -> f a) -> Functor f
<$ :: forall a b. a -> DensityT m b -> DensityT m a
$c<$ :: forall (m :: * -> *) a b.
Functor m =>
a -> DensityT m b -> DensityT m a
fmap :: forall a b. (a -> b) -> DensityT m a -> DensityT m b
$cfmap :: forall (m :: * -> *) a b.
Functor m =>
(a -> b) -> DensityT m a -> DensityT m b
Functor, forall a. a -> DensityT m a
forall a b. DensityT m a -> DensityT m b -> DensityT m a
forall a b. DensityT m a -> DensityT m b -> DensityT m b
forall a b. DensityT m (a -> b) -> DensityT m a -> DensityT m b
forall a b c.
(a -> b -> c) -> DensityT m a -> DensityT m b -> DensityT m c
forall {m :: * -> *}. Monad m => Functor (DensityT m)
forall (m :: * -> *) a. Monad m => a -> DensityT m a
forall (m :: * -> *) a b.
Monad m =>
DensityT m a -> DensityT m b -> DensityT m a
forall (m :: * -> *) a b.
Monad m =>
DensityT m a -> DensityT m b -> DensityT m b
forall (m :: * -> *) a b.
Monad m =>
DensityT m (a -> b) -> DensityT m a -> DensityT m b
forall (m :: * -> *) a b c.
Monad m =>
(a -> b -> c) -> DensityT m a -> DensityT m b -> DensityT m c
forall (f :: * -> *).
Functor f
-> (forall a. a -> f a)
-> (forall a b. f (a -> b) -> f a -> f b)
-> (forall a b c. (a -> b -> c) -> f a -> f b -> f c)
-> (forall a b. f a -> f b -> f b)
-> (forall a b. f a -> f b -> f a)
-> Applicative f
<* :: forall a b. DensityT m a -> DensityT m b -> DensityT m a
$c<* :: forall (m :: * -> *) a b.
Monad m =>
DensityT m a -> DensityT m b -> DensityT m a
*> :: forall a b. DensityT m a -> DensityT m b -> DensityT m b
$c*> :: forall (m :: * -> *) a b.
Monad m =>
DensityT m a -> DensityT m b -> DensityT m b
liftA2 :: forall a b c.
(a -> b -> c) -> DensityT m a -> DensityT m b -> DensityT m c
$cliftA2 :: forall (m :: * -> *) a b c.
Monad m =>
(a -> b -> c) -> DensityT m a -> DensityT m b -> DensityT m c
<*> :: forall a b. DensityT m (a -> b) -> DensityT m a -> DensityT m b
$c<*> :: forall (m :: * -> *) a b.
Monad m =>
DensityT m (a -> b) -> DensityT m a -> DensityT m b
pure :: forall a. a -> DensityT m a
$cpure :: forall (m :: * -> *) a. Monad m => a -> DensityT m a
Applicative, forall a. a -> DensityT m a
forall a b. DensityT m a -> DensityT m b -> DensityT m b
forall a b. DensityT m a -> (a -> DensityT m b) -> DensityT m b
forall (m :: * -> *). Monad m => Applicative (DensityT m)
forall (m :: * -> *) a. Monad m => a -> DensityT m a
forall (m :: * -> *) a b.
Monad m =>
DensityT m a -> DensityT m b -> DensityT m b
forall (m :: * -> *) a b.
Monad m =>
DensityT m a -> (a -> DensityT m b) -> DensityT m b
forall (m :: * -> *).
Applicative m
-> (forall a b. m a -> (a -> m b) -> m b)
-> (forall a b. m a -> m b -> m b)
-> (forall a. a -> m a)
-> Monad m
return :: forall a. a -> DensityT m a
$creturn :: forall (m :: * -> *) a. Monad m => a -> DensityT m a
>> :: forall a b. DensityT m a -> DensityT m b -> DensityT m b
$c>> :: forall (m :: * -> *) a b.
Monad m =>
DensityT m a -> DensityT m b -> DensityT m b
>>= :: forall a b. DensityT m a -> (a -> DensityT m b) -> DensityT m b
$c>>= :: forall (m :: * -> *) a b.
Monad m =>
DensityT m a -> (a -> DensityT m b) -> DensityT m b
Monad)

instance MonadTrans DensityT where
  lift :: forall (m :: * -> *) a. Monad m => m a -> DensityT m a
lift = forall (m :: * -> *) a.
WriterT [Double] (StateT [Double] m) a -> DensityT m a
DensityT forall b c a. (b -> c) -> (a -> b) -> a -> c
. forall (t :: (* -> *) -> * -> *) (m :: * -> *) a.
(MonadTrans t, Monad m) =>
m a -> t m a
lift forall b c a. (b -> c) -> (a -> b) -> a -> c
. forall (t :: (* -> *) -> * -> *) (m :: * -> *) a.
(MonadTrans t, Monad m) =>
m a -> t m a
lift

instance (Monad m) => MonadState [Double] (DensityT m) where
  get :: DensityT m [Double]
get = forall (m :: * -> *) a.
WriterT [Double] (StateT [Double] m) a -> DensityT m a
DensityT forall a b. (a -> b) -> a -> b
$ forall (t :: (* -> *) -> * -> *) (m :: * -> *) a.
(MonadTrans t, Monad m) =>
m a -> t m a
lift forall a b. (a -> b) -> a -> b
$ forall s (m :: * -> *). MonadState s m => m s
get
  put :: [Double] -> DensityT m ()
put = forall (m :: * -> *) a.
WriterT [Double] (StateT [Double] m) a -> DensityT m a
DensityT forall b c a. (b -> c) -> (a -> b) -> a -> c
. forall (t :: (* -> *) -> * -> *) (m :: * -> *) a.
(MonadTrans t, Monad m) =>
m a -> t m a
lift forall b c a. (b -> c) -> (a -> b) -> a -> c
. forall s (m :: * -> *). MonadState s m => s -> m ()
put

instance (Monad m) => MonadWriter [Double] (DensityT m) where
  tell :: [Double] -> DensityT m ()
tell = forall (m :: * -> *) a.
WriterT [Double] (StateT [Double] m) a -> DensityT m a
DensityT forall b c a. (b -> c) -> (a -> b) -> a -> c
. forall w (m :: * -> *). MonadWriter w m => w -> m ()
tell
  listen :: forall a. DensityT m a -> DensityT m (a, [Double])
listen = forall (m :: * -> *) a.
WriterT [Double] (StateT [Double] m) a -> DensityT m a
DensityT forall b c a. (b -> c) -> (a -> b) -> a -> c
. forall w (m :: * -> *) a. MonadWriter w m => m a -> m (a, w)
listen forall b c a. (b -> c) -> (a -> b) -> a -> c
. forall (m :: * -> *) a.
DensityT m a -> WriterT [Double] (StateT [Double] m) a
getDensityT
  pass :: forall a. DensityT m (a, [Double] -> [Double]) -> DensityT m a
pass = forall (m :: * -> *) a.
WriterT [Double] (StateT [Double] m) a -> DensityT m a
DensityT forall b c a. (b -> c) -> (a -> b) -> a -> c
. forall w (m :: * -> *) a. MonadWriter w m => m (a, w -> w) -> m a
pass forall b c a. (b -> c) -> (a -> b) -> a -> c
. forall (m :: * -> *) a.
DensityT m a -> WriterT [Double] (StateT [Double] m) a
getDensityT

instance (MonadDistribution m) => MonadDistribution (DensityT m) where
  random :: DensityT m Double
random = do
    [Double]
trace <- forall s (m :: * -> *). MonadState s m => m s
get
    Double
x <- case [Double]
trace of
      [] -> forall (m :: * -> *). MonadDistribution m => m Double
random
      Double
r : [Double]
xs -> forall s (m :: * -> *). MonadState s m => s -> m ()
put [Double]
xs forall (m :: * -> *) a b. Monad m => m a -> m b -> m b
>> forall (f :: * -> *) a. Applicative f => a -> f a
pure Double
r
    forall w (m :: * -> *). MonadWriter w m => w -> m ()
tell [Double
x]
    forall (f :: * -> *) a. Applicative f => a -> f a
pure Double
x

runDensityT :: (Monad m) => DensityT m b -> [Double] -> m (b, [Double])
runDensityT :: forall (m :: * -> *) b.
Monad m =>
DensityT m b -> [Double] -> m (b, [Double])
runDensityT (DensityT WriterT [Double] (StateT [Double] m) b
m) = forall (m :: * -> *) s a. Monad m => StateT s m a -> s -> m a
evalStateT (forall w (m :: * -> *) a. WriterT w m a -> m (a, w)
runWriterT WriterT [Double] (StateT [Double] m) b
m)