mersenneforum.org  

Go Back   mersenneforum.org > Search Forums

Showing results 1 to 25 of 1000
Search took 0.30 seconds.
Search: Posts Made By: paulunderwood
Forum: Programming 2021-10-15, 19:39
Replies: 41
Views: 981
Posted By paulunderwood
For lengthy numbers you might find this is better...

For lengthy numbers you might find this is better as there is no truncation:


ans=$(echo "print($d*$m)" | gp -q -f)
Forum: Miscellaneous Math 2021-10-14, 22:03
Replies: 2
Views: 99
Posted By paulunderwood
Semi-primes are being stubborn, but when I feed...

Semi-primes are being stubborn, but when I feed in Carmichael numbers counterexamples abound such as [n,a]=[19384289, 8494896]

This is yet another test that shows that X Frobenius tests with X...
Forum: Miscellaneous Math 2021-10-14, 01:56
Replies: 2
Views: 99
Posted By paulunderwood
I can save a few Selfridges by using the weaker...

I can save a few Selfridges by using the weaker form of Fermat's Little Theorem:

{
tst(n,a)=kronecker(a^2-4,n)==-1&&
gcd(a+4,n)==1&&
Mod(a-1,n)^n==a-1&&
Mod(a,n)^n==a&&
Mod(a+1,n)^n==a+1&&...
Forum: Hardware 2021-10-13, 20:21
Replies: 167
Views: 34,934
Posted By paulunderwood
https://thepihut.com/products/argon-one-m-2-raspbe...

https://thepihut.com/products/argon-one-m-2-raspberry-pi-4-case

Boot a pi from an M2 drive + 2 full size HDMI sockets
Forum: Miscellaneous Math 2021-10-13, 19:13
Replies: 2
Views: 99
Posted By paulunderwood
Talking Single Parameter Frobenius test -- 1+1+1+1+2 Selfridges

{
tst(n,a)=kronecker(a^2-4,n)==-1&&
gcd((a^3-a)*(a+4),n)==1&&
Mod(a-1,n)^(n-1)==1&&
Mod(a,n)^(n-1)==1&&
Mod(a+1,n)^(n-1)==1&&
Mod(a+4,n)^(n-1)==1&&
Mod(Mod(x+2,n),x^2-a*x+1)^(n+1)==2*a+5;}
...
Forum: Forum Feedback 2021-10-13, 14:34
Replies: 8
Views: 196
Posted By paulunderwood
I use hightlight then ctrl+shift+c to copy from a...

I use hightlight then ctrl+shift+c to copy from a terminal. Finally pasting with ctrl+v. To paste into a terminal use ctrl+shift+v. HTH :smile:
Forum: Riesel Prime Search 2021-10-08, 23:21
Replies: 706
Views: 76,654
Posted By paulunderwood
Here is another found with the Propper treatment:...

Here is another found with the Propper treatment:

https://primes.utm.edu/primes/page.php?id=132808 (1,935,350 digits)
Forum: Programming 2021-10-07, 02:16
Replies: 10
Views: 409
Posted By paulunderwood
I don't understand how it loses "n:". The first...

I don't understand how it loses "n:". The first line must be n and its value and not an empty line. It could be a better script, dropping "j".
Forum: Programming 2021-10-07, 01:18
Replies: 10
Views: 409
Posted By paulunderwood
#!/bin/bash readarray lines < $1 echo...

#!/bin/bash

readarray lines < $1
echo ${lines[0]}
i=1
j=${#lines[@]}
while ! [[ ${lines[$i]} =~ "skew" ]]; do
i=$(expr $i + 1)
echo ${lines[$(expr $j - $i )]}
done
Forum: Hardware 2021-10-06, 16:12
Replies: 22
Views: 4,275
Posted By paulunderwood
FWIW (not much) I implemented a base 3 Fermat PRP...

FWIW (not much) I implemented a base 3 Fermat PRP multi-limb test on my 32 bit Longan Nano. With a 16 bit model, a mult fits in in 32 bits with room for two 16 bit additions without overflow. I do...
Forum: Miscellaneous Math 2021-10-06, 15:57
Replies: 13
Views: 661
Posted By paulunderwood
David Broadhurst gave me a file (attached to this...

David Broadhurst gave me a file (attached to this message) of 100 small frauds that fool tst2, which took 42 seconds to produce.

He also found this 111 digit counterexample:


{[n,a,t,t2]=[...
Forum: Marin's Mersenne-aries 2021-10-06, 13:50
Replies: 35
Views: 1,276
Posted By paulunderwood
Someone should tell Prof. Chris Caldwell about...

Someone should tell Prof. Chris Caldwell about this so he can update his site/database.
Forum: Hardware 2021-10-06, 11:12
Replies: 22
Views: 4,275
Posted By paulunderwood
Not all the registers are general purpose -- some...

Not all the registers are general purpose -- some seem to be for I/O.

I like the small set of registers and instructions. They are easy to remember and therefore to program. Oh what fun! Takes me...
Forum: Programming 2021-10-05, 19:25
Replies: 10
Views: 409
Posted By paulunderwood
#!/bin/bash while read line; do ...

#!/bin/bash

while read line; do
larray=($line)
echo "n: input_number"
echo "c5: ${larray[0]}"
echo "c4: ${larray[1]}"
echo "c3: ${larray[2]}"
echo "c2: ${larray[3]}"...
Forum: Programming 2021-10-05, 08:00
Replies: 41
Views: 981
Posted By paulunderwood
a="31363920 948583258055568...

a="31363920 948583258055568 -21094061301960794201597 -101781123694367174730692456309"
echo "$a" > testcase
echo "22$a" >> testcase
c="31363920"
m="3"

printf "\nYAFU says $c*$m is "...
Forum: Programming 2021-10-05, 02:24
Replies: 41
Views: 981
Posted By paulunderwood
Just move the input redirection to after done...

Just move the input redirection to after done like this:

while read line
do
larray=($line)
../../yafu/yafu with some elements of the array
done < testcase
Forum: Programming 2021-10-05, 01:45
Replies: 41
Views: 981
Posted By paulunderwood
#!/bin/bash a="31363920 948583258055568...

#!/bin/bash

a="31363920 948583258055568 -21094061301960794201597 -101781123694367174730692456309 1859252514529139843616145918243322920 1347307311985305327578340438707291179185984...
Forum: Programming 2021-10-04, 22:21
Replies: 41
Views: 981
Posted By paulunderwood
One last thing to try: cm4=$(../../yafu/yafu...

One last thing to try:

cm4=$(../../yafu/yafu \"expr($c\*$m)\" -silent)

Then I'll d/l YAFU :smile:
Forum: Programming 2021-10-04, 22:02
Replies: 41
Views: 981
Posted By paulunderwood
And what happens with: cm4=$(../../yafu/yafu...

And what happens with:

cm4=$(../../yafu/yafu \"expr($c*$m)\" -silent)

?
Forum: Programming 2021-10-04, 21:53
Replies: 41
Views: 981
Posted By paulunderwood
I'm not familiar with yafu. Shoud the quotes be...

I'm not familiar with yafu. Shoud the quotes be escaped for yafu's expr in a bash expression?
Forum: Programming 2021-10-04, 21:37
Replies: 41
Views: 981
Posted By paulunderwood
I have c4m=$(../../yafu/yafu $(expr $c4 \* $am)...

I have
c4m=$(../../yafu/yafu $(expr $c4 \* $am) -silent)

You have
cm4=$(../../yafu/yafu "expr($c*$m)" -silent)

Do you see my extra $ for expr?
Forum: Programming 2021-10-04, 21:16
Replies: 41
Views: 981
Posted By paulunderwood
Here is what you should test: am=3 ...

Here is what you should test:

am=3
c4=2346021463111956
echo "c4 is $c4 and am is $am"
c4m=$(echo "expr($c4*$am)")
echo "c4 is $c4 and c4m is $c4m"
...
Forum: Programming 2021-10-04, 18:48
Replies: 41
Views: 981
Posted By paulunderwood
Of course not. I substituted ../../yafu/yafu with...

Of course not. I substituted ../../yafu/yafu with factor. You'll have to rectify it.
Forum: Programming 2021-10-04, 18:16
Replies: 41
Views: 981
Posted By paulunderwood
am=3 c4=2346021463111956 echo "c4 is...

am=3
c4=2346021463111956
echo "c4 is $c4 and am is $am"
c4m=$(echo "expr($c4*$am)")
echo "c4 is $c4 and c4m is $c4m"
c4m=$(factor $(expr $c4 \* $am))
echo "c4 is...
Forum: Programming 2021-10-04, 17:36
Replies: 41
Views: 981
Posted By paulunderwood
Attach the script for us to debug. Specify which...

Attach the script for us to debug. Specify which arguments you might be using with the script.

:smile:
Showing results 1 to 25 of 1000

 
All times are UTC. The time now is 05:11.


Sat Oct 16 05:11:51 UTC 2021 up 84 days, 23:40, 0 users, load averages: 1.03, 1.09, 1.12

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.