def difcost(a,b):
dif=0
k = 0
for i in range(shape(a)[0]):
for j in range(shape(a)[1]):
# Euclidean Distance
k+=1
dif+=pow(a[i,j]-b[i,j],2)
print k,':' ,a[i,j] ,b[i,j] ,dif
if dif!=dif:
return dif
return dif
で実行してみました。結果はk=585で毎回止まるようです。メモリ系のバグでしょうか。
ubuntuもあるのでそちらでも試して見ます。

577 : 0 0.0121184608179 81.0386182152 577 : 0 0.0164147052137 91.5984888312 577 : 0 0.0292418437414 79.7176242734
578 : 0 0.0508183707462 81.041200722 578 : 0 0.0496066791031 91.6009496538 578 : 0 0.0733880158902 79.7230100743
579 : 0 0.0257723893872 81.041864938 579 : 0 0.0274640804479 91.6017039296 579 : 0 0.0441988852923 79.7249636157
580 : 0 0.0193929967783 81.0422410264 580 : 0 0.0304702331206 91.6026323647 580 : 0 0.0250509068254 79.7255911637
581 : 0 0.136879372348 81.0609769889 581 : 0 0.13911625497 91.6219856971 581 : 0 0.155482865404 79.7497660851
582 : 0 0.0558954001978 81.0641012847 582 : 0 0.0732634782209 91.6273532343 582 : 0 0.0667560238416 79.7542224518
583 : 0 0.0431843749592 81.0659661749 583 : 0 0.0296703213474 91.6282335623 583 : 0 0.0374581954731 79.7556255682
584 : 0 0.0525606269812 81.0687287944 584 : 0 0.0382979828685 91.6297002978 584 : 0 0.0666638480736 79.7600696369
585 : 0 nan nan 585 : 0 nan nan 585 : 0 nan nan
>>> >>> >>>