def tile_with_mixed_patterns(config): options = extract_options(config); options["dof_type"] = str(config.get("dof_type", "isotropic")); options["thickness_type"] = str(config.get("thickness_type", "vertex")); networks = load_wires(str(config["wire_list_file"])); guide_mesh = pymesh.load_mesh(config["guide_mesh"]); per_vertex_thickness = options["thickness_type"] == "vertex"; tiler = Tiler(); tiler.tile_with_mixed_patterns(mesh, per_vertex_thickness, options["dof_type"] == "isotropic"); network = tiler.wire_network; inflator = Inflator(network); inflator.inflate(network.get_attribute("thickness").ravel(), per_vertex_thickness); return inflator.mesh;
def tile_with_mixed_patterns(config): options = extract_options(config); options["dof_type"] = str(config.get("dof_type", "isotropic")); options["thickness_type"] = str(config.get("thickness_type", "vertex")); networks = load_wires(str(config["wire_list_file"])); guide_mesh = pymesh.load_mesh(config["guide_mesh"]); per_vertex_thickness = options["thickness_type"] == "vertex"; tiler = Tiler(); tiler.tile_with_mixed_patterns(mesh, per_vertex_thickness, options["dof_type"] == "isotropic"); network = tiler.wire_network; if config.get("trim", False): network.trim(); inflator = Inflator(network); inflator.inflate(network.get_attribute("thickness").ravel(), per_vertex_thickness); return inflator.mesh;
def test_tile_with_mixed_patterns(self): wire_networks = [self.get_brick5(), self.get_cross_3D()] params = [ Parameters(wire_networks[0], 0.1), Parameters(wire_networks[1], 0.0) ] params[0].load_default_isotropic_parameters() params[1].load_default_isotropic_parameters() max_num_dofs = max(params[0].num_dofs, params[1].num_dofs) pattern_id = np.array([0, 1, 1, 0, 1, 0, 0, 1]) mesh = generate_box_mesh(np.zeros(3), np.ones(3), subdiv_order=1, using_simplex=False) mesh.add_attribute("pattern_id") mesh.set_attribute("pattern_id", pattern_id) for i in range(max_num_dofs): dof_name = "dof_{}".format(i) dof = np.array([ params[j].dofs[i] if i < params[j].num_dofs else 0 for j in pattern_id ]) mesh.add_attribute(dof_name) mesh.set_attribute(dof_name, dof) tiler = Tiler(wire_networks) tiler.tile_with_mixed_patterns(mesh) tiled_wire_network = tiler.wire_network self.assert_array_equal(mesh.bbox, tiled_wire_network.bbox) self.assertEqual( (wire_networks[0].num_edges + wire_networks[1].num_edges) * 4, tiled_wire_network.num_edges) self.assertEqual( (wire_networks[0].num_vertices + wire_networks[1].num_vertices) * 4 - 4 * 3, tiled_wire_network.num_vertices)
def test_tile_with_mixed_patterns(self): wire_networks = [ self.get_brick5(), self.get_cross_3D() ]; params = [Parameters(wire_networks[0], 0.1), Parameters(wire_networks[1], 0.0)]; params[0].load_default_isotropic_parameters(); params[1].load_default_isotropic_parameters(); max_num_dofs = max(params[0].num_dofs, params[1].num_dofs); pattern_id = np.array([0, 1, 1, 0, 1, 0, 0, 1]); mesh = generate_box_mesh(np.zeros(3), np.ones(3), subdiv_order=1, using_simplex=False); mesh.add_attribute("pattern_id"); mesh.set_attribute("pattern_id", pattern_id); for i in range(max_num_dofs): dof_name = "dof_{}".format(i); dof = np.array([params[j].dofs[i] if i < params[j].num_dofs else 0 for j in pattern_id ]); mesh.add_attribute(dof_name); mesh.set_attribute(dof_name, dof); tiler = Tiler(wire_networks); tiler.tile_with_mixed_patterns(mesh); tiled_wire_network = tiler.wire_network; self.assert_array_equal( mesh.bbox, tiled_wire_network.bbox); self.assertEqual( (wire_networks[0].num_edges + wire_networks[1].num_edges) * 4, tiled_wire_network.num_edges); self.assertEqual( (wire_networks[0].num_vertices + wire_networks[1].num_vertices) * 4 - 4 * 3, tiled_wire_network.num_vertices);
def tile_with_mixed_patterns(config): options = extract_options(config) options["dof_type"] = str(config.get("dof_type", "isotropic")) options["thickness_type"] = str(config.get("thickness_type", "vertex")) networks = load_wires(str(config["wire_list_file"])) guide_mesh = pymesh.load_mesh(config["guide_mesh"]) per_vertex_thickness = options["thickness_type"] == "vertex" tiler = Tiler() tiler.tile_with_mixed_patterns(mesh, per_vertex_thickness, options["dof_type"] == "isotropic") network = tiler.wire_network if config.get("trim", False): network.trim() if "output_wire_network" in config: network.write_to_file(config["output_wire_network"]) inflator = Inflator(network) inflator.subdivide_order = options["subdiv"] inflator.subdivide_method = options["subdiv_method"] inflator.inflate( network.get_attribute("thickness").ravel(), per_vertex_thickness) return inflator.mesh