Почему эта функция Python прерывается в цикле?

Question

Из этого вопроса

Нахождение лучшей точки компромисса на кривой

Я использую решение Python, и оно прекрасно работает, если я запускаю его сам.

Но у меня есть словарь этих «кривых»:

sse {}

{'00': [4.998436383771836,
  0.0160165895672664,
  0.004512422186993107,
  0.0013171501024742112,
  0.000788783358847752,
  0.0005498425886621068],
'67':[0.13598930783101504,
  0.04717783658889547,
  0.027547125931827038,
  0.021440839841617088,
  0.016775671441391703,
  0.013185864754748117,
  0.010318462898609907],

и так далее ...

, если язапустите это, оно работает отлично:

curve = sse['67']
nPoints = len(curve)
allCoord = np.vstack((range(nPoints), curve)).T
np.array([range(nPoints), curve])
firstPoint = allCoord[0]
lineVec = allCoord[-1] - allCoord[0]
lineVecNorm = lineVec / np.sqrt(np.sum(lineVec**2))
vecFromFirst = allCoord - firstPoint
scalarProduct = np.sum(vecFromFirst * np.matlib.repmat(lineVecNorm, nPoints, 1), axis=1)
vecFromFirstParallel = np.outer(scalarProduct, lineVecNorm)
vecToLine = vecFromFirst - vecFromFirstParallel
distToLine = np.sqrt(np.sum(vecToLine ** 2, axis=1))
idxOfBestPoint = np.argmax(distToLine)

print(idxOfBestPoint)

но когда я пытаюсь сделать это зацикливающей функцией, я получаю странные значения, которые не совпадают, когда запускаются сами по себе:

for k in sse:
    curve = sse[str(k)]
    nPoints = len(curve)
    allCoord = np.vstack((range(nPoints), curve)).T
    np.array([range(nPoints), curve])
    firstPoint = allCoord[0]
    lineVeceVeceVec = allCoord[-1] - allCoord[0]
    lineVecNorm = lineVec / np.sqrt(np.sum(lineVec**2))
    vecFromFirst = allCoord - firstPoint
    scalarProduct = np.sum(vecFromFirst * np.matlib.repmat(lineVecNorm, nPoints, 1), axis=1)
    vecFromFirstParallel = np.outer(scalarProduct, lineVecNorm)
    vecToLine = vecFromFirst - vecFromFirstParallel
    distToLine = np.sqrt(np.sum(vecToLine ** 2, axis=1))
    idxOfBestPoint = np.argmax(distToLine)
    print(idxOfBestPoint, k)

Это будет плеватьчто-то вроде этого:

7 57
98 11
6 45
4 50
98 91
98 00
1 62
98 79
7 48
98 12
98 38
98 23
5 37
98 56
98 67
5 25
7 46
98 22
98 49
2 47
98 41
98 78
98 35
98 68
98 14
98 24
1 0
98 42

Я не могу сказать, есть ли какая-то переменная, которая не сбрасывается, или что могло бы вызвать ее сбой, добавив простой цикл?

Чтобы было понятно, это делаетзапустить, но calc выдает «98» в качестве колена в цикле, но при отдельном запуске будет 7 для списка «67», а не 98.

Answer 1

У вас есть опечатка: lineVeceVeceVec = allCoord[-1] - allCoord[0]. Это должно быть lineVec = allCoord[-1] - allCoord[0]. То, что вы выплевывали, было в последний раз, когда lineVec был назначен результату, и это последнее значение, вероятно, было от того, когда вы запустили не зацикленную версию.