def searchVectorStep(x, r, b): pdf = partD(x, r) n = [-pdf[i] for i in range(len(x))] n.append(1) normN = norm(n) if(normN == 1): normN = normN + 0.001 n[-1] = normN ** 2 - 1 r = ((b - fkt(x, r)) * normN ** 2) / (1 - normN ** 2) d = n / (normN ** 2) return x + r * d[:len(x)]
def searchVectorStep(coords, radi, b): pdf = partD(coords, radi) n = [-pdf[i] for i in range(len(coords))] n.append(1) normN = norm(n) if(normN == 1): normN = normN + 0.001 n[-1] = normN ** 2 - 1 radi = ((b - fkt(coords, radi)) * normN ** 2) / (1 - normN ** 2) d = n / (normN ** 2) return coords + radi * d[:len(coords)]