import itertools
import numpy as np
from .internal import Shape, ShapeDecorator, ShapeAttribute
[docs]@ShapeDecorator
class SphereUnions(Shape):
"""A collection of identical sphere-union bodies in 3D.
A `SphereUnions` object is a union of spheres, each of which has
its own color, radius, and local position. The `SphereUnions`
object can be rigidly rotated and translated via its position and
orientation attributes.
"""
_ATTRIBUTES = list(itertools.starmap(ShapeAttribute, [
('positions', np.float32, (0, 0, 0), 2, True,
'Position of each particle'),
('orientations', np.float32, (1, 0, 0, 0), 2, True,
'Orientation quaternion of each particle'),
('colors', np.float32, (.5, .5, .5, 1), 2, False,
'Color, RGBA, [0, 1] for each sphere in the union'),
('points', np.float32, (0, 0, 0), 2, False,
'Positions in local coordinates for the spheres in the union, to be replicated for each particle'),
('radii', np.float32, .5, 1, True,
'Radius of each sphere in the union')
]))
@property
def diameters(self):
"""Diameter of each particle."""
return 2*self._attributes.get('radii', .5)
@diameters.setter
def diameters(self, value):
self.radii = 0.5*np.atleast_1d(value)