A library for working with phylogenetic and population genetic data.
v0.27.0
tree/function/manipulation.hpp
Go to the documentation of this file.
1 #ifndef GENESIS_TREE_FUNCTION_MANIPULATION_H_
2 #define GENESIS_TREE_FUNCTION_MANIPULATION_H_
3 
4 /*
5  Genesis - A toolkit for working with phylogenetic data.
6  Copyright (C) 2014-2020 Lucas Czech and HITS gGmbH
7 
8  This program is free software: you can redistribute it and/or modify
9  it under the terms of the GNU General Public License as published by
10  the Free Software Foundation, either version 3 of the License, or
11  (at your option) any later version.
12 
13  This program is distributed in the hope that it will be useful,
14  but WITHOUT ANY WARRANTY; without even the implied warranty of
15  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16  GNU General Public License for more details.
17 
18  You should have received a copy of the GNU General Public License
19  along with this program. If not, see <http://www.gnu.org/licenses/>.
20 
21  Contact:
22  Lucas Czech <lucas.czech@h-its.org>
23  Exelixis Lab, Heidelberg Institute for Theoretical Studies
24  Schloss-Wolfsbrunnenweg 35, D-69118 Heidelberg, Germany
25 */
26 
34 #include <cstddef> // size_t
35 #include <functional>
36 
37 #include "genesis/tree/tree.hpp"
39 
40 namespace genesis {
41 namespace tree {
42 
43 // =================================================================================================
44 // Forward Declarations
45 // =================================================================================================
46 
47 class Tree;
48 class TreeNode;
49 class TreeEdge;
50 class TreeLink;
51 class Subtree;
52 
53 // =================================================================================================
54 // Minimal Tree
55 // =================================================================================================
56 
68 
80 template< class NodeDataType = CommonNodeData, class EdgeDataType = CommonEdgeData >
82 {
83  auto tree = minimal_tree_topology();
84  for( size_t i = 0; i < tree.node_count(); ++i ) {
85  tree.node_at(i).reset_data( NodeDataType::create() );
86  }
87  for( size_t i = 0; i < tree.edge_count(); ++i ) {
88  tree.edge_at(i).reset_data( EdgeDataType::create() );
89  }
90  return tree;
91 }
92 
93 // =================================================================================================
94 // Add Nodes
95 // =================================================================================================
96 
131 TreeNode& add_new_node( Tree& tree, TreeNode& target_node );
132 
173 TreeNode& add_new_node(
174  Tree& tree,
175  TreeEdge& target_edge,
176  std::function<void( TreeEdge& target_edge, TreeEdge& new_edge )> adjust_edges = {}
177 );
178 
212 TreeNode& add_new_leaf_node(
213  Tree& tree,
214  TreeEdge& target_edge,
215  std::function<void( TreeEdge& target_edge, TreeEdge& new_edge )> adjust_edges = {}
216 );
217 
218 // =================================================================================================
219 // Delete Nodes
220 // =================================================================================================
221 
235 void delete_node( Tree& tree, TreeNode& target_node );
236 
242 void delete_leaf_node( Tree& tree, TreeNode& target_node );
243 
259 void delete_linear_node(
260  Tree& tree,
261  TreeNode& target_node,
262  std::function<void( TreeEdge& remaining_edge, TreeEdge& deleted_edge )> adjust_edges = {}
263 );
264 
276 void delete_subtree( Tree& tree, Subtree const& subtree );
277 
290 void delete_edge(
291  Tree& tree,
292  TreeEdge& target_edge,
293  std::function<void( TreeNode& remaining_node, TreeNode& deleted_node )> adjust_nodes = {}
294 );
295 
309 void delete_zero_branch_length_edges( Tree& tree, bool include_leaf_edges = false );
310 
311 // =================================================================================================
312 // Rooting, Rerooting, Unrooting
313 // =================================================================================================
314 
332 TreeNode& make_rooted(
333  Tree& tree,
334  TreeEdge& target_edge,
335  std::function<void( TreeEdge& target_edge, TreeEdge& new_edge )> adjust_edges = {}
336 );
337 
351 TreeNode& make_rooted(
352  Tree& tree,
353  std::function<void( TreeEdge& target_edge, TreeEdge& new_edge )> adjust_edges = {}
354 );
355 
365 void make_unrooted(
366  Tree& tree,
367  std::function<void( TreeEdge& remaining_edge, TreeEdge& deleted_edge )> adjust_edges = {}
368 );
369 
398 void change_rooting( Tree& tree, TreeLink& at_link );
399 
407 void change_rooting( Tree& tree, TreeNode& at_node );
408 
409 // =================================================================================================
410 // Ladderize
411 // =================================================================================================
412 
413 enum class LadderizeOrder
414 {
415  kSmallFirst,
417 };
418 
426 void ladderize( Tree& tree, LadderizeOrder order = LadderizeOrder::kSmallFirst );
427 
428 } // namespace tree
429 } // namespace genesis
430 
431 #endif // include guard
genesis::tree::LadderizeOrder::kSmallFirst
@ kSmallFirst
genesis::tree::delete_zero_branch_length_edges
void delete_zero_branch_length_edges(Tree &tree, bool include_leaf_edges)
Delete (contract) all branches of a CommonTree that have branch length zero.
Definition: tree/function/manipulation.cpp:706
genesis::tree::LadderizeOrder
LadderizeOrder
Definition: tree/function/manipulation.hpp:413
tree.hpp
Header of Tree class.
genesis::tree::add_new_node
TreeNode & add_new_node(Tree &tree, TreeNode &target_node)
Add a new Node as a leaf to an existing Node.
Definition: tree/function/manipulation.cpp:109
genesis::tree::add_new_leaf_node
TreeNode & add_new_leaf_node(Tree &tree, TreeEdge &target_edge, std::function< void(TreeEdge &target_edge, TreeEdge &new_edge)> adjust_edges)
Add a new Node as a leaf to an existing Edge, by also adding a new Node in the middle of that Edge.
Definition: tree/function/manipulation.cpp:263
genesis::tree::change_rooting
void change_rooting(Tree &tree, TreeLink &at_link)
"Reroot" the Tree at the given TreeLink.
Definition: tree/function/manipulation.cpp:779
genesis::tree::ladderize
void ladderize(Tree &tree, LadderizeOrder order)
Ladderize a Tree, that is, order its subtrees by size.
Definition: tree/function/manipulation.cpp:840
genesis::tree::Tree
Class for representing phylogenetic trees.
Definition: tree/tree.hpp:97
genesis::tree::make_unrooted
void make_unrooted(Tree &tree, std::function< void(TreeEdge &remaining_edge, TreeEdge &deleted_edge)> adjust_edges)
Unroot a Tree.
Definition: tree/function/manipulation.cpp:768
genesis
Container namespace for all symbols of genesis in order to keep them separate when used as a library.
Definition: placement/formats/edge_color.cpp:42
genesis::tree::LadderizeOrder::kLargeFirst
@ kLargeFirst
tree.hpp
genesis::tree::delete_edge
void delete_edge(Tree &tree, TreeEdge &target_edge, std::function< void(TreeNode &remaining_node, TreeNode &deleted_node)> adjust_nodes)
Delete a TreeEdge from a Tree, that is, contract the two TreeNodes at its ends into one TreeNode.
Definition: tree/function/manipulation.cpp:597
genesis::tree::minimal_tree_topology
Tree minimal_tree_topology()
Create a minimal Tree that can be used with manipulation functions such as add_new_node() or add_new_...
Definition: tree/function/manipulation.cpp:55
genesis::tree::delete_node
void delete_node(Tree &tree, TreeNode &target_node)
Delete a TreeNode from a Tree.
Definition: tree/function/manipulation.cpp:277
genesis::tree::delete_subtree
void delete_subtree(Tree &tree, Subtree const &subtree)
Delete a complete Subtree from a Tree.
Definition: tree/function/manipulation.cpp:518
genesis::tree::minimal_tree
Tree minimal_tree()
Create a minimal Tree that can be used with manipulation functions such as add_new_node() or add_new_...
Definition: tree/function/manipulation.hpp:81
genesis::tree::delete_leaf_node
void delete_leaf_node(Tree &tree, TreeNode &target_node)
Delete a leaf TreeNode.
Definition: tree/function/manipulation.cpp:296
genesis::tree::delete_linear_node
void delete_linear_node(Tree &tree, TreeNode &target_node, std::function< void(TreeEdge &remaining_edge, TreeEdge &deleted_edge)> adjust_edges)
Delete a "linear" TreeNode from a Tree, that is, a node with two neighbours.
Definition: tree/function/manipulation.cpp:385
genesis::tree::make_rooted
TreeNode & make_rooted(Tree &tree, TreeEdge &target_edge, std::function< void(TreeEdge &target_edge, TreeEdge &new_edge)> adjust_edges)
Root a Tree at a given TreeEdge.
Definition: tree/function/manipulation.cpp:747