]> gitweb.michael.orlitzky.com - numerical-analysis.git/log
numerical-analysis.git
10 years agoFix the incorrect definition of "minor".
Michael Orlitzky [Wed, 12 Feb 2014 06:52:35 +0000 (01:52 -0500)]
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.

10 years agoAdd big_K code to FEM.R1.
Michael Orlitzky [Wed, 12 Feb 2014 05:00:27 +0000 (00:00 -0500)]
Add big_K code to FEM.R1.

10 years agoUse a better implementation of backwards_substitute.
Michael Orlitzky [Wed, 12 Feb 2014 01:14:47 +0000 (20:14 -0500)]
Use a better implementation of backwards_substitute.

10 years agoAdd Linear.Matrix.ifoldr2.
Michael Orlitzky [Wed, 12 Feb 2014 00:21:53 +0000 (19:21 -0500)]
Add Linear.Matrix.ifoldr2.

10 years agoImplement forward substitute in terms of a fold.
Michael Orlitzky [Tue, 11 Feb 2014 21:10:53 +0000 (16:10 -0500)]
Implement forward substitute in terms of a fold.

10 years agoAdd Linear.Matrix.set_idx.
Michael Orlitzky [Tue, 11 Feb 2014 21:10:18 +0000 (16:10 -0500)]
Add Linear.Matrix.set_idx.
Loosen constraints on the dot product function.

10 years agoUpdate Integration.Gaussian to use zip2.
Michael Orlitzky [Tue, 11 Feb 2014 20:25:49 +0000 (15:25 -0500)]
Update Integration.Gaussian to use zip2.

10 years agoGeneralize colzip and colzipwith to any matrices.
Michael Orlitzky [Tue, 11 Feb 2014 20:22:35 +0000 (15:22 -0500)]
Generalize colzip and colzipwith to any matrices.
Add zip2three to zip three matrices together.

10 years agoFix backward_substitute.
Michael Orlitzky [Tue, 11 Feb 2014 20:07:53 +0000 (15:07 -0500)]
Fix backward_substitute.
Add 'solve_positive_definite' to Linear.System.
Add tests for everything in Linear.System.

10 years agoAdd reverse2 to Linear.Matrix.
Michael Orlitzky [Tue, 11 Feb 2014 20:07:35 +0000 (15:07 -0500)]
Add reverse2 to Linear.Matrix.

10 years agoAdd Eq constraints throughout Linear.Iteration.
Michael Orlitzky [Tue, 11 Feb 2014 18:16:26 +0000 (13:16 -0500)]
Add Eq constraints throughout Linear.Iteration.

10 years agoAdd another test for 'cholesky'.
Michael Orlitzky [Tue, 11 Feb 2014 18:15:15 +0000 (13:15 -0500)]
Add another test for 'cholesky'.
Improve implementation of upper/lower triangular tests.
Drop some typeclass constraints on upper/lower triangular tests.

10 years agoAdd the imap2 function to Linear.Matrix.
Michael Orlitzky [Sat, 8 Feb 2014 00:33:45 +0000 (19:33 -0500)]
Add the imap2 function to Linear.Matrix.

10 years agoMake safe indexing slower but not dependent on the number of rows/columns.
Michael Orlitzky [Sat, 8 Feb 2014 00:24:14 +0000 (19:24 -0500)]
Make safe indexing slower but not dependent on the number of rows/columns.

10 years agoRemove the "row" function which returned a vector instead of a matrix.
Michael Orlitzky [Sat, 8 Feb 2014 00:10:02 +0000 (19:10 -0500)]
Remove the "row" function which returned a vector instead of a matrix.

10 years agoDrop the 'column' function that returned a vector instead of a matrix.
Michael Orlitzky [Fri, 7 Feb 2014 23:56:21 +0000 (18:56 -0500)]
Drop the 'column' function that returned a vector instead of a matrix.

