Example #1
0
def onecurvepoint(p0):
    tol=1e-4
    while True:
        df = polynomial_spline.df(p0)
        p1 = p0 - polynomial_spline.f(p0)*df/dot(df,df)
        dif = p1-p0
        if sqrt(dot(dif,dif))<tol:
            return p1
        else: p0=p1
Example #2
0
def onecurvepoint(p0):
    tol = 1e-4
    while True:
        df = polynomial_spline.df(p0)
        p1 = p0 - polynomial_spline.f(p0) * df / dot(df, df)
        dif = p1 - p0
        if sqrt(dot(dif, dif)) < tol:
            return p1
        else:
            p0 = p1
Example #3
0
def curvepoint(p0):
    tol=1e-4
    while True:
        df1 = polynomial_spline.df(p0); f1 = polynomial_spline.f(p0)
        df2 = array([2*p0[0],2*p0[1],2*p0[2]]); f2 = p0[0]**2+p0[1]**2+p0[2]**2-Rn**2
        df1df2 = dot(df1,df2); df1df1 = dot(df1,df1); df2df2 = dot(df2,df2)
        beta = (-f1+f2*df1df1/df1df2)/(df1df2-df1df1*df2df2/df1df2)
        alpha = (-f1+f2*df1df2/df2df2)/(df1df1-df1df2*df1df2/df2df2)
        dk = alpha*df1+beta*df2
        p1 = p0+dk
        dif = p1-p0
        if sqrt(dot(dif,dif))<tol:
            return p1
        else: p0=p1
Example #4
0
def curvepoint(p0):
    tol = 1e-4
    while True:
        df1 = polynomial_spline.df(p0)
        f1 = polynomial_spline.f(p0)
        df2 = array([2 * p0[0], 2 * p0[1], 2 * p0[2]])
        f2 = p0[0]**2 + p0[1]**2 + p0[2]**2 - Rn**2
        df1df2 = dot(df1, df2)
        df1df1 = dot(df1, df1)
        df2df2 = dot(df2, df2)
        beta = (-f1 + f2 * df1df1 / df1df2) / (df1df2 -
                                               df1df1 * df2df2 / df1df2)
        alpha = (-f1 + f2 * df1df2 / df2df2) / (df1df1 -
                                                df1df2 * df1df2 / df2df2)
        dk = alpha * df1 + beta * df2
        p1 = p0 + dk
        dif = p1 - p0
        if sqrt(dot(dif, dif)) < tol:
            return p1
        else:
            p0 = p1