]> gitweb.michael.orlitzky.com - numerical-analysis.git/log
numerical-analysis.git
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'.

12 years agoAdd the ODE.IVP module with Euler's method.
Michael Orlitzky [Thu, 6 Dec 2012 02:14:15 +0000 (21:14 -0500)]
Add the ODE.IVP module with Euler's method.
Clean up .ghci a little bit and add the new module.

12 years agoBump test deps for ghc-7.6.
Michael Orlitzky [Thu, 6 Dec 2012 02:13:41 +0000 (21:13 -0500)]
Bump test deps for ghc-7.6.

12 years agoBump deps for ghc-7.6.
Michael Orlitzky [Wed, 5 Dec 2012 16:06:53 +0000 (11:06 -0500)]
Bump deps for ghc-7.6.

12 years agoAdd Integration.Simpson implementing Simpson's rule for numerical integration.
Michael Orlitzky [Fri, 23 Nov 2012 18:40:23 +0000 (13:40 -0500)]
Add Integration.Simpson implementing Simpson's rule for numerical integration.

12 years agoMove the partition function out of Integration.Trapezoid and into Misc.
Michael Orlitzky [Fri, 23 Nov 2012 18:39:42 +0000 (13:39 -0500)]
Move the partition function out of Integration.Trapezoid and into Misc.
Fix some trapezoid doctests.

12 years agoUse filemanip to automatically doctest all source files.
Michael Orlitzky [Fri, 23 Nov 2012 18:39:02 +0000 (13:39 -0500)]
Use filemanip to automatically doctest all source files.

12 years agoAdd the Integration.Trapezoid module; import it in the .ghci file.
Michael Orlitzky [Wed, 21 Nov 2012 03:02:27 +0000 (22:02 -0500)]
Add the Integration.Trapezoid module; import it in the .ghci file.

12 years agoUse the 2-norm instead of the infty-norm for the fixed point error ratios.
Michael Orlitzky [Thu, 18 Oct 2012 15:20:02 +0000 (11:20 -0400)]
Use the 2-norm instead of the infty-norm for the fixed point error ratios.

12 years agoFix a "norm" in Roots.Fast.
Michael Orlitzky [Wed, 17 Oct 2012 21:26:04 +0000 (17:26 -0400)]
Fix a "norm" in Roots.Fast.

12 years agoRename "norm" to "norm_2" in the Vector class.
Michael Orlitzky [Wed, 17 Oct 2012 21:25:31 +0000 (17:25 -0400)]
Rename "norm" to "norm_2" in the Vector class.
Add "norm_infty" to the Vector class.

12 years agoAdd an iteration count to the fixed_point function, rename it, and move it to the...
Michael Orlitzky [Wed, 17 Oct 2012 16:41:12 +0000 (12:41 -0400)]
Add an iteration count to the fixed_point function, rename it, and move it to the Roots.Fast module.
Implement fixed_point and fixed_point_iteration_count in Roots.Simple in terms of the general function.
Add a fixed_point_error_ratios function.

12 years agoRequire Vector instances to be instances of Num as well.
Michael Orlitzky [Wed, 17 Oct 2012 16:40:38 +0000 (12:40 -0400)]
Require Vector instances to be instances of Num as well.

12 years agoImport Aliases and TwoTuple in the .ghci file.
Michael Orlitzky [Mon, 15 Oct 2012 05:15:19 +0000 (01:15 -0400)]
Import Aliases and TwoTuple in the .ghci file.

12 years agoAdd the Aliases module.
Michael Orlitzky [Mon, 15 Oct 2012 05:15:03 +0000 (01:15 -0400)]
Add the Aliases module.

12 years agoAdd a custom Show instance for TwoTuple.
Michael Orlitzky [Mon, 15 Oct 2012 05:14:37 +0000 (01:14 -0400)]
Add a custom Show instance for TwoTuple.

12 years agoChange the type signature of fixed_point to work on Vectors.
Michael Orlitzky [Sun, 14 Oct 2012 18:05:59 +0000 (14:05 -0400)]
Change the type signature of fixed_point to work on Vectors.

12 years agoBump the doctest dependency.
Michael Orlitzky [Sun, 14 Oct 2012 18:05:06 +0000 (14:05 -0400)]
Bump the doctest dependency.

12 years agoAdd a Vector class and make TwoTuple (and several base types) instances of it.
Michael Orlitzky [Sun, 14 Oct 2012 18:04:31 +0000 (14:04 -0400)]
Add a Vector class and make TwoTuple (and several base types) instances of it.

12 years agoBegin work on a simple 2-tuple type.
Michael Orlitzky [Sun, 14 Oct 2012 02:36:49 +0000 (22:36 -0400)]
Begin work on a simple 2-tuple type.

12 years agoMake numbers == 3000.1.* required.
Michael Orlitzky [Fri, 12 Oct 2012 17:16:36 +0000 (13:16 -0400)]
Make numbers == 3000.1.* required.
Import Data.Number.BigFloat in the .ghci file.
Add a doctest for newtons_method using BigFloat.