def test_with_invalid_cuds(self):
# given
cuds = object()
# when/then
with self.assertRaises(TypeError):
cuds2vtk(cuds)
def test_source_from_an_orthorhombic_lattice(self):
# given
lattice = make_orthorhombic_lattice(
'test', (0.5, 0.54, 0.58), (15, 25, 35), (7, 9, 8))
self.add_velocity(lattice)
# when
data_set = cuds2vtk(cuds=lattice)
# then
self.assertEqual(data_set.GetNumberOfPoints(), 15 * 25 * 35)
assert_array_equal(data_set.GetOrigin(), (7.0, 9.0, 8.0))
point_data = data_set.GetPointData()
arrays = {
point_data.GetArray(index).GetName():
vtk_to_numpy(point_data.GetArray(index))
for index in range(point_data.GetNumberOfArrays())}
for node in lattice.iter(item_type=CUBA.NODE):
point_id = data_set.ComputePointId(node.index)
assert_array_equal(
lattice.get_coordinate(node.index),
data_set.GetPoint(point_id))
for key, value in node.data.iteritems():
assert_array_equal(arrays[key.name][point_id], value)
def test_source_from_a_xy_plane_rectangular_lattice(self):
# given
lattice = make_orthorhombic_lattice(
'test', (0.3, 0.35, 0.4), (13, 23, 1), origin=(0.2, -2.7, 0.0))
self.add_velocity(lattice)
# when
data_set = cuds2vtk(cuds=lattice)
# then
self.assertEqual(data_set.GetNumberOfPoints(), 13 * 23)
assert_array_equal(data_set.GetOrigin(), (0.2, -2.7, 0.0))
point_data = data_set.GetPointData()
arrays = {
point_data.GetArray(index).GetName():
vtk_to_numpy(point_data.GetArray(index))
for index in range(point_data.GetNumberOfArrays())}
for node in lattice.iter(item_type=CUBA.NODE):
point_id = data_set.ComputePointId(node.index)
assert_array_equal(
lattice.get_coordinate(node.index),
data_set.GetPoint(point_id))
for key, value in node.data.iteritems():
assert_array_equal(arrays[key.name][point_id], value)
def test_with_empty_cuds_particles(self):
# given
cuds = Particles('test')
# when
data_set = cuds2vtk(cuds)
# then
self.assertEqual(data_set.GetNumberOfPoints(), 0)
self.assertEqual(data_set.GetNumberOfCells(), 0)
def test_source_from_a_xy_plane_hexagonal_lattice(self):
# given
lattice = make_hexagonal_lattice('test', 0.1, 0.2, (5, 4, 1))
self.add_velocity(lattice)
# when
data_set = cuds2vtk(cuds=lattice)
# then
self.assertEqual(data_set.GetNumberOfPoints(), 5 * 4)
points = vtk_to_numpy(data_set.GetPoints().GetData())
for node in lattice.iter(item_type=CUBA.NODE):
position = lattice.get_coordinate(node.index)
point_id = data_set.FindPoint(position)
assert_array_equal(
points[point_id], numpy.asarray(position, dtype=points.dtype))
def test_source_from_a_cubic_lattice(self):
# given
lattice = make_cubic_lattice('test', 0.4, (14, 24, 34), (4, 5, 6))
self.add_velocity(lattice)
# when
data_set = cuds2vtk(cuds=lattice)
# then
self.assertEqual(data_set.GetNumberOfPoints(), 14 * 24 * 34)
assert_array_equal(data_set.GetOrigin(), (4.0, 5.0, 6.0))
point_data = data_set.GetPointData()
arrays = {
point_data.GetArray(index).GetName():
vtk_to_numpy(point_data.GetArray(index))
for index in range(point_data.GetNumberOfArrays())}
for node in lattice.iter(item_type=CUBA.NODE):
point_id = data_set.ComputePointId(node.index)
assert_array_equal(
lattice.get_coordinate(node.index),
data_set.GetPoint(point_id))
for key, value in node.data.iteritems():
assert_array_equal(arrays[key.name][point_id], value)
def test_with_cuds_particles(self):
# given
points = [
[0.0, 0.0, 0.0], [1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]]
bonds = [[0, 1], [0, 3], [1, 3, 2]]
point_temperature = [10., 20., 30., 40.]
bond_temperature = [60., 80., 190.]
cuds = Particles('test')
particle_uids = cuds.add(
Particle(
coordinates=point,
data=DataContainer(
TEMPERATURE=point_temperature[index],
MASS=index))
for index, point in enumerate(points))
cuds.add(
Bond(
particles=[particle_uids[uid] for uid in indices],
data=DataContainer(
TEMPERATURE=bond_temperature[index],
MASS=index))
for index, indices in enumerate(bonds))
# when
data_set = cuds2vtk(cuds)
# then check points
self.assertEqual(data_set.GetNumberOfPoints(), 4)
coordinates = [
particle.coordinates for particle in cuds.iter(
item_type=CUBA.PARTICLE)]
vtk_points = [
data_set.GetPoint(index)
for index in range(len(particle_uids))]
self.assertItemsEqual(vtk_points, coordinates)
point_data = data_set.GetPointData()
self.assertEqual(point_data.GetNumberOfArrays(), 2)
arrays = {
point_data.GetArray(index).GetName(): point_data.GetArray(index)
for index in range(2)}
self.assertItemsEqual(arrays.keys(), ['MASS', 'TEMPERATURE'])
mass = vtk_to_numpy(arrays['MASS'])
self.assertItemsEqual(mass, range(4))
assert_array_equal(
vtk_to_numpy(arrays['TEMPERATURE']),
[point_temperature[int(index)] for index in mass])
# then check bonds
self.assertEqual(data_set.GetNumberOfCells(), 3)
links = [[
particle.coordinates
for particle in cuds.iter(bond.particles)]
for bond in cuds.iter(item_type=CUBA.BOND)]
vtk_lines = data_set.GetLines()
lines = [[
data_set.GetPoint(index) for index in line]
for line in iter_cells(vtk_lines)]
self.assertItemsEqual(lines, links)
cell_data = data_set.GetCellData()
self.assertEqual(cell_data.GetNumberOfArrays(), 2)
arrays = {
cell_data.GetArray(index).GetName(): cell_data.GetArray(index)
for index in range(2)}
self.assertItemsEqual(arrays.keys(), ['MASS', 'TEMPERATURE'])
mass = vtk_to_numpy(arrays['MASS'])
self.assertItemsEqual(mass, range(3))
assert_array_equal(
vtk_to_numpy(arrays['TEMPERATURE']),
[bond_temperature[int(index)] for index in mass])
def test_with_cuds_mesh(self):
# given
points = numpy.array([
[0, 0, 0], [1, 0, 0], [0, 1, 0], [0, 0, 1],
[2, 0, 0], [3, 0, 0], [3, 1, 0], [2, 1, 0],
[2, 0, 1], [3, 0, 1], [3, 1, 1], [2, 1, 1]],
'f')
cells = [
[0, 1, 2, 3], # tetra
[4, 5, 6, 7, 8, 9, 10, 11]] # hex
faces = [[2, 7, 11]]
edges = [[1, 4], [3, 8]]
count = itertools.count()
points = [
Point(coordinates=point, data=DataContainer(TEMPERATURE=index))
for index, point in enumerate(points)]
cuds = Mesh('test')
cuds.add(points)
faces = [
Face(
points=[points[index].uid for index in face],
data=DataContainer(TEMPERATURE=next(count)))
for face in faces]
edges = [
Edge(
points=[points[index].uid for index in edge],
data=DataContainer(TEMPERATURE=next(count)))
for edge in edges]
cells = [
Cell(
points=[points[index].uid for index in cell],
data=DataContainer(TEMPERATURE=next(count)))
for cell in cells]
cuds.add(edges)
cuds.add(faces)
cuds.add(cells)
# when
data_set = cuds2vtk(cuds=cuds)
# then check points
self.assertEqual(data_set.GetNumberOfPoints(), 12)
point_data = data_set.GetPointData()
self.assertEqual(point_data.GetNumberOfArrays(), 1)
temperature = point_data.GetArray(0)
self.assertEqual(temperature.GetName(), 'TEMPERATURE')
self.assertItemsEqual(
vtk_to_numpy(temperature), range(12))
# then check cells
self.assertEqual(data_set.GetNumberOfCells(), 5)
cell_data = data_set.GetCellData()
self.assertEqual(cell_data.GetNumberOfArrays(), 1)
temperature = cell_data.GetArray(0)
self.assertEqual(temperature.GetName(), 'TEMPERATURE')
self.assertItemsEqual(
vtk_to_numpy(temperature), range(5))
# For each cell in the container
# find the corresponding cell in the vtkCellArray and
# verify that they link to points that have the right coordinates.
for cell in itertools.chain(
cuds.iter(item_type=CUBA.EDGE),
cuds.iter(item_type=CUBA.FACE),
cuds.iter(item_type=CUBA.CELL)):
# The temperature value is also the index that the cells
# are expected to have in the vtkCellArray.
value = cell.data[CUBA.TEMPERATURE]
index = numpy.nonzero(vtk_to_numpy(temperature) == value)[0]
vtk_cell = data_set.GetCell(index)
ids = vtk_cell.GetPointIds()
for i, uid in enumerate(cell.points):
cell_point = cuds.get(uid)
assert_array_equal(
data_set.GetPoint(ids.GetId(i)), cell_point.coordinates)