10 years agoRename matmap to map2.
Michael Orlitzky [Fri, 7 Feb 2014 23:46:51 +0000 (18:46 -0500)]
Rename matmap to map2.

10 years agoMore progress on the 1d FEM, the "big F" vector is now computed correctly.
Michael Orlitzky [Fri, 7 Feb 2014 15:20:26 +0000 (10:20 -0500)]
More progress on the 1d FEM, the "big F" vector is now computed correctly.

10 years agoAdd the 'ifoldl2' function to Linear.Matrix.
Michael Orlitzky [Fri, 7 Feb 2014 15:19:53 +0000 (10:19 -0500)]
Add the 'ifoldl2' function to Linear.Matrix.

10 years agoFix a typo.
Michael Orlitzky [Fri, 7 Feb 2014 15:15:39 +0000 (10:15 -0500)]
Fix a typo.

10 years agoAdd a bunch more natural numbers in the Naturals module.
Michael Orlitzky [Fri, 7 Feb 2014 15:15:11 +0000 (10:15 -0500)]
Add a bunch more natural numbers in the Naturals module.

10 years agoAdd the (in progress) FEM.R1 module.
Michael Orlitzky [Thu, 6 Feb 2014 06:36:45 +0000 (01:36 -0500)]
Add the (in progress) FEM.R1 module.

10 years agoFix some docs in Roots.Simple.
Michael Orlitzky [Thu, 6 Feb 2014 06:36:28 +0000 (01:36 -0500)]
Fix some docs in Roots.Simple.

10 years agoImplement Gaussian quadrature and test it.
Michael Orlitzky [Wed, 5 Feb 2014 20:30:47 +0000 (15:30 -0500)]
Implement Gaussian quadrature and test it.

10 years agoRename zipcol to colzip.
Michael Orlitzky [Wed, 5 Feb 2014 20:30:26 +0000 (15:30 -0500)]
Rename zipcol to colzip.
Add a new function, colzipwith.

10 years agoAdd the Integration.Gaussian.nodes_and_weights function.
Michael Orlitzky [Wed, 5 Feb 2014 01:33:31 +0000 (20:33 -0500)]
Add the Integration.Gaussian.nodes_and_weights function.

10 years agoAdd zipcol and matmap functions in Linear.Matrix.
Michael Orlitzky [Wed, 5 Feb 2014 01:33:09 +0000 (20:33 -0500)]
Add zipcol and matmap functions in Linear.Matrix.

10 years agoAdd the Integration.Gaussian module where we begin to implement Gaussian quadrature.
Michael Orlitzky [Tue, 4 Feb 2014 22:41:19 +0000 (17:41 -0500)]
Add the Integration.Gaussian module where we begin to implement Gaussian quadrature.

10 years agoAdd comments to Linear.QR about non-convergence.
Michael Orlitzky [Tue, 4 Feb 2014 22:40:59 +0000 (17:40 -0500)]
Add comments to Linear.QR about non-convergence.

10 years agoFinish the eigenvectors_symmetric tests.
Michael Orlitzky [Tue, 4 Feb 2014 20:03:56 +0000 (15:03 -0500)]
Finish the eigenvectors_symmetric tests.

10 years agoAdd the eigenvectors_symmetric function and begin writing tests for it.
Michael Orlitzky [Tue, 4 Feb 2014 17:16:52 +0000 (12:16 -0500)]
Add the eigenvectors_symmetric function and begin writing tests for it.

10 years agoImplement the QR algorithm for computing eigenvalues.
Michael Orlitzky [Tue, 4 Feb 2014 07:19:04 +0000 (02:19 -0500)]
Implement the QR algorithm for computing eigenvalues.

10 years agoAdd type synonyms for column/row matrices.
Michael Orlitzky [Tue, 4 Feb 2014 07:18:09 +0000 (02:18 -0500)]
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.

10 years agoSmall cleanups in Linear.Vector.
Michael Orlitzky [Mon, 3 Feb 2014 06:03:14 +0000 (01:03 -0500)]
Small cleanups in Linear.Vector.
Add 'diagonal' and 'trace' functions to Linear.Matrix.

