Fix a bug in the relative error.

Remove some commented code from FEM.R1.

Add relative error functions to FEM.R1.

Add natural numbers up to N64.

Speed up the Cholesky factorization with a fold.

Add column synonyms up to Col32.

Add natural numbers up to N32.

Fix bugs in big_K, big_M.

Add Mat6 and Mat7 synonyms.
Move the example FEM into its own FEM.R1.Exampl1 module.
Use the new example FEM module in the FEM.R1 doctests.

Remove solution', because what did it do?

Fix a few warnings.

Commit whatever I was working on last and have forgotten about.

Add a FEM.R1 example (should be moved into the tree after being cleaned up).

Finish the FEM.R1 solution for Dirichlet boundary conditions.

Require a matrix to be nonempty if we're going to QR factorize it.

Generalize the signature of zipwith2.
Implement the dot product using zipwith2.
Implement matrix multiplication using the dot product.
Implement the frobenius_norm with zipwith2.

Use the element_sum2 function to sum the terms in Gaussian integration instead of the incorrect 1-norm.

Fix the type of zip2 to be more general.

Use the Col type synonym in the column Normed instance.

Fix the column matrix 1-norm (it was missing the absolute value).

Add the list_plot' function to Piecewise.

Add 'coefficients' to FEM.R1 which solves the desired system.

Add a preliminary Piecewise module.

Define Col6...Col9 types in Linear.Matrix.

Export everything from Z to N20 from Naturals.

Fix the incorrect definition of "minor".
Add a "preminor" function that does what "minor" used to.
Implement cofactors and a matrix inverse in terms of them.

Add big_K code to FEM.R1.

Use a better implementation of backwards_substitute.

Add Linear.Matrix.ifoldr2.

Implement forward substitute in terms of a fold.

Add Linear.Matrix.set_idx.
Loosen constraints on the dot product function.

Update Integration.Gaussian to use zip2.

Generalize colzip and colzipwith to any matrices.
Add zip2three to zip three matrices together.

Fix backward_substitute.
Add 'solve_positive_definite' to Linear.System.
Add tests for everything in Linear.System.

Add reverse2 to Linear.Matrix.

Add Eq constraints throughout Linear.Iteration.

Add another test for 'cholesky'.
Improve implementation of upper/lower triangular tests.
Drop some typeclass constraints on upper/lower triangular tests.

Add the imap2 function to Linear.Matrix.

Make safe indexing slower but not dependent on the number of rows/columns.

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

Drop the 'column' function that returned a vector instead of a matrix.

Rename matmap to map2.

More progress on the 1d FEM, the "big F" vector is now computed correctly.

Add the 'ifoldl2' function to Linear.Matrix.

Fix a typo.

Add a bunch more natural numbers in the Naturals module.

Add the (in progress) FEM.R1 module.

Fix some docs in Roots.Simple.

Implement Gaussian quadrature and test it.

Rename zipcol to colzip.
Add a new function, colzipwith.

Add the Integration.Gaussian.nodes_and_weights function.

Add zipcol and matmap functions in Linear.Matrix.

Add the Integration.Gaussian module where we begin to implement Gaussian quadrature.

Add comments to Linear.QR about non-convergence.

Finish the eigenvectors_symmetric tests.

Add the eigenvectors_symmetric function and begin writing tests for it.

Implement the QR algorithm for computing eigenvalues.

Add type synonyms for column/row matrices.
Add row' and column' functions that do what you'd expect only they return Mats instead of Vecs like their un-prime counterparts.

Small cleanups in Linear.Vector.
Add 'diagonal' and 'trace' functions to Linear.Matrix.

Clean up imports everywhere.

Add tolerant versions of is_{upper,lower}_triangular.
Fix the QR factorization code.
Add tests for the QR behavior.

Add a Normed instance for (Vec n a).

Switch Linear.Vector to the numeric prelude and add the element_sum function to it.
Use Linear.Vector.element_sum in Linear.Matrix.

Remove old comment from Normed.

Add the Frobenius norm to Linear.Matrix.

Add the Linear.QR module and update the ghci/cabal files.

New function: Linear.Matrix.identity_matrix.

Fix implementation of Linear.Matrix.construct.

Add the Polynomials.Orthogonal module.

Bump dependencies.

Add the rayleigh_quotient function.

Add Gauss-Seidel and SOR iterative methods.

Rename 'diagonal' to 'diagonal_part'.
Add functions to pull out the (strict) lower/upper triangular parts of a matrix.

Fix a bunch of hlint warnings.

Update ghci/cabal files to load the new module.

Add the Linear.Iteration module.

Add another example for forward_substitute.

Add 'diagonal' function.

Constrain Normed instance to column vectors (m-by-1 matrices).

Add the unit_roundoff function.

Update numeric-prelude and fixed-vector.
Use the (!) function from the new fixed-vector.
Add the Fast/Simple trisection method and a test.

Get determinants working.

Replace the whole Matrix implementation with something a little better.
Add the Linear.System module.
Get all doctests working again.
Update .ghci/cabal files.

Add Vec1, Vec4, Vec5 aliases.

Add triangular functions, determinant, minor (all half-baked) to Linear.Matrix.

Add the Linear.System module.

Add the delete function to Vector.hs.

Fix compiler warnings and doctests.

Convert scalars to 1x1 vectors.

A huge pile of crap upon Matrix/Vector.

Update .ghci.

Remove assumptions on the Normed class.

Convert Integration/Simpson.hs and Integration/Trapezoid.hs to numeric-prelude.

Use RealField/RealRing where possible instead of their constituents.
Convert src/ODE/IVP.hs to numeric-prelude.
Add RealRing/RealField instanced to BigFloat.

Rename Aliases.hs to BigFloat.hs, now containing numeric-prelude instances for BigFloats.
Rewrite Normed.hs, Roots/Fast.hs, and Roots/Simple.hs for numeric-prelude.

Make the partition function work with numeric-prelude.

Add numeric-prelude to the deps.

Move the Vector and Matrix modules under Linear.

Bump to fixed-vector-0.2.*.

Replace (fromRational . toRational) with realToFrac everywhere.