townforge/math
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity
5,028 Commits 1 Branch 0 Tags 151 MiB
cc
Commit Graph

2240 Commits

Author SHA1 Message Date
Nick Thompson
ec3412ac28 Better comments, #include guards. 2018-01-12 13:14:34 -06:00
Nick Thompson
358d9d3496 In float precision, make sure that eps < x[j] < 1 - eps, or we hit singularities on the boundary. 2018-01-11 21:57:15 -06:00
Nick Thompson
a13c7e9ab0 Double infinite limits. 2018-01-11 12:08:53 -06:00
Nick Thompson
d5eb6e3995 Support infinite lower and upper bounds. 2018-01-11 11:39:21 -06:00
Nick Thompson
08550d4064 3 template arguments -> 1 template argument. 2018-01-10 14:08:23 -06:00
Nick Thompson
d0e52cf8c9 Policies for monte-carlo 2017-12-30 10:53:53 -07:00
Nick Thompson
0f007c23c6 Reuse temporaries to increase performance. 2017-12-23 17:06:20 -07:00
Nick Thompson
1a29ce99c6 Demonstrate affine invariance of Catmull-Rom splines. 2017-12-22 21:52:52 -07:00
Nick Thompson
4cce07e15e Implement tangent vector computation. Fix index lookup. Close all Catmull-Rom curves and document. 2017-12-22 16:31:36 -07:00
Nick Thompson
cad34ff756 First pass at a Catmull-Rom curve interpolator. 2017-12-21 16:12:24 -07:00
jzmaddock
0dd6df4221 Math.Hypergeometric_1F1: add error handling for large recurrences. 2017-12-02 19:16:58 +00:00
jzmaddock
15a1156b4c Math.Hypergeometric_1F1: Update tests for small random values, fix resulting errors. 2017-12-02 19:10:50 +00:00
Nick Thompson
0da69d2626 Add isfinite to std_real_concept to fix build. 2017-11-29 16:17:08 -07:00
Nick Thompson
d53284370a Use ADL to green some more builds. Add requirements to build on Chebyshev transform. 2017-11-29 15:09:20 -07:00
Nick Thompson
401e353385 Hopefully get a few more builds green; eliminate compiler warnings. 2017-11-28 08:41:23 -06:00
jzmaddock
ef10ba37d2 Merge branch 'lambertw' of https://github.com/boostorg/math into lambertw 2017-11-27 19:08:42 +00:00
jzmaddock
3cebb614a5 Math.LambertW: fix Halley termination condition. 2017-11-27 19:06:53 +00:00
pabristow
66067c3b6b Added include to test_value to ensure that when used standalone with GCC that supports suffix Q, multiprecision float128 is included. 2017-11-24 16:53:18 +00:00
Nick Thompson
57ba888a50 Unit test variance. 2017-11-23 14:18:05 -07:00
Nick Thompson
773076475c Fix exception safety in threads. 2017-11-23 10:53:35 -07:00
Nick Thompson
67578ce12b Naive Monte-Carlo integration. Still needs to handle exceptions in threads. 2017-11-21 18:49:14 -07:00
pabristow
9b9414c256 numerous minor fixes for JM 1st comments on docs 19Nov17 2017-11-20 17:03:29 +00:00
jzmaddock
be326922dc Math: add start of 1F1 implementation. 2017-11-19 11:58:26 +00:00
pabristow
6aa9f286ff Tests OK, including multiprecision. 2017-11-16 17:58:59 +00:00
jzmaddock
587f2d8918 Math.Hypergeometric: Add 2F0 function and tests. 
Also fix a couple of warnings in test_data.hpp, and improve the continued fraction form in 0F1.
 2017-11-15 18:27:35 +00:00
jzmaddock
1783c3a74c Math: Fix inspection report errors. 2017-11-09 19:39:32 +00:00
pabristow
89ae728836 added graphs and updated docs to use them. 2017-11-06 17:18:33 +00:00
Nick Thompson
713f17f76a Fix typo in error message. 2017-11-04 13:50:12 -06:00
pabristow
3c3c91709f w-1 branch improvements and tests OK 2017-11-02 18:47:30 +00:00
jzmaddock
47fa45bee4 Import SOC headers, Test and fix up 1F0 and 0F1. 2017-10-29 19:33:39 +00:00
jzmaddock
81c8a80dd5 Update big_constant.hpp to work standalone. 2017-10-29 19:32:13 +00:00
pabristow
04121d2872 Lambert w for huge z added, and test and more docs 2017-10-28 19:13:58 +01:00
pabristow
f51d987acd added much on W-1 branch (handling tiny z), and more docs. 2017-10-27 18:18:06 +01:00
pabristow
d01d0c4eb8 Working version with C array of precomputed lookup tables, but still inline inside 2017-10-13 09:43:39 +01:00
pabristow
dd06d2beae 100 decimal digit precision values and tests OK. 2017-10-10 15:13:49 +01:00
jzmaddock
97f61f0020 Math.Quaternion: Add workarounds for older compilers and fix up constexpr support. 2017-10-06 19:22:45 +01:00
jzmaddock
cfdd05f158 Math.Quaternion: Add move semantics. 2017-10-05 18:18:21 +01:00
jzmaddock
94e92dc5f4 Math.Quaternion: add swap method. 2017-10-05 18:12:24 +01:00
jzmaddock
68b82322ab Math.Quaternion: remove dead code from header. 2017-10-05 13:06:20 +01:00
jzmaddock
601ee6e00f Math.Quaternion: Rewrite test program, mostly remove valarray usage from header. 2017-10-05 12:50:24 +01:00
pabristow
39846818e1 Expanded docs using snippets. 2017-10-03 15:11:25 +01:00
John Maddock
f7acd35e51 Math.Chebyshev: tidy up and simplify fftw dispatch code with a traits class to abstract the differences. 2017-09-25 19:37:36 +01:00
John Maddock
2eabd87de1 Math.Chebyshev: Fix up chebeshev_transform for more types. 2017-09-21 19:13:36 +01:00
John Maddock
aaa6db8880 Math.Chebyshev: Fix ambiguous calls to pow with gcc-7 2017-09-21 17:56:51 +01:00
jzmaddock
9ee7a93678 Math.Chebyshev: Modify chebyshev_transform to dispatch to correct fftw3 functions. 
Begin to hook up fftw3 tests.
 2017-09-21 13:40:23 +01:00
jzmaddock
a1c18789a0 Math.chebyshev: Hook up conceptual tests and fix resulting errors. 2017-09-21 11:52:59 +01:00
jzmaddock
c7f243239e Math.chebyshev: Backport to C++03. 2017-09-20 19:07:17 +01:00
jzmaddock
82392e3e0c Merge branch 'develop' into chebyshev 2017-09-20 17:59:38 +01:00
Nick Thompson
dd9f2b7033 Since I cannot find a symmetry to reuse function evaluations after Chebyshev grid refinement, just give up for now and simplify the constructor accordingly. 2017-09-19 00:48:20 -06:00
Nick Thompson
956864cbd5 Recurse when grid isn't sufficiently refined. Optimization potential exists by reusing function evaluations. 2017-09-17 21:16:48 -06:00
Nick Thompson
0d85d4d11e Performance improvements. 2017-09-17 19:30:15 -06:00
Nick Thompson
2cb935ec17 Use std::enable_if to obviate the need for C++17 if constexpr. 2017-09-17 16:23:50 -06:00
Nick Thompson
0acb51b78d Linear time and numerically stable evaluation of derivative of Chebyshev series by a modified Clenshaw recurrence. 2017-09-17 14:02:43 -06:00
Nick Thompson
cb470280ba Chebyshev transform 2017-09-13 17:51:50 -06:00
jzmaddock
5755e1d37e Merge pull request #87 from boostorg/gauss 
Merge Gauss Kronrod code into develop.
 2017-09-13 16:57:06 +01:00
jzmaddock
25cc8fdad7 Add Gauss and Gauss-Kronrod docs. 2017-09-09 19:18:08 +01:00
jzmaddock
4c0bec596c Remove dead files. 2017-09-06 19:56:15 +01:00
jzmaddock
d9893f5e01 Quadrature: update conceptual tests for Gauss (Kronrod) quadrature. 2017-09-06 19:53:10 +01:00
jzmaddock
e09ad4c3aa Quadrature: finish off Gauss and Gauss-Kronrod integration code and tests. 2017-09-06 18:03:15 +01:00
pabristow
15568b8d62 Work on precision demo before holiday. OK on VS14.1 and GCC 7.1.0 but docs need much more work. 2017-09-02 12:12:10 +01:00
jzmaddock
7d2002db80 Quadrature: add gauss and gauss-kronrod quadrature. 2017-08-31 19:42:26 +01:00
pabristow
c1b8ffec7c Needed to merge? 2017-08-31 12:50:11 +01:00
pabristow
48435ed996 Commit before merge develop to get new config file 2017-08-30 17:05:34 +01:00
jzmaddock
1527667e20 Next.hpp: Change struct to class to avoid msvc compiler warnings. 2017-08-29 20:00:32 +01:00
pabristow
c408d7ec1c Restored test_out_of_range to original location as is used by other tests 2017-08-29 17:03:06 +01:00
pabristow
34fd59c372 students t allowing df and x infinite 2017-08-29 15:24:45 +01:00
pabristow
0f070dd259 Commit work-in-progress before fixing students t on develop 2017-08-29 08:21:59 +01:00
jzmaddock
e7e915816f Merge branch 'gauss_konrod' of https://github.com/NAThompson/math into gauss 
# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.
 2017-08-24 11:47:26 +01:00
jzmaddock
e100973d53 Merge pull request #81 from lakshayg/noexcept 
make move constructor and operator= noexcept
 2017-08-23 19:02:49 +01:00
jzmaddock
f301e2a197 Merge pull request #72 from boostorg/pr66 
PR66: Moved to here for integration testing.
 2017-08-23 18:47:46 +01:00
Lakshay Garg
4fc8d43b4b make move constructor and operator= noexcept 
This change helps the polynomial class play well with STL
containers.

> If the move constructor for an element type in a container
> is not noexcept then the container will use the copy constructor rather
> than the move constructor -- HIC++ Version 4.0

Benchmarking shows that the number of calls to copy constructor
are reduced.
 2017-08-22 12:41:32 +05:30
jzmaddock
f6495af64c Polynomial: more simplifications and modernizations. 
Also update msvc performance results.
 2017-08-21 13:15:39 +01:00
jzmaddock
d9858963d0 Merge branch 'move' of https://github.com/lakshayg/math into pr79 2017-08-20 09:16:20 +01:00
jzmaddock
b70d7b1359 Add polynomial performance test. 2017-08-20 08:28:41 +01:00
Lakshay Garg
36b272c42b enable move semantics for polynomial 2017-08-17 09:11:40 +05:30
jzmaddock
bc03a4ec79 Remove forward references to gcd/lcm types: 
These are no longer part of Boost.Math, but moved to Integer.
 2017-08-14 19:56:14 +01:00
pabristow
26aea4e7df Big refactor JM small_z and tag_type select code 2017-08-10 17:49:17 +01:00
jzmaddock
68c80b5fc8 Merge branch 'pr66' of https://github.com/NAThompson/math into pr66_78 
# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.
 2017-08-09 18:02:36 +01:00
jzmaddock
3770719aff Quadrature: Remove the minimum levels requirement on tanh_sinh. 2017-08-07 19:04:42 +01:00
jzmaddock
a928ff3df3 Quadrature: Correct logic in finding the first abscissa values which are stored as complements. 2017-08-06 13:14:42 +01:00
jzmaddock
1225612a28 quadrature: Remove unneeded test case, adjust one error rate, and finally silence some warnings. 2017-08-06 12:50:53 +01:00
jzmaddock
70e87cb407 Merge pull request #75 from boostorg/jzmaddock-patch-1 
Remove forward references to gcd/lcm types:
 2017-08-06 08:38:58 +01:00
jzmaddock
623263e726 Disable abs overload for __float128 when building with gcc-7 and later. 2017-08-03 11:38:16 +01:00
jzmaddock
e84195d3a6 Quadrature: Update docs. 
Also add ability to get the # levels consumed from the integrators.
[skip ci]
 2017-07-28 19:26:04 +01:00
jzmaddock
c3eb9e73d9 quadrature: move the tolerance into the integrate method from the constructor. 
Since nothing in the constructor depends on the tolerance, and we may wish to use different integration tolerances with the same integrator object.
Update the tests to match the new interface.  Also add "direct" integration methods which do not adapt the range at all, but assume that we are using the integrators native range.
 2017-07-25 18:26:21 +01:00
jzmaddock
27c1fa591f sinh_sinh: add pre-computed coefficients. 2017-07-23 18:50:38 +01:00
jzmaddock
6232615fd6 Remove forward references to gcd/lcm types: 
These are no longer part of Boost.Math, but moved to Integer.
 2017-07-21 18:19:56 +01:00
jzmaddock
397d8ce9c9 exp_sinh: add support for pre-computed constants. 2017-07-20 19:55:44 +01:00
jzmaddock
9d86bdefa6 Fix atomic namespaces. 2017-07-20 19:54:56 +01:00
jzmaddock
eb18bcd5aa tanh_sinh: add some missing using declarations. 2017-07-18 18:34:51 +01:00
jzmaddock
6259b8c795 tanh_sinh: add support for 2 arg user functions. 
Update tests accordingly and get everything passing locally.
 2017-07-16 19:18:20 +01:00
jzmaddock
3f9ae1db4e tanh_sinh: fix atomic/mutex usage to work on platforms that don't support the std:: versions. 2017-07-15 19:09:23 +01:00
jzmaddock
4334ced6a7 tanh_sinh: make the new version thread safe. 2017-07-15 17:07:59 +01:00
jzmaddock
4984809487 tanh_sinh: Add support for a minimum complement value, 
so we can prune the table size when going to extreme ends of the range isn't required.
 2017-07-15 16:47:55 +01:00
jzmaddock
80125ffb2c tanh_sinh: Hook up infinite limits code to use the available complements of abscissa values where appropriate. 2017-07-12 19:17:29 +01:00
jzmaddock
c986efa8f7 tanh_sinh: add precomputed coefficients for double/long double/__float128. 2017-07-11 18:45:28 +01:00
jzmaddock
2f837ee55a tanh_sinh: add back pre-computed coefficients for float case. 
Also refactored code back into one header/class.
 2017-07-10 19:19:46 +01:00
