1 files changed, 20 insertions, 2 deletions

diff --git a/lib/Numeric/LinearAlgebra/Algorithms.hs b/lib/Numeric/LinearAlgebra/Algorithms.hs
index 1109296..f4b7ee9 100644
--- a/lib/Numeric/LinearAlgebra/Algorithms.hs
+++ b/lib/Numeric/LinearAlgebra/Algorithms.hs
@@ -1,4 +1,4 @@
-{-# OPTIONS_GHC -XFlexibleContexts -XFlexibleInstances #-}
+{-# LANGUAGE FlexibleContexts, FlexibleInstances #-}
 {-# LANGUAGE CPP #-}
 -----------------------------------------------------------------------------
 {- |
@@ -85,7 +85,6 @@ import Numeric.LinearAlgebra.Linear
 import Data.List(foldl1')
 import Data.Array
 
-
 -- | Auxiliary typeclass used to define generic computations for both real and complex matrices.
 class (Normed (Matrix t), Linear Vector t, Linear Matrix t) => Field t where
     svd'         :: Matrix t -> (Matrix t, Vector Double, Matrix t)
@@ -397,6 +396,10 @@ ranksv teps maxdim s = k where
 eps :: Double
 eps =  2.22044604925031e-16
 
+peps :: RealFloat x => x -> x
+peps x = 2.0**(fromIntegral $ 1-floatDigits x)
+
+
 -- | The imaginary unit: @i = 0.0 :+ 1.0@
 i :: Complex Double
 i = 0:+1
@@ -467,6 +470,21 @@ instance Normed (Matrix (Complex Double)) where
     pnorm = pnormCM
 
 -----------------------------------------------------------------------
+-- to be optimized
+
+instance Normed (Vector Float) where
+    pnorm t = pnorm t . double
+
+instance Normed (Vector (Complex Float)) where
+    pnorm t = pnorm t . double
+
+instance Normed (Matrix Float) where
+    pnorm t = pnorm t . double
+
+instance Normed (Matrix (Complex Float)) where
+    pnorm t = pnorm t . double
+
+-----------------------------------------------------------------------
 
 -- | The nullspace of a matrix from its SVD decomposition.
 nullspaceSVD :: Field t