def test_closed_surface_semi_detached(self): orbital_position_container = self.build_system('semi-detached', up.radians(10)) self.assertTrue( testutils.surface_closed( faces=orbital_position_container.primary.faces, points=orbital_position_container.primary.points))
def test_closed_surface_over_contact(self): s = testutils.prepare_binary_system( testutils.BINARY_SYSTEM_PARAMS['over-contact'], spots_primary=testutils.SPOTS_META["primary"], spots_secondary=testutils.SPOTS_META["secondary"]) s.primary.discretization_factor = up.radians(7) s.init() orbital_position_container = testutils.prepare_orbital_position_container( s) orbital_position_container.build_mesh(components_distance=1.0) orbital_position_container.build_faces(components_distance=1.0) for component in ['primary', 'secondary']: instance = getattr(orbital_position_container, component) points = instance.points faces = instance.faces if isinstance(instance.spots, (dict, )): for idx, spot in instance.spots.items(): faces = up.concatenate((faces, spot.faces + len(points)), axis=0) points = up.concatenate((points, spot.points), axis=0) setattr(instance, 'points', points) setattr(instance, 'faces', faces) points = np.concatenate((orbital_position_container.primary.points, orbital_position_container.secondary.points)) faces = np.concatenate( (orbital_position_container.primary.faces, orbital_position_container.secondary.faces + np.shape(orbital_position_container.primary.points)[0])) self.assertTrue(testutils.surface_closed(faces=faces, points=points))
def test_closed_surface_over_contact(self): s = testutils.prepare_binary_system( testutils.BINARY_SYSTEM_PARAMS['over-contact']) s.primary.discretization_factor = up.radians(10) s.init() orbital_position_container = testutils.prepare_orbital_position_container( s) orbital_position_container.build_mesh(components_distance=1.0) orbital_position_container.build_faces(components_distance=1.0) points = np.concatenate((orbital_position_container.primary.points, orbital_position_container.secondary.points)) faces = np.concatenate( (orbital_position_container.primary.faces, orbital_position_container.secondary.faces + np.shape(orbital_position_container.primary.points)[0])) self.assertTrue(testutils.surface_closed(faces=faces, points=points))
def test_closed_surface_squashed(self): position_container = self.build_system('squashed', up.radians(10)) self.assertTrue(testutils.surface_closed(faces=position_container.star.faces, points=position_container.star.points))