The funny thing is that to an extent that I've not measured yet, it might be possible to find families of higher-yield primes!
With a single k value the skew would be trivial to find, but with 12...
Type: Posts; User: FatPhil
The funny thing is that to an extent that I've not measured yet, it might be possible to find families of higher-yield primes!
With a single k value the skew would be trivial to find, but with 12...
I was going to implement it that way. I've not had a chance to work on it yet, as when I finished my other project I fell ill with a flu which has prevented me from thinking straight. I'm psyching...
To counter the flaring that you get a the high ps, you can either band the ps by their logs rather than by their numeric p value, or simply make the bands 10 times wider.
In some ways the log/log...
"""
$ gcc --version
[xxxxxxxx:~/Temporary Items/Prime] music% gcc --version
gcc (GCC) 3.1 20020420 (prerelease)
$ uname -a
Darwin xxxxxx.xxxxx.xxx.xxx 6.3 Darwin Kernel Version 6.3: Sat Dec 14...
Yeah, sorry Nuri, I've been messin' you about with my early morning bug-ridden postings. My second post had the correct requests. As it's early morning again, I'm frightened to make any requests!
...
When I do an 011 build (this weekend) I'll gather together a source archive. I can't distribute the full source, but I can point you to the missing files - a marginally slower version of the prime...
Too early in the morning... the derivative thereof (and the negative of that - third time lucky!)
I.e. it ought to slowly creep below the 1/p line (by a squared logarithmic rate, which to be honest...
Wonderful!
To my eye, it's looking like 1/log(p), which is exactly as predicted. If you want to try to get a straight line, then simply plot #factors removed against log(p). Can you do another...
The program doesn't exit cleanly unless it's finished. The only other way to stop it is to kill it. So detecting that requires signal handling, whichI can do under most Unices, but would know how...
Odd. I decreased the progress reporting frequency by a factor of 8, and in theory didn't touch the batch file frequency (10 mins IIRC). It doesn't write a file when exiting, only when a 10 minute...
Yeah, NbeGon stopped development (or I stopped development on it), 30 seconds before I handed my source code and description to Paul. I realised that the only thing that I could do further was...
It sounds like you're still at my 009 tech. Unrolling by 2, which was my 009-010 speedup, gave another 50%. The C I gave you _should_ translate 1<->1 into assembly, once you've mapped variables onto...
Yes!
These families are deliberately chosen with properties that mean they sieve down to almost nothing. Basically, vast proportions of the candidates have very small divisors, and simply by...
1) Absolutely. While I was soak-testing 010 before release I noticed that had didn;t output (one factor) from Paul's SoBSieve. The maths told me I was right to remove the candidate, so why didn't...
Indeed.
The algorithm that depends on the number of candidates is the one we don't use! (it's effectively just glorified trial-division.) Pretty useless for a situation like this, forget I...
Which didn't work. There's no point in me gnashing if I have no new ideas to chew on, so I've sent the 010 code, explanation and a few musings to Paul. I hope that he can either provide some new...
The n in that is the whole range, not the number of candidates remaining. So it won't speed up anything apart from the file load part of the restart times. The only algorithm that depends on the...
I had an interesting brainwave today - I'll see if I can push you over the limit. I'm hoping I can push Louie's 88000 over too. 12% shouldn't be too hard (famous last words). This idea also splits...
Snap.
My 010 has just passed a very lengthy test, and is available from http://fatphil.org/maths/sierpinski/bin/ . One and a half times the speed of 009 on linux and windows, as hoped.
The...
Thanks Nuutti - I've just worked out what the problem is.
(Literally 30 seconds ago - so yes, I'm not at HUT, I'm debugging instead, but will be at the beer festival - fancy coming? Gallows Bird,...
The first bug report is in! However it is so far only reproducable on alien88's OpenBSD's machine with version 009.
The following I have tested just now and do not have the error spotted:
OpenBSD...
Good suggestion. v010 will be out tomorrow. It promises to be nearly 50% faster (that is a 3:2 ratio) than v009. (I made one lousy typo when doing the optimisation yesterday).
I'm praying that...
I just think it would be stupid for me to spend a day trying to make this optimisation again, and to fail again, so I'm making sure that that with which I am happy is available.
009 is between one...
I think it's cooperation works best when both brains are truly primed for tackling the problem, cobwebs dusted off and all that. I think we're both at ramming speed now, aren't we? The amusing thing...
Have you tried setting the priority to "idle"?
Phil