Я написал это очень плохо оптимизированный код C, который делает простую математическую расчет:оптимизация кода Python (20x медленнее, чем C)
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#define MIN(a, b) (((a) < (b)) ? (a) : (b))
#define MAX(a, b) (((a) > (b)) ? (a) : (b))
unsigned long long int p(int);
float fullCheck(int);
int main(int argc, char **argv){
int i, g, maxNumber;
unsigned long long int diff = 1000;
if(argc < 2){
fprintf(stderr, "Usage: %s maxNumber\n", argv[0]);
return 0;
}
maxNumber = atoi(argv[1]);
for(i = 1; i < maxNumber; i++){
for(g = 1; g < maxNumber; g++){
if(i == g)
continue;
if(p(MAX(i,g)) - p(MIN(i,g)) < diff && fullCheck(p(MAX(i,g)) - p(MIN(i,g))) && fullCheck(p(i) + p(g))){
diff = p(MAX(i,g)) - p(MIN(i,g));
printf("We have a couple %llu %llu with diff %llu\n", p(i), p(g), diff);
}
}
}
return 0;
}
float fullCheck(int number){
float check = (-1 + sqrt(1 + 24 * number))/-6;
float check2 = (-1 - sqrt(1 + 24 * number))/-6;
if(check/1.00 == (int)check)
return check;
if(check2/1.00 == (int)check2)
return check2;
return 0;
}
unsigned long long int p(int n){
return n * (3 * n - 1)/2;
}
И тогда я пытался (просто для удовольствия) портирование под Python посмотреть, как это будет реагировать. Моя первая версия была почти 1: 1 конверсией, которая работает ужасно медленно (120 + сек в Python против < 1сек в C). Я сделал немного оптимизации, и это то, что я получил:
#!/usr/bin/env/python
from cmath import sqrt
import cProfile
from pstats import Stats
def quickCheck(n):
partial_c = (sqrt(1 + 24 * (n)))/-6
c = 1/6 + partial_c
if int(c.real) == c.real:
return True
c = c - 2*partial_c
if int(c.real) == c.real:
return True
return False
def main():
maxNumber = 5000
diff = 1000
for i in range(1, maxNumber):
p_i = i * (3 * i - 1)/2
for g in range(i, maxNumber):
if i == g:
continue
p_g = g * (3 * g - 1)/2
if p_i > p_g:
ma = p_i
mi = p_g
else:
ma = p_g
mi = p_i
if ma - mi < diff and quickCheck(ma - mi):
if quickCheck(ma + mi):
print ('New couple ', ma, mi)
diff = ma - mi
cProfile.run('main()','script_perf')
perf = Stats('script_perf').sort_stats('time', 'calls').print_stats(10)
Это работает примерно 16secs, что лучше, но и почти в 20 раз медленнее, чем C. Теперь, я знаю, C лучше Python для такого рода вычислений, но то, что я хотел бы знать, - это то, что я пропустил (Python-мудрый, как ужасно медленная функция или такой), который мог бы сделать эту функцию быстрее. Обратите внимание, что я использую Python 3.1.1, если это имеет значение.
Пожалуйста, вы можете описать то, что ваша программа делает –
профайлер может замедлить вещи совсем немного –
я вижу, что, как представляется, много дополнительных управляющими операторами.Давным-давно давным-давно в галактике я заметил, что для такого рода алгоритмов, которые уменьшали количество управляющих утверждений, как можно больше, приводили меня к огромным повышениям производительности. – TheJacobTaylor