Col,
Col10,
Mat(..),
- colzipwith,
construct,
fromList,
identity_matrix,
- matmap,
- row',
- transpose )
+ map2,
+ row,
+ transpose,
+ zipwith2 )
import Linear.QR ( eigenvectors_symmetric )
import Normed ( Normed(..) )
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
-- want.
norm_p (1::Int) weighted_values
where
- function_values = matmap f nodes
- weighted_values = colzipwith (*) weights function_values
+ function_values = map2 f nodes
+ weighted_values = zipwith2 (*) weights function_values