fromList,
identity_matrix,
map2,
- row',
+ row,
transpose )
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.