Remove the "row" function which returned a vector instead of a matrix.

[numerical-analysis.git] / src / Integration / Gaussian.hs

diff --git a/src/Integration/Gaussian.hs b/src/Integration/Gaussian.hs
index 8fcca4f710d4fd918e061f475196ad3c9a2af1e5..be92e5c8bdfd85e232cee5db3b1d57b7ce7675b1 100644 (file)
--- a/src/Integration/Gaussian.hs
+++ b/src/Integration/Gaussian.hs
@@ -28,8 +28,8 @@ import Linear.Matrix (
   construct,
   fromList,
   identity_matrix,
-  matmap,
-  row',
+  map2,
+  row,
   transpose )
 import Linear.QR ( eigenvectors_symmetric )
 import Normed ( Normed(..) )
@@ -119,11 +119,11 @@ nodes_and_weights iterations =
     nodes = shifted_nodes + ones -- unshift the nodes
 
     -- Get the first component of each column.
-    first_components = row' vecs 0
+    first_components = row vecs 0
 
     -- Square it and multiply by 2; see the Golub-Welsch paper for
     -- this magic.
-    weights_row = matmap (\x -> (fromInteger 2)*x^2) first_components
+    weights_row = map2 (\x -> (fromInteger 2)*x^2) first_components
 
     weights = transpose $ weights_row
 
@@ -212,5 +212,5 @@ gaussian' f nodes weights =
   -- want.
   norm_p (1::Int) weighted_values
   where
-    function_values = matmap f nodes
+    function_values = map2 f nodes
     weighted_values = colzipwith (*) weights function_values