jzmaddock
4a5871ac71 tanh_sinh: Rework base class to cache refinements on demand. 2017-07-10 17:07:23 +01:00
jzmaddock
18938213e7 tanh_sinh: begin transforming code to calculate and store complements of abscissa values, 
and use these in the range transformations so that there is no cancellation error when transforming to an abscissa value very close to zero.
 2017-07-09 18:22:21 +01:00
jzmaddock
8b40c2cbb1 quadrature: Fix cpp_dec_float tests. 2017-06-26 09:02:09 +01:00
jzmaddock
6f96c8b229 quadrature: Get tests clean with older compilers. 
Also make tanh_sinh constants static members since they're constants.
 2017-06-25 11:44:08 +01:00
jzmaddock
040feb0240 Get the tests hooked up and working with a variety of compilers. 2017-06-24 09:48:32 +01:00
pabristow
bf1b8e8e10 Corrected mismerge to use my new version of constants include exp_minus_one constant. 2017-06-22 18:15:18 +01:00
pabristow
555e0f416b commit merge from develop 2017-06-22 17:05:02 +01:00
jzmaddock
b55167446f Quadrature: Policify the code, 
Make tanh_sinh quadrature arbitrary precision.
Rework tests so the object files are not so large.
Added some more example integrals to the tanh_sinh tests.
 2017-06-21 18:29:46 +01:00
jzmaddock
5e3fd1574e Merge branch 'double_exponential' of https://github.com/NAThompson/math into pr66 2017-06-17 18:10:25 +01:00
jzmaddock
9a8324c989 Merge pull request #71 from boostorg/pr65 
PR based on #65 prepared for CI testing
 2017-06-17 09:02:56 +01:00
jzmaddock
fdd959a32f Merge branch 'cleanup' of https://github.com/NAThompson/math into pr65. 
Also made a couple of  minor modifications.
 2017-06-16 19:35:13 +01:00
jzmaddock
34d413e830 trapezoidal: More fixes prior to merging: 
Add new documentation section for integrators.
Remove isfinite from std_real_concept.hpp - it's not part of our conceptual requirements.
Add workarounds for older compilers that can't handle the default arguments to the function.
Restrict testing test_trapezoidal.cpp to compilers that handle lambdas.
Make the concept checks C++03.
Make test_trapezoidal.cpp compile with VC11's broken lambda support.
 2017-06-16 18:13:42 +01:00
jzmaddock
f435cd3451 Bring trapezoidal.hpp more into line with the rest of Boost.Math: 
Add policy support for error handling and epsilon.
Fix some warnings in the test case (msvc).
 2017-06-15 19:14:31 +01:00
jzmaddock
a8bb78d25b Merge branch 'romberg' of https://github.com/NAThompson/math into pr59 
# Resolved Conflicts:
#	test/Jamfile.v2
 2017-06-15 18:26:48 +01:00
jzmaddock
cf0668eb62 Fix bessel N functions to behave correctly with decimal number types. 2017-06-13 18:54:48 +01:00
jzmaddock
4b90631b84 next.hpp: workarounds for std libs with no max_digist10 
Only invoke numeric_limits<>::max_digits10 for decimal multiprecision types.
Also make sure our traits classes are instantiated on the actual number type, not an expression template.
 2017-06-11 19:49:43 +01:00
jzmaddock
9b2f878103 Update next.hpp and ulp.hpp to work with decimal types. 2017-06-10 11:11:58 +01:00
jzmaddock
5c3624c0b8 Address CI failures: 
* Fix bjam command line in appveyor file.
* Fix assertions triggered in bessel_i0.hpp.
* Remove common factor testing - it's moved to integer.
* Add lots of additional {} in initializers to silence clang warnings.
 2017-06-04 13:12:33 +01:00
jzmaddock
0242c64eec float_distance: Fix cases that span several orders of magnitude. 
Add more tests.
See https://github.com/boostorg/multiprecision/issues/22.
 2017-06-02 19:40:52 +01:00
jzmaddock
37783e519a Bessel_I0: decimal number types shouldn't use numeric_limits::digits to judge precision. 2017-05-21 08:49:06 +01:00
Nick Thompson
23403b61a1 Add concept test and include test. 2017-05-18 19:02:29 -06:00
Nick Thompson
a980bd8af3 Add concept and include test. 2017-05-18 18:23:01 -06:00
Nick
1753d6e9e4 Merge branch 'develop' into romberg 2017-05-18 17:34:52 -06:00
Nick Thompson
7660162d41 Add concept and include test. 2017-05-18 17:33:45 -06:00
jzmaddock
bc37badb91 Change signbit to normalise the result when using std::signbit. 2017-05-15 18:28:32 +01:00
Nick Thompson
78409954f2 Use quadrature namespace. Add example. 2017-05-11 21:51:35 -06:00
Nick Thompson
4801e2d8bf Add 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. 2017-05-11 21:03:00 -06:00
Nick Thompson
3bfeee838d Double exponential quadrature. 2017-05-10 22:03:21 -06:00
pabristow
369ce4312b Halley/Luu version working for 50 decimal digits 'reference' test values. 2017-05-10 18:21:35 +01:00
jzmaddock
5db9e8870b Legendre_stieltjes: hook up tests. 
Add concept checks.
Fix some warnings, and reduce dependencies.
 2017-05-09 19:13:44 +01:00
jzmaddock
83891e2b79 Merge branch 'legendre_stieltjes' of https://github.com/NAThompson/math into stieltjes 2017-05-09 18:36:53 +01:00
jzmaddock
27ded3276d Remove C++11 isms from, legendre_prime header and tests - we don't want to break existing C++03 code using that header. 2017-05-09 12:18:43 +01:00
jzmaddock
97a4184df9 Hook up concept checks for legendre_p_prime. 
Fix errors from above, add domain-error handling to the function.
 2017-05-08 19:47:02 +01:00
jzmaddock
150a0b5afa Merge branch 'legendre_derivatives' of https://github.com/NAThompson/math into legrende_prime 
# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.
 2017-05-08 19:14:16 +01:00
jzmaddock
1d0880a9ef Remove dependencies on format and type_index. 2017-05-06 13:16:31 +01:00
jzmaddock
87ab68456d Improve interpolator documentation and examples. 2017-05-06 12:26:37 +01:00
jzmaddock
6e8687812f Add concept checks for barycentric_rational. 2017-05-05 13:12:25 +01:00
jzmaddock
80e7f7027c Merge branch 'barycentric' of https://github.com/NAThompson/math into interpolate_rat 
Fix some msvc failures.
Suppress some compiler warnings.
Added iterator based constructor.
Moved docs to better position.

# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.
 2017-05-05 13:01:56 +01:00
jzmaddock
d300f98fdc Fix msvc compiler errors in barycentric code. 2017-05-04 13:27:41 +01:00
jzmaddock
adc2234d95 Merge branch 'barycentric' of https://github.com/NAThompson/math into interpolate_rat 
# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.
 2017-05-04 13:18:05 +01:00
Nick Thompson
cf35d745d1 Remove 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. 2017-05-03 10:13:18 -06:00
jzmaddock
b9362cb000 Add concept checks for cubic-spline code, and fix resulting errors. 2017-05-03 13:32:53 +01:00
jzmaddock
127178518d Add iterator based construction to cubic B spline. 2017-05-03 13:07:54 +01:00
Nick Thompson
45226c7ef1 Numerical differentiation by finite differences and the complex step derivative. 2017-05-02 18:29:03 -06:00
jzmaddock
3141e075ed Fix min/max usage. 2017-05-02 20:06:53 +01:00
jzmaddock
9dd532e01a Merge branch 'develop' of https://github.com/NAThompson/math into interpolate 2017-05-02 19:30:31 +01:00
jzmaddock
6407bedb75 Bring polynomial gcd into line with new Boost.Integer code, and improve tests. 2017-05-02 19:12:18 +01:00
Nick Thompson
48a0cf714b Put license information in each file. Do not throw when condition number is approaching an undefined value; use mollification. Improve documentation. 2017-05-02 12:02:02 -06:00
Nick Thompson
2bbadc69fb Place copyright notices in file. 2017-05-01 18:36:26 -06:00
Nick Thompson
4b4c4e9364 Legendre-Stieltjes polynomials. 2017-04-30 15:58:10 -06:00
jzmaddock
ceae76a526 Switch to using Boost.Integer gcd/lcm. 2017-04-24 19:26:23 +01:00
Nick Thompson
a918d0f331 Return all zeros of the Legendre polynomial at once; there seems to be no benefit to having the user calculate them one at a time. 2017-04-22 15:47:22 -05:00
jzmaddock
1c10933943 polynomial.hpp should not be including private Boost.Config header. 2017-04-20 19:23:02 +01:00
jzmaddock
fafda3176c Fix gcd tests and code for clang on Linux 2017-04-20 14:43:28 +01:00
jzmaddock
82e3949024 GCD: Don't enable variadics unless the compiler supports them. 
Update docs.
Escape from variadic GCD early when the result hits unity.
 2017-04-15 18:31:01 +01:00
jzmaddock
c1ab95b6cc Add 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.
 2017-04-15 18:13:30 +01:00
Nick Thompson
4110a69416 Zeros 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. 2017-04-14 15:51:39 -05:00
Nick Thompson
21bcf34a51 Derivative of Legendre polynomials of the first kind. Used to generate the weights of Gaussian quadratures. 2017-04-12 14:16:26 -05:00
jzmaddock
70e82419d5 Fix for expression template types. 2017-04-08 20:12:23 +01:00
jzmaddock
f74da283aa Add variadic gcd and lcm variations, update tests to match. 2017-04-08 19:42:33 +01:00
jzmaddock
e71629f1c0 Make gcd/lcm functions constexpr and noexcept where applicable. 2017-04-08 19:13:08 +01:00
jzmaddock
8ee33649c7 Remove dead code orphaned by previous commits. 2017-04-07 09:55:31 +01:00
jzmaddock
5ca0e2ead8 Tighten up polynomial-GCD algorithm selection. 
Disallow GCD's of polynomials of bounded integers as overflow is too likely to occur.
 2017-04-07 09:44:48 +01:00
jzmaddock
8e98f30cb2 Remove SFINAE-logic and replace with simpler numeric_limits based algorithm selection. 2017-04-07 09:43:40 +01:00
jzmaddock
39006ecaca Reinstate gcd of rational numbers. 2017-04-06 19:22:45 +01:00
jzmaddock
7795ad3fd1 Merge branch 'gcd_polynomial' into gcd_cpp14 2017-04-06 11:28:32 +01:00
jzmaddock
9d0dac1df9 gcd: add the function objects back. 
Improve the gcd/lcm tests.
 2017-04-04 18:29:45 +01:00
jzmaddock
e4c0fdb0f2 Make gcd/lcm independent of type_traits and mpl. 
Later we'll try and reconcile this version with that in Boost.Integer.
 2017-04-02 20:10:46 +01:00
jzmaddock
8c9709b317 Fix cases where wchar_t is not an intrinsic type. 
Fixes: https://svn.boost.org/trac/boost/ticket/12944
 2017-03-30 19:35:59 +01:00
Nick Thompson
fff12aa16b Improve documentation for trapezoidal quadrature, explaining the role of the L1 norm calculation. Throw correct error type. 2017-03-28 13:56:18 -05:00
Nick Thompson
9f58f0d337 Add 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. 2017-03-26 13:45:50 -05:00
pabristow
c5ee664a51 Refactored with policies and passes tests and timing. 2017-03-23 14:15:37 +00:00
pabristow
bdee37743d Bug in values near to -exp(-1) corrected. 2017-03-13 10:03:33 +00:00
pabristow
0048d41dd6 Added new constants expminusone and sixth 2017-03-10 13:06:44 +00:00
pabristow
cea2465e3e Need expminusone constant 2017-03-06 18:23:54 +00:00
pabristow
94d3cf4043 refactored to use local test_value.hpp 2017-03-06 18:10:52 +00:00
Nick Thompson
4efe710855 Add 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. 2017-03-06 11:12:51 -06:00
Nick Thompson
136e7411f5 Adaptive Trapezoidal Quadrature 
This routine estimates the definite integral of a function f.
Assuming that f is periodic, it can be shown that this routine converges exponentially fast.
In fact, the test cases given exhibit exponential convergence with decreasing stepsize.

A potential improvement is using the Bulirsch sequence rather than the Romberg sequence to schedule the refinements.
However, the convergence is so rapid for functions of the class specified above that there seems to be no need at present.

This code is cppcheck clean, and runs successfully under AddressSanitizer and UndefinedBehaviorSanitizer.
 2017-03-05 19:05:41 -06:00
Nick Thompson
446f8c415d Cite Kress's book 'Numerical Analysis' in MLA format. 2017-03-05 10:49:22 -06:00
Nick Thompson
2506d7ed1f Use PIMPL so that the interpolant can be passed as a functor. 2017-03-02 17:16:21 -06:00
Nick Thompson
2b9d88d352 Fixes suggested by code review of PR #57. 2017-03-02 11:54:55 -06:00
Nick Thompson
b157403fd9 Do not use Kahan summation to compute average; use update procedure that cannot overflow recommended by Knuth. 2017-02-27 20:58:42 -06:00
Nick Thompson
fee20ab932 Given 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.

The interpolant is constructed in O(N) time and can be evaluated in constant time.
Its error is O(h^4), and obeys the interpolating condition s(x_j) = f(x_j) for all samples.
In addition, f' can be estimated from s', albeit with lower accuracy.

This routine is cppcheck clean, and is clean under AddressSanitizer and MemorySanitizer.
 2017-02-23 18:21:06 -06:00
jzmaddock
4c19a1ec34 Remove unneeded parameters in new beesel J/K code. 
Also fix some multiprecision-conceptual errors.
 2017-01-29 19:18:36 +00:00
jzmaddock
38eb9d3e38 Update docs for new Bessel K0 and K1 approximation code. 2017-01-24 18:35:54 +00:00
jzmaddock
fe66909b23 Reorder some declarations to keep GCC happy. 
Remove some dead code.
 2017-01-23 19:34:34 +00:00
jzmaddock
7922883e7b Add improved approximations for K0 and K1. 
Based on http://www.advanpix.com/2016/01/05/rational-approximations-for-the-modified-bessel-function-of-the-second-kind-k1-for-computations-with-double-precision/
 2017-01-23 17:45:51 +00:00
jzmaddock
2dc47a631a Improve constant initialization for I0 and I1. 2017-01-23 17:45:00 +00:00
jzmaddock
cc8e300bae Update I0 and I1 to use better approximations from Pavel Holoborodko 
See: http://www.advanpix.com/2015/11/11/rational-approximations-for-the-modified-bessel-function-of-the-first-kind-i0-computations-double-precision/
 2017-01-18 18:56:27 +00:00
