X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FIntegration%2FGaussian.hs;fp=src%2FIntegration%2FGaussian.hs;h=be92e5c8bdfd85e232cee5db3b1d57b7ce7675b1;hb=9ef91ad4ec3a5c0966f0850d40310722b6c38b68;hp=38214a0a5be7468b0b061e484ca0bb16d01556ff;hpb=a96be0c7818e0408b4aac9076df131c5b3ec6ff4;p=numerical-analysis.git diff --git a/src/Integration/Gaussian.hs b/src/Integration/Gaussian.hs index 38214a0..be92e5c 100644 --- a/src/Integration/Gaussian.hs +++ b/src/Integration/Gaussian.hs @@ -29,7 +29,7 @@ import Linear.Matrix ( fromList, identity_matrix, map2, - row', + row, transpose ) import Linear.QR ( eigenvectors_symmetric ) import Normed ( Normed(..) ) @@ -119,7 +119,7 @@ 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.