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