10 years agoClean up imports everywhere.
Michael Orlitzky [Mon, 3 Feb 2014 02:10:01 +0000 (21:10 -0500)]
Clean up imports everywhere.

10 years agoAdd tolerant versions of is_{upper,lower}_triangular.
Michael Orlitzky [Mon, 3 Feb 2014 01:21:33 +0000 (20:21 -0500)]
Add tolerant versions of is_{upper,lower}_triangular.
Fix the QR factorization code.
Add tests for the QR behavior.

10 years agoAdd a Normed instance for (Vec n a).
Michael Orlitzky [Sun, 2 Feb 2014 23:46:48 +0000 (18:46 -0500)]
Add a Normed instance for (Vec n a).

10 years agoSwitch Linear.Vector to the numeric prelude and add the element_sum function to it.
Michael Orlitzky [Sun, 2 Feb 2014 23:46:05 +0000 (18:46 -0500)]
Switch Linear.Vector to the numeric prelude and add the element_sum function to it.
Use Linear.Vector.element_sum in Linear.Matrix.

10 years agoRemove old comment from Normed.
Michael Orlitzky [Sun, 2 Feb 2014 20:41:58 +0000 (15:41 -0500)]
Remove old comment from Normed.

10 years agoAdd the Frobenius norm to Linear.Matrix.
Michael Orlitzky [Sun, 2 Feb 2014 17:34:42 +0000 (12:34 -0500)]
Add the Frobenius norm to Linear.Matrix.

10 years agoAdd the Linear.QR module and update the ghci/cabal files.
Michael Orlitzky [Sun, 2 Feb 2014 17:12:09 +0000 (12:12 -0500)]
Add the Linear.QR module and update the ghci/cabal files.

10 years agoNew function: Linear.Matrix.identity_matrix.
Michael Orlitzky [Sun, 2 Feb 2014 16:33:03 +0000 (11:33 -0500)]
New function: Linear.Matrix.identity_matrix.

10 years agoFix implementation of Linear.Matrix.construct.
Michael Orlitzky [Sun, 2 Feb 2014 04:16:01 +0000 (23:16 -0500)]
Fix implementation of Linear.Matrix.construct.

10 years agoAdd the Polynomials.Orthogonal module.
Michael Orlitzky [Thu, 30 Jan 2014 21:02:38 +0000 (16:02 -0500)]
Add the Polynomials.Orthogonal module.

10 years agoBump dependencies.
Michael Orlitzky [Thu, 30 Jan 2014 14:19:17 +0000 (09:19 -0500)]
Bump dependencies.

11 years agoAdd the rayleigh_quotient function.
Michael Orlitzky [Thu, 15 Aug 2013 17:11:43 +0000 (13:11 -0400)]
Add the rayleigh_quotient function.

11 years agoAdd Gauss-Seidel and SOR iterative methods.
Michael Orlitzky [Sun, 21 Jul 2013 21:21:12 +0000 (17:21 -0400)]
Add Gauss-Seidel and SOR iterative methods.

11 years agoRename 'diagonal' to 'diagonal_part'.
Michael Orlitzky [Sun, 21 Jul 2013 21:19:35 +0000 (17:19 -0400)]
Rename 'diagonal' to 'diagonal_part'.
Add functions to pull out the (strict) lower/upper triangular parts of a matrix.

11 years agoFix a bunch of hlint warnings.
Michael Orlitzky [Sun, 21 Jul 2013 18:25:30 +0000 (14:25 -0400)]
Fix a bunch of hlint warnings.

11 years agoUpdate ghci/cabal files to load the new module.
Michael Orlitzky [Sun, 21 Jul 2013 14:48:47 +0000 (10:48 -0400)]
Update ghci/cabal files to load the new module.

11 years agoAdd the Linear.Iteration module.
Michael Orlitzky [Sun, 21 Jul 2013 14:48:23 +0000 (10:48 -0400)]
Add the Linear.Iteration module.

