X-Git-Url: http://gitweb.michael.orlitzky.com/?p=numerical-analysis.git;a=blobdiff_plain;f=src%2FLinear%2FSystem.hs;h=2d75f611de1dc85825677ce59e72fe008bdbd21d;hp=85a1dd4b4a619affd6401a5d0e890e19ac007efc;hb=ae914d13235a4582077a5cb2b1edd630d9c6ad62;hpb=6b6bae4206bab66823617e2ba77cdf3e8d3fb752 diff --git a/src/Linear/System.hs b/src/Linear/System.hs index 85a1dd4..2d75f61 100644 --- a/src/Linear/System.hs +++ b/src/Linear/System.hs @@ -2,16 +2,17 @@ {-# LANGUAGE ScopedTypeVariables #-} {-# LANGUAGE TypeFamilies #-} -module Linear.System +module Linear.System ( + backward_substitute, + forward_substitute ) where -import Data.Vector.Fixed (Arity, N1) +import Data.Vector.Fixed ( Arity, N1 ) +import NumericPrelude hiding ( (*), abs ) +import qualified NumericPrelude as NP ( (*) ) +import qualified Algebra.Field as Field ( C ) -import Linear.Matrix - -import NumericPrelude hiding ((*), abs) -import qualified NumericPrelude as NP ((*)) -import qualified Algebra.Field as Field +import Linear.Matrix ( Mat(..), (!!!), construct, transpose ) -- | Solve the system m' * x = b', where m' is upper-triangular. Will @@ -19,6 +20,8 @@ import qualified Algebra.Field as Field -- -- Examples: -- +-- >>> import Linear.Matrix ( Mat2, Mat3, fromList, vec2d, vec3d ) +-- -- >>> let identity = fromList [[1,0,0],[0,1,0],[0,0,1]] :: Mat3 Double -- >>> let b = vec3d (1, 2, 3::Double) -- >>> forward_substitute identity b @@ -69,6 +72,8 @@ forward_substitute m' b' = x' -- -- Examples: -- +-- >>> import Linear.Matrix ( Mat3, fromList, vec3d ) +-- -- >>> let identity = fromList [[1,0,0],[0,1,0],[0,0,1]] :: Mat3 Double -- >>> let b = vec3d (1, 2, 3::Double) -- >>> backward_substitute identity b