pabristow
ffb025ca2c First 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. 2016-12-22 18:30:27 +00:00
pabristow
78320adb7a return 0 to T(0) in float_distance changed because fixed-point requires explicit construction 2016-12-05 16:41:18 +00:00
jzmaddock
1844134d07 Changes to call correct error handlers in C99 annex F cases. 2016-11-30 18:57:11 +00:00
zerotypos-found
3871f6a7e8 Fix issues in the inspection report 
- Remove Non-ASCII characters (i.e. replace dashes with hyphens).
- Protect min from macro expansion.
 2016-10-25 11:56:09 +09:00
jzmaddock
37fe39f2d9 Merge pull request #50 from kartikmohta/fix/boost-zero-polynomial-evaluate 
Fix boost::math::polynomial::evaluate for the zero polynomial
 2016-10-19 19:33:24 +01:00
jzmaddock
3370ea2f79 Add missing #include. 
Fixes https://svn.boost.org/trac/boost/ticket/12537.
 2016-10-19 19:27:37 +01:00
jzmaddock
21ab379ee9 Merge branch 'develop' of https://github.com/boostorg/math into develop 2016-10-19 19:22:51 +01:00
jzmaddock
4787d870c6 Suppress some msvc warnings. 2016-10-19 19:22:37 +01:00
Kartik Mohta
9dad3e2b62
Fix boost::math::polynomial::evaluate for the zero polynomial 2016-10-17 16:13:25 -04:00
zerotypos-found
4b9823768b Replace tabs to whitespace 2016-10-11 12:44:32 +09:00
jzmaddock
50e7ea3863 Sub-resultant GCD is only for integer types. 2016-08-26 18:40:11 +01:00
jzmaddock
accfe2f7eb Fix msvc compile failure, and add a using declaration so we overload boost::math::gcd. 2016-08-26 18:33:10 +01:00
jzmaddock
f888dcb129 Including polynomial.hpp should get the gcd overloads as well. 2016-08-26 18:32:27 +01:00
Jeremy W. Murphy
64c887a3c9 Add overload of gcd() for polynomials with multi-precision integral coeffs. 
And add a unit test to show that it works.
 2016-08-25 23:02:03 +10:00
Jeremy W. Murphy
d817274f3f Update some doxygen. 2016-08-25 23:00:59 +10:00
Jeremy W. Murphy
783dc6af46 Merge branch 'develop' into polynomial_ufd_gcd 2016-08-25 22:08:48 +10:00
inkstink
c9d9f446c5 fixed typos for real this time 2016-08-19 20:42:14 -04:00
inkstink
950d3990e4 fixed typos 2016-08-19 20:33:24 -04:00
pabristow
7786c7d5a8 Added Thomas Luu references 2016-08-09 16:05:20 +01:00
jzmaddock
53403a58f7 minor Intel fix. 2016-07-29 11:41:28 +01:00
Jeremy W. Murphy
3348f83fb7 Assert the loop invariant, degree(u) >= degree(v). 2016-07-27 17:47:03 +10:00
Jeremy W. Murphy
3b1f362a37 Let's not provide the subresultant algorithm as a default. 2016-07-26 19:07:51 +10:00
Jeremy W. Murphy
af290b272d The algorithm appears to assume that u is not smaller than v. 2016-07-26 19:07:21 +10:00
jzmaddock
d3a570e91a Fix Intel C++ test failures 2016-07-22 09:27:15 +01:00
Jeremy W. Murphy
2e024bd02d Merge branch 'polynomial_ufd_gcd' into polynomial_ufd_gcd_2 
# Conflicts:
#	test/test_polynomial.cpp
 2016-07-15 18:46:35 +10:00
Jeremy W. Murphy
579d6681bc Merge branch 'develop' into HEAD 
# Conflicts:
#	test/test_polynomial.cpp
 2016-07-15 18:42:53 +10:00
Jeremy W. Murphy
ce1d534414 Add case for zero to leading_coefficient and unit test. 2016-07-15 18:12:51 +10:00
jzmaddock
2165b95f52 Two more GCC 4.4 workarounds 2016-07-15 02:15:30 +01:00
jzmaddock
9127a2f707 Merge branch 'develop' of https://github.com/boostorg/math into develop 2016-07-15 01:31:17 +01:00
jzmaddock
12feba86fa Remove references to removed gcd/lcd code. 2016-07-15 18:36:55 +01:00
jzmaddock
cbb0cce89b Disable initializer_list support for GCC-4.4 as it doesn't fully support the new syntax and there appears to be no workaround. 2016-07-15 01:24:18 +01:00
jzmaddock
5d125381f6 Enable SSE2 code for more cases. 2016-07-14 13:24:09 +01:00
jzmaddock
b043126170 Add long double specializations where appropriate. 2016-07-14 13:23:44 +01:00
jzmaddock
3891523510 Fix tgamma_delta_ratio for cases where argument is very large compared to delta. 2016-07-14 13:22:16 +01:00
Jeremy W. Murphy
efa83033fa Merge branch 'develop' into polynomial_ufd_gcd 2016-07-12 21:39:49 +10:00
jzmaddock
5e54ce2427 Fix thread_local code: 
1) Use new Boost.Config macro.
2) Fix constants to cache value and initialize on startup when used with variable precision types.
3) Fix factorials to initialize on startup with UDT types.
4) Streamline some other variable precision code.
 2016-07-03 18:27:49 +01:00
jzmaddock
44a105aaf0 Tidy up clang/thread_local detection a little. 2016-06-27 19:42:09 +01:00
jzmaddock
3f0eaf7c46 Fix clang thread_local support detection 2016-06-27 18:54:38 +01:00
Jeremy W. Murphy
882dbaddd2 [polynomial] No need to call primitive_part() on v if it equals 1. 2016-06-23 11:43:56 +10:00
Jeremy W. Murphy
6f05d647ac [polynomial] Tidy up headers. 2016-06-23 11:28:23 +10:00
Jeremy W. Murphy
6f5537d273 [polynomial] Move gcd code into a new file. 2016-06-20 09:58:39 +10:00
Jeremy W. Murphy
95d4976b86 [polynomial] Rename gcd_ufd to subresultant_gcd. 2016-06-19 15:37:42 +10:00
Jeremy W. Murphy
ebc14aacd3 [polynomial] Correct the requirements and handling of zero for gcd_ufd. 2016-06-18 21:44:28 +10:00
Jeremy W. Murphy
21438ba03a gcd_ufd: neaten code and improve documentation. 2016-06-17 22:52:43 +10:00
Jeremy W. Murphy
fc109276f2 Be more specific about kinds of polynomials. 2016-06-16 23:50:43 +10:00
Jeremy W. Murphy
66740e77df Update code to reflect operator bool() now exists; overload primitive_part. 2016-06-16 23:49:17 +10:00
Jeremy W. Murphy
48e856cf34 Merge branch 'develop' into polynomial_ufd_gcd 2016-06-16 23:21:48 +10:00
jzmaddock
a96792d309 Merge branch 'develop' of https://github.com/boostorg/math into develop 2016-06-05 18:29:14 +01:00
jzmaddock
afb8e89642 Change selection of thread_local to be more conservative. 2016-06-05 18:28:56 +01:00
jzmaddock
e4cd194686 Fix misc clang issues 2016-06-04 19:26:59 +01:00
jzmaddock
5ce39f2673 squash warnings and forward policies 2016-06-03 15:11:34 +01:00
jzmaddock
e946c45e9e Fixes for cases where number type has no infinity and divide by zero is a hard error. 2016-06-01 19:25:03 +01:00
jzmaddock
5966b186f9 Add needed typecast to factorial-from-string multiprecision code 2016-06-01 00:09:49 +01:00
Daniel James
cd496c9ace Fix typo 2016-05-31 14:04:39 +01:00
jzmaddock
1cb146a001 Merge branch 'develop' of https://github.com/boostorg/math into develop 2016-05-31 13:36:01 +01:00
jzmaddock
7499075a01 Recalculate cached constants if the precision changes in zeta.hpp. 2016-05-31 13:33:51 +01:00
jzmaddock
f02c9cba13 Calculate lgamma near 1 and 2 via series expansion in terms of zeta functions of integers. 
Likewise for tgamma1pm1.
 2016-05-31 13:33:08 +01:00
jzmaddock
609e8373b8 We can calculate the limit as an int, no need to invoke FP code. 2016-05-31 13:31:52 +01:00
jzmaddock
6961f14eb3 Don't initialize erf's cached constants for variable-precision types. 2016-05-31 13:30:54 +01:00
jzmaddock
bd39b86ca1 Warning suppression. 2016-05-31 10:12:20 +01:00
jzmaddock
fb8528e47b Allow our cached data to be recalculated if the precision changes in unchecked_factorial.hpp. 2016-05-31 10:11:35 +01:00
jzmaddock
8b59c26e73 Warning disable. 2016-05-31 10:10:51 +01:00
jzmaddock
a97fedee55 Don't cache values for variable-precision types in precision.hpp 2016-05-31 10:10:07 +01:00
jzmaddock
cde7271817 Update bernoulli/polygamma code to avoid caching values if the precision has changed. 2016-05-31 10:09:15 +01:00
jzmaddock
04378a1a87 Warning disable in error_handling.hpp 2016-05-31 10:07:24 +01:00
jzmaddock
dd4db9a28c Update constants to work with types whose precision varies. 2016-05-31 10:07:01 +01:00
Jeremy W. Murphy
d922f95fd5 Merge branch 'develop' into polynomial_ufd_gcd 
# Conflicts:
#	include/boost/math/tools/polynomial.hpp
#	test/test_polynomial.cpp
 2016-05-16 09:48:15 +10:00
Jeremy W. Murphy
94526e54c7 Merge branch 'develop' into polynomial_zero 
# Conflicts:
#	test/test_polynomial.cpp
 2016-05-15 21:03:12 +10:00
Jeremy W. Murphy
e608dc6657 Rename clear() to set_zero(). 2016-05-13 18:27:31 +10:00
Jeremy W. Murphy
18a49862a9 Use typedef hack to prevent accidental conversion to an arithmetic type. 2016-05-12 21:02:34 +10:00
Kolya Matteo
51575f6474 Don't try pow to negative powers even for constants 
This failed when pow(T,int) cannot be found. It was failing for
boost::multiprecision types.
 2016-05-11 15:02:38 -04:00
Kolya Matteo
4c1a2eacef Use default policy for error handling 2016-05-11 13:50:49 -04:00
Kolya Matteo
15e163b2a5 Throw an exception, instead of asserting 2016-05-11 13:50:49 -04:00
Kolya Matteo
f5fc788b4f Fix merge 2016-05-11 13:50:49 -04:00
Kolya Matteo
141ef67b6a Add pow() function for polynomials and integer powers. 2016-05-11 13:50:49 -04:00
Jeremy W. Murphy
a737a3f5c2 operator bool() and clear(): test for zero and set to zero respectively. 
Checking for zero or non-zero and setting to zero are common operations,
so fast implementations
 2016-05-11 14:02:16 +10:00
jzmaddock
bea587f3fd Merge branch 'gcd_n' of https://github.com/jeremy-murphy/math into develop 
# Resolved Conflicts:
#	include/boost/math/common_factor_rt.hpp
#	test/test_gcd.cpp
 2016-05-09 18:59:17 +01:00
jzmaddock
c9bf1df944 Merge branch 'polynomial_shift' of https://github.com/jeremy-murphy/math into develop 
# Resolved Conflicts:
#	test/test_polynomial.cpp
 2016-05-09 18:29:22 +01:00
jzmaddock
cace9b59de Merge branch 'gcd' into develop 
# Resolved Conflicts:
#	test/test_polynomial.cpp
 2016-05-09 18:15:54 +01:00
jzmaddock
993cdcb42e Reduce dependencies by removing all use of Boost.Format. 
We made only cursory use of the lib, so it's just as easy to do with it.
Also fixed a number of missing #includes which were hidden by including format.hpp.
 2016-05-09 11:38:13 +01:00
Jeremy W. Murphy
a56745b479 Merge branch 'develop' into polynomial_ufd_gcd 
# Conflicts:
#	test/test_polynomial.cpp
 2016-05-07 18:32:42 +10:00
Jeremy W. Murphy
a822400320 Merge branch 'develop' into polynomial_shift 
# Conflicts:
#	test/test_polynomial.cpp
 2016-05-07 18:27:46 +10:00
Jeremy W. Murphy
e3101bc0eb Return pair of <T, I> that includes the termination position. 2016-05-07 14:33:35 +10:00
Jeremy W. Murphy
df8c6584dc Rename gcd_n to gcd_range. 2016-05-07 10:58:43 +10:00
Kolya Matteo
8336375c78 Fix corner cases for raise_error 
This enables raise_error to work when no %1% appears in the `function` or
`message` parameters and for types without numeric_limits specializations.
Also removes misformat when BOOST_NO_RTTI is defined.
 2016-05-05 18:10:43 -04:00
Kolya Matteo
430239841d Add missing normalization in polynomial initializer_list assignment 2016-05-02 14:08:25 -04:00
Kolya Matteo
589a18f63f Don't assume default-constructibility 2016-05-02 09:52:41 -04:00
Kolya Matteo
627a4cedfd Multiply into a new vector and swap 2016-05-01 23:55:28 -04:00
Kolya Matteo
9c8e3f2d4b Fix self-multiply-assign by backing down 2016-05-01 23:18:50 -04:00
Kolya Matteo
ca07aa8a0a Fix indefinite hang with p *= p 2016-05-01 20:47:54 -04:00
Jeremy W. Murphy
2168585e98 Compute delta as late as possible. 2016-04-27 12:13:49 +10:00
Jeremy W. Murphy
8e25e0a487 Improve comments, remove superfluous using directive. 2016-04-27 12:12:20 +10:00
Jeremy W. Murphy
60d6a05d2c Prefer immediate return to goto. 2016-04-22 13:32:34 +10:00
Jeremy W. Murphy
d484fa657c Remove inferior generalized gcd; include integral size promotion. 2016-04-20 22:53:56 +10:00
Jeremy W. Murphy
dc81bc8e6e Add 4.6.1E, generalized gcd. 
Should be useful for cross-validating 4.6.1C.
 2016-04-20 01:38:05 +10:00