11 years agoAdd another example for forward_substitute.
Michael Orlitzky [Sun, 21 Jul 2013 14:14:19 +0000 (10:14 -0400)]
Add another example for forward_substitute.

11 years agoAdd 'diagonal' function.
Michael Orlitzky [Sun, 21 Jul 2013 14:13:50 +0000 (10:13 -0400)]
Add 'diagonal' function.

11 years agoConstrain Normed instance to column vectors (m-by-1 matrices).
Michael Orlitzky [Sun, 21 Jul 2013 14:13:24 +0000 (10:13 -0400)]
Constrain Normed instance to column vectors (m-by-1 matrices).

11 years agoAdd the unit_roundoff function.
Michael Orlitzky [Mon, 1 Jul 2013 20:00:59 +0000 (16:00 -0400)]
Add the unit_roundoff function.

11 years agoUpdate numeric-prelude and fixed-vector.
Michael Orlitzky [Sat, 8 Jun 2013 21:19:36 +0000 (17:19 -0400)]
Update numeric-prelude and fixed-vector.
Use the (!) function from the new fixed-vector.
Add the Fast/Simple trisection method and a test.

11 years agoGet determinants working.
Michael Orlitzky [Sun, 24 Feb 2013 14:55:10 +0000 (09:55 -0500)]
Get determinants working.

11 years agoReplace the whole Matrix implementation with something a little better.
Michael Orlitzky [Sun, 24 Feb 2013 07:41:03 +0000 (02:41 -0500)]
Replace the whole Matrix implementation with something a little better.
Add the Linear.System module.
Get all doctests working again.
Update .ghci/cabal files.

11 years agoAdd Vec1, Vec4, Vec5 aliases.
Michael Orlitzky [Sun, 24 Feb 2013 07:01:24 +0000 (02:01 -0500)]
Add Vec1, Vec4, Vec5 aliases.

11 years agoAdd triangular functions, determinant, minor (all half-baked) to Linear.Matrix.
Michael Orlitzky [Sun, 24 Feb 2013 00:05:26 +0000 (19:05 -0500)]
Add triangular functions, determinant, minor (all half-baked) to Linear.Matrix.

11 years agoAdd the Linear.System module.
Michael Orlitzky [Sun, 24 Feb 2013 00:05:10 +0000 (19:05 -0500)]
Add the Linear.System module.

11 years agoAdd the delete function to Vector.hs.
Michael Orlitzky [Sun, 24 Feb 2013 00:04:51 +0000 (19:04 -0500)]
Add the delete function to Vector.hs.

11 years agoFix compiler warnings and doctests.
Michael Orlitzky [Thu, 21 Feb 2013 04:08:56 +0000 (23:08 -0500)]
Fix compiler warnings and doctests.

11 years agoConvert scalars to 1x1 vectors.
Michael Orlitzky [Thu, 21 Feb 2013 03:34:29 +0000 (22:34 -0500)]
Convert scalars to 1x1 vectors.

11 years agoA huge pile of crap upon Matrix/Vector.
Michael Orlitzky [Thu, 21 Feb 2013 02:47:50 +0000 (21:47 -0500)]
A huge pile of crap upon Matrix/Vector.

11 years agoUpdate .ghci.
Michael Orlitzky [Thu, 21 Feb 2013 02:47:15 +0000 (21:47 -0500)]
Update .ghci.

11 years agoRemove assumptions on the Normed class.
Michael Orlitzky [Thu, 21 Feb 2013 02:46:56 +0000 (21:46 -0500)]
Remove assumptions on the Normed class.

11 years agoConvert Integration/Simpson.hs and Integration/Trapezoid.hs to numeric-prelude.
Michael Orlitzky [Thu, 21 Feb 2013 00:06:10 +0000 (19:06 -0500)]
Convert Integration/Simpson.hs and Integration/Trapezoid.hs to numeric-prelude.

