mersenneforum.org  

Go Back   mersenneforum.org > Great Internet Mersenne Prime Search > Software

Reply
 
Thread Tools
Old 2021-02-09, 22:55   #540
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

2×3,217 Posts
Default

I have released 2.2.0:

Code:
   framework:
      Updated OpenCL on Windows.  See makefile for details.
      Updated primesieve to 7.6.
   
   psieve, psievecl:  version 1.4
      Some refactoring to support OpenCL worker.
      First release of psievecl.
      Verify factors from -I input file

   srsieve2, srsieve2cl:  version 1.5.1
      Fixed bug that was introduced in the refactoring of 1.5 that impacts generic sieving
      while using multiple threads.
      Added -R to remove sequences.  Use -Rk*b^n+c format to remove a single sequence or
      use -R with a file that has multiple sequences.  This is not tested yet.
psievecl is about 20x faster than psieve. The main slowdown is factor validation, which is less noticeable as factors become more sparse.

One odd behavior is that I noticed that srsieve2cl is slower than the previous release, but I do not know why. Even after reverting the framework changes for the release it was slower. I thought it was the update to OpenCL or primesieve, but reverting to the older versions of those made no difference. I'm likely missing something, but I don't know what.

I'm hoping that someone is willing to give the -R option with srsieve2 a spin.
rogue is online now   Reply With Quote
Old 2021-02-10, 02:51   #541
Happy5214
 
Happy5214's Avatar
 
"Alexander"
Nov 2008
The Alamo City

3·11·23 Posts
Default

Quote:
Originally Posted by rogue View Post
I'm hoping that someone is willing to give the -R option with srsieve2 a spin.
No luck:

Code:
$ ./srsieve2 -i b2_n.abcd -o b2_n.abcd -R "658687*2^n-1"
srsieve2 v1.5.1, a program to find factors of k*b^n+c numbers for fixed b and variable k and n
2018 terms for sequence 658687*2^n-1 have been removed
Must use generic sieving logic because there is more than one sequence
Sieving with generic logic for p >= 982453051
Fatal Error:  Expected 62636 terms when building sequences, but counted only 60618
Input file attached (with .txt extension added).
Attached Files
File Type: txt b2_n.abcd.txt (219.7 KB, 77 views)
Happy5214 is offline   Reply With Quote
Old 2021-02-10, 04:16   #542
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

11001001000102 Posts
Default

Quote:
Originally Posted by Happy5214 View Post
No luck:

Code:
$ ./srsieve2 -i b2_n.abcd -o b2_n.abcd -R "658687*2^n-1"
srsieve2 v1.5.1, a program to find factors of k*b^n+c numbers for fixed b and variable k and n
2018 terms for sequence 658687*2^n-1 have been removed
Must use generic sieving logic because there is more than one sequence
Sieving with generic logic for p >= 982453051
Fatal Error:  Expected 62636 terms when building sequences, but counted only 60618
Input file attached (with .txt extension added).
That will be easy to fix. When the sequence is removed the code isn't decrementing il_TermCount.
rogue is online now   Reply With Quote
Old 2021-02-23, 17:15   #543
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

643410 Posts
Default

I have released 2.2.1. Here are the changes:

[quote]
framework: no changes

gfndsieve, gfndsievecl: version 2.0
Moved GFN divisor testing to GFNDivisorTester class so that GFNDivisorApp is smaller
and so that future support of non-x86 is easier since GFNDivisorTester
calls a number of x86 asm methods directly.

For gfndsievecl do not report any terms with factors < 50. This reduces the size
needed for the buffer that is used to report factors.

Added -r and -R options to support functionality similar to ppsieve.
-r will not generate a bitmap for tracking terms. It will only generate an
output file of factors.
-R is used with -r. If a term has a factor below 32767 (the default value),
then the program will not output any factors for the term.
-r and -x are mutually exclusive with -r overriding -x.

Added various speed improvements.

srsieve2, srsieve2cl: version 1.5.2
Fixed issue with CisOne logic as it tries to rebuild sequences when there are
multiple sequences as that is not yet supported.
[quote]

