|
| 1 | +import pytest |
| 2 | + |
| 3 | +import numpy as np |
| 4 | + |
| 5 | +from ipygany import PolyMesh, Data, Component |
| 6 | + |
| 7 | + |
| 8 | +vertices = np.array([ |
| 9 | + [0., 0., 0.], |
| 10 | + [0., 0., 0.], |
| 11 | + [0., 0., 0.], |
| 12 | +]) |
| 13 | + |
| 14 | +triangles = np.array([ |
| 15 | + [0, 1, 2], |
| 16 | +]) |
| 17 | + |
| 18 | +data_1d = Data(name='1d', components=[Component('x', np.array([0., 0., 0.]))]) |
| 19 | +data_3d = Data('3d', [ |
| 20 | + Component(name='x', array=np.array([1., 1., 1.])), |
| 21 | + Component('y', np.array([2., 2., 2.])), |
| 22 | + Component('z', np.array([3., 3., 3.])), |
| 23 | +]) |
| 24 | + |
| 25 | + |
| 26 | +def test_data_access(): |
| 27 | + poly = PolyMesh(vertices=vertices, triangle_indices=triangles, data=[data_1d, data_3d]) |
| 28 | + |
| 29 | + assert np.all(np.equal(poly['1d', 'x'].array, np.array([0., 0., 0.]))) |
| 30 | + assert np.all(np.equal(poly['3d', 'y'].array, np.array([2., 2., 2.]))) |
| 31 | + |
| 32 | + |
| 33 | +def test_mesh_data_creation(): |
| 34 | + poly = PolyMesh(vertices=vertices, triangle_indices=triangles, data={ |
| 35 | + '1d': [Component('x', np.array([0., 0., 0.]))], |
| 36 | + '2d': [Component('x', np.array([1., 1., 1.])), Component('y', np.array([2., 2., 2.]))] |
| 37 | + }) |
| 38 | + |
| 39 | + assert np.all(np.equal(poly['1d', 'x'].array, np.array([0., 0., 0.]))) |
| 40 | + assert np.all(np.equal(poly['2d', 'x'].array, np.array([1., 1., 1.]))) |
| 41 | + assert np.all(np.equal(poly['2d', 'y'].array, np.array([2., 2., 2.]))) |
| 42 | + |
| 43 | + poly = PolyMesh(vertices=vertices, triangle_indices=triangles, data={ |
| 44 | + '1d': {'x': np.array([0., 0., 0.])}, |
| 45 | + '2d': {'x': np.array([1., 1., 1.]), 'y': np.array([2., 2., 2.])} |
| 46 | + }) |
| 47 | + |
| 48 | + assert np.all(np.equal(poly['1d', 'x'].array, np.array([0., 0., 0.]))) |
| 49 | + assert np.all(np.equal(poly['2d', 'x'].array, np.array([1., 1., 1.]))) |
| 50 | + assert np.all(np.equal(poly['2d', 'y'].array, np.array([2., 2., 2.]))) |
| 51 | + |
| 52 | + |
| 53 | +def test_component_creation(): |
| 54 | + comp = Component('z', np.array([1., 2., 3.])) |
| 55 | + |
| 56 | + assert comp.name == 'z' |
| 57 | + assert comp.min == 1. |
| 58 | + assert comp.max == 3. |
0 commit comments