11 years agoUse RealField/RealRing where possible instead of their constituents.
Michael Orlitzky [Wed, 20 Feb 2013 23:56:13 +0000 (18:56 -0500)]
Use RealField/RealRing where possible instead of their constituents.
Convert src/ODE/IVP.hs to numeric-prelude.
Add RealRing/RealField instanced to BigFloat.

11 years agoRename Aliases.hs to BigFloat.hs, now containing numeric-prelude instances for BigFloats.
Michael Orlitzky [Wed, 20 Feb 2013 23:35:05 +0000 (18:35 -0500)]
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.

11 years agoMake the partition function work with numeric-prelude.
Michael Orlitzky [Wed, 20 Feb 2013 19:31:41 +0000 (14:31 -0500)]
Make the partition function work with numeric-prelude.

11 years agoAdd numeric-prelude to the deps.
Michael Orlitzky [Wed, 20 Feb 2013 19:31:33 +0000 (14:31 -0500)]
Add numeric-prelude to the deps.

11 years agoMove the Vector and Matrix modules under Linear.
Michael Orlitzky [Tue, 12 Feb 2013 17:16:39 +0000 (12:16 -0500)]
Move the Vector and Matrix modules under Linear.

11 years agoBump to fixed-vector-0.2.*.
Michael Orlitzky [Mon, 11 Feb 2013 14:34:30 +0000 (09:34 -0500)]
Bump to fixed-vector-0.2.*.

11 years agoReplace (fromRational . toRational) with realToFrac everywhere.
Michael Orlitzky [Tue, 5 Feb 2013 03:17:11 +0000 (22:17 -0500)]
Replace (fromRational . toRational) with realToFrac everywhere.

11 years agoUpdate cabal file for the separate project/repo.
Michael Orlitzky [Tue, 5 Feb 2013 02:29:36 +0000 (21:29 -0500)]
Update cabal file for the separate project/repo.

11 years agoFix import and module list names.
Michael Orlitzky [Tue, 5 Feb 2013 01:42:10 +0000 (20:42 -0500)]
Fix import and module list names.

11 years agoAdd a 2D fixed-point doctest (example) to Vector.hs.
Michael Orlitzky [Tue, 5 Feb 2013 01:41:36 +0000 (20:41 -0500)]
Add a 2D fixed-point doctest (example) to Vector.hs.
Remove unused imports.

11 years agoAdd a Normed instance for Float.
Michael Orlitzky [Tue, 5 Feb 2013 01:41:09 +0000 (20:41 -0500)]
Add a Normed instance for Float.

11 years agoRemove non-fixed Matrix module.
Michael Orlitzky [Tue, 5 Feb 2013 01:25:14 +0000 (20:25 -0500)]
Remove non-fixed Matrix module.
Drop the "Fixed" prefix for Vector/Matrix modules.
Update ghci includes.

11 years agoClean up imports in FixedMatrix.hs.
Michael Orlitzky [Tue, 5 Feb 2013 01:23:33 +0000 (20:23 -0500)]
Clean up imports in FixedMatrix.hs.
Add matrix multiplication and Cholesky factorization to FixedMatrix.hs.
Add some doctests to FixedMatrix.hs.

11 years agoClean up imports in FixedVector.hs.
Michael Orlitzky [Tue, 5 Feb 2013 01:22:47 +0000 (20:22 -0500)]
Clean up imports in FixedVector.hs.

11 years agoRemove reimplemented vector functions from FixedVector.hs.
Michael Orlitzky [Tue, 5 Feb 2013 00:01:00 +0000 (19:01 -0500)]
Remove reimplemented vector functions from FixedVector.hs.
Replace custom Vn functions with vector equivalents.

11 years agoAdd a Vector instance for Vn.
Michael Orlitzky [Mon, 4 Feb 2013 23:49:47 +0000 (18:49 -0500)]
Add a Vector instance for Vn.

