Beispiel #1
0
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)]
Beispiel #2
0
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)]