jzmaddock
1ef4bb6490 Add mixed binary algorithm plus gcd_traits class with bitscan support. 2016-04-18 13:50:01 +01:00
Jeremy W. Murphy
aa19132d74 4.6.1C: Greatest common divisor over a unique factorization domain. 2016-04-17 14:21:01 +10:00
Jeremy W. Murphy
5bfe2a7b15 Efficient conversion to bool. 
This will save a lot of (x == zero_element<...>(...)) syntax that is not
only unwieldy but inefficient.
 2016-04-17 13:42:50 +10:00
Jeremy W. Murphy
d6dd40def6 Contextual comment and zero case for pp(). 2016-04-15 14:17:42 +10:00
Jeremy W. Murphy
573b9fc399 content and primitive_part of a polynomial over a ufd. 2016-04-15 13:14:18 +10:00
Jeremy W. Murphy
b3d398636c gcd_n: Further comments. 2016-04-15 07:29:38 +10:00
Jeremy W. Murphy
af28716e31 Algorithm description. 2016-04-14 22:58:46 +10:00
Jeremy W. Murphy
5b8c448ce0 gcd of n values. 2016-04-14 14:22:12 +10:00
jzmaddock
e9bafa3236 We need some using std::swap declarations. 2016-04-12 18:35:01 +01:00
Jeremy W. Murphy
e963c67af2 Normalize after left-shift assignment to prevent denormalized zeroes. 2016-04-10 01:27:36 +10:00
jzmaddock
07f5e1666b Get the tests passing with Jeremy's new gcd. 2016-04-06 19:16:47 +01:00
jzmaddock
9f5d39c5e1 Merge branch 'gcd_revamp' of https://github.com/jeremy-murphy/math into gcd 
# Fixed Conflicts:
#	reporting/performance/test_gcd.cpp
#	test/test_gcd.cpp
 2016-04-06 18:09:45 +01:00
jzmaddock
11500b787b Merge branch 'polynomial_shift' of https://github.com/jeremy-murphy/math into gcd 
# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.
 2016-04-06 17:58:21 +01:00
Jeremy W. Murphy
473cc9f209 Use std::numeric_limits rather than Boost's built-in type_traits. 
The Boost type traits are too strict about what qualifies as a signed
number.
 2016-04-05 09:48:15 +10:00
Jeremy W. Murphy
fd11b9a4da Use <<= rather than <<. 2016-04-05 09:44:28 +10:00
Jeremy W. Murphy
232a64424f Use unqualified swap to allow ADL. 2016-04-05 09:41:59 +10:00
Jeremy W. Murphy
11e502b4ac New gcd implementation. 2016-04-04 18:36:53 +10:00
Jeremy W. Murphy
cff8705458 Consistent white space. 2016-03-28 16:08:04 +11:00
Jeremy W. Murphy
238cd6078b polynomial: odd/even predicates. 
We use the definition of odd and even consistent with considering x as the
smallest prime factor. That is, an even polynomial has a constant of zero.
 2016-03-28 16:04:27 +11:00
Jeremy W. Murphy
e9da8a10bf polynomial: Left and right shift operators. 
Shifting adds or removes a factor of x in the same way that shifting adds
or removes a factor of 2 to integers.
 2016-03-27 22:03:33 +11:00
Jeremy W. Murphy
521ee8a0af Use the template parameter N and container's size_type. 2016-03-08 19:24:45 +11:00
jzmaddock
65f7440c2f Minor fixes for polynomial division: 
* Use Boost.Config macros not direct checks on C++ version.
* Use dedicated integer power function - fixes VC10 build failure, also prevents loss of digits when integer type is wider than double.
 2016-03-02 19:37:46 +00:00
jzmaddock
2ed1361e73 Merge branch 'polynomial_initializer_list' of https://github.com/jeremy-murphy/math into develop 2016-03-02 18:47:11 +00:00
jzmaddock
7cecd8ddf4 Fix up mpfr support 2016-02-05 13:59:45 +00:00
Jeremy W. Murphy
7969de6fae Add C++11 initializer_list constructor and assignment to polynomial class. 2016-01-31 19:04:37 +11:00
Jeremy W. Murphy
0fae5d590a Explicitly specify the return type for Boost lambda on * and /. 
Appears that return type deduction failed for multiprecision types.
 2016-01-16 14:51:35 +11:00
jzmaddock
8e7660e909 Merge branch 'develop' of https://github.com/boostorg/math into develop 2016-01-10 11:51:38 +00:00
jzmaddock
1adea9929f Add missing inline declarations. 
See https://svn.boost.org/trac/boost/ticket/11892.
 2016-01-10 11:51:24 +00:00
Kohei Takahashi
fbd1393858 Qualify std:: for isnan in some situation. 
Because isnan is implemented as a macro and libstdc++ undef it within
<cmath> (at least FreeBSD 10).
 2016-01-06 19:39:55 +09:00
jzmaddock
74ff2db959 Change <quadmath.h> config to not use it at all if we don't have __has_include as GCC may be configured with --disable-libquadmath but still signal that it supports __float128 2016-01-01 16:49:48 +00:00
jzmaddock
587eb660aa Merge branch 'develop' of https://github.com/boostorg/math into develop 2015-12-14 18:03:22 +00:00
jzmaddock
b04e78b79e Fix msvc warnings. 2015-12-14 18:03:05 +00:00
pabristow
7a7d3321fe Moved polystuff to tools, added JWM as author and numerous minor edits 2015-12-14 16:13:05 +00:00
jzmaddock
c36506f852 Merge pull request #17 from jeremy-murphy/polynomial_division 
Polynomial division
 2015-12-13 12:06:00 +00:00
Jeremy W. Murphy
25ed0ffbca Negation depends only on negation, not on multiplication. 2015-12-13 12:15:09 +11:00
jzmaddock
56f74cbfe9 Fix msvc warning. 2015-11-27 18:52:54 +00:00
Jeremy W. Murphy
09dd190e5a Move identity transformation into detail namespace. 2015-11-21 14:28:39 +11:00
Jeremy W. Murphy
6773baf318 Move arithmetic implementation into private functions. 
Primarily so that they can be used inside loops where normalization is not
required or is even erroneous.
 2015-11-21 14:26:00 +11:00
Jeremy W. Murphy
f8d6220e4e More documentation updates and parameter name consistency. 2015-11-20 12:06:13 +11:00
Jeremy W. Murphy
292bfc6752 Denormalize zero in multiplication to fit algorithmic requirements. 
Fixes bug introduced by normalizing zero to the empty set.
 2015-11-20 11:33:00 +11:00
Jeremy W. Murphy
e6e1e9440b Make the single-argument constructor explicit. 
Prevent unexpected conversion to polynomial!
 2015-11-18 12:15:35 +11:00
Jeremy W. Murphy
37e8063420 Add missing heterogeneous division operators. 2015-11-18 12:14:45 +11:00
jzmaddock
477ab19672 Restrict use of __float128 to cases where quadmath.h is available 2015-11-16 18:35:16 +00:00
Jeremy W. Murphy
87268a735a Division of polynomials over integers. 
Test int and double throughout most of test suite.
 2015-11-16 00:41:37 +11:00
Jeremy W. Murphy
b981995cc0 Replace synthetic division with algorithm D from Knuth. 
Simpler, specific to division of polynomials over a field.
 2015-11-15 14:08:40 +11:00
jzmaddock
e8c74af768 Fix error in cyl_neumann(2, x) introduced by previous reorganisation. 2015-11-11 19:19:11 +00:00
jzmaddock
3f342bfa60 More MacOS fixes: apply fixes from main template to <double> specialization. 2015-11-11 17:49:50 +00:00
jzmaddock
7a84f91de1 Merge branch 'develop' of https://github.com/boostorg/math into develop 2015-11-10 12:47:26 +00:00
jzmaddock
824c28f09f Another try at getting MacOS tests passing. 2015-11-10 12:47:06 +00:00
Jeremy W. Murphy
94e961ec7c Division and modulus by a number. 
Test that multiplication by 0.5 is equal to division by 2.
 2015-11-09 00:43:29 +11:00
Jeremy W. Murphy
3c147eb77a Revert "is_field_like checks that type T is not integral, nor complex, and has /=." 
This reverts commit 9c6d10644b.
 2015-11-09 00:15:48 +11:00
Jeremy W. Murphy
9c6d10644b is_field_like checks that type T is not integral, nor complex, and has /=. 2015-11-07 22:36:38 +11:00
Jeremy W. Murphy
115585481c Merge remote-tracking branch 'refs/remotes/jeremy/polynomial_division' into polynomial_division 2015-11-05 00:02:31 +11:00
Jeremy W. Murphy
ec913ea3df Multiplicative identity element. 2015-11-05 00:02:24 +11:00
Jeremy Murphy
4510ecaa8a Disable division and remainder for integer types. 
Use std::numeric_limits<T>::is_integer.
Also removed the "unchecked_" prefix from synthetic_division.
 2015-11-04 20:17:06 +11:00
Jeremy W. Murphy
cb36a289e3 Simplify addition: there is actually no need to check for zero. 
The likely frequency of zero makes it silly to check for it.
 2015-11-04 12:24:31 +11:00
Jeremy W. Murphy
72f35334d3 Unary minus (negate) operator. 2015-11-03 23:57:11 +11:00
Jeremy W. Murphy
babd53e857 Multiplication bugfixes for zero. 2015-11-03 23:56:53 +11:00
Jeremy W. Murphy
cba0dc6e7e Throw logic_error for degree() if size() == 0. 
Use size() in place of degree() internally.
 2015-11-03 23:55:37 +11:00
Jeremy W. Murphy
066885dc5a Whoops, need sign function for addition of constant values too. 2015-11-03 14:43:03 +11:00
Jeremy W. Murphy
b39391cb02 Change exception to assertion (according to my design philosophy). 2015-11-02 23:42:13 +11:00
Jeremy W. Murphy
fed9a626be Rename linear_translation to addition. 2015-11-02 23:24:06 +11:00
Jeremy W. Murphy
964c49c70e Generalize addition and subtraction of constant values. 
Simplify the implementation by not checking for infrequent cases.
 2015-11-02 23:23:47 +11:00
Jeremy W. Murphy
38edca726d Rename remove_high_degree_zeroes to normalize. 2015-11-02 23:04:32 +11:00
Jeremy W. Murphy
eaa54b1bab remove_high_degree_zeroes(): trim zeroes from the right of m_data. 2015-11-01 16:01:26 +11:00
Jeremy W. Murphy
2479c2ae90 Fixed some installation problems, now lambda works. 2015-11-01 13:25:46 +11:00
Jeremy W. Murphy
f6d014d15e Tidying up and trying out algorithmic forms of loops. 2015-11-01 12:49:24 +11:00
Jeremy W. Murphy
739c056f2b Enforce representation of zero polynomial as {}, not {0}. 
Simplify addition and subtraction into one function.
Add short-circuit logic for zero to arithmetic.
 2015-11-01 03:38:13 +11:00
Jeremy W. Murphy
6080afee77 Fix quotient_remainder to use the correct ordering of coefficients. 
Add equality operator and zero_element() for multiplication.
 2015-11-01 03:28:59 +11:00
pabristow
57a71ba5f8 #11768: Skewness formula for triangular distribution corrected, tests added and docs updated. 2015-10-29 18:19:46 +00:00
Jeremy W. Murphy
86a92e2602 Use Boost.Operators for dividable and modable. 2015-10-27 08:53:51 +11:00
jzmaddock
d6e6ec366e Move some code around and simplify things a bit to improve performance. 2015-10-26 18:16:14 +00:00
jzmaddock
3eed7acedf Fix argument names. 2015-10-26 18:14:59 +00:00
Jeremy W. Murphy
1fd285721c Add some post-condition assertions. 2015-10-25 14:55:39 +11:00
Jeremy W. Murphy
240e4e6e0f Separate out into unchecked_synthetic_division(). 2015-10-25 12:05:32 +11:00
Jeremy W. Murphy
35041b318c Add / and %. 2015-10-25 10:49:47 +11:00
Jeremy W. Murphy
cf840c4ce2 First draft of working quotient_remainder. 2015-10-25 01:46:20 +11:00
jzmaddock
a7a0117859 Tentative VC7.1 fix. 2015-10-22 10:51:04 +01:00
jzmaddock
5d42feb5c1 Reduce conceptual requirements on log_min_value/log_max_value. 2015-10-22 10:12:55 +01:00
jzmaddock
1690851386 Fix for __float128 support in fpclassify. 2015-09-24 10:18:29 +01:00
jzmaddock
2f31042470 Merge branch 'develop' of https://github.com/boostorg/math into develop 2015-09-11 12:02:29 +01:00
jzmaddock
3287c922ea Fix some corner cases in the beta functions. 
Add better tests for the incomplete beta derivative.
Update docs to match.
 2015-09-11 12:02:04 +01:00
Adam Wulkiewicz
00714c1a6d Fix type passed into BOOST_MATH_NOEXCEPT() in specialization of epsilon(). 2015-09-09 15:06:29 +02:00
jzmaddock
d719e4b567 Remove no longer needed special_functions sub-directory. 
Also fix lots of inspect-report failures.
 2015-09-07 19:35:29 +01:00
