Dann eine Weg suchen zu parallelisieren. 
Code: Alles auswählen
psy@BrAIN1ak> ./letzte_primzahl.py 100000
99991 ( 0.0576159954071 secs)
psy@BrAIN1ak> ./letzte_primzahl.py 1000000
999983 ( 0.548065900803 secs)
psy@BrAIN1ak> ./letzte_primzahl.py 10000000
9999991 ( 5.59474992752 secs)
Code: Alles auswählen
psy@BrAIN1ak> ./letzte_primzahl_psyco.py 100000
99991 ( 0.00420188903809 secs)
psy@BrAIN1ak> ./letzte_primzahl_psyco.py 1000000
999983 ( 0.0389289855957 secs)
psy@BrAIN1ak> ./letzte_primzahl_psyco.py 10000000
9999991 ( 0.379332065582 secs)
Code: Alles auswählen
#!/usr/bin/python
from sys import argv
from time import time
import psyco
psyco.full()
stop = int(argv[1])
def quer(x):
x = [int(n) for n in list(str(x))]
return sum(x)
def is_prim(x):
if str(x).endswith("5"): return False
elif int(str(x)[-1])%2 == 0: return False
elif x>10:
if quer(x)%3 == 0: return False
n = 2
while n<x:
if x%n == 0: break
else: n+=1
else:
return True
return False
stamp = time()
while not is_prim(stop): stop -= 1
stamp=time()-stamp
print stop, '(', stamp, 'secs)'
Code: Alles auswählen
psy@BrAIN1ak> ./a.out
99991 (0 clockcycles) # <1 secs
999983 (30000 clockcycles) # 0.03 secs
9999991 (270000 clockcycles) # 0.27 secs
Code: Alles auswählen
/*
* Gibt die letzte Primzahl vor der per Argument uebergebenen Zahl aus.
*/
#include <stdio.h>
#include <time.h>
int quer(int x)
{
int wertigkeit, hilf, summe = 0;
hilf=x;
wertigkeit = 1;
while((hilf=hilf/10) != 0) wertigkeit=wertigkeit * 10;
do {
summe = summe + x / wertigkeit;
x = x % wertigkeit;
wertigkeit = wertigkeit / 10;
} while(wertigkeit > 0);
return summe;
}
int is_prim(int x)
{
int n = 2;
if (x%5 == 0) return 0;
if (x%2 == 0) return 0;
if (x>10) {
if (quer(x)%3 == 0) return 0;
}
while (n<x) {
if (x%n == 0) return 0;
else n++;
}
return 1;
}
int main()
{
int stop;
clock_t stamp;
stop = 100000;
stamp = clock();
while (!is_prim(stop)) stop--;
stamp = clock() - stamp;
printf("%d (%d clockcycles)\n", stop, stamp);
stop = 1000000;
stamp = clock();
while (!is_prim(stop)) stop--;
stamp = clock() - stamp;
printf("%d (%d clockcycles)\n", stop, stamp);
stop = 10000000;
stamp = clock();
while (!is_prim(stop)) stop--;
stamp = clock() - stamp;
printf("%d (%d clockcycles)\n", stop, stamp);
return 0;
}
ist der lang 
Code: Alles auswählen
$ ./a.out
99991 (0 clockcycles)
999983 (0 clockcycles)
9999991 (0 clockcycles)Code: Alles auswählen
#include <math.h>
#include <stdio.h>
#include <time.h>
int is_prim(int x)
{
int n;
if (x == 2) return 1;
if (!(x & 1)) return 0;
for (n = 3; n < sqrt(x); n += 2) {
if (!(x % n)) return 0;
}
return 1;
}
int main()
{
int i, stop;
clock_t stamp;
for (i = 0; i < 3; ++i) {
stop = 100000 * pow(10, i);
stamp = clock();
while (!is_prim(stop)) stop--;
stamp = clock() - stamp;
printf("%d (%d clockcycles)\n", stop, (int) stamp);
}
return 0;
}
Code: Alles auswählen
#include <stdio.h>
#include <math.h>
int main(int argc, char *argv)
{
printf("> %f <\n", sqrt(pow(5, 2)));
return 0;
}
Code: Alles auswählen
#!/usr/bin/python
from sys import argv
from time import time
import psyco
psyco.full()
stop = int(argv[1])
def is_prim(x):
if x==2: return True
if not (x & 1): return False
n = 3
while n < (x**0.5):
if not x%n: return False
n += 2
return True
stamp = time()
while not is_prim(stop): stop -= 1
stamp=time()-stamp
print stop, '(', stamp, 'secs)'
Code: Alles auswählen
$ gcc -O3 -Wall -ansi -pedantic test.c
/tmp/ccGLyIIq.o: In function `is_prim':test.c:(.text+0x47): undefined reference to `sqrt'
/tmp/ccGLyIIq.o: In function `main':test.c:(.text+0xcd): undefined reference to `pow'
:test.c:(.text+0x1a1): undefined reference to `sqrt'
collect2: ld returned 1 exit statusCode: Alles auswählen
$ gcc -O3 -Wall -ansi -pedantic -lm test.c
Da hast du wohl recht... Allerdings ist es IMHO z.Z. so, das ein Programmiere keine Möglichkeiten hat sein Programm zu parallelisieren. Der einzige Weg ist es wohl ein Python Interpreter pro Prozessor/Core zu starten.BlackJack hat geschrieben:Ich weiss nicht ob sich das jeder der Diskutanten klar macht: Ein Programm kann nicht auf magische Weise mehrere Prozessoren/Kerne benutzen. Da muss der Programmierer immer etwas für tun.
)Das relevante ist aber das Paket python-mpi welches Teil von Scientific python ist (ganz recht, da kommt ein 404, leider), daneben gibt es noch pythonMPI (dito), pympi und pypar.Rebecca hat geschrieben:(mpichpython: Meinten Sie teppichpython? lampython: Meinten Sie baumpython?