4801e2d8bfAdd namespace boost::math::quadrature. Remove throw when condition number of summation exceeds precision of type; how to properly mollify the condition number is not clear and should be done consistently rather than ad-hoc.
Nick Thompson
2017-05-11 21:03:00 -0600
3bfeee838dDouble exponential quadrature.
Nick Thompson
2017-05-10 22:03:21 -0600
369ce4312bHalley/Luu version working for 50 decimal digits 'reference' test values.
pabristow
2017-05-10 18:21:35 +0100
135d05fe73legendre_stieltjes: add docs to math.qbk.
jzmaddock
2017-05-10 18:12:20 +0100
47a126c938Legendre-Stieltjes: Disable tests in C++03 mode which would fail.
jzmaddock
2017-05-10 18:11:37 +0100
5db9e8870bLegendre_stieltjes: hook up tests. Add concept checks. Fix some warnings, and reduce dependencies.
jzmaddock
2017-05-09 19:13:44 +0100
0e4c9ce418Merge branch 'legendre_prime' into develop
jzmaddock
2017-05-09 12:37:52 +0100
27ded3276dRemove C++11 isms from, legendre_prime header and tests - we don't want to break existing C++03 code using that header.
jzmaddock
2017-05-09 12:18:43 +0100
97a4184df9Hook up concept checks for legendre_p_prime. Fix errors from above, add domain-error handling to the function.
jzmaddock
2017-05-08 19:47:02 +0100
6afcd60c81Bring gcd/lcm tests into line with the inteter versions.
John Maddock
2017-05-04 11:22:18 +0100
4a569bf1c4Fix header order include so that ADL works
John Maddock
2017-05-04 11:17:11 +0100
7fab3cbb0echange namespace for gcd/lcm test
John Maddock
2017-05-04 11:11:12 +0100
bcff4eae78Remove unused variables, variables set but unused, set floating point number that would flush to zero to identically zero.
Nick Thompson
2017-05-03 15:30:24 -0600
cf35d745d1Remove unused typedef. Fix dangling reference to policy in airy. Pipe unused variables in compile test to /dev/null to suppress compiler warnings. Remove comment within comment.
Nick Thompson
2017-05-03 10:13:18 -0600
6407bedb75Bring polynomial gcd into line with new Boost.Integer code, and improve tests.
jzmaddock
2017-05-02 19:12:18 +0100
48a0cf714bPut license information in each file. Do not throw when condition number is approaching an undefined value; use mollification. Improve documentation.
Nick Thompson
2017-05-02 12:02:02 -0600
8b55166363Bring tests into line with new Boost.Integer gcd/lcm.
jzmaddock
2017-05-02 18:50:00 +0100
4632cb1169Remove typos in documentation. Remove hacks for float128 in unit tests.
Nick Thompson
2017-05-02 10:33:36 -0600
2bbadc69fbPlace copyright notices in file.
Nick Thompson
2017-05-01 18:36:26 -0600
4b4c4e9364Legendre-Stieltjes polynomials.
Nick Thompson
2017-04-30 15:58:10 -0600
69579ae210Update gcd/lcm docs to point to Boost.Integer.
jzmaddock
2017-04-24 19:42:22 +0100
ceae76a526Switch to using Boost.Integer gcd/lcm.
jzmaddock
2017-04-24 19:26:23 +0100
a918d0f331Return all zeros of the Legendre polynomial at once; there seems to be no benefit to having the user calculate them one at a time.
Nick Thompson
2017-04-22 15:47:22 -0500
1c10933943polynomial.hpp should not be including private Boost.Config header.
jzmaddock
2017-04-20 19:23:02 +0100
060f9e4da5Merge branch 'gcd_cpp14' into develop
jzmaddock
2017-04-20 19:18:41 +0100
fafda3176cFix gcd tests and code for clang on Linux
jzmaddock
2017-04-20 14:43:28 +0100
82e3949024GCD: Don't enable variadics unless the compiler supports them. Update docs. Escape from variadic GCD early when the result hits unity.
jzmaddock
2017-04-15 18:31:01 +0100
c1ab95b6ccAdd rational and multiprecision types to common_factor_test.cpp. Disable rational overloads if there is no numeric_limits support as boost::rational won't compile in that case anyway.
jzmaddock
2017-04-15 18:13:30 +0100
4110a69416Zeros of Legendre polynomials. This uses a root bracketing given by Szego with an asymptotic by Tricomi to get a domain and an initial guess for the root, then refines it via Newton's method.
Nick Thompson
2017-04-14 15:51:39 -0500
21bcf34a51Derivative of Legendre polynomials of the first kind. Used to generate the weights of Gaussian quadratures.
Nick Thompson
2017-04-12 14:16:26 -0500
e731cc8c3dUpdate performance test code to work with latest gcd.
jzmaddock
2017-04-10 12:28:57 +0100
89d28ff323Fix up documentation for adaptive trapezoidal quadrature.
Nick Thompson
2017-04-08 19:48:52 -0500
76a5423666Finally learned how to use quickbook. The documentation for the cubic b spline looks good now.
Nick Thompson
2017-04-08 18:47:42 -0500
70e82419d5Fix for expression template types.
jzmaddock
2017-04-08 20:12:23 +0100
f74da283aaAdd variadic gcd and lcm variations, update tests to match.
jzmaddock
2017-04-08 19:42:33 +0100
e71629f1c0Make gcd/lcm functions constexpr and noexcept where applicable.
jzmaddock
2017-04-08 19:13:08 +0100
8ee33649c7Remove dead code orphaned by previous commits.
jzmaddock
2017-04-07 09:55:31 +0100
5ca0e2ead8Tighten up polynomial-GCD algorithm selection. Disallow GCD's of polynomials of bounded integers as overflow is too likely to occur.
jzmaddock
2017-04-07 09:44:48 +0100
8e98f30cb2Remove SFINAE-logic and replace with simpler numeric_limits based algorithm selection.
jzmaddock
2017-04-07 09:43:40 +0100
39006ecacaReinstate gcd of rational numbers.
jzmaddock
2017-04-06 19:22:45 +0100
141fcf0b15Disable test case which invokes undefined behaviour
jzmaddock
2017-04-06 11:34:08 +0100
7795ad3fd1Merge branch 'gcd_polynomial' into gcd_cpp14
jzmaddock
2017-04-06 11:28:32 +0100
9d0dac1df9gcd: add the function objects back. Improve the gcd/lcm tests.
jzmaddock
2017-04-04 18:29:45 +0100
e4c0fdb0f2Make gcd/lcm independent of type_traits and mpl. Later we'll try and reconcile this version with that in Boost.Integer.
jzmaddock
2017-04-02 20:10:46 +0100
fff12aa16bImprove documentation for trapezoidal quadrature, explaining the role of the L1 norm calculation. Throw correct error type.
Nick Thompson
2017-03-28 13:56:18 -0500
9f58f0d337Add facility to check L1 norm and hence condition number of the quadrature sum. Throw when condition number exceeds the number of significant digits of the type.
Nick Thompson
2017-03-26 13:45:50 -0500
c5ee664a51Refactored with policies and passes tests and timing.
pabristow
2017-03-23 14:15:37 +0000
bdee37743dBug in values near to -exp(-1) corrected.
pabristow
2017-03-13 10:03:33 +0000
6734b363ddInclude spline documentation in math.qbk so it gets added to the build.
Nick Thompson
2017-03-11 15:56:34 -0600
0048d41dd6Added new constants expminusone and sixth
pabristow
2017-03-10 13:06:44 +0000
94d3cf4043refactored to use local test_value.hpp
pabristow
2017-03-06 18:10:52 +0000
36159a4426Discuss how the interpolant extrapolates.
Nick
2017-03-06 11:57:25 -0600
4efe710855Add tests detailing the use cases which make it difficult to use std::unique_ptr for pimpl idiom. Also, test that evaluating the spline outside the interpolation region succeeds as best it can.
Nick Thompson
2017-03-06 11:12:51 -0600
136e7411f5Adaptive Trapezoidal Quadrature
Nick Thompson
2017-03-05 19:05:41 -0600
446f8c415dCite Kress's book 'Numerical Analysis' in MLA format.
Nick Thompson
2017-03-05 10:49:02 -0600
7dbd074343Add population example.
Nick Thompson
2017-03-03 08:58:18 -0600
2506d7ed1fUse PIMPL so that the interpolant can be passed as a functor.
Nick Thompson
2017-03-02 17:16:21 -0600
47714b4e72Doc update of TODO list.
jzmaddock
2017-03-02 18:56:57 +0000
c42e8feebeRemove all traces of the linguistically incorrect word 'routine'.
Nick Thompson
2017-03-02 12:30:12 -0600
2b9d88d352Fixes suggested by code review of PR #57.
Nick Thompson
2017-03-02 11:54:55 -0600
b157403fd9Do not use Kahan summation to compute average; use update procedure that cannot overflow recommended by Knuth.
Nick Thompson
2017-02-27 20:58:42 -0600
f35a3551b0Merge pull request #56 from NAThompson/patch-2
jzmaddock
2017-02-27 18:12:02 +0000
b0b0810701Update history and regenerate docs for 1.64.
jzmaddock
2017-02-25 18:48:39 +0000
fee20ab932Given a function f, known at evenly spaced samples y_j = f(a + jh), this function constructs an interpolant using compactly supported cubic b splines. The advantage of using splines of compact support over traditional cubic splines is that compact support makes the splines well-conditioned.
Nick Thompson
2017-02-23 18:21:06 -0600
fa4dd6ff4cRemove finished TODO list entries
Nick
2017-02-22 14:08:22 -0600
4c19a1ec34Remove unneeded parameters in new beesel J/K code. Also fix some multiprecision-conceptual errors.
jzmaddock
2017-01-29 19:18:36 +0000
38eb9d3e38Update docs for new Bessel K0 and K1 approximation code.
jzmaddock
2017-01-24 18:35:54 +0000
fe66909b23Reorder some declarations to keep GCC happy. Remove some dead code.
jzmaddock
2017-01-23 19:34:34 +0000
Added missing header path, gmp and mpfr libraries
Michael Krasnyk
2016-12-29 23:37:40 +0100
ffb025ca2cFirst very rough prototype of Lambert W function, example of calculating diode current versus voltage, and some tests, including multiprecision and fixed_point types. Not yet using policies and trouble near the singularity at z=-exp(-1) and large z.
pabristow
2016-12-22 18:30:27 +0000
78320adb7areturn 0 to T(0) in float_distance changed because fixed-point requires explicit construction
pabristow
2016-12-05 16:41:18 +0000
1844134d07Changes to call correct error handlers in C99 annex F cases.
jzmaddock
2016-11-30 18:57:11 +0000