コード例 #1
0
def ns_implicit_mid():
    global x_nsimp
    x_nsimp = np.array([x0, xp0])
    traj = x_nsimp

    for i in range(1, steps + 1):
        x_new = fs(NSIMP, x_nsimp)
        traj = np.vstack([traj, x_new])
        x_nsimp = x_new

    return traj
コード例 #2
0
def getpoint_spline(dist, tck):
    x = fs(getlengthfromspline,0.0,(tck,dist))[0]
    return x,interpolate.splev(x,tck,der=0)
コード例 #3
0
def getpoint_spline(dist,tck):
    x = fs(getlengthfromspline,0.0,(tck,dist))[0] #fs is a solver. It solves for dist_spline  =
    #req_dist
    return x,interpolate.splev(x,tck,der=0).item()
コード例 #4
0
 def d_to_x_old(
         self, length, start_x
 ):  #Bad idea! The solver is calling splev more than required.
     x_l = fs(lambda x: self.x_to_d(start_x, x) - length, self.road_x[0])
     return x_l
コード例 #5
0
def calcu():
    ##########初始化坐标##########
    global R1
    R1 = float(e_CIRCLE_1.get())
    global R2
    R2 = float(e_CIRCLE_2.get())
    global R3
    R3 = float(e_CIRCLE_3.get())
    global R4
    R4 = float(e_CIRCLE_4.get())

    global CIRCLE_1
    CIRCLE_1 = Circle(0, 3, R1)
    global CIRCLE_4
    CIRCLE_4 = Circle(21, 43, R4)

    global POINT_A
    global POINT_P
    global POINT_B
    global POINT_C
    POINT_A = Point(float(e_POINT_A_X.get()), float(e_POINT_A_Y.get()))
    POINT_P = Point(float(e_POINT_P_X.get()), float(e_POINT_P_Y.get()))
    POINT_B = Point(float(e_POINT_B_X.get()), float(e_POINT_B_Y.get()))
    POINT_C = Point(float(e_POINT_C_X.get()), float(e_POINT_C_Y.get()))

    ##########得到M点的坐标##########
    global POINT_M
    POINT_M = fs(qiedian, [1, 1])
    POINT_M = Point(POINT_M[0], POINT_M[1])
    ##########求LINE1和LINE2##########
    global LINE_1  ## 全局变量必须设置在方法体内部
    LINE_1 = lineBy2Point(POINT_M, POINT_C, 1)
    LINE_1.sayHello(1)

    global LINE_2
    LINE_2 = lineBy2Point(POINT_C, POINT_B, 2)
    LINE_2.sayHello(2)
    ##########求CIRCLE2和CIRCLE3##########
    R2_cicle = fs(
        solve_circle_2,
        [1, 1],
    )
    global CIRCLE_2
    CIRCLE_2 = Circle(R2_cicle[0], R2_cicle[1], 15)
    CIRCLE_2.sayHello(2)

    R3_cicle = fs(solve_circle_3, [1, 1])
    global CIRCLE_3
    CIRCLE_3 = Circle(R3_cicle[0], R3_cicle[1], 3)
    CIRCLE_3.sayHello(3)
    ################# 交点 #################
    global POINT_E
    POINT_E = Point(0, 3)
    POINT_E.sayHello("E")
    POINT_M.sayHello("M")

    global POINT_N
    POINT_N = fs(line1ToCircle2, [1, 1])
    POINT_N = Point(POINT_N[0], POINT_N[1])
    POINT_N.sayHello("N")

    global POINT_L
    POINT_L = fs(line2ToCircle2, [1, 1])
    POINT_L = Point(POINT_L[0], POINT_L[1])
    POINT_L.sayHello("L")

    global POINT_Y
    POINT_Y = fs(line2ToCircle3, [1, 1])
    POINT_Y = Point(POINT_Y[0], POINT_Y[1])
    POINT_Y.sayHello("Y")

    global POINT_X
    POINT_X = fs(line3ToCircle3, [1, 1])
    POINT_X = Point(POINT_X[0], POINT_X[1])
    POINT_X.sayHello("X")

    ################# 临界面积 #################
    print "\n现在开始求积分:"
    SUIWEI = e_suiwei.get()
    area = client(float(SUIWEI))

    lable1 = Label(root, text="面积是:" + str(float(area)))
    lable1.pack()
コード例 #6
0
def getpoint_spline(dist, tck):
    x = fs(getlengthfromspline, 0.0, (tck, dist))[0]
    return x, interpolate.splev(x, tck, der=0)
コード例 #7
0
	def d_to_x_old(self,length,start_x): #Bad idea! The solver is calling splev more than required.
 		x_l = fs(lambda x: self.x_to_d(start_x,x) - length,self.road_x[0])
		return x_l