Last active
August 10, 2018 10:58
-
-
Save johnmyleswhite/5556201 to your computer and use it in GitHub Desktop.
Naively counting Pythagorean triples in Python and Julia
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
total = 0 | |
N = 300 | |
start_time = time() | |
for a in 0:(N - 1) | |
for b in 0:(N - 1) | |
for c in 0:(N - 1) | |
if a^2 + b^2 == c^2 | |
total = total + 1 | |
end | |
end | |
end | |
end | |
end_time = time() | |
end_time - start_time | |
# Repeat now that JIT has done its work | |
total = 0 | |
N = 300 | |
start_time = time() | |
for a in 0:(N - 1) | |
for b in 0:(N - 1) | |
for c in 0:(N - 1) | |
if a^2 + b^2 == c^2 | |
total = total + 1 | |
end | |
end | |
end | |
end | |
end_time = time() | |
println(end_time - start_time) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
function loop(N::Integer) | |
total = 0 | |
start_time = time() | |
for a in 0:(N - 1) | |
for b in 0:(N - 1) | |
for c in 0:(N - 1) | |
if a^2 + b^2 == c^2 | |
total = total + 1 | |
end | |
end | |
end | |
end | |
end_time = time() | |
return end_time - start_time | |
end | |
loop(300) | |
# Repeat now that JIT has done its work | |
println(loop(300)) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
total = 0 | |
N = 300 | |
from time import time | |
start = time() | |
for a in range(N): | |
for b in range(N): | |
for c in range(N): | |
if a**2 + b**2 == c**2: | |
total = total + 1 | |
end = time() | |
print(end - start) |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
If
N
is defined asconst
in JuliaGlobals, the time changes from 660 ms to 34 ms on my machine