jzmaddock
c6593b149b Fix GCC test failures when building with -fno-exceptions. 2015-09-06 19:21:17 +01:00
jzmaddock
ba365b5aaa Enable running of tests with exception handling off. 
And fix a few bugs that showed up in the process.
Currently msvc-only (not GCC as Boost.Test doesn't build).
 2015-09-05 19:08:28 +01:00
jzmaddock
a0f0b96921 Fix case where non-centrality parameter for chi-squared is zero. 
See https://svn.boost.org/trac/boost/ticket/11557.
 2015-09-04 09:40:08 +01:00
jzmaddock
9d3d2e96a0 Fix code to match documentation in hypergeometric distro. 2015-09-03 19:23:10 +01:00
jzmaddock
789131a224 Add some needed typecasts to args of std::max. 2015-09-03 18:43:36 +01:00
jzmaddock
c2b706da15 Add __float128 support. 2015-09-03 18:43:05 +01:00
jzmaddock
e4ec61699a Fix tests for FLT_EVAL_TYPE. 
See https://svn.boost.org/trac/boost/ticket/11429.
 2015-09-03 18:13:20 +01:00
jzmaddock
1d1d400762 Add some needed typecasts. 2015-09-03 10:29:28 +01:00
jzmaddock
429825535b Fix type of boolean values. 2015-09-03 09:34:42 +01:00
jzmaddock
7dcabda78b Merge branch 'accuracy' into develop 
Conflicts:
	test/float128/test_polygamma.cpp
	test/powm1_sqrtp1m1_test.hpp
 2015-09-03 09:12:11 +01:00
jzmaddock
673c229a6c More constexpr and noexcept specifiers. 2015-09-01 17:34:50 +01:00
jzmaddock
146f0471ec Fix previous commit! 2015-08-29 11:21:55 +01:00
jzmaddock
3d8d58262d Add workaround for early/broken constexpr implementations. 2015-08-29 09:20:28 +01:00
jzmaddock
57219acee9 Make more functions constexpr and noexcept. 2015-08-28 12:27:41 +01:00
jzmaddock
06cff2b46f Merge branch 'accuracy' of https://github.com/boostorg/math into accuracy 2015-08-27 10:43:25 +01:00
jzmaddock
4c278283a7 More noexcept specs. 2015-08-27 10:42:45 +01:00
jzmaddock
dac39d03e1 Change noexcept to be less greedy in instantiating numeric_limits: fixes some multiprecision failures 2015-08-25 09:04:29 +01:00
jzmaddock
bf29ef6628 Begin marking up functions as constexpr and noexcept where possible. 2015-08-23 17:50:56 +01:00
jzmaddock
8f252f78a0 Building with hidden scope doesn't work on Oracle. 
Also added needed typecast for Oracle 64-bit build.
 2015-08-20 08:53:00 +01:00
jzmaddock
ba98215010 Oracle compiler needs GNU align syntax. 2015-08-17 19:07:37 +01:00
jzmaddock
277dbd4f5c Add compiler comparison option. 2015-08-17 13:05:50 +01:00
jzmaddock
8a29241cc3 Fix bug in root bracketing code 2015-08-11 18:17:56 +01:00
jzmaddock
3e67c80eeb Improve quantile calculation so as to avoid underflow causing lots of PDF calls and poor performance. 2015-08-09 17:29:34 +01:00
jzmaddock
c3793dbff2 One more try to get powm1 accuracy down. 2015-08-02 18:26:47 +01:00
jzmaddock
801a739595 Fix GCC warnings. 2015-08-01 19:13:22 +01:00
jzmaddock
258a2aeb4c Big warning suppression fix for VC14. 2015-08-01 18:00:41 +01:00
jzmaddock
a18661ad40 Tweak powm1 one more time near the boundary conditions. 2015-07-31 19:19:47 +01:00
jzmaddock
da958ee36a Some more tweaks to powm1. 
Plus some more tests.
 2015-07-31 09:25:31 +01:00
jzmaddock
ce953afcda Yet another "better powm1". 2015-07-28 19:23:53 +01:00
jzmaddock
1177477903 Remove deprecated type_traits usage. 2015-07-25 11:45:49 +01:00
jzmaddock
e973194d55 Fix some corner cases in the non-central distros. 2015-07-18 16:41:28 +01:00
jzmaddock
541f0d2839 Fix method selection condition in powm1. 2015-07-13 10:39:35 +01:00
jzmaddock
74b8427084 Fix domain_error for k==1. 2015-07-11 09:49:17 +01:00
jzmaddock
869bf00304 Improve powm1 accuracy slightly. 2015-07-11 09:16:55 +01:00
jzmaddock
bab82f2eb1 Correct overflow handling in bessel derivatives. 2015-07-10 19:50:05 +01:00
jzmaddock
43de5afcf2 Oops correct previous commit. 2015-07-07 18:45:45 +01:00
jzmaddock
faacf0520d Fix function template specialization definition. 2015-07-07 18:42:01 +01:00
jzmaddock
d9dc112aee Fix bug in ellint_2 for very small angles. 
With thanks to Rocco Romeo.
 2015-07-05 19:42:31 +01:00
jzmaddock
32b31bbb72 Add explicit overloads for std::pair to handle C++03 compilers. 2015-06-19 10:05:12 +01:00
jzmaddock
aedcfb355b Fix template name used only on Darwin. 2015-06-18 18:35:15 +01:00
jzmaddock
2a94e603e5 Add relative_difference tests. 
And fix resulting issues.
 2015-05-26 17:17:47 +01:00
jzmaddock
afc1237939 Stop using modfl on Mingw - it crashes! 2015-05-25 19:20:09 +01:00
jzmaddock
69b2524d13 Document ulp function. 
Improve tests, and fix issues that show up.
 2015-05-25 08:47:34 +01:00
jzmaddock
94b8195aad Add ulp function. 
Update tests for relative_difference etc.
Add new functions to math_fwd.hpp.
Add new functions to concept tests.
 2015-05-24 12:01:06 +01:00
jzmaddock
73c27a1a6d Rewrite float-comparison docs to use new public interfaces. 2015-05-23 12:06:36 +01:00
jzmaddock
e7b5c6a0be Add relative_difference.hpp 2015-05-18 18:24:26 +01:00
jzmaddock
74a42f98ea Fix various typos. 
Change Schroeder to Schroder.
 2015-05-18 18:23:47 +01:00
jzmaddock
8aaf48ecaa Add test case for root-finding. 
Also fix missing using declaration in roots.hpp.
 2015-05-14 19:02:40 +01:00
jzmaddock
66296fc6bd Fix bug in newton_raphson_iterate. 
Which may cause switching to bisection too early.
 2015-05-12 19:32:59 +01:00
jzmaddock
0b166f37ab Merge branch 'develop' into rooting 2015-05-11 19:08:00 +01:00
jzmaddock
764738d596 Change root finding code to be better boiler-plated. 
Re-run tests accordingly.
 2015-05-10 10:06:08 +01:00
jzmaddock
308f39c78a Fix iteration count housekeeping so we don't accidentally miss one if we get an exact result because f(x) is zero. 2015-05-05 13:15:23 +01:00
jzmaddock
c8eb3c9cf4 Add default constructor for eps_tolerance. 2015-04-29 18:47:09 +01:00
jzmaddock
da33ba4425 Fix heuristic in schroeder_iterate. 
Brings schroeder_iterate more into line with halley_iterate.
 2015-04-28 19:16:33 +01:00
Jessica Hamilton
d4240296b5 Haiku: disable float128 and long double support 2015-04-09 09:38:27 +12:00
jzmaddock
2f09aca3cd Merge branch 'develop' of https://github.com/boostorg/math into develop 2015-04-07 18:01:54 +01:00
jzmaddock
0cbed0315d Fix digit count for 128-bit quad-floart types. 
As per discussion here: http://www.exploringbinary.com/7-bits-are-not-enough-for-2-digit-accuracy/comment-page-1/#comment-9358
 2015-04-07 17:57:07 +01:00
Lauri Nurmi
0c273d09a9 Fix spelling of "occurr*". 2015-03-28 22:11:16 +02:00
pabristow
1f0552daea merged develop into rooting and rebuilt html 2015-03-12 12:31:01 +00:00
jzmaddock
c6e152efd4 Disable __float128 on Solaris - it segfaults when used at program startup (tested solaris 11). 2015-03-08 18:58:15 +00:00
jzmaddock
b423b350a8 Add needed member. 2015-02-28 19:42:39 +00:00
jzmaddock
97cd1ea0ea Change not to use cstdfloat as that header may not actually define any types! 2015-02-28 19:31:05 +00:00
jzmaddock
5e92ce7ae6 Solaris doesn't have fpclassify 2015-02-26 11:55:23 +00:00
jzmaddock
58edb99ca8 Workaround for bug in Solaris compiler. 2015-02-23 17:15:23 +00:00
jzmaddock
1ef4e9c1c2 Tentative fix for solaris/oracle compiler issue https://svn.boost.org/trac/boost/ticket/10925 2015-02-22 18:48:09 +00:00
jzmaddock
ec7a31e982 Merge branch 'carlson' into develop 2015-02-20 13:30:24 +00:00
jzmaddock
7fa11a133f Add Heuman Lambda Function. 2015-02-19 19:33:56 +00:00
jzmaddock
072d6ef395 Add Jacobi Zeta Function. 2015-02-19 16:49:03 +00:00
jzmaddock
026d5d2ee2 Hook up a full set of tests for ellint_d. 2015-02-18 18:04:02 +00:00
jzmaddock
c1200bf7cd Add initial version of Elliptic Integral D. 2015-02-17 19:03:50 +00:00
jzmaddock
bf3f5528eb Merge branch 'carlson' into develop 2015-02-17 13:08:21 +00:00
jzmaddock
4705e25f70 Fix some compiler warnings 2015-02-17 12:43:54 +00:00
jzmaddock
61cc79bebd Fix test failures. 2015-02-17 12:43:40 +00:00
jzmaddock
5eb3bbbfcc Apply Carlson's latest relations for the Legendre Elliptic integrals. 2015-02-16 17:49:56 +00:00
jzmaddock
f7c55524be Add missing #include. 2015-02-15 18:57:14 +00:00
jzmaddock
110a5327b2 Apply more special cases to RD and RG plus tests. 2015-02-15 18:48:33 +00:00
jzmaddock
cbb1d9cb5f Add RF(x, y, 0) as a special case. 2015-02-15 11:34:44 +00:00
jzmaddock
545961746c Fix for issue https://svn.boost.org/trac/boost/ticket/10938 2015-02-05 19:02:52 +00:00
jzmaddock
20965d162e Tidy up dead code, make comment a little clearer. See https://svn.boost.org/trac/boost/ticket/10985 2015-02-05 18:34:49 +00:00
jzmaddock
529190e194 Merge branch 'master' of https://github.com/alexhenrie/math into develop 2015-02-05 18:10:30 +00:00
Alex Henrie
a4664e00d5 Correct spelling of "suppress" 2015-01-30 15:28:30 -07:00
jzmaddock
63525295e4 [Carlson RG] Fix test failures and concept-checking failures. 2015-01-21 17:46:14 +00:00
jzmaddock
b358c793d5 [Carlson Elliptic Integrals] Add doc for RG. 2015-01-08 18:15:18 +00:00
jzmaddock
11d4d7eb11 [Carlson elliptic Integrals] Add ellint_rg plus tests. 2015-01-08 13:43:37 +00:00
jzmaddock
2ba0b25c16 [Carlson Elliptic Integrals] Fix #includes and conceptual errors. Tweak RD convergence and error rates. 2015-01-08 09:47:20 +00:00
jzmaddock
64bf08fce4 [Carlson Elliptic Integrals] Add special cases for RF. 2015-01-07 19:00:25 +00:00
jzmaddock
ef2312dfd5 [Carlson elliptic integrals] Fix missing return statement. 2015-01-07 18:38:11 +00:00
jzmaddock
89fcff4266 [Carlson Elliptic Integrals] Update RD to use Carlson's latest methods. Add more RD test cases. Update copyright info. 2015-01-07 18:24:08 +00:00
jzmaddock
9561623736 [Beta/Gamma]Fix concept failures caused by previous commit. 2015-01-06 18:06:10 +00:00
jzmaddock
08a365a1f5 [legendre polys] Disable policy based overloads when the final argument is not in fact a policy. See https://svn.boost.org/trac/boost/ticket/10902 2015-01-05 19:08:56 +00:00
jzmaddock
b6e56ccd22 [constants] Fix namespace scope constants so they are constexpr, see https://svn.boost.org/trac/boost/ticket/10901 2015-01-05 18:32:30 +00:00
jzmaddock
5bd2b77ec4 [gamma_p_derivative] Fix spurious underflow spotted by Rocco Romeo. 2015-01-05 18:16:51 +00:00
jzmaddock
bc3d8ed2e3 [Incomplete Beta] Fix some more spurious underflows spotted by Rocco Romeo. 2015-01-05 17:39:13 +00:00
jzmaddock
b3d11e9545 [Ellint Carlson] Update RJ to carlson's latest algorithm. Add missing include to ellint_rc.hpp. Add more test cases for RJ special cases. 2015-01-05 11:46:34 +00:00
jzmaddock
fe6d89db43 [tgamma] Apply fix from https://svn.boost.org/trac/boost/ticket/10905. 2014-12-31 18:18:37 +00:00
jzmaddock
669214556a [ellint RC] rewrite to use elementary functions only. 2014-12-30 13:24:41 +00:00
jzmaddock
04abb9d194 [ellint RF] rewrite to use carlson's latest methods. 2014-12-30 13:24:08 +00:00
jzmaddock
1def3942bf [ellint 3] Move k==0 case to earlier in the logic 2014-12-30 13:23:11 +00:00
jzmaddock
4c9624786c [Ellint Pi] Change v < 0 case to use Rocco Romeo's simplified expressions. 2014-12-27 18:57:32 +00:00
jzmaddock
e6b86ba93d [Ellint Pi] Fix spurious underflow in sqrt(-v*N). 2014-12-27 13:30:51 +00:00
jzmaddock
47c2f9254c [Ellint Pi] Add some more special case handling, plus tests. 2014-12-26 17:21:10 +00:00
jzmaddock
18dd27295d [Ellint Pi] Reorganize which reduction occurs in which order, and fix numerical stability issues identified by Rocco Romeo. 2014-12-25 16:52:37 +00:00
jzmaddock
83155a7a67 [Elliptic Pi] Change expression to not suffer spurious under/over flow, add test case. 2014-12-22 16:55:06 +00:00
jzmaddock
2ab808f547 [Ellint Pi] Try and actually use the correct equation for negative v in Pi[v, k] ! 2014-12-21 19:08:18 +00:00
jzmaddock
b885aa757a [bernoulli] suppress some compiler warnings. 2014-12-21 13:10:16 +00:00
jzmaddock
57b38b0dfc Add binomial_coefficient to math_fwd.hpp and test. 2014-12-21 13:09:40 +00:00
jzmaddock
5af7e8fdb8 [incomplete gamma] Fix corner cases identified by Rocco Romeo. 2014-12-20 19:18:25 +00:00
jzmaddock
0c1a603e17 [incomplete beta] Fix corner cases identified by Rocco Romeo. 2014-12-20 18:18:16 +00:00
jzmaddock
160beaa477 [Ellint Pi] Fix corner cases identified by Rocco Romeo and add tests for them. 2014-12-20 09:47:06 +00:00
jzmaddock
50bffe8005 Update conceptual requirements for root finding algorithms, and add concept test for them. 2014-12-14 09:44:03 +00:00
pabristow
ee8d6338fc Merge branch 'develop' into arcsine_dist 
Conflicts:
	doc/html/backgrounders.html
	doc/html/constants.html
	doc/html/cstdfloat.html
	doc/html/dist.html
	doc/html/extern_c.html
	doc/html/gcd_lcm.html
	doc/html/index.html
	doc/html/indexes.html
	doc/html/indexes/s01.html
	doc/html/indexes/s02.html
	doc/html/indexes/s03.html
	doc/html/indexes/s04.html
	doc/html/indexes/s05.html
	doc/html/inverse_complex.html
	doc/html/math_toolkit/acknowledgement.html
	doc/html/math_toolkit/acknowledgements.html
	doc/html/math_toolkit/acos.html
	doc/html/math_toolkit/acosh.html
	doc/html/math_toolkit/airy.html
	doc/html/math_toolkit/airy/ai.html
	doc/html/math_toolkit/airy/aip.html
	doc/html/math_toolkit/airy/bi.html
	doc/html/math_toolkit/airy/bip.html
	doc/html/math_toolkit/archetypes.html
	doc/html/math_toolkit/asin.html
	doc/html/math_toolkit/asinh.html
	doc/html/math_toolkit/atan.html
	doc/html/math_toolkit/atanh.html
	doc/html/math_toolkit/bessel.html
	doc/html/math_toolkit/bessel/bessel_derivatives.html
	doc/html/math_toolkit/bessel/bessel_first.html
	doc/html/math_toolkit/bessel/bessel_over.html
	doc/html/math_toolkit/bessel/bessel_root.html
	doc/html/math_toolkit/bessel/mbessel.html
	doc/html/math_toolkit/bessel/sph_bessel.html
	doc/html/math_toolkit/building.html
	doc/html/math_toolkit/c99.html
	doc/html/math_toolkit/comp_compilers.html
	doc/html/math_toolkit/comparisons.html
	doc/html/math_toolkit/compile_time.html
	doc/html/math_toolkit/compilers_overview.html
	doc/html/math_toolkit/complex_history.html
	doc/html/math_toolkit/complex_implementation.html
	doc/html/math_toolkit/concepts.html
	doc/html/math_toolkit/config_macros.html
	doc/html/math_toolkit/constants.html
	doc/html/math_toolkit/constants_faq.html
	doc/html/math_toolkit/constants_intro.html
	doc/html/math_toolkit/contact.html
	doc/html/math_toolkit/conventions.html
	doc/html/math_toolkit/create.html
	doc/html/math_toolkit/credits.html
	doc/html/math_toolkit/demo.html
	doc/html/math_toolkit/directories.html
	doc/html/math_toolkit/dist_concept.html
	doc/html/math_toolkit/dist_ref.html
	doc/html/math_toolkit/dist_ref/dist_algorithms.html
	doc/html/math_toolkit/dist_ref/dists.html
	doc/html/math_toolkit/dist_ref/dists/bernoulli_dist.html
	doc/html/math_toolkit/dist_ref/dists/beta_dist.html
	doc/html/math_toolkit/dist_ref/dists/binomial_dist.html
	doc/html/math_toolkit/dist_ref/dists/cauchy_dist.html
	doc/html/math_toolkit/dist_ref/dists/chi_squared_dist.html
	doc/html/math_toolkit/dist_ref/dists/exp_dist.html
	doc/html/math_toolkit/dist_ref/dists/extreme_dist.html
	doc/html/math_toolkit/dist_ref/dists/f_dist.html
	doc/html/math_toolkit/dist_ref/dists/gamma_dist.html
	doc/html/math_toolkit/dist_ref/dists/geometric_dist.html
	doc/html/math_toolkit/dist_ref/dists/hyperexponential_dist.html
	doc/html/math_toolkit/dist_ref/dists/hypergeometric_dist.html
	doc/html/math_toolkit/dist_ref/dists/inverse_chi_squared_dist.html
	doc/html/math_toolkit/dist_ref/dists/inverse_gamma_dist.html
	doc/html/math_toolkit/dist_ref/dists/inverse_gaussian_dist.html
	doc/html/math_toolkit/dist_ref/dists/laplace_dist.html
	doc/html/math_toolkit/dist_ref/dists/logistic_dist.html
	doc/html/math_toolkit/dist_ref/dists/lognormal_dist.html
	doc/html/math_toolkit/dist_ref/dists/nc_beta_dist.html
	doc/html/math_toolkit/dist_ref/dists/nc_chi_squared_dist.html
	doc/html/math_toolkit/dist_ref/dists/nc_f_dist.html
	doc/html/math_toolkit/dist_ref/dists/nc_t_dist.html
	doc/html/math_toolkit/dist_ref/dists/negative_binomial_dist.html
	doc/html/math_toolkit/dist_ref/dists/normal_dist.html
	doc/html/math_toolkit/dist_ref/dists/pareto.html
	doc/html/math_toolkit/dist_ref/dists/poisson_dist.html
	doc/html/math_toolkit/dist_ref/dists/rayleigh.html
	doc/html/math_toolkit/dist_ref/dists/skew_normal_dist.html
	doc/html/math_toolkit/dist_ref/dists/students_t_dist.html
	doc/html/math_toolkit/dist_ref/dists/triangular_dist.html
	doc/html/math_toolkit/dist_ref/dists/uniform_dist.html
	doc/html/math_toolkit/dist_ref/dists/weibull_dist.html
	doc/html/math_toolkit/dist_ref/nmp.html
	doc/html/math_toolkit/ellint.html
	doc/html/math_toolkit/ellint/ellint_1.html
	doc/html/math_toolkit/ellint/ellint_2.html
	doc/html/math_toolkit/ellint/ellint_3.html
	doc/html/math_toolkit/ellint/ellint_carlson.html
	doc/html/math_toolkit/ellint/ellint_intro.html
	doc/html/math_toolkit/error_handling.html
	doc/html/math_toolkit/exact_typdefs.html
	doc/html/math_toolkit/examples.html
	doc/html/math_toolkit/exp.html
	doc/html/math_toolkit/expint.html
	doc/html/math_toolkit/expint/expint_i.html
	doc/html/math_toolkit/expint/expint_n.html
	doc/html/math_toolkit/factorials.html
	doc/html/math_toolkit/factorials/sf_binomial.html
	doc/html/math_toolkit/factorials/sf_double_factorial.html
	doc/html/math_toolkit/factorials/sf_factorial.html
	doc/html/math_toolkit/factorials/sf_falling_factorial.html
	doc/html/math_toolkit/factorials/sf_rising_factorial.html
	doc/html/math_toolkit/fastest_typdefs.html
	doc/html/math_toolkit/float128.html
	doc/html/math_toolkit/float128/exp_function.html
	doc/html/math_toolkit/float128/overloading_template_functions_w.html
	doc/html/math_toolkit/float128/typeinfo.html
	doc/html/math_toolkit/fp_facets.html
	doc/html/math_toolkit/fp_facets/examples.html
	doc/html/math_toolkit/fp_facets/facets_intro.html
	doc/html/math_toolkit/fp_facets/portability.html
	doc/html/math_toolkit/fp_facets/rationale.html
	doc/html/math_toolkit/fp_facets/reference.html
	doc/html/math_toolkit/fpclass.html
	doc/html/math_toolkit/future.html
	doc/html/math_toolkit/gcd_credits.html
	doc/html/math_toolkit/gcd_function_object.html
	doc/html/math_toolkit/gcd_header.html
	doc/html/math_toolkit/gcd_history.html
	doc/html/math_toolkit/getting_best.html
	doc/html/math_toolkit/greatest_typdefs.html
	doc/html/math_toolkit/hankel.html
	doc/html/math_toolkit/hankel/cyl_hankel.html
	doc/html/math_toolkit/hankel/sph_hankel.html
	doc/html/math_toolkit/high_precision.html
	doc/html/math_toolkit/high_precision/e_float.html
	doc/html/math_toolkit/high_precision/float128.html
	doc/html/math_toolkit/high_precision/use_mpfr.html
	doc/html/math_toolkit/high_precision/use_multiprecision.html
	doc/html/math_toolkit/high_precision/use_ntl.html
	doc/html/math_toolkit/high_precision/using_test.html
	doc/html/math_toolkit/high_precision/why_high_precision.html
	doc/html/math_toolkit/hints.html
	doc/html/math_toolkit/history1.html
	doc/html/math_toolkit/history2.html
	doc/html/math_toolkit/internals1.html
	doc/html/math_toolkit/internals1/cf.html
	doc/html/math_toolkit/internals1/minima.html
	doc/html/math_toolkit/internals1/rational.html
	doc/html/math_toolkit/internals1/roots.html
	doc/html/math_toolkit/internals1/roots2.html
	doc/html/math_toolkit/internals1/series_evaluation.html
	doc/html/math_toolkit/internals1/tuples.html
	doc/html/math_toolkit/internals2.html
	doc/html/math_toolkit/internals2/error_test.html
	doc/html/math_toolkit/internals2/minimax.html
	doc/html/math_toolkit/internals2/polynomials.html
	doc/html/math_toolkit/internals2/test_data.html
	doc/html/math_toolkit/internals_overview.html
	doc/html/math_toolkit/interp.html
	doc/html/math_toolkit/intro_pol_overview.html
	doc/html/math_toolkit/introduction.html
	doc/html/math_toolkit/inv_hyper.html
	doc/html/math_toolkit/inv_hyper/acosh.html
	doc/html/math_toolkit/inv_hyper/asinh.html
	doc/html/math_toolkit/inv_hyper/atanh.html
	doc/html/math_toolkit/inv_hyper/inv_hyper_over.html
	doc/html/math_toolkit/issues.html
	doc/html/math_toolkit/jacobi.html
	doc/html/math_toolkit/jacobi/jac_over.html
	doc/html/math_toolkit/jacobi/jacobi_cd.html
	doc/html/math_toolkit/jacobi/jacobi_cn.html
	doc/html/math_toolkit/jacobi/jacobi_cs.html
	doc/html/math_toolkit/jacobi/jacobi_dc.html
	doc/html/math_toolkit/jacobi/jacobi_dn.html
	doc/html/math_toolkit/jacobi/jacobi_ds.html
	doc/html/math_toolkit/jacobi/jacobi_elliptic.html
	doc/html/math_toolkit/jacobi/jacobi_nc.html
	doc/html/math_toolkit/jacobi/jacobi_nd.html
	doc/html/math_toolkit/jacobi/jacobi_ns.html
	doc/html/math_toolkit/jacobi/jacobi_sc.html
	doc/html/math_toolkit/jacobi/jacobi_sd.html
	doc/html/math_toolkit/jacobi/jacobi_sn.html
	doc/html/math_toolkit/lanczos.html
	doc/html/math_toolkit/lcm_function_object.html
	doc/html/math_toolkit/macros.html
	doc/html/math_toolkit/main_faq.html
	doc/html/math_toolkit/main_intro.html
	doc/html/math_toolkit/main_tr1.html
	doc/html/math_toolkit/mem_typedef.html
	doc/html/math_toolkit/minimum_typdefs.html
	doc/html/math_toolkit/namespaces.html
	doc/html/math_toolkit/navigation.html
	doc/html/math_toolkit/new_const.html
	doc/html/math_toolkit/next_float.html
	doc/html/math_toolkit/next_float/float_advance.html
	doc/html/math_toolkit/next_float/float_distance.html
	doc/html/math_toolkit/next_float/float_next.html
	doc/html/math_toolkit/next_float/float_prior.html
	doc/html/math_toolkit/next_float/nextafter.html
	doc/html/math_toolkit/number_series.html
	doc/html/math_toolkit/number_series/bernoulli_numbers.html
	doc/html/math_toolkit/number_series/primes.html
	doc/html/math_toolkit/number_series/tangent_numbers.html
	doc/html/math_toolkit/oct_create.html
	doc/html/math_toolkit/oct_header.html
	doc/html/math_toolkit/oct_history.html
	doc/html/math_toolkit/oct_mem_fun.html
	doc/html/math_toolkit/oct_non_mem.html
	doc/html/math_toolkit/oct_overview.html
	doc/html/math_toolkit/oct_specialization.html
	doc/html/math_toolkit/oct_synopsis.html
	doc/html/math_toolkit/oct_tests.html
	doc/html/math_toolkit/oct_todo.html
	doc/html/math_toolkit/oct_trans.html
	doc/html/math_toolkit/oct_typedefs.html
	doc/html/math_toolkit/oct_value_ops.html
	doc/html/math_toolkit/octonion.html
	doc/html/math_toolkit/overview.html
	doc/html/math_toolkit/overview_tr1.html
	doc/html/math_toolkit/owens_t.html
	doc/html/math_toolkit/perf_over1.html
	doc/html/math_toolkit/perf_over2.html
	doc/html/math_toolkit/perf_test_app.html
	doc/html/math_toolkit/pol_overview.html
	doc/html/math_toolkit/pol_ref.html
	doc/html/math_toolkit/pol_ref/assert_undefined.html
	doc/html/math_toolkit/pol_ref/discrete_quant_ref.html
	doc/html/math_toolkit/pol_ref/error_handling_policies.html
	doc/html/math_toolkit/pol_ref/internal_promotion.html
	doc/html/math_toolkit/pol_ref/iteration_pol.html
	doc/html/math_toolkit/pol_ref/namespace_pol.html
	doc/html/math_toolkit/pol_ref/pol_ref_ref.html
	doc/html/math_toolkit/pol_ref/policy_defaults.html
	doc/html/math_toolkit/pol_ref/precision_pol.html
	doc/html/math_toolkit/pol_tutorial.html
	doc/html/math_toolkit/pol_tutorial/ad_hoc_dist_policies.html
	doc/html/math_toolkit/pol_tutorial/ad_hoc_sf_policies.html
	doc/html/math_toolkit/pol_tutorial/changing_policy_defaults.html
	doc/html/math_toolkit/pol_tutorial/namespace_policies.html
	doc/html/math_toolkit/pol_tutorial/policy_tut_defaults.html
	doc/html/math_toolkit/pol_tutorial/policy_usage.html
	doc/html/math_toolkit/pol_tutorial/understand_dis_quant.html
	doc/html/math_toolkit/pol_tutorial/user_def_err_pol.html
	doc/html/math_toolkit/pol_tutorial/what_is_a_policy.html
	doc/html/math_toolkit/powers.html
	doc/html/math_toolkit/powers/cbrt.html
	doc/html/math_toolkit/powers/cos_pi.html
	doc/html/math_toolkit/powers/ct_pow.html
	doc/html/math_toolkit/powers/expm1.html
	doc/html/math_toolkit/powers/hypot.html
	doc/html/math_toolkit/powers/log1p.html
	doc/html/math_toolkit/powers/powm1.html
	doc/html/math_toolkit/powers/sin_pi.html
	doc/html/math_toolkit/powers/sqrt1pm1.html
	doc/html/math_toolkit/quat.html
	doc/html/math_toolkit/quat_header.html
	doc/html/math_toolkit/quat_history.html
	doc/html/math_toolkit/quat_mem_fun.html
	doc/html/math_toolkit/quat_non_mem.html
	doc/html/math_toolkit/quat_overview.html
	doc/html/math_toolkit/quat_synopsis.html
	doc/html/math_toolkit/quat_tests.html
	doc/html/math_toolkit/quat_todo.html
	doc/html/math_toolkit/rationale.html
	doc/html/math_toolkit/rationale0.html
	doc/html/math_toolkit/refs.html
	doc/html/math_toolkit/relative_error.html
	doc/html/math_toolkit/remez.html
	doc/html/math_toolkit/result_type.html
	doc/html/math_toolkit/rounding.html
	doc/html/math_toolkit/rounding/modf.html
	doc/html/math_toolkit/rounding/round.html
	doc/html/math_toolkit/rounding/trunc.html
	doc/html/math_toolkit/run_time.html
	doc/html/math_toolkit/sf_beta.html
	doc/html/math_toolkit/sf_beta/beta_derivative.html
	doc/html/math_toolkit/sf_beta/beta_function.html
	doc/html/math_toolkit/sf_beta/ibeta_function.html
	doc/html/math_toolkit/sf_beta/ibeta_inv_function.html
	doc/html/math_toolkit/sf_erf.html
	doc/html/math_toolkit/sf_erf/error_function.html
	doc/html/math_toolkit/sf_erf/error_inv.html
	doc/html/math_toolkit/sf_gamma.html
	doc/html/math_toolkit/sf_gamma/digamma.html
	doc/html/math_toolkit/sf_gamma/gamma_derivatives.html
	doc/html/math_toolkit/sf_gamma/gamma_ratios.html
	doc/html/math_toolkit/sf_gamma/igamma.html
	doc/html/math_toolkit/sf_gamma/igamma_inv.html
	doc/html/math_toolkit/sf_gamma/lgamma.html
	doc/html/math_toolkit/sf_gamma/tgamma.html
	doc/html/math_toolkit/sf_implementation.html
	doc/html/math_toolkit/sf_poly.html
	doc/html/math_toolkit/sf_poly/hermite.html
	doc/html/math_toolkit/sf_poly/laguerre.html
	doc/html/math_toolkit/sf_poly/legendre.html
	doc/html/math_toolkit/sf_poly/sph_harm.html
	doc/html/math_toolkit/sign_functions.html
	doc/html/math_toolkit/sinc.html
	doc/html/math_toolkit/sinc/sinc_overview.html
	doc/html/math_toolkit/sinc/sinc_pi.html
	doc/html/math_toolkit/sinc/sinhc_pi.html
	doc/html/math_toolkit/spec.html
	doc/html/math_toolkit/special_tut.html
	doc/html/math_toolkit/special_tut/special_tut_impl.html
	doc/html/math_toolkit/special_tut/special_tut_test.html
	doc/html/math_toolkit/stat_tut.html
	doc/html/math_toolkit/stat_tut/dist_params.html
	doc/html/math_toolkit/stat_tut/overview.html
	doc/html/math_toolkit/stat_tut/overview/complements.html
	doc/html/math_toolkit/stat_tut/overview/generic.html
	doc/html/math_toolkit/stat_tut/overview/headers.html
	doc/html/math_toolkit/stat_tut/overview/objects.html
	doc/html/math_toolkit/stat_tut/overview/parameters.html
	doc/html/math_toolkit/stat_tut/overview/summary.html
	doc/html/math_toolkit/stat_tut/variates.html
	doc/html/math_toolkit/stat_tut/weg.html
	doc/html/math_toolkit/stat_tut/weg/binom_eg.html
	doc/html/math_toolkit/stat_tut/weg/binom_eg/binom_conf.html
	doc/html/math_toolkit/stat_tut/weg/binom_eg/binom_size_eg.html
	doc/html/math_toolkit/stat_tut/weg/binom_eg/binomial_coinflip_example.html
	doc/html/math_toolkit/stat_tut/weg/binom_eg/binomial_quiz_example.html
	doc/html/math_toolkit/stat_tut/weg/c_sharp.html
	doc/html/math_toolkit/stat_tut/weg/cs_eg.html
	doc/html/math_toolkit/stat_tut/weg/cs_eg/chi_sq_intervals.html
	doc/html/math_toolkit/stat_tut/weg/cs_eg/chi_sq_size.html
	doc/html/math_toolkit/stat_tut/weg/cs_eg/chi_sq_test.html
	doc/html/math_toolkit/stat_tut/weg/dist_construct_eg.html
	doc/html/math_toolkit/stat_tut/weg/error_eg.html
	doc/html/math_toolkit/stat_tut/weg/f_eg.html
	doc/html/math_toolkit/stat_tut/weg/find_eg.html
	doc/html/math_toolkit/stat_tut/weg/find_eg/find_location_eg.html
	doc/html/math_toolkit/stat_tut/weg/find_eg/find_mean_and_sd_eg.html
	doc/html/math_toolkit/stat_tut/weg/find_eg/find_scale_eg.html
	doc/html/math_toolkit/stat_tut/weg/geometric_eg.html
	doc/html/math_toolkit/stat_tut/weg/inverse_chi_squared_eg.html
	doc/html/math_toolkit/stat_tut/weg/nag_library.html
	doc/html/math_toolkit/stat_tut/weg/nccs_eg.html
	doc/html/math_toolkit/stat_tut/weg/nccs_eg/nccs_power_eg.html
	doc/html/math_toolkit/stat_tut/weg/neg_binom_eg.html
	doc/html/math_toolkit/stat_tut/weg/neg_binom_eg/neg_binom_conf.html
	doc/html/math_toolkit/stat_tut/weg/neg_binom_eg/neg_binom_size_eg.html
	doc/html/math_toolkit/stat_tut/weg/neg_binom_eg/negative_binomial_example1.html
	doc/html/math_toolkit/stat_tut/weg/neg_binom_eg/negative_binomial_example2.html
	doc/html/math_toolkit/stat_tut/weg/normal_example.html
	doc/html/math_toolkit/stat_tut/weg/normal_example/normal_misc.html
	doc/html/math_toolkit/stat_tut/weg/st_eg.html
	doc/html/math_toolkit/stat_tut/weg/st_eg/paired_st.html
	doc/html/math_toolkit/stat_tut/weg/st_eg/tut_mean_intervals.html
	doc/html/math_toolkit/stat_tut/weg/st_eg/tut_mean_size.html
	doc/html/math_toolkit/stat_tut/weg/st_eg/tut_mean_test.html
	doc/html/math_toolkit/stat_tut/weg/st_eg/two_sample_students_t.html
	doc/html/math_toolkit/synopsis.html
	doc/html/math_toolkit/threads.html
	doc/html/math_toolkit/tr1_ref.html
	doc/html/math_toolkit/trans.html
	doc/html/math_toolkit/tuning.html
	doc/html/math_toolkit/tutorial.html
	doc/html/math_toolkit/tutorial/non_templ.html
	doc/html/math_toolkit/tutorial/templ.html
	doc/html/math_toolkit/tutorial/user_def.html
	doc/html/math_toolkit/value_op.html
	doc/html/math_toolkit/zetas.html
	doc/html/math_toolkit/zetas/zeta.html
	doc/html/octonions.html
	doc/html/overview.html
	doc/html/perf.html
	doc/html/policy.html
	doc/html/quaternions.html
	doc/html/special.html
	doc/html/status.html
	doc/html/toolkit.html
	doc/html/using_udt.html
	doc/html/utils.html
 2014-12-08 12:45:58 +00:00
jzmaddock
5b9d72131f [polygamma] Fix integer arithmetic bug. 2014-12-01 13:16:41 +00:00
pabristow
f9b2addcdc Improved docs and tests and graphs. 2014-11-27 16:44:42 +00:00
jzmaddock
f064f2e36b [signbit] Fix call to work when signbit is a macro. 2014-11-26 19:15:24 +00:00
jzmaddock
b4111f8312 [polygamma] Fix signed/unsigned comparison warning 2014-11-26 19:06:10 +00:00
jzmaddock
85b59d348b [Polygamma] Fix a few corner cases, and bail out early in reflection formula if n is too large. 2014-11-26 18:22:30 +00:00
jzmaddock
e21bca611f [polygamma] Remove series acceleration from polygamma_nearzero as it appears to make no real difference to the number of terms required. 2014-11-26 17:42:30 +00:00
jzmaddock
4194fea6d0 [polygamma] Fix integer overflow in branch check case, and bail out of transition code is the number of series evaluations is too high. 2014-11-24 18:08:31 +00:00
jzmaddock
a663232f38 [polygamma] Add two special cases for x = 1 and x = 0.5. 2014-11-22 18:37:31 +00:00
jzmaddock
e2cd2e72dc [polygamma] Document new method for negative x in code comments, simply some code, change table to coefficients to store only non-zero values. 2014-11-20 09:56:21 +00:00
jzmaddock
075f9f893b [polygamma] Change low order negative-x cases to use new polynomial method. Change expected error rates to match. 2014-11-18 18:27:21 +00:00
jzmaddock
495d4655b4 [polygamma] remove dead code. 2014-11-17 18:23:30 +00:00
jzmaddock
f4b86c2996 [polygamma] Switch to new better method for handling negative arguments. 2014-11-17 18:12:25 +00:00
jzmaddock
26f063eca4 Fix accuracy issue in cos_pi. 
Add tests for cos_pi and sin_pi.
Doc regen.
 2014-11-16 17:18:17 +00:00
jzmaddock
9d6df5d835 [polygamma] Remove some dead code and tidy up initialization and policy usage. 2014-11-14 17:05:41 +00:00
jzmaddock
2bcdfed1a4 [Polygamma] Remove dead code. 2014-11-08 11:03:45 +00:00
jzmaddock
dd92123098 [digamma] fix reference. 2014-11-07 16:57:18 +00:00
jzmaddock
50778a40c4 [digamma]Oops fix missing return statement on half integer special case. 2014-11-07 16:56:15 +00:00
jzmaddock
0e98fd90c2 [Polygamma] Add optimized trigamma version that uses rational approximations. 2014-11-07 16:32:41 +00:00
jzmaddock
8477de53b4 [polygamma]Fix spurious overflow for very large x. 2014-11-04 18:27:48 +00:00
jzmaddock
3478ee9c46 [polygamma] Fix GCC failures, a few bugs and forward policy arguments. Fix test data to be 128-bit safe. 2014-11-03 17:55:21 +00:00
pabristow
04c43aba75 Cosmetic tidy. 2014-11-03 12:11:50 +00:00
jzmaddock
d50cf889f6 Add basic thread safety and initialization code. 2014-11-03 12:01:11 +00:00
jzmaddock
eccec791c7 [polygamma] Fix spurious overflow in transition zone, update cot derivative to work to arbitrary level. 2014-11-03 11:50:10 +00:00
jzmaddock
e70f53e670 [Polygamma] Fix issues with small negative arguments for x. 
Tidy up a couple of cosmetic coding issues.
 2014-10-31 18:51:14 +00:00
jzmaddock
8438b8a84b [Polygamma] Extend derivatives-of-cot method to larger orders. 
Extend test cases for negative x.
Fix missing return statement on error condition.
 2014-10-31 18:07:01 +00:00
jzmaddock
31c345108c [polygamma] Fix gcc-11 test failures. 2014-10-30 17:13:39 +00:00
pabristow
397a3cf4ff Merge branch 'arcsine_dist' of https://github.com/boostorg/math into arcsine_dist 2014-10-30 12:15:13 +00:00
pabristow
2de39afe7e Fixed silly mistake causing real_concept tests to fail 2014-10-30 12:14:33 +00:00
jzmaddock
6bf9d2bffc [arcsine]Fix concept failures 2014-10-30 11:46:05 +00:00
jzmaddock
812122d56f Fix include order to be alphabetical, update copyright. 2014-10-30 11:13:31 +00:00
jzmaddock
3f58ab9009 Fix copyright info. 2014-10-30 11:11:50 +00:00
jzmaddock
081abc704a [zeta/polygamma]Add cache of odd-integer zeta values. 
Fix even integer case to call fast closed-form code more often.
Fix initialization code to initializes new cache of even-integer values.
 2014-10-30 11:08:36 +00:00
jzmaddock
8a36ba07d1 {Polygamma]Fix polygamma_nearzero to call correct zeta function overload. 2014-10-30 11:06:28 +00:00
jzmaddock
232e7af69b [polygamma] short circuit series evaluation when it can make no difference to the result. 2014-10-29 17:38:38 +00:00
jzmaddock
8362ab2a08 [polygamma] 
Replace small-x expansion with one that's accelerated and avoids spurious overflow.
Also tweak selection logic, and reuse variables better.
 2014-10-28 19:39:35 +00:00
pabristow
477e194490 Corrected previous failures using real_concept. tested locally OK on MSVC and gcc 4.8 2014-10-28 16:37:43 +00:00
pabristow
0a932fa52b Expanded tests from double, but errors - see email 27 Oct14 2014-10-27 15:51:54 +00:00
jzmaddock
a3bed60089 More test failure fixes. 2014-10-26 16:25:19 +00:00
jzmaddock
22bc7212e7 [Polygamma] 
Hook up concept checks.
Move test instantiations into instances lib.
Fix a few test failures.
Add Nikhar's trigamma tests.
 2014-10-26 09:42:11 +00:00
jzmaddock
dffdb444ee Expand test cases for negative x. 
Add comment on the calculation of derivatives of cot(x).
 2014-10-25 17:28:09 +01:00
jzmaddock
bda3331ed3 [Polygamma] 
Fix digits_base10 to actually return the correct value!
Add support for negative x.
Add missing using declaration to digamma.
Update tests with negative x tests.
 2014-10-24 18:50:37 +01:00
jzmaddock
f2ae6940ca [Polygamma] 
Optimize zeta function for integer arguments.
 2014-10-23 17:08:19 +01:00
jzmaddock
29e86c8fcb [Polygamma] 
Extend existing digamma implementation to true multiprecision.
Add integer and half-integer test cases for digamma.
Change polygamma to call digamma when order is zero.
 2014-10-23 11:42:57 +01:00
jzmaddock
d9ea302cff Hook up basic multiprecision digamma. 2014-10-23 09:59:08 +01:00
jzmaddock
4bc3b6076c [Polygamma] 
Fix real_concept compilation and runtime.
Add digits_base10 support function to policies.
 2014-10-21 19:03:26 +01:00
jzmaddock
5f89e70efd [Polygamma] 
Add test cases.
Rewrite polygamma_atinfinityplus to avoid spurious underflow/overflow.
Rewrite polygamma_attransitionplus to call polygamma_atinfinityplus rather than have it's own routine.
Change condition which selects when polygamma_atinfinityplus can be called.
 2014-10-21 13:12:26 +01:00
pabristow
716e7eff75 First draft of working docs. Need more tests. 2014-10-20 14:36:03 +01:00
jzmaddock
eb2525b516 Merge branch 'develop' into bernoulli-tgamma 
Conflicts:
	doc/equations/bernoulli_numbers2.mml
	doc/equations/bernoulli_numbers2.png
	doc/equations/bernoulli_numbers2.svg
	doc/equations/generate.sh
	doc/equations/tangent_numbers.png
	doc/equations/tangent_numbers.svg
	doc/math.qbk
	include/boost/cstdfloat.hpp
	include/boost/math/policies/error_handling.hpp
	include/boost/math/special_functions/detail/bernoulli_details.hpp
 2014-10-18 17:49:49 +01:00
pabristow
7d303aaab0 Initial rough drafts of arcsine (only committed so can pick up better diagnostics in error_handling from inf_dist branch) 2014-10-03 10:12:05 +01:00
jzmaddock
bd807c74e1 Fix more spurious overflows detected by Rocco Romeo. 2014-09-24 18:02:00 +01:00
jzmaddock
b471723bd1 Fix bug cases reported by Rocco Romeo: 
* Add another approximation for x very small.
* Add code to handle cases where tgamma(a) would overflow.
* Add tests for the above.
 2014-09-14 17:05:15 +01:00
jzmaddock
d1d6d5b063 Fix derivative of incomplete beta in special cases where a = b = 0.5 and b = 1. 
Add additional tests to non-central F and T distributions to verify above fix.
Fixes https://svn.boost.org/trac/boost/ticket/10480.
 2014-09-12 16:29:29 +01:00
jzmaddock
321d8f2eda Merge branch 'hyper_exponential' into develop 2014-09-08 18:12:52 +01:00
jzmaddock
17e4f40fce Revert "[distributions] Hyper-Exponential: changed return type of 'probabilities' and 'rates' from 'vector<...>' to 'vector<...> const\&' to allow users to do things like this: 'std::copy(he.probabilities().begin(), he.probabilities().end(), std::ostream_iterator<double>(std::cout, " "));'." 
This reverts commit a8efc9a637.
 2014-09-08 18:10:53 +01:00
jzmaddock
5d37faf873 Merge branch 'feature/hyperexponential_dist' of https://github.com/sguazt/math into hyper_exponential 2014-09-08 17:14:23 +01:00
sguazt
a8efc9a637 [distributions] Hyper-Exponential: changed return type of 'probabilities' and 'rates' from 'vector<...>' to 'vector<...> const\&' to allow users to do things like this: 'std::copy(he.probabilities().begin(), he.probabilities().end(), std::ostream_iterator<double>(std::cout, " "));'. 2014-09-07 10:29:09 +02:00
jzmaddock
f73ecae7f5 Add some tests for the error handlers. 2014-09-05 18:56:49 +01:00
jzmaddock
bc77596f23 Improve calculation of phase probabilities when normalizing so they really should add up to 1. 
Allow slight tolerance in checking for sum of phase probabilities.
Add missing call to normalize in iterator constructor.
Remove asserts that duplicated the error handling code - better to provide consistent behavior via the error handlers.
Add comment to examples about other containers working equally well.
 2014-09-03 19:25:20 +01:00
sguazt
43de825473 [distributions] Hyper-Exponential: BOOST_MSVC -> _MSC_VER (as _MSC_VER is always defined for MS VS, whereas BOOST_MSVC isn’t defined until *after* the first call to boost/config.hpp). Also, added MSVC pragma 4389. 2014-09-01 21:25:11 +02:00
jzmaddock
4a6cbc1135 Initializer-list constructors don't need to be templates - indeed it is better if they're not. 2014-09-01 18:14:35 +01:00
sguazt
f5cec49a60 [distributions] Hyper-Exponential: added include of 'iterator' for 'std::distance' and aesthetic changes. 2014-09-01 09:50:16 +02:00
sguazt
34e42a36c5 [distributions] Hyper-Exponential: added include file for macro BOOST_NO_CXX11_HDR_INITIALIZER_LIST and removed 'const\&' from 'std::initializer_list' (as from (Stroustrup,TCPL4e). 2014-08-31 23:22:15 +02:00
jzmaddock
60a6730fe2 Modify Marco's 1 and 2 arg constructors to use SFINAE to differentiate between iterators and ranges. 
Add a 1 arg initializer list constructor.
Update tests accordingly.
 2014-08-30 17:03:26 +01:00
sguazt
62ee849eab [distributions] Hyper-Exponential: added more constructors and class construction test case. 2014-08-28 10:12:41 +02:00
jzmaddock
08308ab945 Merge branch 'feature/hyperexponential_dist' of https://github.com/sguazt/math into hyper_exponential 2014-08-26 18:57:25 +01:00
sguazt
fa198a0b04 [distributions] Hyper-Exponential: clean-up. 2014-08-26 17:08:05 +02:00
jzmaddock
2e7eede435 Add comment linking to discussion. 2014-08-25 17:04:52 +01:00
jzmaddock
d8c97bec27 Add initializer_list support and tests. 2014-08-25 13:11:51 +01:00
sguazt
c693510588 [distributions] Hyper-Exponential: use exponential's PDF in hyperexponential PDF. 2014-08-25 10:45:08 +02:00
jzmaddock
d424e2ef23 Fix some more warnings and enable tests for all the test types. 2014-08-24 18:40:21 +01:00
jzmaddock
b925b0c0e8 Add hyperexponential_distribution to concept checks, fix for warnings/bugs. 2014-08-24 15:21:36 +01:00
jzmaddock
ffc483f907 Up precision of tests. 
Fix some conceptual compiler errors in hyperexponential_distribution.
 2014-08-24 09:29:11 +01:00
Robin Eckert
fb835eb029 fix rounding errors 
- the direct predecessor/successor for 0.5/-0.5 was formerly wrongly
   rounded to 1.0/-1.0 rather than 0.0

 - integral numbers of the format 1xxxxxx1 filling mantissa completely
   were formerly rounded to the next integral away from zero
 2014-08-23 19:42:02 +02:00
sguazt
0735d8177e [distributions] Hyper-exponential: range arguments to contructor should be passed by const-ref. 2014-08-22 18:32:22 +02:00
sguazt
3621b083e2 [distributions] Hyper-Exponential: registered tests and fixed some test-related issue. 2014-08-22 16:37:35 +02:00
sguazt
144db22c2f [distributions] Hyper-Exponential: solved 'unused parameter' warning in 'mode' function. 2014-08-22 15:06:12 +02:00
sguazt
412bba6a24 [distributions] Added Hyper-Exponential distribution. 2014-08-22 15:01:23 +02:00
Kohei Takahashi
3c07a49aac Don't check long double capability for BOOST_CSTDFLOAT_xx. 
For example IA-32, those macros suggest long double as 80-bit floating
  point type, however some math runtime implementations are broken for
  long doublle even if supported by compiler (e.g. FreeBSD 9, ...).

Signed-off-by: Kohei Takahashi <flast@flast.jp>
 2014-08-18 05:13:00 +09:00
jzmaddock
97cdf3ca3a Disable some msvc warnings, see: d02ae84866 (commitcomment-7257997) 2014-08-04 19:11:53 +01:00
jzmaddock
6e8145a076 Add code to profile changes to toms748_solve.hpp to see if their effect is beneficial. 
Change the heuristics used when the initial guess is really bad.
Change example to return a fail-code when an error occurs.
 2014-07-31 12:48:23 +01:00
Marcel Raad
7621f1e931 Remove redundant variable declarations 
... to silence compiler warnings about variable shadowing.
 2014-07-12 12:44:45 +02:00
Rogier van Dalen
d9ba6ae8b0 Fix missing include 2014-07-10 15:45:49 +01:00
jzmaddock
866492cb02 remove dependency to tr1 - we don't really need it anymore. 2014-06-13 13:32:08 +01:00
jzmaddock
2f65e6cb65 Fix typos. 2014-06-13 08:20:14 +01:00
jzmaddock
ee8edd4c33 Fix up functions that were returning tools::max_value to call overflow_error instead. 
See https://svn.boost.org/trac/boost/ticket/10111.
 2014-06-12 13:50:30 +01:00
jzmaddock
e20b44e8a9 Move headers that are mostly used for internal maintenance and testing into an "include_private" directory - this is to reduce dependencies on other Boost modules when using this library. 2014-06-09 16:08:35 +01:00
jzmaddock
09fceb2142 Merge branch 'develop' of https://github.com/boostorg/math into develop 2014-05-28 13:25:10 +01:00
jzmaddock
7cb3316d06 Fix for 106-bit "double double" types. 2014-05-28 13:24:40 +01:00
Kohei Takahashi
c50b5a8bcc Fix missing header 
Signed-off-by: Kohei Takahashi <flast@flast.jp>
 2014-05-27 11:07:23 +09:00
jzmaddock
a29a74bdaf Add needed include to config.hpp, fix runtime checks under valgrind in erf_inv initialization, see https://svn.boost.org/trac/boost/ticket/10005 2014-05-26 19:14:49 +01:00
jzmaddock
ab6cc524dd Fix Halley iteration to work when derivative is zero, even though the second derivative is not, see https://svn.boost.org/trac/boost/ticket/10046. 2014-05-25 15:35:08 +01:00
jzmaddock
e6996e1188 Fix round and trunc functions to work with integer arguments, see https://svn.boost.org/trac/boost/ticket/10066. 2014-05-24 09:47:21 +01:00
jzmaddock
5a9688e4ef Apply fix from https://svn.boost.org/trac/boost/ticket/9884 2014-04-14 15:49:43 +01:00
jzmaddock
0c01f682eb Fix some corner cases in negative_binomial_distribution, see https://svn.boost.org/trac/boost/ticket/9834 2014-04-06 11:59:10 +01:00
jzmaddock
88b1695b58 More tentative apple fixes. 2014-04-04 17:20:09 +01:00
jzmaddock
1c1ae42e6e Tentative fix for more clang/apple failures. 2014-04-04 17:15:19 +01:00
jzmaddock
9199b09304 Fix missing policy parameter in some SF forwarding calls. See https://svn.boost.org/trac/boost/ticket/9833. 2014-04-04 15:08:27 +01:00
jzmaddock
2d0d36d865 Provide a further customization point for constants initialization (to allow for explicit construction). 2014-04-02 16:47:36 +01:00
jzmaddock
41fb7e1c91 Tentative fix for clang/libc++ on Apple issue. 2014-03-31 18:47:13 +01:00
jzmaddock
a0ad3d8764 Add some static casts to avoid C++11 errors in initializer lists. 2014-03-28 18:15:59 +00:00
jzmaddock
fa1896fbc4 Make sure the forward declarations get included before the real definitions in all headers. 2014-03-27 09:17:33 +00:00
jzmaddock
ef123b3504 Fix GCC release mode failures in test_next.cpp by explicitly testing SSE2 settings in next.hpp. 2014-03-22 19:16:36 +00:00
jzmaddock
4087c44efd Fix random number usage in the post-TR1 world. 2014-03-21 12:11:49 +00:00
jzmaddock
200f384c14 Tentative fix for gcc-4.7 and earlier failures. 2014-03-21 11:36:57 +00:00
jzmaddock
6bddb5be85 Merge branch 'cstdfloat' into develop 
Conflicts:
	doc/overview/roadmap.qbk
	example/cstdfloat_example.cpp
	include/boost/cstdfloat.hpp
	test/test_gamma.cpp
	test/test_gamma.hpp
 2014-03-19 13:38:42 +00:00
jzmaddock
296173dbe2 Some mostly cosmetic changes to the cstdfloat docs. 2014-03-17 16:58:58 +00:00
jzmaddock
cae085c873 More fixes for factorials and gamma function ratios from Rocco Romeo. 2014-03-16 09:59:11 +00:00
jzmaddock
3e338ec355 Fix some corner cases in tgamma_delta_ratio, with thanks to Rocco Romeo. 2014-03-14 18:18:22 +00:00
jzmaddock
c3e530716a Merge branch 'cstdfloat' of https://github.com/boostorg/math into cstdfloat 2014-03-13 18:43:28 +00:00
jzmaddock
80894e9025 Fix errors in corner cases of rising/falling factorials and gamma ratios raised by Rocco Romeo. 2014-03-13 16:30:10 +00:00
ckormanyos
849853f337 In <cstdfloat>, correct (again) the handling of cancellation errors in 128-bit sinh(x). 2014-03-12 21:56:08 +01:00
ckormanyos
aae19f5f08 In <cstdfloat>, correct the handling of cancellation errors in 128-bit sinh(x). 2014-03-12 21:23:44 +01:00
ckormanyos
9876cecd8f In <cstdfloat>, properly handle cancellation errors in 128-bit sinh(x). 2014-03-12 20:26:13 +01:00
jzmaddock
c9d65cfa1f Add basic tests for __float128 functions. 
Fix one complex function.
Add factorials test.
 2014-03-10 18:25:08 +00:00
jzmaddock
9f8ffee4b7 Disable some native operators when --fast-math or the Intel compiler is in use as they simply don't work in that environment, see https://svn.boost.org/trac/boost/ticket/9748#comment:3 2014-03-10 09:40:51 +00:00
jzmaddock
b4e1fe2375 Fix some compile errors in instrumentation code and extend a little into the Bessel functions. 2014-03-10 09:39:38 +00:00
jzmaddock
52a2b9edc3 Make sure all the error handlers return a result and don't just rely on an exception being thrown. 2014-03-05 09:13:26 +00:00
ckormanyos
8fd67e5edf In <boost/cstdfloat.hpp>, improve clarity of a compiler switch and make better comments. 2014-03-01 23:45:15 +01:00
jzmaddock
4fef819a00 Fix debug build. 2014-02-28 19:07:51 +00:00
jzmaddock
d3780ea58a Change to use C++11 <random> when available. 2014-02-28 13:40:10 +00:00
ckormanyos
3eca6efc55 In <boost/cstdfloat.hpp> and its quadmath part, do style modifications and add some informative comments. Undefine the internal BOOST_NO_FLOAT128_T from public availability. 2014-02-27 21:03:44 +01:00
jzmaddock
049f707a67 Fix Linux GCC failures 2014-02-27 17:07:25 +00:00
jzmaddock
b4759e553b Enable broken <cmath> function overloads only when needed. 2014-02-26 19:15:23 +00:00
jzmaddock
8459a54c90 Add missing return statements. 2014-02-26 19:09:29 +00:00