X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=forester%2Fjava%2Fsrc%2Forg%2Fforester%2Fsdi%2FSDIR.java;h=a21fd6b6655c09764f158c34d587bdeeb808605a;hb=8aadbec1b0627ebd71e57e60e06621d4038bb79a;hp=c9224daa3a10ee347886e3ff24195fc3b07625be;hpb=ecfa899d355303acc65e93601e00922aad822883;p=jalview.git
diff --git a/forester/java/src/org/forester/sdi/SDIR.java b/forester/java/src/org/forester/sdi/SDIR.java
index c9224da..a21fd6b 100644
--- a/forester/java/src/org/forester/sdi/SDIR.java
+++ b/forester/java/src/org/forester/sdi/SDIR.java
@@ -23,7 +23,7 @@
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA
//
// Contact: phylosoft @ gmail . com
-// WWW: www.phylosoft.org/forester
+// WWW: https://sites.google.com/site/cmzmasek/home/software/forester
package org.forester.sdi;
@@ -45,11 +45,11 @@ import org.forester.phylogeny.iterators.PhylogenyNodeIterator;
*
Phylogeny height - which is the largest distance from
* root to external node (minimizing of which is the same as "midpoint rooting")
*
- *
+ *
* @see SDIse
- *
+ *
* @see SDI
- *
+ *
* @author Christian M. Zmasek
*/
public class SDIR {
@@ -76,7 +76,7 @@ public class SDIR {
/**
* Returns the number of differently rooted trees which minimize the
* (rooting) "criterion" - as determined by method "infer".
- *
+ *
* @see #infer(Phylogeny,Phylogeny,boolean,boolean,boolean,boolean,int,boolean)
* @return number of differently rooted trees which minimized the criterion
*/
@@ -97,7 +97,7 @@ public class SDIR {
* not necessarily zero.
*
* (Last modified: 01/22/00)
- *
+ *
* @see #infer(Phylogeny,Phylogeny,boolean,boolean,boolean,boolean,int,boolean)
* @return the minimal difference in tree heights -- IF calculated by
* "infer"
@@ -113,7 +113,7 @@ public class SDIR {
* IMPORTANT : If the tree is not rooted by minimizing the sum of
* duplications or the mapping cost L, then this number is NOT NECESSARILY
* the MINIMAL number of duplications.
- *
+ *
* @see #infer(Phylogeny,Phylogeny,boolean,boolean,boolean,boolean,int,boolean)
* @return (minimal) number of duplications
*/
@@ -126,7 +126,7 @@ public class SDIR {
* minimize_mapping_cost is set to true.
*
* (Last modified: 11/07/00)
- *
+ *
* @see #infer(Phylogeny,Phylogeny,boolean,boolean,boolean,boolean,int,boolean)
* @return the minimal mapping cost "L" -- IF calculated by "infer"
*/
@@ -142,7 +142,7 @@ public class SDIR {
* first criterion.
*
* (Last modified: 01/12/00)
- *
+ *
* @see #infer(Phylogeny,Phylogeny,boolean,boolean,boolean,boolean,int,boolean)
* @return the minimal tree height -- IF calculated by "infer"
*/
@@ -153,7 +153,7 @@ public class SDIR {
/**
* Returns the sum of times (in ms) needed to run method infer of class SDI.
* Final variable TIME needs to be set to true.
- *
+ *
* @return sum of times (in ms) needed to run method infer of class SDI
*/
public long getTimeSumSDI() {
@@ -188,7 +188,7 @@ public class SDIR {
*
*
* (Last modified: 10/01/01)
- *
+ *
* @param gene_tree
* a binary (except deepest node) gene Phylogeny
* @param species_tree
@@ -212,7 +212,7 @@ public class SDIR {
* Array) must be no lower than 1
* @return array of rooted Trees with duplication vs. speciation assigned if
* return_trees is set to true, null otherwise
- * @throws SDIException
+ * @throws SDIException
*/
public Phylogeny[] infer( final Phylogeny gene_tree,
final Phylogeny species_tree,
@@ -222,7 +222,7 @@ public class SDIR {
final boolean return_trees,
int max_trees_to_return ) throws SDIException {
init();
- SDIse sdise = null;
+ SDI sdise = null;
final ArrayList trees = new ArrayList();
Phylogeny[] tree_array = null;
List branches = null;
@@ -266,26 +266,23 @@ public class SDIR {
final PhylogenyNode n = iter.next();
if ( n.isRoot() ) {
if ( ( n.getNumberOfDescendants() != 2 ) && ( n.getNumberOfDescendants() != 3 ) ) {
- throw new IllegalArgumentException( "attempt to run SDI on gene tree with "
- + n.getNumberOfDescendants() + " child nodes at its root" );
+ throw new SDIException( "gene tree has " + n.getNumberOfDescendants() + " descendents at its root" );
}
}
else if ( !n.isExternal() && ( n.getNumberOfDescendants() != 2 ) ) {
- throw new IllegalArgumentException( "attempt to run SDI on gene tree which is not completely binary [found node with "
- + n.getNumberOfDescendants() + " child nodes]" );
+ throw new SDIException( "gene tree is not completely binary" );
}
}
for( final PhylogenyNodeIterator iter = species_tree.iteratorPostorder(); iter.hasNext(); ) {
final PhylogenyNode n = iter.next();
if ( !n.isExternal() && ( n.getNumberOfDescendants() != 2 ) ) {
- throw new IllegalArgumentException( "attempt to run SDI with a species tree which is not completely binary (after stripping) [found node with "
- + n.getNumberOfDescendants() + " child nodes]" );
+ throw new SDIException( "species tree (after stripping) is not completely binary" );
}
}
g.reRoot( g.getFirstExternalNode() );
branches = SDIR.getBranchesInPreorder( g );
if ( minimize_mapping_cost || minimize_sum_of_dup ) {
- sdise = new SDIse( g, species_tree );
+ sdise = new SDI( g, species_tree );
duplications = sdise.getDuplicationsSum();
}
final Set used_root_placements = new HashSet();
@@ -295,7 +292,7 @@ public class SDIR {
prev_root_c2 = prev_root.getChildNode2();
prev_root_was_dup = prev_root.isDuplication();
final PhylogenyBranch current_branch = branches.get( j );
- g.reRoot( current_branch );
+ GSDIR.reRoot( current_branch, g );
if ( minimize_mapping_cost || minimize_sum_of_dup ) {
duplications = sdise.updateM( prev_root_was_dup, prev_root_c1, prev_root_c2 );
}
@@ -421,7 +418,7 @@ public class SDIR {
height = height__diff[ 0 ];
diff = height__diff[ 1 ];
if ( Math.abs( diff ) < SDIR.ZERO_DIFF ) {
- sdise = new SDIse( g, species_tree );
+ sdise = new SDI( g, species_tree );
min_duplications = sdise.getDuplicationsSum();
min_height = height;
min_diff = Math.abs( diff );
@@ -497,8 +494,8 @@ public class SDIR {
branches.add( new PhylogenyBranch( t.getRoot().getChildNode1(), t.getRoot().getChildNode2() ) );
return branches;
}
- final Set one = new HashSet();
- final Set two = new HashSet();
+ final Set one = new HashSet();
+ final Set two = new HashSet();
PhylogenyNode node = t.getRoot();
while ( !node.isRoot() || !two.contains( node.getId() ) ) {
if ( !node.isExternal() && !two.contains( node.getId() ) ) {