In what limited performance testing I have done it appears that gnfdsievecl with -r is about 5x faster than the OpenCL version of ppsieve. I think some of that is due to using a much higher value for -S in gnfdsievecl than what ppsievecl does, but that doesn't explain the entirely of the speed gain. The CPU only version is about 3x slower than ppsieve, but that is due to a lot of fine-tuned assembler in the ppsieve CPU code. I don't expect anyone to use it for that, but at the same time gfndsieve should be about 50% faster for typical gfn divisor sieving.
rogue is online now   Reply With Quote
Old 2021-02-27, 18:43   #544
YaoPlaysMC
 
Nov 2020

2×5 Posts
Default

Can someone make a program for sieving sequences of the form k1*b1^n+k2*b2^n+c with variable n?

Last fiddled with by YaoPlaysMC on 2021-02-27 at 18:43
YaoPlaysMC is offline   Reply With Quote
Old 2021-08-12, 19:53   #545
mathwiz
 
Mar 2019

3108 Posts
Default

There appears to be something broken with gfndsieve at head. This is from the latest svn build:

Code:
./gfndsieve -P 1e11 -W 36 -n 10000 -N 100000 -k 1e5 -K 1e6 -o gfn.txt

gfndsieve v2.0, a program to find factors of k*2^n+1 numbers for variable k and n
Sieve started: 3 < p < 1e11 with 40500450000 terms (100001 <= k <= 999999, 10000 <= n <= 100000, k*2^n+1) (expecting 38743756770 factors)
Fatal Error:  Invalid factor: 100007*2^99457+1 mod 100271 = 13484
mathwiz is offline   Reply With Quote
Old 2021-08-12, 20:09   #546
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

2·3,217 Posts
Default

Quote:
Originally Posted by mathwiz View Post
There appears to be something broken with gfndsieve at head. This is from the latest svn build:

Code:
./gfndsieve -P 1e11 -W 36 -n 10000 -N 100000 -k 1e5 -K 1e6 -o gfn.txt

gfndsieve v2.0, a program to find factors of k*2^n+1 numbers for variable k and n
Sieve started: 3 < p < 1e11 with 40500450000 terms (100001 <= k <= 999999, 10000 <= n <= 100000, k*2^n+1) (expecting 38743756770 factors)
Fatal Error:  Invalid factor: 100007*2^99457+1 mod 100271 = 13484
I get the same error. I'll take a look.
rogue is online now   Reply With Quote
Old 2021-08-13, 05:49   #547
houding
 
houding's Avatar
 
"Adolf"
Nov 2013
South Africa

2·31 Posts
Default

Quote:
Originally Posted by rogue View Post
I get the same error. I'll take a look.

I had the same the same issue with a previous build.


You found the problem back then.


https://www.mersenneforum.org/showthread.php?t=22890
houding is offline   Reply With Quote
Old 2021-08-13, 12:03   #548
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

192216 Posts
Default

Quote:
Originally Posted by houding View Post
I had the same the same issue with a previous build.


You found the problem back then.


https://www.mersenneforum.org/showthread.php?t=22890
It might or not be the same cause. I haven't had a chance to look at it yet.
rogue is online now   Reply With Quote
Old 2021-08-13, 15:03   #549
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

2·3,217 Posts
Default

The problem with gfnsieve is fixed in sourceforge.
rogue is online now   Reply With Quote
Old 2021-08-13, 17:16   #550
Plutie
 
"Evan"
Dec 2020
Montreal

7110 Posts
Default

Similar? issue to what I reported on page 46 - latest svn

Code:
/Math/mtsieve$ ./srsieve2cl -n 100000 -N 500000 -o out.txt -s "(88*10^n-7)/9"
srsieve2cl v1.5.3, a program to find factors of k*b^n+c numbers for fixed b and variable k and n
Sieving with generic logic for p >= 3
Creating CPU worker to use until p >= 1000000
GPU primes per worker is 25600
Sieve started: 3 < p < 2^62 with 400001 terms (100000 < n < 500000, k*10^n+c)
Fatal Error:  Invalid factor: (88*10^100000-7)/9 mod 3 = 18446744073709551610
Plutie is offline   Reply With Quote
Reply

Thread Tools


All times are UTC. The time now is 17:10.


Sun Oct 24 17:10:24 UTC 2021 up 93 days, 11:39, 0 users, load averages: 1.10, 1.32, 1.31

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.