def _run_parallel_or_opposite(self, sgn, gt):
n = 64
ps = 1 * q.um
y = np.linspace(0, sgn * n, num=10)
x = z = np.zeros(y.shape)
traj = Trajectory(zip(x, y, z) * ps, velocity=ps / q.s, pixel_size=ps)
mb = MetaBall(traj, n * ps / 16, orientation=geom.Y_AX)
mb.move(0 * q.s)
np.testing.assert_almost_equal(gt, mb.transform_matrix)
def create_metaballs(params, pixel_size):
x, y, z, r = list(zip(*params))
objects = b""
metaballs = []
for i in range(len(params)):
c_points = [(x[i], y[i], z[i])] * q.um
trajectory = Trajectory(c_points, pixel_size=pixel_size)
metaball = MetaBall(trajectory, r[i] * q.um)
metaball.move(0 * q.s)
metaballs.append(metaball)
objects += metaball.pack()
return metaballs, objects
def create_metaballs(params):
x, y, z, r = zip(*params)
objects = ""
metaballs = []
for i in range(len(params)):
c_points = [(x[i], y[i], z[i])] * q.um
trajectory = Trajectory(c_points)
metaball = MetaBall(trajectory, r[i] * q.um)
metaball.move(0 * q.s)
metaballs.append(metaball)
objects += metaball.pack()
return metaballs, objects
