1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:
46:
47:
48:
49:
50:
51:
52:
53:
54:
55:
56:
57:
58:
59:
60:
61:
62:
63:
64:
65:
66:
67:
68:
69:
70:
71:
72:
73:
74:
75:
76:
77:
78:
79:
80:
81:
82:
83:
84:
85:
86:
87:
88:
89:
90:
91:
92:
93:
94:
95:
96:
97:
98:
99:
100:
101:
102:
103:
104:
105:
106:
107:
108:
109:
110:
111:
112:
113:
114:
115:
116:
117:
118:
119:
120:
121:
122:
123:
124:
125:
126:
127:
128:
129:
130:
131:
132:
133:
134:
135:
136:
137:
138:
139:
140:
141:
142:
143:
144:
145:
146:
147:
148:
149:
150:
151:
152:
153:

Function IsPrime(Number: DWORD): Boolean;
Asm
  PUSH  EBX
  PUSH   ESI
  XOR    EBX, EBX
  MOV    ESI, EAX
  CMP    ESI, $00000002
  JB    @@IsPrime_NoPrime
  JZ    @@IsPrime_Prime

  TEST  ESI, $00000001
  JZ    @@IsPrime_NoPrime
  
  MOV    ECX, $03
  JMP    @@IsPrime_CheckIfLess

@@IsPrime_Loop:
  MOV    EAX, ESI
  XOR    EDX, EDX
  DIV   ECX
  TEST  EDX, EDX
  JZ    @@IsPrime_NoPrime

  ADD   ECX, $02

@@IsPrime_CheckIfLess:
  MOV    EAX, ECX
  MUL    ECX
  CMP    EAX, ESI
  JA    @@IsPrime_Loop
  JZ    @@IsPrime_NoPrime
  
@@IsPrime_Prime:
  MOV    BL, $01
@@IsPrime_NoPrime:
  MOV    EAX, EBX
@@Ret:
  POP   ESI
  POP   EBX
End;

Function IsPrime2(Number: DWORD): Boolean;
Var
  Teiler: DWORD;
Begin
  Result := False;
  If Number < 2 Then
    Exit
  Else If Number = 2 Then
  Begin
    Result := True;
    Exit;
  End;

  If Number Mod 2 = 0 Then
    Exit;

  Teiler := 3;
  While Teiler * Teiler <= Number Do
  Begin
    If Teiler * Teiler <= Number Then
      Exit;

    If Number Mod Teiler = 0 Then
      Exit;

    Inc(Teiler, 2);
  End;

  Result := True;
End;

Procedure TForm1.TestClick(Sender: TObject);
Var
  T1, T2, T3, T4, T5: Int64;
  X, Y, Z: DWORD;
  B: Boolean;
Const
  NumTests = 500;
  ScaleBy = 50000;
  ScaleByY = 25;
  MaxNumber = 50000000;
  PerPixel = $7F;
Begin
  QueryPerformanceFrequency(T1);
  Image1.Height := 500;
  Image1.Width := 1 + MaxNumber Div ScaleBy;
  Image1.Picture.Bitmap.Height := Image1.Height;
  Image1.Picture.Bitmap.Width := Image1.Width;
  Image1.Picture.Bitmap.Canvas.Brush.Color := clBlack;
  Image1.Picture.Bitmap.Canvas.FillRect(Image1.Picture.Bitmap.Canvas.ClipRect);
  QueryPerformanceCounter(T2);
  QueryPerformanceCounter(T4);
  T4 := T4 - T2;
  With Image1.Picture.Bitmap, Canvas Do
  Begin
    T5 := 0;
    For Y := 1 To MaxNumber Do
    Begin
      B := True;
      QueryPerformanceCounter(T2);
      For X := 1 To NumTests Do
        B := IsPrime(Y) And B;
      QueryPerformanceCounter(T3);

      T3 := T3 - T2 - T4;
      If T3 < NumTests Then
        T3 := NumTests;
      T5 := T5 + T3;

      X := Height - Round(ScaleByY * Ln(T3 / NumTests)) - 1;
      Z := Pixels[Y Div ScaleBy, X];
      If Z and $000000FF < $100-PerPixel Then
      Z := Z +PerPixel;
      Pixels[Y Div ScaleBy, X] := Z;

      If Y Mod ScaleBy = 0 Then
      Begin
//        X := Height - Round(ScaleByY * Ln(T5 / NumTests / ScaleBy)) - 1;
//        Pixels[Y Div ScaleBy, X] := Pixels[Y Div ScaleBy, X] Or $000000FF;
//        T5 := 0;
        Update;
      End;
    End;
    T5 := 0;
    For Y := 1 To MaxNumber Do
    Begin
      B := True;
      QueryPerformanceCounter(T2);
      For X := 1 To NumTests Do
        B := IsPrime2(Y) And B;
      QueryPerformanceCounter(T3);
      T3 := T3 - T2 - T4;
      If T3 < NumTests Then
        T3 := NumTests;
      T5 := T5 + T3;

      X := Height - Round(ScaleByY * Ln(T3 / NumTests)) - 1;
      Z := Pixels[Y Div ScaleBy, X];
      If Z and $0000FF00 < $100 * ($100-PerPixel) Then
      Z := Z +$100*PerPixel;
      Pixels[Y Div ScaleBy, X] := Z;

      If Y Mod ScaleBy = 0 Then
      Begin
//                X := Height - Round(ScaleByY * Ln(T5 / NumTests / ScaleBy)) - 1;
//                Pixels[Y Div ScaleBy, X] := Pixels[Y Div ScaleBy, X] Or $0000FF00;
//                T5 := 0;
        Update;
      End;
    End;
  End;
End;