[libmath-prime-util-perl] 17/25: Give the M adjustment a steeper slope

[Partha P. Mukherjee]

This is an automated email from the git hooks/post-receive script.

ppm-guest pushed a commit to annotated tag v0.35
in repository libmath-prime-util-perl.

commit 6d0b318a4309778fb0ec66e95bb4bf043f64fd49
Author: Dana Jacobsen <dana at acm.org>
Date:   Fri Dec 6 12:11:21 2013 -0800

    Give the M adjustment a steeper slope
---
 lmo.c | 22 +++++++++++-----------
 1 file changed, 11 insertions(+), 11 deletions(-)

diff --git a/lmo.c b/lmo.c
index 0a9707e..cd54252 100644
--- a/lmo.c
+++ b/lmo.c
@@ -35,22 +35,22 @@
  *
  *  |   n   |  Meissel |  Lehmer  |   LMOS   |    LMO    |
  *  +-------+----------+----------+----------+-----------+
- *  | 10^19 |          |          |          | 3707.43   |
- *  | 10^18 |          |          |          |  778.70   |
- *  | 10^17 |          |          |          |  172.12   |
- *  | 10^16 | 1410.2   | 1043.6   | 1058.9   |   37.351  |
- *  | 10^15 |  137.1   |  137.3   |  142.7   |    8.072  |
- *  | 10^14 |   26.18  |   21.74  |   21.29  |    1.906  |
- *  | 10^13 |    5.155 |    3.947 |    3.353 |    0.455  |
- *  | 10^12 |    1.059 |    0.700 |    0.626 |    0.111  |
- *  | 10^11 |    0.227 |    0.138 |    0.124 |    0.0269 |
- *  | 10^10 |    0.0509|    0.0309|    0.0286|    0.00702|
+ *  | 10^19 |          |          |          | 3765.02   |
+ *  | 10^18 |          |          |          |  778.21   |
+ *  | 10^17 |          |          | 8844.2   |  163.77   |
+ *  | 10^16 | 1410.2   | 1043.6   | 1058.9   |   34.81   |
+ *  | 10^15 |  137.1   |  137.3   |  142.7   |    7.905  |
+ *  | 10^14 |   26.18  |   21.74  |   21.29  |    1.726  |
+ *  | 10^13 |    5.155 |    3.947 |    3.353 |    0.405  |
+ *  | 10^12 |    1.059 |    0.700 |    0.626 |    0.0936 |
+ *  | 10^11 |    0.227 |    0.138 |    0.124 |    0.0227 |
+ *  | 10^10 |    0.0509|    0.0309|    0.0286|    0.00589|
  */
 
 /* Below this size, just sieve (with table speedup). */
 #define SIEVE_LIMIT  60000000
 /* Adjust to get best performance */
-#define M_FACTOR(n)     (UV) ((double)n * (log(n)/log(6)))
+#define M_FACTOR(n)     (UV) ((double)n * (log(n)/log(5.2)) * (log(log(n))-1.4))
 /* Size of segment used for previous primes, must be >= 21 */
 #define PREV_SIEVE_SIZE 512
 /* Phi sieve multiplier, adjust for best performance and memory use. */

-- 
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/pkg-perl/packages/libmath-prime-util-perl.git