+ final String g2_0_ = "(([&&NHX:S=a1],[&&NHX:S=a2]),([&&NHX:S=o2],[&&NHX:S=o4]))";
+ final Phylogeny g2_0p = TestGSDI.createPhylogeny( g2_0_ );
+ g2_0.setRooted( true );
+ final GSDI sdi2_0p = new GSDI( g2_0p, s2, false, false, false );
+ if ( sdi2_0p.getDuplicationsSum() != 0 ) {
+ return false;
+ }
+ //--
+ final Phylogeny tol_143_ = ParserUtils.readPhylogenies( PATH_TO_TEST_DATA + "tol_143.xml" )[ 0 ];
+ final Phylogeny gene_tree_tax_code_4_ = ParserUtils.readPhylogenies( PATH_TO_TEST_DATA
+ + "gene_tree_tax_code_4.xml" )[ 0 ];
+ final GSDI gsdi_143_4_1 = new GSDI( gene_tree_tax_code_4_.copy(), tol_143_.copy(), false, true, true );
+ if ( gsdi_143_4_1.getDuplicationsSum() != 21 ) {
+ return false;
+ }
+ if ( gsdi_143_4_1.getSpeciationsSum() != 28 ) {
+ return false;
+ }
+ if ( gsdi_143_4_1.getSpeciationOrDuplicationEventsSum() != 6 ) {
+ return false;
+ }
+ //--
+ final Phylogeny gsdi_test_gene_tree_sn_wnt = ParserUtils.readPhylogenies( PATH_TO_TEST_DATA
+ + "gsdi_test_gene_tree_sn_wnt.xml" )[ 0 ];
+ gsdi_test_gene_tree_sn_wnt.setRooted( true );
+ final GSDI a = new GSDI( gsdi_test_gene_tree_sn_wnt.copy(), tol_143_.copy(), false, true, true );
+ if ( a.getDuplicationsSum() != 33 ) {
+ return false;
+ }
+ if ( a.getSpeciationsSum() != 31 ) {
+ return false;
+ }
+ if ( a.getSpeciationOrDuplicationEventsSum() != 0 ) {
+ return false;
+ }
+ if ( a.getTaxCompBase() != TaxonomyComparisonBase.SCIENTIFIC_NAME ) {
+ return false;
+ }
+ if ( a.getMappedExternalSpeciesTreeNodes().size() != 26 ) {
+ return false;
+ }
+ if ( a.getReMappedScientificNamesFromGeneTree().size() != 0 ) {
+ return false;
+ }
+ //--
+ final Phylogeny gsdi_test_species_tree_sn_xml = ParserUtils.readPhylogenies( PATH_TO_TEST_DATA
+ + "gsdi_test_species_tree_sn.xml" )[ 0 ];
+ final GSDI b = new GSDI( gsdi_test_gene_tree_sn_wnt.copy(),
+ gsdi_test_species_tree_sn_xml.copy(),
+ false,
+ true,
+ true );
+ if ( b.getDuplicationsSum() != 8 ) {
+ return false;
+ }
+ if ( b.getSpeciationsSum() != 2 ) {
+ return false;
+ }
+ if ( b.getSpeciationOrDuplicationEventsSum() != 0 ) {
+ return false;
+ }
+ if ( b.getTaxCompBase() != TaxonomyComparisonBase.SCIENTIFIC_NAME ) {
+ return false;
+ }
+ if ( b.getMappedExternalSpeciesTreeNodes().size() != 2 ) {
+ return false;
+ }
+ if ( b.getReMappedScientificNamesFromGeneTree().size() != 0 ) {
+ return false;
+ }
+ if ( b.getStrippedExternalGeneTreeNodes().size() != 87 ) {
+ return false;
+ }
+ if ( b.getStrippedSpeciesTreeNodes().size() != 17 ) {
+ return false;
+ }
+ //--
+ final Phylogeny gsdi_test_species_tree_sn_nh = TestGSDI
+ .createPhylogeny( "((((((('Homo sapiens','Mus musculus')Euarchontoglires,'Petromyzon marinus')Vertebrata,'Nematostella vectensis')'Bilateria Cnidaria',(('Mycosphaerella graminicola','Mycosphaerella pini')Mycosphaerella,'Saccharomyces cerevisiae')'Pezizomycotina Saccharomycetales')Opisthokonta,('Plasmodium chabaudi','Plasmodium falciparum','Plasmodium yoelii yoelii')Plasmodium)Eukaryota,'Pyrococcus horikoshii')Neomura,(('Kineococcus radiotolerans','Kocuria rhizophila','Streptomyces coelicolor','Thermobifida fusca','Microlunatus phosphovorus'),'Bacteroides thetaiotaomicron'))'cellular organisms';" );
+ PhylogenyMethods.transferNodeNameToField( gsdi_test_species_tree_sn_nh,
+ PhylogenyMethods.PhylogenyNodeField.TAXONOMY_SCIENTIFIC_NAME,
+ true );
+ final GSDI c = new GSDI( gsdi_test_gene_tree_sn_wnt.copy(),
+ gsdi_test_species_tree_sn_nh.copy(),
+ false,
+ true,
+ true );
+ if ( c.getDuplicationsSum() != 8 ) {
+ return false;
+ }
+ if ( c.getSpeciationsSum() != 2 ) {
+ return false;
+ }
+ if ( c.getSpeciationOrDuplicationEventsSum() != 0 ) {
+ return false;
+ }
+ if ( c.getTaxCompBase() != TaxonomyComparisonBase.SCIENTIFIC_NAME ) {
+ return false;
+ }
+ if ( c.getMappedExternalSpeciesTreeNodes().size() != 2 ) {
+ return false;
+ }
+ if ( c.getReMappedScientificNamesFromGeneTree().size() != 0 ) {
+ return false;
+ }
+ if ( c.getStrippedExternalGeneTreeNodes().size() != 87 ) {
+ return false;
+ }
+ if ( c.getStrippedSpeciesTreeNodes().size() != 15 ) {
+ return false;
+ }
+ //--
+ final Phylogeny gsdi_test_gene_tree_codes_xml = ParserUtils.readPhylogenies( PATH_TO_TEST_DATA
+ + "gsdi_test_gene_tree_codes.xml" )[ 0 ];
+ final Phylogeny gsdi_test_species_tree_codes_xml = ParserUtils.readPhylogenies( PATH_TO_TEST_DATA
+ + "gsdi_test_species_tree_codes.xml" )[ 0 ];
+ final GSDI d = new GSDI( gsdi_test_gene_tree_codes_xml.copy(),
+ gsdi_test_species_tree_codes_xml.copy(),
+ false,
+ true,
+ true );
+ if ( d.getDuplicationsSum() != 21 ) {
+ return false;
+ }
+ if ( d.getSpeciationsSum() != 28 ) {
+ return false;
+ }
+ if ( d.getSpeciationOrDuplicationEventsSum() != 6 ) {
+ return false;
+ }
+ if ( d.getTaxCompBase() != TaxonomyComparisonBase.CODE ) {
+ return false;
+ }
+ if ( d.getMappedExternalSpeciesTreeNodes().size() != 17 ) {
+ return false;
+ }
+ if ( d.getReMappedScientificNamesFromGeneTree().size() != 0 ) {
+ return false;
+ }
+ if ( d.getStrippedExternalGeneTreeNodes().size() != 12 ) {
+ return false;
+ }
+ if ( d.getStrippedSpeciesTreeNodes().size() != 3 ) {
+ return false;
+ }
+ //--
+ final Phylogeny gsdi_test_gene_tree_sn_xml = ParserUtils.readPhylogenies( PATH_TO_TEST_DATA
+ + "gsdi_test_gene_tree_sn.xml" )[ 0 ];
+ final GSDI e = new GSDI( gsdi_test_gene_tree_sn_xml.copy(),
+ gsdi_test_species_tree_sn_xml.copy(),
+ false,
+ true,
+ true );
+ if ( e.getDuplicationsSum() != 7 ) {
+ return false;
+ }
+ if ( e.getSpeciationsSum() != 9 ) {
+ return false;
+ }
+ if ( e.getSpeciationOrDuplicationEventsSum() != 1 ) {
+ return false;
+ }
+ if ( e.getTaxCompBase() != TaxonomyComparisonBase.SCIENTIFIC_NAME ) {
+ return false;
+ }
+ if ( e.getMappedExternalSpeciesTreeNodes().size() != 12 ) {
+ return false;
+ }
+ if ( e.getReMappedScientificNamesFromGeneTree().size() != 8 ) {
+ return false;
+ }
+ if ( e.getStrippedExternalGeneTreeNodes().size() != 3 ) {
+ return false;
+ }
+ if ( e.getStrippedSpeciesTreeNodes().size() != 7 ) {
+ return false;
+ }
+ }
+ catch ( final Exception e ) {
+ e.printStackTrace( System.out );
+ return false;
+ }
+ return true;
+ }
+
+ private static boolean testGSDIR_general() {
+ try {
+ final PhylogenyFactory factory0 = ParserBasedPhylogenyFactory.getInstance();
+ final String s0 = "([&&NHX:S=1]);";
+ final String gene_0_str = "([&&NHX:S=1],[&&NHX:S=1]);";
+ final Phylogeny s_0 = factory0.create( s0, new NHXParser() )[ 0 ];
+ final Phylogeny gene_0 = factory0.create( gene_0_str, new NHXParser() )[ 0 ];
+ s_0.setRooted( true );
+ gene_0.setRooted( true );
+ final GSDIR sdi0 = new GSDIR( gene_0, s_0, true, true );
+ if ( sdi0.getSpeciationsSum() != 0 ) {
+ return false;
+ }
+ if ( sdi0.getMinDuplicationsSum() != 1 ) {
+ return false;
+ }
+ //
+ final PhylogenyFactory factory00 = ParserBasedPhylogenyFactory.getInstance();
+ final String s00 = "[&&NHX:S=1];";
+ final String gene_00_str = "([&&NHX:S=1],[&&NHX:S=1]);";
+ final Phylogeny s_00 = factory00.create( s00, new NHXParser() )[ 0 ];
+ final Phylogeny gene_00 = factory00.create( gene_00_str, new NHXParser() )[ 0 ];
+ s_00.setRooted( true );
+ gene_00.setRooted( true );
+ final GSDIR sdi00 = new GSDIR( gene_00, s_00, true, true );
+ if ( sdi00.getSpeciationsSum() != 0 ) {
+ return false;
+ }
+ if ( sdi00.getMinDuplicationsSum() != 1 ) {
+ return false;
+ }
+ //
+ final String s1str = "(((([&&NHX:S=HUMAN],([&&NHX:S=MOUSE],[&&NHX:S=RAT])),([&&NHX:S=CAEEL],[&&NHX:S=CAEBR])),[&&NHX:S=YEAST]),[&&NHX:S=ARATH])";
+ final Phylogeny s1 = ParserBasedPhylogenyFactory.getInstance().create( s1str, new NHXParser() )[ 0 ];
+ s1.setRooted( true );
+ final Phylogeny g1 = TestGSDI
+ .createPhylogeny( "(HUMAN[&&NHX:S=HUMAN],(RAT[&&NHX:S=RAT],(CAEEL[&&NHX:T=:S=CAEEL],YEAST[&&NHX:S=YEAST])))" );
+ final GSDIR sdi1 = new GSDIR( g1.copy(), s1.copy(), false, false );
+ if ( sdi1.getMinDuplicationsSum() != 0 ) {
+ return false;
+ }
+ final Phylogeny g2 = TestGSDI
+ .createPhylogeny( "(((HUMAN[&&NHX:S=HUMAN],RAT[&&NHX:S=RAT]),CAEEL[&&NHX:T=:S=CAEEL]),YEAST[&&NHX:S=YEAST])" );
+ final GSDIR sdi2 = new GSDIR( g2.copy(), s1.copy(), false, false );
+ if ( sdi2.getMinDuplicationsSum() != 0 ) {
+ return false;
+ }
+ final Phylogeny g3 = TestGSDI
+ .createPhylogeny( "(RAT[&&NHX:S=RAT],HUMAN[&&NHX:S=HUMAN],(YEAST[&&NHX:S=YEAST],CAEEL[&&NHX:T=:S=CAEEL]))" );
+ final GSDIR sdi3 = new GSDIR( g3.copy(), s1.copy(), false, false );
+ if ( sdi3.getMinDuplicationsSum() != 0 ) {
+ return false;
+ }
+ final Phylogeny g4 = TestGSDI
+ .createPhylogeny( "(((((MOUSE[&&NHX:S=MOUSE],[&&NHX:S=RAT]),[&&NHX:S=HUMAN]),([&&NHX:S=ARATH],[&&NHX:S=YEAST])),[&&NHX:S=CAEEL]),[&&NHX:S=CAEBR])" );
+ final GSDIR sdi4 = new GSDIR( g4.copy(), s1.copy(), false, false );
+ if ( sdi4.getMinDuplicationsSum() != 0 ) {
+ return false;
+ }