Delphi-Quellcode:
for k:=0 to MaxPrime div PrimeBits do begin
FillChar(Primes[0], PrimeLen, 0);
for i:=0 to Trunc(Sqrt((k+1)*PrimeBits)/30) do
Sqrt(k+1); was wäre wenn k von Anfang an schon von 0 bis Sqrt(MaxPrime div PrimeBits) laufen würde ?
Die Schleife I würde dann von 0 bis (k * PrimeBits) div 30 laufen, richtig ?
Ergo: in der äußeren Schleife j würde der Endwert durch den Compiler nur einmalig berechnet und in der Schleife j entfällt die langsamme Fließkommaarithmetik komplett und wird durch viel schnellere Integer Arithmetik ersetzt.
Natürlich kannst du hier auch wieder einige Mul's einsparen.
Schau dir mal meinen Source genauer an. Er macht ja exakt das gleiche wie der deinige, kommt aber ohne Fließkommazahlen aus und benutzt weniger Divisionen.
1.) Fließkomma weg
2.) Divisionen raus
3.) Multiplikationen durch Shift+Adds ersetzen
So sollte deine Optimierungsstrategie sein.
Gruß Hagen