![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
| 言語 | 処理時間 | 備考 |
| Python 3.7.4 | 324.276秒 | コードは(9-5) |
| Ruby 2.6.4 | 98.800秒 | コードは(9-6) |
| ES-BASIC (32bit) | 2.906秒 | コードは(9-3)、JITコンパイルに要した時間も含む |
| gcc(C言語)(32bit) | 1.954秒 | コードは(9-2)、もちろん最適化レベル最大、コンパイル時間含まず |
| ES-BASIC (32bit、レジスタ変数利用) | 1.820秒 | コードは(9-7)、JITコンパイルに要した時間も含む |
| gcc(C言語)(64bit) | 1.195秒 | コードは(9-4)、もちろん最適化レベル最大、コンパイル時間含まず |
#include <stdio.h>
#include <time.h>
#define MAX 447
#define STEP (14.0 / 0x100000)
#define DN 4
int main()
{
int x, y, n, sx, sy, sn, h = 0;
double zx, zy, cx, cy, xx, yy;
for (y = 0; y < 384; y++) {
for (x = 0; x < 512; x++) {
sn = 0;
for (sx = 0; sx < DN; sx++) {
cx = (double) 0x4750 / 0x10000 + (x * DN + sx) * (STEP / DN);
for (sy = 0; sy < DN; sy++) {
cy = (double) 0x1e8 / 0x10000 + (y * DN + sy) * (STEP / DN);
zx = zy = 0;
for (n = 0; n < MAX; n++) {
xx = zx * zx;
yy = zy * zy;
if (xx + yy > 4) break;
zy = zx * zy * 2 - cy;
zx = xx - yy + cx;
}
sn += n;
}
}
// n = sn / (DN * DN);
n = sn >> 4;
h += n; // 表示の代わりに集計
}
}
printf("%d\n", h);
printf("(%.3f[sec])\n", clock() / (double) CLOCKS_PER_SEC);
return 0;
}1 H=0 2 FOR Y=0,383 3 FOR X=0,511 4 SN=0 5 FOR SX=0,3 6 CX=(X*4+SX)*56+4673536 7 FOR SY=0,3 8 CY=(Y*4+SY)*56+124928 9 ZX=0 10 ZY=0 11 FOR N=0,446 12 XX=ZX*>>ZX,24; // これはZXとZXをIMULで計算して64bitの結果を生成した後で、SHRD命令で24ビット右シフトして下位32bitを結果とする命令 13 YY=ZY*>>ZY,24 14 T=XX+YY 15 IF T>0X4000000 GOTO SKIP 16 ZY=ZX*>>ZY,23 17 ZX=XX-YY+CX 18 ZY=ZY-CY 19 NEXT 20 LABEL SKIP 21 SN=SN+N 22 NEXT 23 NEXT 24 N=SN>>4; // N=SN/16 (N=0...447) 25 H=H+N; // 表示の代わりに集計 26 NEXT 27 NEXT 28 PRINT H
#include <stdio.h>
#include <time.h>
#define MAX 447
#define DN 4
int main()
{
long long x, y, sx, sy, zx, zy, cx, cy, xx, yy;
int n, sn, h = 0;
for (y = 0; y < 384; y++) {
for (x = 0; x < 512; x++) {
sn = 0;
for (sx = 0; sx < DN; sx++) {
cx = (x * 4 + sx) * 56 + 4673536;
for (sy = 0; sy < DN; sy++) {
cy = (y * 4 + sy) * 56 + 124928;
zx = zy = 0;
for (n = 0; n < MAX; n++) {
xx = (zx * zx) >> 24;
yy = (zy * zy) >> 24;
if (xx + yy > 0x4000000) break;
zy = (zx * zy) >> 23;
zx = xx - yy + cx;
zy = zy - cy;
}
sn += n;
}
}
// n = sn / (DN * DN);
n = sn >> 4;
h += n; // 表示の代わりに集計
}
}
printf("%d\n", h);
printf("(%.3f[sec])\n", clock() / (double) CLOCKS_PER_SEC);
return 0;
}import time
start = time.time()
h = 0
CX0 = 0x4750 / 0x10000
CY0 = 0x1e8 / 0x10000
STP = 14.0 / 0x100000 / 4
for y in range(384):
for x in range(512):
sn = 0
for sx in range(4):
cx = CX0 + (x * 4 + sx) * STP
for sy in range(4):
cy = CY0 + (y * 4 + sy) * STP
zx = zy = 0
n = 0
while n < 447:
xx = zx * zx
yy = zy * zy
if xx + yy > 4: break
zy = zx * zy * 2 - cy
zx = xx - yy + cx
n = n + 1
sn += n
n = sn >> 4
h += n
print(h)
print(time.time() - start)start_time = Time.now
h = 0
CX0 = 0x4750 / 65536.0
CY0 = 0x1e8 / 65536.0
STP = 14.0 / 0x100000 / 4
y = 0
while y < 384 do
x = 0
while x < 512 do
sn = 0
sx = 0
while sx < 4 do
cx = CX0 + (x * 4 + sx) * STP
sy = 0
while sy < 4 do
cy = CY0 + (y * 4 + sy) * STP
zx = zy = 0
n = 0
while n < 447 do
xx = zx * zx
yy = zy * zy
if xx + yy > 4 then
break
end
zy = zx * zy * 2 - cy
zx = xx - yy + cx
n = n + 1
end
sn += n
sy = sy + 1
end
sx = sx + 1
end
n = sn >> 4
h += n
x = x + 1
end
y = y + 1
end
p h
p Time.now - start_time1 H=0 2 $EDI = EDI; // EDIレジスタを自由に使いたいときの前処理 3 FOR Y=0,383 4 FOR X=0,511 5 SN=0 6 FOR SX=0,3 7 CX=(X*4+SX)*56+4673536 8 FOR SY=0,3 9 CY=(Y*4+SY)*56+124928 10 ALIAS ZX:ECX, ZY:EBX, XX:ESI, YY:EDI 11 ZX=0 12 ZY=0 13 FOR N=0,446 14 XX=ZX*>>ZX,24 15 YY=ZY*>>ZY,24; // 結果はEAXにも入っている. 16 EAX=EAX+XX 17 IF EAX>0X4000000 GOTO SKIP 18 ZY=ZX*>>ZY,23 19 XX=XX-YY+CX 20 ZX=XX 21 ZY=ZY-CY 22 NEXT 23 LABEL SKIP 24 SN=SN+N 25 NEXT 26 NEXT 27 N=SN>>4; // N=SN/16 (N=0...447) 28 H=H+N; // 表示の代わりに集計 29 NEXT 30 NEXT 31 PRINT H
| コメント | お名前 | NameLink | |
