mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Software (https://www.mersenneforum.org/forumdisplay.php?f=10)
-   -   Software for IBM AIX (https://www.mersenneforum.org/showthread.php?t=13801)

pacionet 2010-08-30 09:18

Software for IBM AIX
 
I can run some software on some IBM AIX servers.
Until now the only program I can build and run successfully is glucas.
All the others have compile problems, etc.

Are there any other software that can be compiled and run (in a simple way !!!) for IBM AIX (powerpc) to do some prime search ?

I am looking for some software for sieve and/or search of primes; also "small" primes would be good.

ldesnogu 2010-08-30 10:03

"All the others" is rather fuzzy :smile: Which one did you try exactly?

Did you try [url=http://www.hogranch.com/mayer/README.html]mlucas[/url]?

pacionet 2010-08-30 10:13

Yes, I remember some difficulties to compile it.
Anyway I am searching some software for a different kind of primes, for Mersenne I am using glucas.

I didn't think that compile and build on this platform was so annoying !

mdettweiler 2010-08-30 15:27

You may want to check out [url=http://www.mersenneforum.org/showthread.php?p=202851#post202851]Phrot[/url], which is based on Glucas and can do:

-PRP tests on all k*b^n+-c numbers
-Proth tests (a full primality test) on k*2^n+1
-PRP tests on Generalized Fermat numbers

Note that any positive results returned by the PRP tests will need to be proven prime separately; aside from k*2^n+1, which Phrot can prove directly, you'll need to prove the primes with another program such as LLR or PFGW on an x86 machine. However, for most searches that would only be a small fraction of the overall computing time.

pacionet 2010-08-30 19:39

After some time spending trying to install Phrot i find some errors during compiling (it seems that it misses libm functions (sqrt, etc.))

It seems that I should install a particular libm which is not in standard AIX 5.3 distribution ... quite frustrating !!!

xilman 2010-08-30 19:53

[quote=pacionet;227758]After some time spending trying to install Phrot i find some errors during compiling (it seems that it misses libm functions (sqrt, etc.))

It seems that I should install a particular libm which is not in standard AIX 5.3 distribution ... quite frustrating !!![/quote]Sounds most unlikely to me. Over 12 years ago, when I last played with AIX, libm came as standard and had the all usual functions available. Are you sure that you are linking with "-lm" and that the linker can find the libm library? Try "find / -name libm.\*" and see what turns up.

Paul

rogue 2010-08-31 01:19

[QUOTE=pacionet;227758]After some time spending trying to install Phrot i find some errors during compiling (it seems that it misses libm functions (sqrt, etc.))

It seems that I should install a particular libm which is not in standard AIX 5.3 distribution ... quite frustrating !!![/QUOTE]

It could be a 32-bit vs 64-bit issue.

pacionet 2010-08-31 10:09

After some modifications:

1) export OBJECT_MODE=64
2) Adding -maix64 to compile flags
3) Commentig this line of mathtypes.h, because of conflicting definition

//typedef long long int64;

I got the following error:

[CODE]gcc -maix64 -lm -Wall \
-ffast-math \
-o phrot.g5 \
phrot.c \
-O3 -DY_MEM_THRESHOLD=8192 -DY_KILL_BRANCHES -DY_VECTORIZE \
-I../glucas ../glucas/libyeafft.a
In file included from phrot.c:57:
mathtypes.h:7: warning: ignoring #pragma warn
phrot.c: In function 'doPRPTest':
phrot.c:886: warning: format '%lli' expects type 'long long int', but argument 3 has type 'int64'
phrot.c:995: warning: format '%lli' expects type 'long long int', but argument 3 has type 'int64'
[B]ld: 0711-317 ERROR: Undefined symbol: Y_NTHREADS
ld: 0711-317 ERROR: Undefined symbol: Y_2NN
ld: 0711-317 ERROR: Undefined symbol: Y_BS
ld: 0711-317 ERROR: Undefined symbol: Y_DI
ld: 0711-317 ERROR: Undefined symbol: Y_DJ
ld: 0711-317 ERROR: Undefined symbol: Y_IR
ld: 0711-317 ERROR: Undefined symbol: Y_NC
ld: 0711-317 ERROR: Undefined symbol: Y_2N0[/B]
ld: 0711-345 Use the -bloadmap or -bnoquiet option to obtain more information.
collect2: ld returned 8 exit status
make: *** [phrot.g5] Error 1[/CODE]

Some libraries still miss ?

pacionet 2010-08-31 10:26

Solved with :

1) ar rc libyeafft.a *.o

instead of

ar rc libyeafft.a [drty]*.o

as in Proth documentation

and changing Makefile

gcc [B]-maix64 -lm -lpthread[/B] ${GCC_MARCH} -Wall

Still got this warnings

[B]phrot.c:886: warning: format '%lli' expects type 'long long int', but argument 3 has type 'int64'
phrot.c:995: warning: format '%lli' expects type 'long long int', but argument 3 has type 'int64'[/B]

Now I try to run some self test (if any) to test Proth working

ldesnogu 2010-08-31 10:56

[quote=pacionet;227868]gcc [B]-maix64 -lm -lpthread[/B] ${GCC_MARCH} -Wall

Still got this warnings

[B]phrot.c:886: warning: format '%lli' expects type 'long long int', but argument 3 has type 'int64'
phrot.c:995: warning: format '%lli' expects type 'long long int', but argument 3 has type 'int64'[/B][/quote]
To fix this warning the best solution is to use PRId64 (from <inttypes.h>. The less portable one is to replace %lli with %li.

pacionet 2010-08-31 11:44

Mmmmm...
Now I builded successfully Phrot (no errors or warning) but I got the following result

bash-3.2$ [B]phrot.g5 -q 5028*10^83982+1[/B]
Phil Carmody's Phrot (0.72)
Input 5028*10^83982+1 : Actually testing 5028*1000000^13997+1 (witness=3 13998/28672 limbs)
5028*10^83982+1 is composite LLR64=0000000000000001. (e=0.00000 (0.0824785~0@0.000) t=133.91s)

On this proven prime number:

[url]http://primes.utm.edu/primes/page.php?id=79257[/url]


I fear that the problem is in maththypes.h

[B]//typedef long long int64;[/B]

while in <inttypes.h>

[B]#ifdef __64BIT__
typedef long int64;
#else /* _ILP32 */
#if defined(_LONG_LONG)
typedef signed long long int64;
#endif
#endif[/B]

If i don't comment the line

[B]//typedef long long int64;[/B]

I got the following error:

In file included from phrot.c:57:
[B]mathtypes.h:9: error: conflicting types for 'int64'[/B]
/usr/include/inttypes.h:622: error: previous declaration of 'int64' was here


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

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