def getPara2(self):
runout = self.runout
r = self.r
para = [[]] * len(r)
for i in range(int(len(r) / 4)):
vector_down = pre.pfit(pre.gcfl(r[i * 4], runout[i * 4], 0))
vector_up = pre.pfit(pre.gcfl(r[i * 4 + 1], runout[i * 4 + 1], 0))
args_1 = list(zip([vector_down] * 36,
[-i * 10
for i in range(36)])) # 需对相位取负结果才和之前一样，原因未知？？
args_2 = list(zip([vector_up] * 36, [-i * 10 for i in range(36)]))
para[i * 4] = [i for i in map(self.spin_vec, args_1)]
para[i * 4 + 1] = [i for i in map(self.spin_vec, args_2)]
center_down = pre.circ(pre.gcra(r[i * 4 + 2], runout[i * 4 + 2],
0))
center_up = pre.circ(pre.gcra(r[i * 4 + 3], runout[i * 4 + 3], 0))
args_3 = list(zip([center_down] * 36, [i * 10 for i in range(36)]))
args_4 = list(zip([center_up] * 36, [i * 10 for i in range(36)]))
para[i * 4 + 2] = [i for i in map(self.spin_cen, args_3)]
para[i * 4 + 3] = [i for i in map(self.spin_cen, args_4)]
return para
def vec(self, args):
runout = args[0]
r = args[1]
theta = args[2]
res = pre.pfit(pre.gcfl(r, pre.spin(runout, theta), 0))
return res
ad_part1_runout = 'G:/190708-190712data/MNJ-HPC-002ZP1-20190709.csv'
ad_part2_runout = 'G:/190708-190712data/MNJ-HPT-001ZP1-20190708.csv'
ad_inp_data = "G:/19.12.25whole model (less grids)/Step2----inp_all_node_data.csv"
ad_part1_flat_node_down = "G:/19.12.25whole model (less grids)/Step3----part1_flat_node_data_down.csv"
ad_part1_flat_node_up = "G:/19.12.25whole model (less grids)/Step3----part1_flat_node_data_up.csv"
ad_part2_flat_node_down = "G:/19.12.25whole model (less grids)/Step3----part2_flat_node_data_down.csv"
ad_part2_flat_node_up = "G:/19.12.25whole model (less grids)/Step3----part2_flat_node_data_up.csv"
ad_part1_radial_node_down = "G:/19.12.25whole model (less grids)/Step3----part1_radial_node_data_down.csv"
ad_part1_radial_node_up = "G:/19.12.25whole model (less grids)/Step3----part1_radial_node_data_up.csv"
ad_part2_radial_node_down = "G:/19.12.25whole model (less grids)/Step3----part2_radial_node_data_down.csv"
ad_part2_radial_node_up = "G:/19.12.25whole model (less grids)/Step3----part2_radial_node_data_up.csv"
ad_part2_node_number = "G:/19.12.25whole model (less grids)/Step4----part2_node_number.csv"
"""============================================================flat============================================================"""
#part1frure2
part1_flat_runout_down = -pre.idata(ad_part1_runout, [3])
part1_flat_runout_coordinate_down = pre.gcfl(80 / 2, part1_flat_runout_down, 0)
part1_normal_vector_down = pre.planefit(part1_flat_runout_coordinate_down)
part1_flat_node_data_down = pre.idata(ad_part1_flat_node_down, [0, 1, 2, 3])
new_part1_flat_node_data_down = rep.frure(80 / 2, part1_normal_vector_down[0],
part1_normal_vector_down[1],
part1_normal_vector_down[2],
part1_normal_vector_down[3],
part1_flat_runout_down,
part1_flat_node_data_down, -H1)
part1_flat_runout_up = pre.idata(ad_part1_runout, [1])
part1_flat_runout_coordinate_up = pre.gcfl(206.5 / 2, part1_flat_runout_up, 0)
part1_normal_vector_up = pre.planefit(part1_flat_runout_coordinate_up)
part1_flat_node_data_up = pre.idata(ad_part1_flat_node_up, [0, 1, 2, 3])
new_part1_flat_node_data_up = rep.frure(206.5 / 2, part1_normal_vector_up[0],
part1_normal_vector_up[1],
