mersenneforum.org  

Go Back   mersenneforum.org > Prime Search Projects > Sierpinski/Riesel Base 5

Reply
 
Thread Tools
Old 2006-11-19, 08:14   #1
michaf
 
michaf's Avatar
 
Jan 2005

479 Posts
Default Request from Phil

I just read on the primenumbers forum the following request:
Can we provide him with such a list?

Cheers, Micha

I'm trying to resuscitate some old code for testing proth-like
numbers (in part so that people with non-x86 architectures can
participate in such searches), and would very much like a _complete_
list of all of the primes found in the Riesel/Sierpinski-5 searches
so that I can test my code very thoroughly.

(Mark R has already been testing old versions of the code and and
not-so-merrily finding several unpleasant, but in retrospect doomed
to fail, error cases on his G5 - many many thanks Mark!)

The Mersenne forum pages do seem to contain most of the information,
but it's all smooshed out across multiple threads and bogillions of
posts, and alas the most promising posts contained dead links, which
is a tad inconvenient.

Who's actually running the project? I assumed it was Robert as I
thought it was his braichild, but the forum posts imply that perhaps
he's not, but maybe I misunderstand.

Phil
michaf is offline   Reply With Quote
Old 2006-11-19, 09:39   #2
axn
 
axn's Avatar
 
Jun 2003

111318 Posts
Default

http://geocities.com/base5_sierpinski_riesel/

The "primes.txt" under Riesel & Sierpinski directories has all the project primes. (It is not yet updated with the latest 5 or 6 primes).
axn is offline   Reply With Quote
Old 2006-11-19, 11:50   #3
fatphil
 
fatphil's Avatar
 
May 2003

3·7·11 Posts
Default Many thanks!

I don't intend to test with all of the largest primes, as I don't believe size is the troublesome issue, so I'll test with tens of thousands of small primes first. If that all goes well, I'll get back to you. Mark has found a failure case on PPC, but so far everything works OK on x86. Does anyone else apart from Mark and I have non-x86 machines?
fatphil is offline   Reply With Quote
Old 2006-11-19, 11:54   #4
fatphil
 
fatphil's Avatar
 
May 2003

3×7×11 Posts
Default

Oh bollocks:

"The GeoCities web site you were trying to view has temporarily exceeded its data transfer limit. Please try again later."

I think I grabbed the whole riesel file, not didn't get the sierpinski one. Fortunately, I'm only looking at riesel ones currently.
fatphil is offline   Reply With Quote
Old 2006-11-19, 12:27   #5
axn
 
axn's Avatar
 
Jun 2003

7×11×61 Posts
Default

Quote:
Originally Posted by fatphil View Post
Oh bollocks:

"The GeoCities web site you were trying to view has temporarily exceeded its data transfer limit. Please try again later."

I think I grabbed the whole riesel file, not didn't get the sierpinski one. Fortunately, I'm only looking at riesel ones currently.
PM me your email-id and I can send you the both the files.
axn is offline   Reply With Quote
Old 2006-11-19, 13:05   #6
axn
 
axn's Avatar
 
Jun 2003

7·11·61 Posts
Default

Quote:
Originally Posted by fatphil View Post
I think I grabbed the whole riesel file, not didn't get the sierpinski one. Fortunately, I'm only looking at riesel ones currently.
YGM
axn is offline   Reply With Quote
Old 2006-11-20, 23:15   #7
fatphil
 
fatphil's Avatar
 
May 2003

3·7·11 Posts
Default What kinds of speeds is LLR running at?

I'm reticent to hit asm as I'd like to keep my code portable.
Can anyone tell me roughly what times a 2GHz AthlonXP takes for the following, for example. I don't believe more than a couple of percent improvement would be possible on a register starved architecture like the x86. I guess I could actually download LLR at some point...

