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