added generic konst

1 files changed, 4 insertions, 6 deletions

diff --git a/lib/Data/Packed/Random.hs b/lib/Data/Packed/Random.hs
index b30b299..c34f212 100644
--- a/lib/Data/Packed/Random.hs
+++ b/lib/Data/Packed/Random.hs
@@ -22,10 +22,8 @@ module Data.Packed.Random (
 import Numeric.GSL.Vector
 import Data.Packed
 import Numeric.Container
-import Data.Packed.Internal(constantD)
 import Numeric.LinearAlgebra.Algorithms
 
-constant k v = constantD k v
 
 -- | Obtains a matrix whose rows are pseudorandom samples from a multivariate
 -- Gaussian distribution.
@@ -36,7 +34,7 @@ gaussianSample :: Int -- ^ seed
                -> Matrix Double -- ^ result
 gaussianSample seed n med cov = m where
     c = dim med
-    meds = constant 1 n `outer` med
+    meds = konst 1 n `outer` med
     rs = reshape c $ randomVector seed Gaussian (c * n)
     m = rs `mXm` cholSH cov `add` meds
 
@@ -52,7 +50,7 @@ uniformSample seed n rgs = m where
     cs = zipWith subtract as bs
     d = dim a
     dat = toRows $ reshape n $ randomVector seed Uniform (n*d)
-    am = constant 1 n `outer` a
+    am = konst 1 n `outer` a
     m = fromColumns (zipWith scale cs dat) `add` am
 
 ------------ utilities -------------------------------
@@ -62,7 +60,7 @@ meanCov :: Matrix Double -> (Vector Double, Matrix Double)
 meanCov x = (med,cov) where
     r    = rows x
     k    = 1 / fromIntegral r
-    med  = constant k r `vXm` x
-    meds = constant 1 r `outer` med
+    med  = konst k r `vXm` x
+    meds = konst 1 r `outer` med
     xc   = x `sub` meds
     cov  = flip scale (trans xc `mXm` xc) (recip (fromIntegral (r-1)))