本文整理汇总了Python中bspline.Bspline类的典型用法代码示例。如果您正苦于以下问题:Python Bspline类的具体用法?Python Bspline怎么用?Python Bspline使用的例子?那么, 这里精选的类代码示例或许可以为您提供帮助。
在下文中一共展示了Bspline类的13个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: __init__
def __init__(self, l = 2.0, A = 2.0, m = -0.2, ncp = 6, type = 1, p = 3):
self._l = l
self._A = A
self._m = m
self._ncp = ncp
self._type = type
self._xp = self._spacing()
self._Pig = self._init_cp()
self._P = self._cp()
Bspline.__init__(self, self._P, p)
示例2: __init__
def __init__(self, cp_s, cp_p, p = 3, pitch = 0.0, height = 0.0 ):
self._p = p
self._pitch = pitch
self._height = height
if height == 0.0 and pitch == 0.0:
self._cp_s = cp_s
self._cp_p = cp_p
else:
self._cp_s = self._redefCP(cp_s)
self._cp_p = self._redefCP(cp_p)
self._prof_s = Bspline(self._cp_s, p = p)
self._prof_p = Bspline(self._cp_p, p = p)
示例3: __init__
def __init__(self,upper_points,lower_points,center_iine_controls,thickness_controls,name='shell'):
self.Po = upper_points
self.Pi = lower_points
self.Po_bar = upper_points.copy()
self.Pi_bar = lower_points.copy()
self.name = name
self.Cc = center_iine_controls
self.Ct = thickness_controls
self.bsc_o = Bspline(self.Cc,upper_points)
self.bsc_i = Bspline(self.Cc,lower_points)
self.bst_o = Bspline(self.Ct,upper_points)
self.bst_i = Bspline(self.Ct,lower_points)
self.r_mag = np.average(upper_points[:,1])
#for revolution of 2-d profile
self.n_theta = 20
self.Theta = np.linspace(0,2*np.pi,self.n_theta)
self.ones = np.ones(self.n_theta)
self.sin_theta = np.sin(self.Theta)
self.cos_theta = np.cos(self.Theta)
#calculate derivatives
#in polar coordinates
self.dPo_bar_xqdCc = self.bsc_o.B.flatten()
self.dPo_bar_rqdCc = self.r_mag*self.bsc_o.B.flatten()
self.dPi_bar_xqdCc = self.bsc_i.B.flatten()
self.dPi_bar_rqdCc = self.r_mag*self.bsc_i.B.flatten()
self.dPo_bar_rqdCt = self.r_mag*self.bst_o.B.flatten()
self.dPi_bar_rqdCt = -1*self.r_mag*self.bst_i.B.flatten()
#Project Polar derivatives into revolved cartisian coordinates
self.dXoqdCc = np.outer(self.dPo_bar_xqdCc,self.ones)
self.dYoqdCc = np.outer(self.dPo_bar_rqdCc,self.sin_theta)
self.dZoqdCc = np.outer(self.dPo_bar_rqdCc,self.cos_theta)
self.dXiqdCc = np.outer(self.dPi_bar_xqdCc,self.ones)
self.dYiqdCc = np.outer(self.dPi_bar_rqdCc,self.sin_theta)
self.dZiqdCc = np.outer(self.dPi_bar_rqdCc,self.cos_theta)
self.dYoqdCt = np.outer(self.dPo_bar_rqdCt,self.sin_theta)
self.dZoqdCt = np.outer(self.dPo_bar_rqdCt,self.cos_theta)
self.dYiqdCt = np.outer(self.dPi_bar_rqdCt,self.sin_theta)
self.dZiqdCt = np.outer(self.dPi_bar_rqdCt,self.cos_theta)
示例4: Corner
class Corner(object):
"""
we define the needle to be inserted into the divergent part of the nozzle
"""
def __init__(self, needleDef):
"""
Constructor
l-> lenght
A-> inlet Area
m-> inlet derivative
ncp-> number of Control Point
type -> stretching function
"""
self._needleDef = needleDef
self._P = needleDef.getCPc()
self._corner = Bspline(self._P, p = 2)
def __call__(self, u):
return self._corner(u)
def plot(self):
return self._corner.plot()
示例5: __init__
class NozzleComp:
def __init__(self, cdDef, p = 3):
self._cdDef = cdDef
self._P = cdDef.getCP()
self._p = p
self._convDiv = Bspline(self._P, p = p)
def __call__(self, u):
return self._convDiv(u)
def plot(self):
self._convDiv.plot()
def getCP(self):
return self._P
示例6: _err
def _err(self):
xy_p, s = self._xy_p, self._s
x = np.ones(len(s))
y = np.ones(len(s))
err = np.ones(len(s))
for i in xrange(len(s)):
x[i], y[i], nul1 = Bspline.__call__(self, s[i])
err[i] = abs(xy_p[0,i]- x[i]) + abs(xy_p[1,i]- y[i]) #it may be vectorized
# err = abs(xy_p[0,:])
return err
示例7: __init__
def __init__(self, needleDef):
"""
Constructor
l-> lenght
A-> inlet Area
m-> inlet derivative
ncp-> number of Control Point
type -> stretching function
"""
self._needleDef = needleDef
self._P = needleDef.getCPt()
self._tip = Bspline(self._P, p = 2)
示例8: range
## Escribo bien los nodos y los pesos para las integrales en r
r_nodos = np.array([]);
wr_pesos = np.array([]);
for i in range(N_intervalos_r+1):
aux_x = 0.5*(knots_r[i+grado+1]-knots_r[i+grado])*x + 0.5*(knots_r[i+grado+1]+knots_r[i+grado]);
aux_w = 0.5*(knots_r[i+grado+1]-knots_r[i+grado])*w;
r_nodos = np.hstack((r_nodos, aux_x));
wr_pesos = np.hstack((wr_pesos, aux_w));
w_pesos = wr_pesos
wr_pesos = np.tile(wr_pesos, (N_splines_r-1, 1));
r_nodos2 = np.tile(r_nodos, (N_splines_r-1, 1));
## B-splines en la coordenada r
basis = Bspline(knots_r, grado);
## Calculo la matriz con los bsplines y sus derivadas en r
bsr = [basis._Bspline__basis(i, basis.p) for i in r_nodos]; # evaluo los bsplines
dbsr = [basis.d(i) for i in r_nodos]; # evaluo las derivadas de los bsplines
bsr = np.array(bsr)
bsr = bsr[:,0:N_splines_r-1]
dbsr = np.array(dbsr)
dbsr = dbsr[:,0:N_splines_r-1]
## Matriz de solapamiento en r
Sr = np.dot(np.transpose(bsr), (np.transpose(r_nodos2)*np.transpose(wr_pesos)*bsr));
# Sr = np.array([[Sr[i][j] for i in range(N_splines_r-1)] for j in range(N_splines_r-1)]);
## Matriz de energia cinetica en r
示例9: __init__
def __init__(self, stl, controls, name="body", r_ref=None, x_ref=None):
"""stl must be an STL object"""
self.stl = stl
geom_points = stl.points
self.coords = Coordinates(geom_points, cartesian=True)
self.P = self.coords.cylindrical
self.P_cart = self.coords.cartesian
self.P_bar = geom_points.copy() # just initialization
if isinstance(controls, int):
X = geom_points[:, 0]
x_max = np.max(X)
x_min = np.min(X)
C_x = np.linspace(x_min, x_max, controls)
C_r = np.zeros((controls,))
control_points = np.array(zip(C_x, C_r))
self.C = control_points
self.n_controls = controls
else:
self.C = controls
self.n_controls = len(controls)
self.C_bar = self.C.copy()
self.delta_C = np.zeros(self.C.shape)
self.bs = Bspline(self.C, geom_points)
self.name = name
if x_ref is not None:
self.x_mag = float(x_ref)
else:
self.x_mag = 10 ** np.floor(np.log10(np.average(geom_points[:, 0])))
if r_ref is not None:
self.r_mag = float(r_ref)
else:
indecies = np.logical_and(abs(geom_points[:, 2]) < 2e-6, geom_points[:, 1] > 0)
points = geom_points[indecies]
self.r_mag = 10 ** np.floor(
np.log10(np.average(points[:, 1]))
) # grab the order of magnitude of the average
# for revolution of 2-d profile
# self.n_theta = 20
# sgrab the theta values from the points
self.Theta = self.P[:, 2]
# this is too complex. shouldn't need to tile, then flatten later.
self.sin_theta = np.tile(np.sin(self.Theta), (self.n_controls, 1)).T.flatten()
self.cos_theta = np.tile(np.cos(self.Theta), (self.n_controls, 1)).T.flatten()
# self.sin_theta = np.tile(np.sin(self.Theta),self.n_controls)
# self.cos_theta = np.tile(np.cos(self.Theta),self.n_controls)
# calculate derivatives
# in polar coordinates
self.dP_bar_xqdC = np.array(self.x_mag * self.bs.B.flatten())
self.dP_bar_rqdC = np.array(self.r_mag * self.bs.B.flatten())
# Project Polar derivatives into revolved cartisian coordinates
self.dXqdC = self.dP_bar_xqdC.reshape(-1, self.n_controls)
self.dYqdC = (self.dP_bar_rqdC * self.sin_theta).reshape(-1, self.n_controls)
self.dZqdC = (self.dP_bar_rqdC * self.cos_theta).reshape(-1, self.n_controls)
示例10: Body
class Body(object):
"""FFD class for solid bodies which only have one surface"""
def __init__(self, stl, controls, name="body", r_ref=None, x_ref=None):
"""stl must be an STL object"""
self.stl = stl
geom_points = stl.points
self.coords = Coordinates(geom_points, cartesian=True)
self.P = self.coords.cylindrical
self.P_cart = self.coords.cartesian
self.P_bar = geom_points.copy() # just initialization
if isinstance(controls, int):
X = geom_points[:, 0]
x_max = np.max(X)
x_min = np.min(X)
C_x = np.linspace(x_min, x_max, controls)
C_r = np.zeros((controls,))
control_points = np.array(zip(C_x, C_r))
self.C = control_points
self.n_controls = controls
else:
self.C = controls
self.n_controls = len(controls)
self.C_bar = self.C.copy()
self.delta_C = np.zeros(self.C.shape)
self.bs = Bspline(self.C, geom_points)
self.name = name
if x_ref is not None:
self.x_mag = float(x_ref)
else:
self.x_mag = 10 ** np.floor(np.log10(np.average(geom_points[:, 0])))
if r_ref is not None:
self.r_mag = float(r_ref)
else:
indecies = np.logical_and(abs(geom_points[:, 2]) < 2e-6, geom_points[:, 1] > 0)
points = geom_points[indecies]
self.r_mag = 10 ** np.floor(
np.log10(np.average(points[:, 1]))
) # grab the order of magnitude of the average
# for revolution of 2-d profile
# self.n_theta = 20
# sgrab the theta values from the points
self.Theta = self.P[:, 2]
# this is too complex. shouldn't need to tile, then flatten later.
self.sin_theta = np.tile(np.sin(self.Theta), (self.n_controls, 1)).T.flatten()
self.cos_theta = np.tile(np.cos(self.Theta), (self.n_controls, 1)).T.flatten()
# self.sin_theta = np.tile(np.sin(self.Theta),self.n_controls)
# self.cos_theta = np.tile(np.cos(self.Theta),self.n_controls)
# calculate derivatives
# in polar coordinates
self.dP_bar_xqdC = np.array(self.x_mag * self.bs.B.flatten())
self.dP_bar_rqdC = np.array(self.r_mag * self.bs.B.flatten())
# Project Polar derivatives into revolved cartisian coordinates
self.dXqdC = self.dP_bar_xqdC.reshape(-1, self.n_controls)
self.dYqdC = (self.dP_bar_rqdC * self.sin_theta).reshape(-1, self.n_controls)
self.dZqdC = (self.dP_bar_rqdC * self.cos_theta).reshape(-1, self.n_controls)
def copy(self):
return copy.deepcopy(self)
def deform(self, delta_C):
"""returns new point locations for the given motion of the control
points"""
self.delta_C = delta_C
self.delta_C[:, 0] = self.delta_C[:, 0] * self.x_mag
self.C_bar = self.C + self.delta_C
delta_P = self.bs.calc(self.C_bar)
self.P_bar = self.P.copy()
self.P_bar[:, 0] = delta_P[:, 0]
self.P_bar[:, 1] = self.P[:, 1] + self.r_mag * delta_P[:, 1]
# transform to cartesian coordinates
self.coords = Coordinates(self.P_bar, cartesian=False)
self.P_bar_cart = self.coords.cartesian
self.Xo = self.P_bar_cart[:, 0]
self.Yo = self.P_bar_cart[:, 1]
self.Zo = self.P_bar_cart[:, 2]
self.stl.update_points(self.P_bar_cart)
return self.P_bar
示例11: Shell
class Shell(object):
"""FFD class for shell bodies which have two connected surfaces"""
def __init__(
self, outer_stl, inner_stl, center_line_controls, thickness_controls, name="shell", r_ref=None, x_ref=None
):
self.outer_stl = outer_stl
self.inner_stl = inner_stl
outer_points = outer_stl.points
inner_points = inner_stl.points
self.n_outer = len(outer_points)
self.n_inner = len(inner_points)
self.outer_coords = Coordinates(outer_points, cartesian=True)
self.inner_coords = Coordinates(inner_points, cartesian=True)
self.Po = self.outer_coords.cylindrical
self.Pi = self.inner_coords.cylindrical
self.Po_cart = self.outer_coords.cartesian
self.Pi_cart = self.inner_coords.cartesian
# just initialization for array size
self.Po_bar = outer_points.copy()
self.Pi_bar = inner_points.copy()
self.name = name
if isinstance(center_line_controls, int):
X = outer_points[:, 0]
x_max = np.max(X)
x_min = np.min(X)
C_x = np.linspace(x_min, x_max, center_line_controls)
C_r = np.zeros((center_line_controls,))
control_points = np.array(zip(C_x, C_r))
self.Cc = control_points
self.n_c_controls = center_line_controls
else:
self.Cc = center_line_controls
self.n_c_controls = len(center_line_controls)
self.Cc_bar = self.Cc.copy()
self.delta_Cc = np.zeros(self.Cc.shape)
if isinstance(thickness_controls, int):
X = inner_points[:, 0]
x_max = np.max(X)
x_min = np.min(X)
C_x = np.linspace(x_min, x_max, thickness_controls)
C_r = np.zeros((thickness_controls,))
control_points = np.array(zip(C_x, C_r))
self.Ct = control_points
self.n_t_controls = thickness_controls
else:
self.Ct = thickness_controls
self.n_t_controls = len(thickness_controls)
self.Ct_bar = self.Ct.copy()
self.delta_Ct = np.zeros(self.Ct.shape)
self.bsc_o = Bspline(self.Cc, outer_points)
self.bsc_i = Bspline(self.Cc, inner_points)
self.bst_o = Bspline(self.Ct, outer_points)
self.bst_i = Bspline(self.Ct, inner_points)
self.name = name
if x_ref is not None:
self.x_mag = float(x_ref)
else:
self.x_mag = 10 ** np.floor(np.log10(np.average(outer_points[:, 0])))
if r_ref is not None:
self.r_mag = float(r_ref)
else:
indecies = np.logical_and(abs(outer_points[:, 2]) < 2e-6, outer_points[:, 1] > 0)
points = outer_points[indecies]
self.r_mag = 10 ** np.floor(
np.log10(np.average(points[:, 1]))
) # grab the order of magnitude of the average
self.outer_theta = self.Po[:, 2]
self.sin_outer_c_theta = np.tile(np.sin(self.outer_theta), (self.n_c_controls, 1)).T.flatten()
self.cos_outer_c_theta = np.tile(np.cos(self.outer_theta), (self.n_c_controls, 1)).T.flatten()
self.sin_outer_t_theta = np.tile(np.sin(self.outer_theta), (self.n_t_controls, 1)).T.flatten()
self.cos_outer_t_theta = np.tile(np.cos(self.outer_theta), (self.n_t_controls, 1)).T.flatten()
self.inner_theta = self.Pi[:, 2]
self.sin_inner_c_theta = np.tile(np.sin(self.inner_theta), (self.n_c_controls, 1)).T.flatten()
self.cos_inner_c_theta = np.tile(np.cos(self.inner_theta), (self.n_c_controls, 1)).T.flatten()
self.sin_inner_t_theta = np.tile(np.sin(self.inner_theta), (self.n_t_controls, 1)).T.flatten()
self.cos_inner_t_theta = np.tile(np.cos(self.inner_theta), (self.n_t_controls, 1)).T.flatten()
# calculate derivatives
# in polar coordinates
self.dPo_bar_xqdCc = np.array(self.x_mag * self.bsc_o.B.flatten())
self.dPo_bar_rqdCc = np.array(self.r_mag * self.bsc_o.B.flatten())
self.dPi_bar_xqdCc = np.array(self.x_mag * self.bsc_i.B.flatten())
self.dPi_bar_rqdCc = np.array(self.r_mag * self.bsc_i.B.flatten())
#.........这里部分代码省略.........
示例12: Body
class Body(object):
"""FFD class for solid bodyies which only have one surface"""
def __init__(self,geom_points,control_points,name="body"):
self.P = geom_points
self.P_bar = geom_points.copy()
self.C = control_points
self.bs = Bspline(control_points,geom_points)
self.name = name
self.r_mag = np.average(geom_points[:,1])
#for revolution of 2-d profile
self.n_theta = 20
self.Theta = np.linspace(0,2*np.pi,self.n_theta)
self.ones = np.ones(self.n_theta)
self.sin_theta = np.sin(self.Theta)
self.cos_theta = np.cos(self.Theta)
#calculate derivatives
#in polar coordinates
self.dP_bar_xqdC = self.bs.B.flatten()
self.dP_bar_rqdC = self.r_mag*self.bs.B.flatten()
#Project Polar derivatives into revolved cartisian coordinates
self.dXqdC = np.outer(self.dP_bar_xqdC,self.ones)
self.dYqdC = np.outer(self.dP_bar_rqdC,self.sin_theta)
self.dZqdC = np.outer(self.dP_bar_rqdC,self.cos_theta)
def deform(self,delta_C):
"""returns new point locations for the given motion of the control
points"""
self.C_bar = self.C+delta_C
delta_P = self.bs.calc(self.C_bar)
self.P_bar = self.P.copy()
self.P_bar[:,0] = delta_P[:,0]
self.P_bar[:,1] = self.P[:,1]+self.r_mag*delta_P[:,1]
#Perform axial roation of 2-d polar coordiantes
P_bar_r = self.P_bar[:,1]
self.Xo = np.outer(self.P_bar[:,0],self.ones)
self.Yo = np.outer(P_bar_r,self.sin_theta)
self.Zo = np.outer(P_bar_r,self.cos_theta)
return self.P_bar
def plot_spline(self,ax,point_color='r',line_color='b'):
map_points = self.bs(np.linspace(0,1,100))
ax.plot(map_points[:,0],map_points[:,1],c=line_color,label="%s b-spline"%self.name)
ax.scatter(self.C_bar[:,0],self.C_bar[:,1],c=point_color,label="%s control points"%self.name,s=50)
def plot_geom(self,ax,initial_color='g',ffd_color='k'):
if initial_color:
ax.scatter(self.P[:,0],self.P[:,1],c=initial_color,s=50,label="%s initial geom"%self.name)
ax.plot(self.P[:,0],self.P[:,1],c=initial_color)
if ffd_color:
ax.scatter(self.P_bar[:,0],self.P_bar[:,1],c=ffd_color,s=50,label="%s ffd geom"%self.name)
ax.plot(self.P_bar[:,0],self.P_bar[:,1],c=ffd_color)
示例13: Shell
class Shell(object):
"""FFD class for shell bodies which have two connected surfaces"""
def __init__(self,upper_points,lower_points,center_iine_controls,thickness_controls,name='shell'):
self.Po = upper_points
self.Pi = lower_points
self.Po_bar = upper_points.copy()
self.Pi_bar = lower_points.copy()
self.name = name
self.Cc = center_iine_controls
self.Ct = thickness_controls
self.bsc_o = Bspline(self.Cc,upper_points)
self.bsc_i = Bspline(self.Cc,lower_points)
self.bst_o = Bspline(self.Ct,upper_points)
self.bst_i = Bspline(self.Ct,lower_points)
self.r_mag = np.average(upper_points[:,1])
#for revolution of 2-d profile
self.n_theta = 20
self.Theta = np.linspace(0,2*np.pi,self.n_theta)
self.ones = np.ones(self.n_theta)
self.sin_theta = np.sin(self.Theta)
self.cos_theta = np.cos(self.Theta)
#calculate derivatives
#in polar coordinates
self.dPo_bar_xqdCc = self.bsc_o.B.flatten()
self.dPo_bar_rqdCc = self.r_mag*self.bsc_o.B.flatten()
self.dPi_bar_xqdCc = self.bsc_i.B.flatten()
self.dPi_bar_rqdCc = self.r_mag*self.bsc_i.B.flatten()
self.dPo_bar_rqdCt = self.r_mag*self.bst_o.B.flatten()
self.dPi_bar_rqdCt = -1*self.r_mag*self.bst_i.B.flatten()
#Project Polar derivatives into revolved cartisian coordinates
self.dXoqdCc = np.outer(self.dPo_bar_xqdCc,self.ones)
self.dYoqdCc = np.outer(self.dPo_bar_rqdCc,self.sin_theta)
self.dZoqdCc = np.outer(self.dPo_bar_rqdCc,self.cos_theta)
self.dXiqdCc = np.outer(self.dPi_bar_xqdCc,self.ones)
self.dYiqdCc = np.outer(self.dPi_bar_rqdCc,self.sin_theta)
self.dZiqdCc = np.outer(self.dPi_bar_rqdCc,self.cos_theta)
self.dYoqdCt = np.outer(self.dPo_bar_rqdCt,self.sin_theta)
self.dZoqdCt = np.outer(self.dPo_bar_rqdCt,self.cos_theta)
self.dYiqdCt = np.outer(self.dPi_bar_rqdCt,self.sin_theta)
self.dZiqdCt = np.outer(self.dPi_bar_rqdCt,self.cos_theta)
def plot_geom(self,ax,initial_color='g',ffd_color='k'):
if initial_color:
ax.scatter(self.Po[:,0],self.Po[:,1],c=initial_color,s=50,label="%s initial geom"%self.name)
ax.scatter(self.Pi[:,0],self.Pi[:,1],c=initial_color,s=50)
ax.plot(self.Po[:,0],self.Po[:,1],c=initial_color)
ax.plot(self.Pi[:,0],self.Pi[:,1],c=initial_color)
if ffd_color:
ax.scatter(self.Po_bar[:,0],self.Po_bar[:,1],c=ffd_color,s=50,label="%s ffd geom"%self.name)
ax.scatter(self.Pi_bar[:,0],self.Pi_bar[:,1],c=ffd_color,s=50)
ax.plot(self.Po_bar[:,0],self.Po_bar[:,1],c=ffd_color)
ax.plot(self.Pi_bar[:,0],self.Pi_bar[:,1],c=ffd_color)
def plot_centerline_spline(self,ax,point_color='r',line_color='b'):
ax.scatter(self.Cc_bar[:,0],self.Cc_bar[:,1],c=point_color,s=50,label="%s Centerline Control Points"%self.name)
map_points = self.bsc_o(np.linspace(0,1,100))
ax.plot(map_points[:,0],map_points[:,1],label="Centerline b-spline Curve",c=line_color)
def plot_thickness_spline(self,ax,point_color='r',line_color='b'):
ax.scatter(self.Ct_bar[:,0],self.Ct_bar[:,1],c=point_color,s=50,label="%s Thickness Control Points"%self.name)
map_points = self.bsc_o(np.linspace(0,1,100))
ax.plot(map_points[:,0],map_points[:,1],label="Thickness b-spline Curve",c=line_color)
def deform(self,delta_Cc,delta_Ct):
"""returns new point locations for the given motion of the control
points for center-line and thickness"""
self.Cc_bar = self.Cc+delta_Cc
delta_Pc_o = self.bsc_o.calc(self.Cc_bar)
delta_Pc_i = self.bsc_i.calc(self.Cc_bar)
self.Ct_bar = self.Ct+delta_Ct
delta_Pt_o = self.bst_o.calc(self.Ct_bar)
delta_Pt_i = self.bst_i.calc(self.Ct_bar)
self.Po_bar = self.Po.copy()
self.Pi_bar = self.Pi.copy()
self.Po_bar[:,0] = delta_Pc_o[:,0]
self.Po_bar[:,1] = self.Po[:,1]+self.r_mag*(delta_Pc_o[:,1]+delta_Pt_o[:,1])
self.Pi_bar[:,0] = delta_Pc_i[:,0]
self.Pi_bar[:,1] = self.Pi[:,1]+self.r_mag*(delta_Pc_i[:,1]-delta_Pt_i[:,1])
#.........这里部分代码省略.........