Actually testing 194196250*1953125^11117-1 (11119/24576 limbs)
310714*5^100057-1 [] is PRP. (t=364.70s)
Actually testing 5718593750*1953125^11264-1 (11266/24576 limbs)
73198*5^101383-1 [] is PRP. (t=376.42s)
Actually testing 34058000*1953125^11296-1 (11298/24576 limbs)
272464*5^101667-1 [] is PRP. (t=379.82s)
Actually testing 97943750*1953125^11394-1 (11396/24576 limbs)
31342*5^102551-1 [] is PRP. (t=384.77s)
Actually testing 172598750*1953125^11440-1 (11442/24576 limbs)
276158*5^102964-1 [] is PRP. (t=369.83s)
Actually testing 1876900*1953125^11509-1 (11510/24576 limbs)
75076*5^103583-1 [] is PRP. (t=377.25s)
Actually testing 197946250*1953125^11871-1 (11873/24576 limbs)
316714*5^106843-1 [] is PRP. (t=395.33s)
Actually testing 41354500*1953125^11925-1 (11927/24576 limbs)
330836*5^107328-1 [] is PRP. (t=399.21s)

(No erronious results on x86 - from about 15000 tests - looks good!)

Phil
fatphil is offline   Reply With Quote
Old 2006-11-20, 23:55   #8
axn
 
axn's Avatar
 
Jun 2003

469710 Posts
Default

Quote:
Originally Posted by fatphil View Post
Can anyone tell me roughly what times a 2GHz AthlonXP takes for the following, for example.
Don't have any data on Athlons, but on my 2.8GHz Pentium D, tests of this size were taking appr. 460 sec. So I would guess that an Athlon XP 2GHz might take something in the vicinity of 800 sec?
axn is offline   Reply With Quote
Old 2006-11-21, 00:25   #9
tnerual
 
tnerual's Avatar
 
Oct 2006

7×37 Posts
Default

i will start a test on a athlon xp 2500+ @1980mhz.

is it good for your benchmark ?

i can't get closer

edit:
310714*5^100057-1 is a probable prime. Time: 1068.439 sec.

Last fiddled with by tnerual on 2006-11-21 at 00:58
tnerual is offline   Reply With Quote
Old 2006-11-21, 08:14   #10
michaf
 
michaf's Avatar
 
Jan 2005

479 Posts
Default

With an Athlon 64 X2 3800+ running at 2260 MHz I get about 3.78 ms per bit on number of size 5^138000-ish
michaf is offline   Reply With Quote
Old 2006-11-21, 08:56   #11
fatphil
 
fatphil's Avatar
 
May 2003

E716 Posts
Default Redux

My stake in the ground (hmmmm, steak...) is 360s for 5^100000 at 2.083GHz, say. And the responses vary from:

* 1068s on 1980MHz -> 1015s
* ~460s on 2.8GHz P4 -> ~800s
* bigger number on Ath64 -> ~620s

I, so far, have only built my code for win/P3 and linux/ath. In windows, I use cygwin. The code is GPL, and so people can build their own (and are invited to, tweak the compiler options - tweak the code - or anything to make it better!). However, it requires building YEAFFT (from GLucas), and last time I looked, SSE2 didn't work in YEAFFT on cygwin. (I patched it last year, but the patches never got into the main line.)

I will try to build for other architectures when I can. I hope that linux/P4 and linux/Ath64 will make good use of SSE2.

So, those who are familiar with the concept of how I code - stuff's _always_ work in progress (ask Mark - he's had 5 code drops in as many days!) - are invited to volunteer to help me test this on real-world candidates. Please state chip type and OS, so we can work out the best way to get an executable to you (I build or you build).
fatphil is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
Request pinhodecarlos Lounge 3 2017-10-26 18:58
Bug/request Dubslow YAFU 4 2012-03-31 03:07
Odd request? Xyzzy Lounge 23 2011-03-08 17:50
GMP-ECM Request rogue GMP-ECM 4 2009-11-23 15:07
A little request 10metreh Aliquot Sequences 8 2009-10-07 18:13

All times are UTC. The time now is 07:16.

Sun Sep 27 07:16:57 UTC 2020 up 17 days, 4:27, 0 users, load averages: 2.23, 1.73, 1.58

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2020, 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.