11 years agoBegin the fixed-size matrix implementation.
Michael Orlitzky [Mon, 4 Feb 2013 20:16:12 +0000 (15:16 -0500)]
Begin the fixed-size matrix implementation.

11 years agoEncode the vector type in Vn.
Michael Orlitzky [Mon, 4 Feb 2013 20:14:38 +0000 (15:14 -0500)]
Encode the vector type in Vn.
Add FixedVector.length.
Add FixedVector.map.
Include a special 4D type.

11 years agoAdd toList and fromList functions in FixedVector.hs.
Michael Orlitzky [Mon, 4 Feb 2013 17:31:48 +0000 (12:31 -0500)]
Add toList and fromList functions in FixedVector.hs.

11 years agoRemove TwoTuple.hs.
Michael Orlitzky [Mon, 4 Feb 2013 17:31:32 +0000 (12:31 -0500)]
Remove TwoTuple.hs.

11 years agoAdd 2d/3d wrapper types for easy pattern matching.
Michael Orlitzky [Sun, 3 Feb 2013 22:20:38 +0000 (17:20 -0500)]
Add 2d/3d wrapper types for easy pattern matching.
Implement both safe and unsafe indexing.

11 years agoReplace the custom 'Vector' typeclass with 'Normed' everywhere.
Michael Orlitzky [Sun, 3 Feb 2013 20:25:53 +0000 (15:25 -0500)]
Replace the custom 'Vector' typeclass with 'Normed' everywhere.

11 years agoAdd cleaner Show instance for Vn.
Michael Orlitzky [Sun, 3 Feb 2013 20:09:30 +0000 (15:09 -0500)]
Add cleaner Show instance for Vn.
Fix compiler warnings in FixedVector.hs.
Add more doctests to FixedVector.hs.

11 years agoAdd FixedVector and Matrix to the compiled modules.
Michael Orlitzky [Sun, 3 Feb 2013 20:08:38 +0000 (15:08 -0500)]
Add FixedVector and Matrix to the compiled modules.
Update deps for the new modules.

11 years agoFix compiler warnings in Matrix.hs.
Michael Orlitzky [Sun, 3 Feb 2013 20:01:19 +0000 (15:01 -0500)]
Fix compiler warnings in Matrix.hs.

11 years agoBegin implementation of normed, fixed-length vectors.
Michael Orlitzky [Sun, 3 Feb 2013 17:16:43 +0000 (12:16 -0500)]
Begin implementation of normed, fixed-length vectors.

11 years agoFix two compiler warnings.
Michael Orlitzky [Sun, 3 Feb 2013 02:43:16 +0000 (21:43 -0500)]
Fix two compiler warnings.

11 years agoRemove unused dependencies.
Michael Orlitzky [Sun, 3 Feb 2013 02:34:00 +0000 (21:34 -0500)]
Remove unused dependencies.
Set looser base version bounds.
Expose (and thus compile) more modules.

11 years agoAdd a simple unsafe matrix class to test the Cholesky algorithm.
Michael Orlitzky [Sun, 3 Feb 2013 01:38:13 +0000 (20:38 -0500)]
Add a simple unsafe matrix class to test the Cholesky algorithm.

11 years agoAdd an alternate book, and a gitignore for the numerical-analysis/dist directory.
Michael Orlitzky [Sun, 23 Dec 2012 05:51:23 +0000 (00:51 -0500)]
Add an alternate book, and a gitignore for the numerical-analysis/dist directory.

11 years agoFix an off-by-one error in Euler's method.
Michael Orlitzky [Sun, 9 Dec 2012 09:22:37 +0000 (04:22 -0500)]
Fix an off-by-one error in Euler's method.

11 years agoAdd an implementation of Euler's method that takes the step size 'h' as an argument...
Michael Orlitzky [Sun, 9 Dec 2012 08:22:25 +0000 (03:22 -0500)]
Add an implementation of Euler's method that takes the step size 'h' as an argument instead of the number of subintervals 'n'.