근데.... 처음에 만들었던 파이썬 루틴보다 3배 가까이 시간이 걸립니다... (파이선 9분대, C 25분대)
이리저리 연구를 좀 해봤는데 파이썬에서는 '편하게' Numpy array 를 사용해서 벡터에 대고서 그냥 계산을 한 반면에, C 루틴에서는 모든 값들을 for loop 두개 돌려서 훑으며 계산한게 느린거 같습니다.
아마, for loop 두개 돌리는게 훨씬 느린거라 생각되는데 원래 이런건가요...?
1) Python
f = np.zeros(len(delnu))
f1 = np.zeros(len(delnu))
for j in range(len(hitran_data)):
f1 = inten[j] * ((1/np.pi)*(gamma[j]/(gamma[j]**2 + (delnu - (nu[j] + delar[j]*prs/1013))**2)))
f = f1 + f
참고 delnu 사이즈 : 1900000, hitran 사이즈 : 19615
2) C
for(i=0;i<lines;i++)
for(j=0;j<delnu_size;j++)
{
f[j] += inten[i] * ((1.0/pi) * (gamma[i] / (pow(gamma[i],2.0) + pow(delnu[j] - (nu[i] + del_air[i]*prs/1013.0),2.0))));
}
참고 delnu_size 사이즈 : 1900000, lines 사이즈 : 19615
혹시 제곱승 나타내는 pow() 함수가 느린걸까요?
잡담 | 1337명이 읽었어요. 3.143.0.134
