我需要使用Paraview创建一个可编程过滤器。 我们的想法是创建一个名为Speed的向量,其等于非旋转部分中的速度等于旋转速度+旋转速度。
问题在于我无法接受每个单元格中的速度值。
input0 = inputs[0]
radius=3
Speed1=input0.PointData["U"]
K=vtk.vtkDoubleArray()
X=input0.PointData["X"]
Y=input0.PointData["Y"]
Z=input0.PointData["Z"]
pdi = self.GetInput()
numPts = pdi.GetNumberOfPoints()
for i in range(0, numPts):
if X.getvalue(i)^2+Y.getvalue(i)^2<radius:
temp=U.getvalue(i)
else:
temp=U.getvalue(i)+rot
Speed.InsertNextValue(1)
output.PointData.append(Speed, "Speed")
问题是X.getvalue(i)
无效。