Fusible (//)
[darcs-mirrors/vector.git] / Data / Vector / MVector.hs
index d519e61..38f0e8d 100644 (file)
 -- Generic interface to mutable vectors
 --
 
+#include "phases.h"
+
 module Data.Vector.MVector (
-  MVector(..),
+  MVectorPure(..), MVector(..),
 
-  slice, new, newWith, read, write, copy, grow, unstream
+  slice, new, newWith, read, write, copy, grow, unstream, update
 ) where
 
 import qualified Data.Vector.Stream      as Stream
@@ -33,56 +35,47 @@ import Prelude hiding ( length, read )
 gROWTH_FACTOR :: Double
 gROWTH_FACTOR = 1.5
 
--- | Class of mutable vectors. The type @m@ is the monad in which the mutable
--- vector can be transformed and @a@ is the type of elements. A vector does
--- not necessarily have to be generic in either of them (indeed, it would be
--- unusual for a vector to be generic in the monad). Use GADTs if this is the
--- case. For instance, regular boxed vectors are defined as
---
--- > data Vector m a where
--- >   Vector :: !Int -> !Int -> MutableArray# s a -> Vector (ST s) a
---
--- This is a bit clumsy but I haven't been able to find a better solution. In
--- particular, using a type function for the monad triggers
--- <http://hackage.haskell.org/trac/ghc/ticket/2440> and is probably less
--- portable.
---
-class Monad m => MVector v m a where
+-- | Basic pure functions on mutable vectors
+class MVectorPure v a where
   -- | Length of the mutable vector
-  length           :: v a -> Int
+  length           :: v a -> Int
 
   -- | Yield a part of the mutable vector without copying it. No range checks!
-  unsafeSlice      :: v m a -> Int  -- ^ starting index
-                            -> Int  -- ^ length of the slice
-                            -> v m a
+  unsafeSlice      :: v a -> Int  -- ^ starting index
+                          -> Int  -- ^ length of the slice
+                          -> v a
+
+  -- Check whether two vectors overlap.
+  overlaps         :: v a -> v a -> Bool
 
+-- | Class of mutable vectors. The type @m@ is the monad in which the mutable
+-- vector can be transformed and @a@ is the type of elements.
+--
+class (Monad m, MVectorPure v a) => MVector v m a where
   -- | Create a mutable vector of the given length. Length is not checked!
-  unsafeNew        :: Int -> m (v a)
+  unsafeNew        :: Int -> m (v a)
 
   -- | Create a mutable vector of the given length and fill it with an
   -- initial value. Length is not checked!
-  unsafeNewWith    :: Int -> a -> m (v a)
+  unsafeNewWith    :: Int -> a -> m (v a)
 
   -- | Yield the element at the given position. Index is not checked!
-  unsafeRead       :: v a -> Int -> m a
+  unsafeRead       :: v a -> Int -> m a
 
   -- | Replace the element at the given position. Index is not checked!
-  unsafeWrite      :: v a -> Int -> a -> m ()
+  unsafeWrite      :: v a -> Int -> a -> m ()
 
   -- | Write the value at each position.
-  set              :: v a -> a -> m ()
+  set              :: v a -> a -> m ()
 
   -- | Copy a vector. The two vectors may not overlap. This is not checked!
-  unsafeCopy       :: v a   -- ^ target
-                   -> v a   -- ^ source
+  unsafeCopy       :: v a   -- ^ target
+                   -> v a   -- ^ source
                    -> m ()
 
   -- | Grow a vector by the given number of elements. The length is not
   -- checked!
-  unsafeGrow       :: v m a -> Int -> m (v m a)
-
-  -- Check whether two vectors overlap.
-  overlaps         :: v m a -> v m a -> Bool
+  unsafeGrow       :: v a -> Int -> m (v a)
 
   {-# INLINE unsafeNewWith #-}
   unsafeNewWith n x = do
@@ -120,50 +113,50 @@ class Monad m => MVector v m a where
       n = length v
 
 -- | Test whether the index is valid for the vector
-inBounds :: MVector v m a => v m a -> Int -> Bool
+inBounds :: MVectorPure v a => v a -> Int -> Bool
 {-# INLINE inBounds #-}
 inBounds v i = i >= 0 && i < length v
 
 -- | Yield a part of the mutable vector without copying it. Safer version of
 -- 'unsafeSlice'.
-slice :: MVector v m a => v m a -> Int -> Int -> v m a
+slice :: MVectorPure v a => v a -> Int -> Int -> v a
 {-# INLINE slice #-}
 slice v i n = assert (i >=0 && n >= 0 && i+n <= length v)
             $ unsafeSlice v i n
 
 -- | Create a mutable vector of the given length. Safer version of
 -- 'unsafeNew'.
-new :: MVector v m a => Int -> m (v a)
+new :: MVector v m a => Int -> m (v a)
 {-# INLINE new #-}
 new n = assert (n >= 0) $ unsafeNew n
 
 -- | Create a mutable vector of the given length and fill it with an
 -- initial value. Safer version of 'unsafeNewWith'.
-newWith :: MVector v m a => Int -> a -> m (v a)
+newWith :: MVector v m a => Int -> a -> m (v a)
 {-# INLINE newWith #-}
 newWith n x = assert (n >= 0) $ unsafeNewWith n x
 
 -- | Yield the element at the given position. Safer version of 'unsafeRead'.
-read :: MVector v m a => v a -> Int -> m a
+read :: MVector v m a => v a -> Int -> m a
 {-# INLINE read #-}
 read v i = assert (inBounds v i) $ unsafeRead v i
 
 -- | Replace the element at the given position. Safer version of
 -- 'unsafeWrite'.
-write :: MVector v m a => v a -> Int -> a -> m ()
+write :: MVector v m a => v a -> Int -> a -> m ()
 {-# INLINE write #-}
 write v i x = assert (inBounds v i) $ unsafeWrite v i x
 
 -- | Copy a vector. The two vectors may not overlap. Safer version of
 -- 'unsafeCopy'.
-copy :: MVector v m a => v m a -> v m a -> m ()
+copy :: MVector v m a => v a -> v a -> m ()
 {-# INLINE copy #-}
 copy dst src = assert (not (dst `overlaps` src) && length dst == length src)
              $ unsafeCopy dst src
 
 -- | Grow a vector by the given number of elements. Safer version of
 -- 'unsafeGrow'.
-grow :: MVector v m a => v m a -> Int -> m (v m a)
+grow :: MVector v m a => v a -> Int -> m (v a)
 {-# INLINE grow #-}
 grow v by = assert (by >= 0)
           $ unsafeGrow v by
@@ -172,13 +165,13 @@ grow v by = assert (by >= 0)
 -- | Create a new mutable vector and fill it with elements from the 'Stream'.
 -- The vector will grow logarithmically if the 'Size' hint of the 'Stream' is
 -- inexact.
-unstream :: MVector v m a => Stream a -> m (v a)
-{-# INLINE unstream #-}
+unstream :: MVector v m a => Stream a -> m (v a)
+{-# INLINE_STREAM unstream #-}
 unstream s = case upperBound (Stream.size s) of
                Just n  -> unstreamMax     s n
                Nothing -> unstreamUnknown s
 
-unstreamMax :: MVector v m a => Stream a -> Int -> m (v a)
+unstreamMax :: MVector v m a => Stream a -> Int -> m (v a)
 {-# INLINE unstreamMax #-}
 unstreamMax s n
   = do
@@ -187,7 +180,7 @@ unstreamMax s n
       n' <- Stream.foldM put 0 s
       return $ slice v 0 n'
 
-unstreamUnknown :: MVector v m a => Stream a -> m (v a)
+unstreamUnknown :: MVector v m a => Stream a -> m (v a)
 {-# INLINE unstreamUnknown #-}
 unstreamUnknown s
   = do
@@ -208,3 +201,11 @@ unstreamUnknown s
                                  . double2Int
                                  $ int2Double (length v) * gROWTH_FACTOR
 
+
+update :: MVector v m a => v a -> Stream (Int, a) -> m ()
+{-# INLINE update #-}
+update v s = Stream.mapM_ put s
+  where
+    {-# INLINE put #-}
+    put (i, x) = write v i x
+