academic/PhyML: Updated for version 20131112.

Signed-off-by: Willy Sudiarto Raharjo <willysr@slackbuilds.org>

1 files changed, 1591 insertions, 0 deletions

diff --git a/academic/PhyML/20131112.patch b/academic/PhyML/20131112.patch
new file mode 100644
index 0000000000..e7eb2670b6
--- /dev/null
+++ b/academic/PhyML/20131112.patch
@@ -0,0 +1,1591 @@
+unchanged:
+--- ./phyml-20120412/src/help.c	2011-09-22 14:51:57.000000000 +1200
++++ ./phyml-20120412-new-bak/src/help.c	2013-02-18 21:51:51.448662440 +1300
+@@ -96,7 +96,7 @@
+   PhyML_Printf("\n");
+ 
+   
+-  PhyML_Printf("%s\n\t-q (or --sequential)\n",BOLD);
++  PhyML_Printf("%s\n\t-q (or --sequential)%s\n",BOLD,FLAT);
+   PhyML_Printf("%s\t\tChanges interleaved format (default) to sequential format.\n",FLAT);
+   PhyML_Printf("\n");
+ 
+unchanged:
+--- ./phyml-20120412/src/lk.c	2012-04-05 14:47:37.000000000 +1200
++++ ./phyml-20120412-new-bak/src/lk.c	2013-02-19 10:22:22.665284529 +1300
+@@ -472,7 +472,7 @@
+   phydbl site_lk_cat, site_lk;
+   int sum_scale_left, sum_scale_rght;
+   int fact_sum_scale;
+-  phydbl max_sum_scale;
++  phydbl max_sum_scale,min_sum_scale;
+   phydbl sum;
+   int ambiguity_check,state;
+   int catg,ns,k,l,site;
+@@ -481,11 +481,12 @@
+   phydbl multiplier;
+   int exponent, piecewise_exponent;
+   phydbl tmp;
+-  phydbl logbig;
++  phydbl logbig,logsmall;
+   phydbl inv_site_lk;
+ 
+ 
+-  logbig = LOG((phydbl)BIG);
++  logbig   = LOG((phydbl)BIG);
++  logsmall = LOG((phydbl)SMALL);
+ 
+   dim1 = tree->mod->n_catg * tree->mod->ns;
+   dim2 = tree->mod->ns;
+@@ -555,7 +556,7 @@
+                             b->Pij_rr[catg*dim3+k*dim2+l] *
+                             b->p_lk_left[site*dim1+catg*dim2+l];
+                         }
+-
++		      
+                       site_lk_cat +=
+                         sum *
+                         tree->mod->pi->v[k] *
+@@ -578,7 +579,7 @@
+                         b->Pij_rr[catg*dim3+k*dim2+l] *
+                         b->p_lk_left[site*dim1+catg*dim2+l];
+                     }
+-
++		  
+                   site_lk_cat +=
+                     sum *
+                     tree->mod->pi->v[k] *
+@@ -588,8 +589,9 @@
+ 	}
+       tree->site_lk_cat[catg] = site_lk_cat;
+     }
+-    
++
+   max_sum_scale =  (phydbl)BIG;
++  min_sum_scale =  (phydbl)SMALL;
+   For(catg,tree->mod->n_catg)
+     {
+       sum_scale_left_cat[catg] =
+@@ -606,18 +608,30 @@
+ 
+       if(sum < .0)
+ 	{
+-	  PhyML_Printf("\n. sum = %G",sum);
+-	  PhyML_Printf("\n. Err in file %s at line %d\n\n",__FILE__,__LINE__);
++	  PhyML_Printf("\n== sum = %G",sum);
++	  PhyML_Printf("\n== Err in file %s at line %d\n\n",__FILE__,__LINE__);
+ 	  Warn_And_Exit("\n");
+ 	}
+ 
+       tmp = sum + (logbig - LOG(tree->site_lk_cat[catg]))/(phydbl)LOG2;
+       if(tmp < max_sum_scale) max_sum_scale = tmp; /* min of the maxs */
++
++      tmp = sum + (logsmall - LOG(tree->site_lk_cat[catg]))/(phydbl)LOG2;
++      if(tmp > min_sum_scale) min_sum_scale = tmp; /* max of the mins */
+     }
+ 
+-/*   fact_sum_scale = (int)((max_sum_scale + min_sum_scale) / 2); */
++  if(min_sum_scale > max_sum_scale)
++    {
++      /* PhyML_Printf("\n== Numerical precision issue alert."); */
++      /* PhyML_Printf("\n== min_sum_scale = %G max_sum_scale = %G",min_sum_scale,max_sum_scale); */
++      /* PhyML_Printf("\n== Err in file %s at line %d\n\n",__FILE__,__LINE__); */
++      /* Warn_And_Exit("\n"); */
++      min_sum_scale = max_sum_scale;
++    }
++
++  fact_sum_scale = (int)((max_sum_scale + min_sum_scale) / 2);
+ 
+-  fact_sum_scale = (int)(max_sum_scale / 2);
++  /* fact_sum_scale = (int)(max_sum_scale / 2); */
+ 
+   /* Apply scaling factors */
+   For(catg,tree->mod->n_catg)
+@@ -651,12 +665,13 @@
+ 
+       if(isinf(site_lk_cat))
+ 	{
+-	  PhyML_Printf("\n+ site=%4d cat=%4d site_lk_cat=%G sum_scale=%d max=%G fact=%d expo=%d dbl=%G",
++	  PhyML_Printf("\n+ site=%4d cat=%4d site_lk_cat=%G sum_scale=%d max=%G min=%G fact=%d expo=%d dbl=%G",
+ 		       tree->curr_site,
+ 		       catg,
+ 		       tree->site_lk_cat[catg],
+ 		       sum_scale_left_cat[catg]+sum_scale_rght_cat[catg],
+ 		       max_sum_scale,
++		       min_sum_scale,
+ 		       fact_sum_scale,
+ 		       -(sum_scale_left_cat[catg]+sum_scale_rght_cat[catg])+fact_sum_scale,
+ 		       (double)tree->site_lk_cat[catg] * pow(2.,-(sum_scale_left_cat[catg]+sum_scale_rght_cat[catg])+fact_sum_scale));
+@@ -678,13 +693,11 @@
+ 			   fact_sum_scale,
+ 			   -(sum_scale_left_cat[catg]+sum_scale_rght_cat[catg])+fact_sum_scale,
+ 			   (double)tree->site_lk_cat[catg] * pow(2.,-(sum_scale_left_cat[catg]+sum_scale_rght_cat[catg])+fact_sum_scale));
+-
++	      
+ 	      Exit("\n");
+ 	    }
+-
+ 	  site_lk_cat = .0;
+ 	}
+-
+       tree->site_lk_cat[catg] = site_lk_cat;
+     }
+ 
+@@ -694,6 +707,24 @@
+       site_lk += tree->site_lk_cat[catg] * tree->mod->gamma_r_proba->v[catg];
+     }
+ 
++  if(isinf(site_lk))
++    {
++      PhyML_Printf("\n= Site = %d",site);
++      PhyML_Printf("\n= invar = %d",tree->data->invar[site]);
++      PhyML_Printf("\n= scale_left = %d scale_rght = %d",sum_scale_left,sum_scale_rght);
++      PhyML_Printf("\n= lk = %G log(lk) = %f < %G",site_lk,log_site_lk,-BIG);
++      PhyML_Printf("\n= fact_sum_scale = %d",fact_sum_scale);
++      For(catg,tree->mod->n_catg) PhyML_Printf("\n. rr=%f p=%f lk=%G expo=%d",
++					       tree->mod->gamma_rr->v[catg],
++					       tree->mod->gamma_r_proba->v[catg],
++					       tree->site_lk_cat[catg],
++					       -(sum_scale_left_cat[catg]+sum_scale_rght_cat[catg])+fact_sum_scale);
++      PhyML_Printf("\n= pinv = %G",tree->mod->pinvar->v);
++      PhyML_Printf("\n= bl mult = %G",tree->mod->br_len_multiplier->v);
++      PhyML_Printf("\n. Err in file %s at line %d",__FILE__,__LINE__);
++      Warn_And_Exit("\n");
++    }
++
+ 
+   inv_site_lk = 0.;
+   
+@@ -741,19 +772,22 @@
+ 
+   log_site_lk = LOG(site_lk) - (phydbl)LOG2 * fact_sum_scale;
+ 
+-
+   For(catg,tree->mod->n_catg) tree->log_site_lk_cat[catg][site] = LOG(tree->site_lk_cat[catg]) - (phydbl)LOG2 * fact_sum_scale;
+   
+   if(isinf(log_site_lk) || isnan(log_site_lk))
+     {
+-      PhyML_Printf("\n. Site = %d",site);
+-      PhyML_Printf("\n. invar = %d",tree->data->invar[site]);
+-      PhyML_Printf("\n. scale_left = %d scale_rght = %d",sum_scale_left,sum_scale_rght);
+-      PhyML_Printf("\n. inv_site_lk = %f",inv_site_lk);
+-      PhyML_Printf("\n. Lk = %G LOG(Lk) = %f < %G",site_lk,log_site_lk,-BIG);
+-      For(catg,tree->mod->n_catg) PhyML_Printf("\n. rr=%f p=%f",tree->mod->gamma_rr->v[catg],tree->mod->gamma_r_proba->v[catg]);
+-      PhyML_Printf("\n. pinv = %G",tree->mod->pinvar->v);
+-      PhyML_Printf("\n. bl mult = %G",tree->mod->br_len_multiplier->v);
++      PhyML_Printf("\n= Site = %d",site);
++      PhyML_Printf("\n= invar = %d",tree->data->invar[site]);
++      PhyML_Printf("\n= scale_left = %d scale_rght = %d",sum_scale_left,sum_scale_rght);
++      PhyML_Printf("\n= inv_site_lk = %f",inv_site_lk);
++      PhyML_Printf("\n= lk = %G log(lk) = %f < %G",site_lk,log_site_lk,-BIG);
++      For(catg,tree->mod->n_catg) PhyML_Printf("\n. rr=%f p=%f lk=%G",
++					       tree->mod->gamma_rr->v[catg],
++					       tree->mod->gamma_r_proba->v[catg],
++					       tree->site_lk_cat[catg]);
++      PhyML_Printf("\n= pinv = %G",tree->mod->pinvar->v);
++      PhyML_Printf("\n= bl mult = %G",tree->mod->br_len_multiplier->v);
++      PhyML_Printf("\n= fact_sum_scale = %d",fact_sum_scale);
+ 
+       /* int i; */
+       /* For(i,2*tree->n_otu-3) */
+@@ -802,6 +836,7 @@
+     Alias_One_Subpatt((d==b->left)?(b->rght):(b->left),d,tree);
+ 
+   if(d->tax) return;
++
+   
+   if(tree->mod->use_m4mod == NO)
+     {
+@@ -813,6 +848,10 @@
+ 	{
+ 	  Update_P_Lk_AA(tree,b,d);
+ 	}
++      else
++        {
++	  Update_P_Lk_Generic(tree,b,d);
++        }
+     }
+   else
+     {
+@@ -1044,7 +1083,7 @@
+ 		  
+ 		  p_lk[site*NsNg+catg*Ns+i] = p1_lk1 * p2_lk2;
+ 		  
+-		  /* 	      PhyML_Printf("\n+ %G",p_lk[site*NsNg+catg*Ns+i]); */
++                  /* PhyML_Printf("\n+ %G",p_lk[site*NsNg+catg*Ns+i]); */
+ 		  
+ 		  if(p_lk[site*NsNg+catg*Ns+i] < smallest_p_lk) smallest_p_lk = p_lk[site*NsNg+catg*Ns+i] ;
+ 		}
+unchanged:
+--- ./phyml-20120412/src/Makefile	1970-01-01 12:00:00.000000000 +1200
++++ ./phyml-20120412-new-bak/src/Makefile	2013-02-18 21:54:46.957806553 +1300
+@@ -0,0 +1,1082 @@
++# Makefile.in generated by automake 1.11.1 from Makefile.am.
++# src/Makefile.  Generated from Makefile.in by configure.
++
++# Copyright (C) 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002,
++# 2003, 2004, 2005, 2006, 2007, 2008, 2009  Free Software Foundation,
++# Inc.
++# This Makefile.in is free software; the Free Software Foundation
++# gives unlimited permission to copy and/or distribute it,
++# with or without modifications, as long as this notice is preserved.
++
++# This program is distributed in the hope that it will be useful,
++# but WITHOUT ANY WARRANTY, to the extent permitted by law; without
++# even the implied warranty of MERCHANTABILITY or FITNESS FOR A
++# PARTICULAR PURPOSE.
++
++
++
++
++pkgdatadir = $(datadir)/phyml
++pkgincludedir = $(includedir)/phyml
++pkglibdir = $(libdir)/phyml
++pkglibexecdir = $(libexecdir)/phyml
++am__cd = CDPATH="$${ZSH_VERSION+.}$(PATH_SEPARATOR)" && cd
++install_sh_DATA = $(install_sh) -c -m 644
++install_sh_PROGRAM = $(install_sh) -c
++install_sh_SCRIPT = $(install_sh) -c
++INSTALL_HEADER = $(INSTALL_DATA)
++transform = $(program_transform_name)
++NORMAL_INSTALL = :
++PRE_INSTALL = :
++POST_INSTALL = :
++NORMAL_UNINSTALL = :
++PRE_UNINSTALL = :
++POST_UNINSTALL = :
++build_triplet = x86_64-unknown-linux-gnu
++host_triplet = x86_64-unknown-linux-gnu
++bin_PROGRAMS = phyml$(EXEEXT)
++#bin_PROGRAMS = test$(EXEEXT)
++#bin_PROGRAMS = phyml-mpi$(EXEEXT)
++#bin_PROGRAMS = rf$(EXEEXT)
++#bin_PROGRAMS = m4$(EXEEXT)
++#bin_PROGRAMS = tiporder$(EXEEXT)
++#bin_PROGRAMS = part$(EXEEXT)
++#bin_PROGRAMS =  \
++#	phycont$(EXEEXT)
++#bin_PROGRAMS = phytime$(EXEEXT)
++subdir = src
++DIST_COMMON = $(srcdir)/Makefile.am $(srcdir)/Makefile.in
++ACLOCAL_M4 = $(top_srcdir)/aclocal.m4
++am__aclocal_m4_deps = $(top_srcdir)/configure.ac
++am__configure_deps = $(am__aclocal_m4_deps) $(CONFIGURE_DEPENDENCIES) \
++	$(ACLOCAL_M4)
++mkinstalldirs = $(install_sh) -d
++CONFIG_HEADER = $(top_builddir)/config.h
++CONFIG_CLEAN_FILES =
++CONFIG_CLEAN_VPATH_FILES =
++am__installdirs = "$(DESTDIR)$(bindir)"
++PROGRAMS = $(bin_PROGRAMS)
++am__m4_SOURCES_DIST = main.c utilities.c utilities.h optimiz.c \
++	optimiz.h lk.c lk.h bionj.c bionj.h models.c models.h free.c \
++	free.h help.c help.h simu.c simu.h eigen.c eigen.h pars.c \
++	pars.h alrt.c alrt.h interface.c interface.h cl.c cl.h mg.c \
++	mg.h times.c times.h mcmc.c mcmc.h rates.c rates.h spr.c spr.h \
++	draw.c draw.h stats.c stats.h tiporder.c tiporder.h m4.c m4.h
++#am_m4_OBJECTS = main.$(OBJEXT) \
++#	utilities.$(OBJEXT) \
++#	optimiz.$(OBJEXT) \
++#	lk.$(OBJEXT) \
++#	bionj.$(OBJEXT) \
++#	models.$(OBJEXT) \
++#	free.$(OBJEXT) \
++#	help.$(OBJEXT) \
++#	simu.$(OBJEXT) \
++#	eigen.$(OBJEXT) \
++#	pars.$(OBJEXT) \
++#	alrt.$(OBJEXT) \
++#	interface.$(OBJEXT) \
++#	cl.$(OBJEXT) \
++#	mg.$(OBJEXT) \
++#	times.$(OBJEXT) \
++#	mcmc.$(OBJEXT) \
++#	rates.$(OBJEXT) \
++#	spr.$(OBJEXT) \
++#	draw.$(OBJEXT) \
++#	stats.$(OBJEXT) \
++#	tiporder.$(OBJEXT) \
++#	m4.$(OBJEXT)
++m4_OBJECTS = $(am_m4_OBJECTS)
++m4_DEPENDENCIES =
++am__part_SOURCES_DIST = main.c utilities.c utilities.h optimiz.c \
++	optimiz.h lk.c lk.h bionj.c bionj.h models.c models.h free.c \
++	free.h help.c help.h simu.c simu.h eigen.c eigen.h pars.c \
++	pars.h alrt.c alrt.h interface.c interface.h cl.c cl.h mg.c \
++	mg.h spr.c spr.h m4.c m4.h draw.c draw.h stats.c stats.h \
++	tiporder.c tiporder.h
++#am_part_OBJECTS = main.$(OBJEXT) \
++#	utilities.$(OBJEXT) \
++#	optimiz.$(OBJEXT) \
++#	lk.$(OBJEXT) \
++#	bionj.$(OBJEXT) \
++#	models.$(OBJEXT) \
++#	free.$(OBJEXT) \
++#	help.$(OBJEXT) \
++#	simu.$(OBJEXT) \
++#	eigen.$(OBJEXT) \
++#	pars.$(OBJEXT) \
++#	alrt.$(OBJEXT) \
++#	interface.$(OBJEXT) \
++#	cl.$(OBJEXT) \
++#	mg.$(OBJEXT) \
++#	spr.$(OBJEXT) \
++#	m4.$(OBJEXT) \
++#	draw.$(OBJEXT) \
++#	stats.$(OBJEXT) \
++#	tiporder.$(OBJEXT)
++part_OBJECTS = $(am_part_OBJECTS)
++part_DEPENDENCIES =
++am__phycont_SOURCES_DIST = main.c utilities.c utilities.h optimiz.c \
++	optimiz.h lk.c lk.h bionj.c bionj.h models.c models.h free.c \
++	free.h help.c help.h simu.c simu.h eigen.c eigen.h pars.c \
++	pars.h alrt.c alrt.h interface.c interface.h cl.c cl.h spr.c \
++	spr.h times.c times.h m4.c m4.h draw.c draw.h rates.c rates.h \
++	mcmc.c mcmc.h stats.c stats.h mg.c mg.h tiporder.c tiporder.h
++#am_phycont_OBJECTS =  \
++#	main.$(OBJEXT) \
++#	utilities.$(OBJEXT) \
++#	optimiz.$(OBJEXT) \
++#	lk.$(OBJEXT) \
++#	bionj.$(OBJEXT) \
++#	models.$(OBJEXT) \
++#	free.$(OBJEXT) \
++#	help.$(OBJEXT) \
++#	simu.$(OBJEXT) \
++#	eigen.$(OBJEXT) \
++#	pars.$(OBJEXT) \
++#	alrt.$(OBJEXT) \
++#	interface.$(OBJEXT) \
++#	cl.$(OBJEXT) \
++#	spr.$(OBJEXT) \
++#	times.$(OBJEXT) \
++#	m4.$(OBJEXT) \
++#	draw.$(OBJEXT) \
++#	rates.$(OBJEXT) \
++#	mcmc.$(OBJEXT) \
++#	stats.$(OBJEXT) \
++#	mg.$(OBJEXT) \
++#	tiporder.$(OBJEXT)
++phycont_OBJECTS = $(am_phycont_OBJECTS)
++phycont_DEPENDENCIES =
++am__phyml_SOURCES_DIST = main.c utilities.c utilities.h optimiz.c \
++	optimiz.h lk.c lk.h bionj.c bionj.h models.c models.h free.c \
++	free.h help.c help.h simu.c simu.h eigen.c eigen.h pars.c \
++	pars.h alrt.c alrt.h interface.c interface.h cl.c cl.h spr.c \
++	spr.h draw.c draw.h stats.c stats.h rates.c rates.h mcmc.c \
++	mcmc.h times.c times.h tiporder.c tiporder.h mg.c mg.h m4.c \
++	m4.h
++am_phyml_OBJECTS = main.$(OBJEXT) \
++	utilities.$(OBJEXT) \
++	optimiz.$(OBJEXT) \
++	lk.$(OBJEXT) \
++	bionj.$(OBJEXT) \
++	models.$(OBJEXT) \
++	free.$(OBJEXT) \
++	help.$(OBJEXT) \
++	simu.$(OBJEXT) \
++	eigen.$(OBJEXT) \
++	pars.$(OBJEXT) \
++	alrt.$(OBJEXT) \
++	interface.$(OBJEXT) \
++	cl.$(OBJEXT) \
++	spr.$(OBJEXT) \
++	draw.$(OBJEXT) \
++	stats.$(OBJEXT) \
++	rates.$(OBJEXT) \
++	mcmc.$(OBJEXT) \
++	times.$(OBJEXT) \
++	tiporder.$(OBJEXT) \
++	mg.$(OBJEXT) \
++	m4.$(OBJEXT)
++phyml_OBJECTS = $(am_phyml_OBJECTS)
++phyml_DEPENDENCIES =
++am__phyml_mpi_SOURCES_DIST = main.c utilities.c utilities.h optimiz.c \
++	optimiz.h lk.c lk.h bionj.c bionj.h models.c models.h free.c \
++	free.h help.c help.h simu.c simu.h eigen.c eigen.h pars.c \
++	pars.h alrt.c alrt.h interface.c interface.h cl.c cl.h spr.c \
++	spr.h draw.c draw.h mpi_boot.c mpi_boot.h stats.c stats.h \
++	rates.c rates.h mcmc.c mcmc.h times.c times.h tiporder.c \
++	tiporder.h mg.c mg.h m4.c m4.h xml.c xml.h
++#am_phyml_mpi_OBJECTS = main.$(OBJEXT) \
++#	utilities.$(OBJEXT) \
++#	optimiz.$(OBJEXT) \
++#	lk.$(OBJEXT) \
++#	bionj.$(OBJEXT) \
++#	models.$(OBJEXT) \
++#	free.$(OBJEXT) \
++#	help.$(OBJEXT) \
++#	simu.$(OBJEXT) \
++#	eigen.$(OBJEXT) \
++#	pars.$(OBJEXT) \
++#	alrt.$(OBJEXT) \
++#	interface.$(OBJEXT) \
++#	cl.$(OBJEXT) \
++#	spr.$(OBJEXT) \
++#	draw.$(OBJEXT) \
++#	mpi_boot.$(OBJEXT) \
++#	stats.$(OBJEXT) \
++#	rates.$(OBJEXT) \
++#	mcmc.$(OBJEXT) \
++#	times.$(OBJEXT) \
++#	tiporder.$(OBJEXT) \
++#	mg.$(OBJEXT) \
++#	m4.$(OBJEXT) \
++#	xml.$(OBJEXT)
++phyml_mpi_OBJECTS = $(am_phyml_mpi_OBJECTS)
++phyml_mpi_DEPENDENCIES =
++am__phytime_SOURCES_DIST = main.c utilities.c utilities.h optimiz.c \
++	optimiz.h lk.c lk.h bionj.c bionj.h models.c models.h free.c \
++	free.h help.c help.h simu.c simu.h eigen.c eigen.h pars.c \
++	pars.h alrt.c alrt.h interface.c interface.h cl.c cl.h spr.c \
++	spr.h times.c times.h m4.c m4.h draw.c draw.h rates.c rates.h \
++	mcmc.c mcmc.h stats.c stats.h mg.c mg.h tiporder.c tiporder.h
++#am_phytime_OBJECTS = main.$(OBJEXT) \
++#	utilities.$(OBJEXT) optimiz.$(OBJEXT) \
++#	lk.$(OBJEXT) bionj.$(OBJEXT) \
++#	models.$(OBJEXT) free.$(OBJEXT) \
++#	help.$(OBJEXT) simu.$(OBJEXT) \
++#	eigen.$(OBJEXT) pars.$(OBJEXT) \
++#	alrt.$(OBJEXT) interface.$(OBJEXT) \
++#	cl.$(OBJEXT) spr.$(OBJEXT) times.$(OBJEXT) \
++#	m4.$(OBJEXT) draw.$(OBJEXT) rates.$(OBJEXT) \
++#	mcmc.$(OBJEXT) stats.$(OBJEXT) mg.$(OBJEXT) \
++#	tiporder.$(OBJEXT)
++phytime_OBJECTS = $(am_phytime_OBJECTS)
++phytime_DEPENDENCIES =
++am__rf_SOURCES_DIST = main.c utilities.c utilities.h optimiz.c \
++	optimiz.h lk.c lk.h bionj.c bionj.h models.c models.h free.c \
++	free.h help.c help.h simu.c simu.h eigen.c eigen.h pars.c \
++	pars.h alrt.c alrt.h interface.c interface.h cl.c cl.h mg.c \
++	mg.h times.c times.h mcmc.c mcmc.h rates.c rates.h spr.c spr.h \
++	draw.c draw.h stats.c stats.h tiporder.c tiporder.h m4.c m4.h
++#am_rf_OBJECTS = main.$(OBJEXT) \
++#	utilities.$(OBJEXT) \
++#	optimiz.$(OBJEXT) \
++#	lk.$(OBJEXT) \
++#	bionj.$(OBJEXT) \
++#	models.$(OBJEXT) \
++#	free.$(OBJEXT) \
++#	help.$(OBJEXT) \
++#	simu.$(OBJEXT) \
++#	eigen.$(OBJEXT) \
++#	pars.$(OBJEXT) \
++#	alrt.$(OBJEXT) \
++#	interface.$(OBJEXT) \
++#	cl.$(OBJEXT) \
++#	mg.$(OBJEXT) \
++#	times.$(OBJEXT) \
++#	mcmc.$(OBJEXT) \
++#	rates.$(OBJEXT) \
++#	spr.$(OBJEXT) \
++#	draw.$(OBJEXT) \
++#	stats.$(OBJEXT) \
++#	tiporder.$(OBJEXT) \
++#	m4.$(OBJEXT)
++rf_OBJECTS = $(am_rf_OBJECTS)
++rf_DEPENDENCIES =
++am__test_SOURCES_DIST = main.c utilities.c utilities.h optimiz.c \
++	optimiz.h lk.c lk.h bionj.c bionj.h models.c models.h free.c \
++	free.h help.c help.h simu.c simu.h eigen.c eigen.h pars.c \
++	pars.h alrt.c alrt.h interface.c interface.h cl.c cl.h spr.c \
++	spr.h draw.c draw.h stats.c stats.h rates.c rates.h mcmc.c \
++	mcmc.h times.c times.h tiporder.c tiporder.h mg.c mg.h m4.c \
++	m4.h xml.c xml.h
++#am_test_OBJECTS = main.$(OBJEXT) \
++#	utilities.$(OBJEXT) \
++#	optimiz.$(OBJEXT) \
++#	lk.$(OBJEXT) \
++#	bionj.$(OBJEXT) \
++#	models.$(OBJEXT) \
++#	free.$(OBJEXT) \
++#	help.$(OBJEXT) \
++#	simu.$(OBJEXT) \
++#	eigen.$(OBJEXT) \
++#	pars.$(OBJEXT) \
++#	alrt.$(OBJEXT) \
++#	interface.$(OBJEXT) \
++#	cl.$(OBJEXT) \
++#	spr.$(OBJEXT) \
++#	draw.$(OBJEXT) \
++#	stats.$(OBJEXT) \
++#	rates.$(OBJEXT) \
++#	mcmc.$(OBJEXT) \
++#	times.$(OBJEXT) \
++#	tiporder.$(OBJEXT) \
++#	mg.$(OBJEXT) \
++#	m4.$(OBJEXT) \
++#	xml.$(OBJEXT)
++test_OBJECTS = $(am_test_OBJECTS)
++test_DEPENDENCIES =
++am__tiporder_SOURCES_DIST = main.c utilities.c utilities.h optimiz.c \
++	optimiz.h lk.c lk.h bionj.c bionj.h models.c models.h free.c \
++	free.h help.c help.h simu.c simu.h eigen.c eigen.h pars.c \
++	pars.h alrt.c alrt.h interface.c interface.h cl.c cl.h mg.c \
++	mg.h times.c times.h mcmc.c mcmc.h rates.c rates.h spr.c spr.h \
++	m4.c m4.h draw.c draw.h stats.c stats.h tiporder.c tiporder.h
++#am_tiporder_OBJECTS = main.$(OBJEXT) \
++#	utilities.$(OBJEXT) \
++#	optimiz.$(OBJEXT) \
++#	lk.$(OBJEXT) \
++#	bionj.$(OBJEXT) \
++#	models.$(OBJEXT) \
++#	free.$(OBJEXT) \
++#	help.$(OBJEXT) \
++#	simu.$(OBJEXT) \
++#	eigen.$(OBJEXT) \
++#	pars.$(OBJEXT) \
++#	alrt.$(OBJEXT) \
++#	interface.$(OBJEXT) \
++#	cl.$(OBJEXT) \
++#	mg.$(OBJEXT) \
++#	times.$(OBJEXT) \
++#	mcmc.$(OBJEXT) \
++#	rates.$(OBJEXT) \
++#	spr.$(OBJEXT) \
++#	m4.$(OBJEXT) \
++#	draw.$(OBJEXT) \
++#	stats.$(OBJEXT) \
++#	tiporder.$(OBJEXT)
++tiporder_OBJECTS = $(am_tiporder_OBJECTS)
++tiporder_DEPENDENCIES =
++DEFAULT_INCLUDES = -I. -I$(top_builddir)
++depcomp = $(SHELL) $(top_srcdir)/depcomp
++am__depfiles_maybe = depfiles
++am__mv = mv -f
++COMPILE = $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) \
++	$(CPPFLAGS) $(AM_CFLAGS) $(CFLAGS)
++AM_V_CC = $(am__v_CC_$(V))
++am__v_CC_ = $(am__v_CC_$(AM_DEFAULT_VERBOSITY))
++am__v_CC_0 = @echo "  CC    " $@;
++AM_V_at = $(am__v_at_$(V))
++am__v_at_ = $(am__v_at_$(AM_DEFAULT_VERBOSITY))
++am__v_at_0 = @
++CCLD = $(CC)
++LINK = $(CCLD) $(AM_CFLAGS) $(CFLAGS) $(AM_LDFLAGS) $(LDFLAGS) -o $@
++AM_V_CCLD = $(am__v_CCLD_$(V))
++am__v_CCLD_ = $(am__v_CCLD_$(AM_DEFAULT_VERBOSITY))
++am__v_CCLD_0 = @echo "  CCLD  " $@;
++AM_V_GEN = $(am__v_GEN_$(V))
++am__v_GEN_ = $(am__v_GEN_$(AM_DEFAULT_VERBOSITY))
++am__v_GEN_0 = @echo "  GEN   " $@;
++SOURCES = $(m4_SOURCES) $(part_SOURCES) $(phycont_SOURCES) \
++	$(phyml_SOURCES) $(phyml_mpi_SOURCES) $(phytime_SOURCES) \
++	$(rf_SOURCES) $(test_SOURCES) $(tiporder_SOURCES)
++DIST_SOURCES = $(am__m4_SOURCES_DIST) $(am__part_SOURCES_DIST) \
++	$(am__phycont_SOURCES_DIST) $(am__phyml_SOURCES_DIST) \
++	$(am__phyml_mpi_SOURCES_DIST) $(am__phytime_SOURCES_DIST) \
++	$(am__rf_SOURCES_DIST) $(am__test_SOURCES_DIST) \
++	$(am__tiporder_SOURCES_DIST)
++ETAGS = etags
++CTAGS = ctags
++DISTFILES = $(DIST_COMMON) $(DIST_SOURCES) $(TEXINFOS) $(EXTRA_DIST)
++ACLOCAL = ${SHELL} /home/guindon/Downloads/phyml-20120412/missing --run aclocal-1.11
++AMTAR = ${SHELL} /home/guindon/Downloads/phyml-20120412/missing --run tar
++AM_DEFAULT_VERBOSITY = 1
++AUTOCONF = ${SHELL} /home/guindon/Downloads/phyml-20120412/missing --run autoconf
++AUTOHEADER = ${SHELL} /home/guindon/Downloads/phyml-20120412/missing --run autoheader
++AUTOMAKE = ${SHELL} /home/guindon/Downloads/phyml-20120412/missing --run automake-1.11
++AWK = gawk
++CC = gcc
++CCDEPMODE = depmode=gcc3
++CFLAGS = -ansi -pedantic -Wall -std=c99 -O0 -g
++CPP = gcc -E
++CPPFLAGS = 
++CYGPATH_W = echo
++DEFS = $(REVISION)
++DEPDIR = .deps
++ECHO_C = 
++ECHO_N = -n
++ECHO_T = 
++EGREP = /usr/bin/grep -E
++EXEEXT = 
++GREP = /usr/bin/grep
++INSTALL = /usr/bin/install -c
++INSTALL_DATA = ${INSTALL} -m 644
++INSTALL_PROGRAM = ${INSTALL}
++INSTALL_SCRIPT = ${INSTALL}
++INSTALL_STRIP_PROGRAM = $(install_sh) -c -s
++LDFLAGS = 
++LIBOBJS = 
++LIBS = -lm 
++LTLIBOBJS = 
++MAKEINFO = ${SHELL} /home/guindon/Downloads/phyml-20120412/missing --run makeinfo
++MKDIR_P = /usr/bin/mkdir -p
++OBJEXT = o
++PACKAGE = phyml
++PACKAGE_BUGREPORT = s.guindon@auckland.ac.nz
++PACKAGE_NAME = PhyML
++PACKAGE_STRING = PhyML 20120412
++PACKAGE_TARNAME = phyml
++PACKAGE_URL = 
++PACKAGE_VERSION = 20120412
++PATH_SEPARATOR = :
++SET_MAKE = 
++SHELL = /bin/sh
++STRIP = 
++VERSION = 20120412
++abs_builddir = /home/guindon/Downloads/phyml-20120412/src
++abs_srcdir = /home/guindon/Downloads/phyml-20120412/src
++abs_top_builddir = /home/guindon/Downloads/phyml-20120412
++abs_top_srcdir = /home/guindon/Downloads/phyml-20120412
++ac_ct_CC = gcc
++am__include = include
++am__leading_dot = .
++am__quote = 
++am__tar = ${AMTAR} chof - "$$tardir"
++am__untar = ${AMTAR} xf -
++bindir = ${exec_prefix}/bin
++build = x86_64-unknown-linux-gnu
++build_alias = 
++build_cpu = x86_64
++build_os = linux-gnu
++build_vendor = unknown
++builddir = .
++datadir = ${datarootdir}
++datarootdir = ${prefix}/share
++docdir = ${datarootdir}/doc/${PACKAGE_TARNAME}
++dvidir = ${docdir}
++exec_prefix = ${prefix}
++host = x86_64-unknown-linux-gnu
++host_alias = 
++host_cpu = x86_64
++host_os = linux-gnu
++host_vendor = unknown
++htmldir = ${docdir}
++includedir = ${prefix}/include
++infodir = ${datarootdir}/info
++install_sh = ${SHELL} /home/guindon/Downloads/phyml-20120412/install-sh
++libdir = ${exec_prefix}/lib
++libexecdir = ${exec_prefix}/libexec
++localedir = ${datarootdir}/locale
++localstatedir = ${prefix}/var
++mandir = ${datarootdir}/man
++mkdir_p = /usr/bin/mkdir -p
++oldincludedir = /usr/include
++pdfdir = ${docdir}
++prefix = /usr/local
++program_transform_name = s,x,x,
++psdir = ${docdir}
++sbindir = ${exec_prefix}/sbin
++sharedstatedir = ${prefix}/com
++srcdir = .
++sysconfdir = ${prefix}/etc
++target_alias = 
++top_build_prefix = ../
++top_builddir = ..
++top_srcdir = ..
++PROG = PHYML
++#PROG = TEST
++#PROG = PHYML
++#PROG = RF
++#PROG = M4
++#PROG = TIPORDER
++#PROG = RWRAP
++#PROG = PART
++#PROG = PHYCONT
++#PROG = PHYTIME
++#phytime_SOURCES = main.c \
++#utilities.c utilities.h\
++#optimiz.c optimiz.h\
++#lk.c lk.h\
++#bionj.c bionj.h\
++#models.c models.h\
++#free.c free.h\
++#help.c help.h\
++#simu.c simu.h\
++#eigen.c eigen.h\
++#pars.c pars.h\
++#alrt.c alrt.h\
++#interface.c interface.h\
++#cl.c cl.h\
++#spr.c spr.h\
++#times.c times.h\
++#m4.c m4.h\
++#draw.c draw.h\
++#rates.c rates.h\
++#mcmc.c mcmc.h\
++#stats.c stats.h\
++#mg.c mg.h\
++#tiporder.c tiporder.h
++
++#phytime_LDADD = -lm
++#phycont_SOURCES = main.c \
++#utilities.c utilities.h\
++#optimiz.c optimiz.h\
++#lk.c lk.h\
++#bionj.c bionj.h\
++#models.c models.h\
++#free.c free.h\
++#help.c help.h\
++#simu.c simu.h\
++#eigen.c eigen.h\
++#pars.c pars.h\
++#alrt.c alrt.h\
++#interface.c interface.h\
++#cl.c cl.h\
++#spr.c spr.h\
++#times.c times.h\
++#m4.c m4.h\
++#draw.c draw.h\
++#rates.c rates.h\
++#mcmc.c mcmc.h\
++#stats.c stats.h\
++#mg.c mg.h\
++#tiporder.c tiporder.h
++
++# continuous.c continuous.h
++#phycont_LDADD = -lm
++# if WANT_RWRAP
++# lib_LTLIBRARIES = librwrap.la
++# librwrap_la_SOURCES = main.c \
++# utilities.c utilities.h\
++# optimiz.c optimiz.h\
++# lk.c lk.h\
++# bionj.c bionj.h\
++# models.c models.h\
++# free.c free.h\
++# help.c help.h\
++# simu.c simu.h\
++# eigen.c eigen.h\
++# pars.c pars.h\
++# alrt.c alrt.h\
++# interface.c interface.h\
++# cl.c cl.h\
++# mg.c mg.h\
++# spr.c spr.h\
++# m4.c m4.h\
++# draw.c draw.h\
++# mcmc.c mcmc.h\
++# stats.c stats.h\
++# times.c times.h\
++# tiporder.c tiporder.h\
++# rates.c rates.h\
++# rwrapper.c rwrapper.h
++# librwrap_la_LIBADD = -lm
++# librwrap_la_LDFLAGS = -I/Library/Frameworks/R.framework/Resources/include -I/Library/Frameworks/R.framework/Resources/include/i386 -I/usr/local/include -flat_namespace -undefined suppress -shared -module
++# librwrap_la_CFLAGS=-std=gnu99 -fPIC
++# else 
++#part_SOURCES = main.c \
++#utilities.c utilities.h\
++#optimiz.c optimiz.h\
++#lk.c lk.h\
++#bionj.c bionj.h\
++#models.c models.h\
++#free.c free.h\
++#help.c help.h\
++#simu.c simu.h\
++#eigen.c eigen.h\
++#pars.c pars.h\
++#alrt.c alrt.h\
++#interface.c interface.h\
++#cl.c cl.h\
++#mg.c mg.h\
++#spr.c spr.h\
++#m4.c m4.h\
++#draw.c draw.h\
++#stats.c stats.h\
++#tiporder.c tiporder.h
++
++#part_LDADD = -lm
++#tiporder_SOURCES = main.c \
++#utilities.c utilities.h\
++#optimiz.c optimiz.h\
++#lk.c lk.h\
++#bionj.c bionj.h\
++#models.c models.h\
++#free.c free.h\
++#help.c help.h\
++#simu.c simu.h\
++#eigen.c eigen.h\
++#pars.c pars.h\
++#alrt.c alrt.h\
++#interface.c interface.h\
++#cl.c cl.h\
++#mg.c mg.h\
++#times.c times.h\
++#mcmc.c mcmc.h\
++#rates.c rates.h\
++#spr.c spr.h\
++#m4.c m4.h\
++#draw.c draw.h\
++#stats.c stats.h\
++#tiporder.c tiporder.h
++
++#tiporder_LDADD = -lm
++#m4_SOURCES = main.c \
++#utilities.c utilities.h\
++#optimiz.c optimiz.h\
++#lk.c lk.h\
++#bionj.c bionj.h\
++#models.c models.h\
++#free.c free.h\
++#help.c help.h\
++#simu.c simu.h\
++#eigen.c eigen.h\
++#pars.c pars.h\
++#alrt.c alrt.h\
++#interface.c interface.h\
++#cl.c cl.h\
++#mg.c mg.h\
++#times.c times.h\
++#mcmc.c mcmc.h\
++#rates.c rates.h\
++#spr.c spr.h\
++#draw.c draw.h\
++#stats.c stats.h\
++#tiporder.c tiporder.h\
++#m4.c m4.h
++
++#m4_LDADD = -lm
++#rf_SOURCES = main.c \
++#utilities.c utilities.h\
++#optimiz.c optimiz.h\
++#lk.c lk.h\
++#bionj.c bionj.h\
++#models.c models.h\
++#free.c free.h\
++#help.c help.h\
++#simu.c simu.h\
++#eigen.c eigen.h\
++#pars.c pars.h\
++#alrt.c alrt.h\
++#interface.c interface.h\
++#cl.c cl.h\
++#mg.c mg.h\
++#times.c times.h\
++#mcmc.c mcmc.h\
++#rates.c rates.h\
++#spr.c spr.h\
++#draw.c draw.h\
++#stats.c stats.h\
++#tiporder.c tiporder.h\
++#m4.c m4.h
++
++#rf_LDADD = -lm
++#phyml_mpi_SOURCES = main.c \
++#utilities.c  utilities.h\
++#optimiz.c  optimiz.h\
++#lk.c  lk.h\
++#bionj.c bionj.h\
++#models.c  models.h\
++#free.c  free.h\
++#help.c  help.h\
++#simu.c  simu.h\
++#eigen.c  eigen.h\
++#pars.c  pars.h\
++#alrt.c  alrt.h\
++#interface.c  interface.h\
++#cl.c  cl.h\
++#spr.c spr.h\
++#draw.c  draw.h\
++#mpi_boot.c  mpi_boot.h\
++#stats.c stats.h\
++#rates.c rates.h\
++#mcmc.c mcmc.h\
++#times.c times.h\
++#tiporder.c tiporder.h\
++#mg.c mg.h\
++#m4.c m4.h\
++#xml.c xml.h
++
++#phyml_mpi_LDADD = -lm
++#test_SOURCES = main.c \
++#utilities.c  utilities.h\
++#optimiz.c  optimiz.h\
++#lk.c  lk.h\
++#bionj.c bionj.h\
++#models.c  models.h\
++#free.c  free.h\
++#help.c  help.h\
++#simu.c  simu.h\
++#eigen.c  eigen.h\
++#pars.c  pars.h\
++#alrt.c  alrt.h\
++#interface.c  interface.h\
++#cl.c  cl.h\
++#spr.c spr.h\
++#draw.c  draw.h\
++#stats.c stats.h\
++#rates.c rates.h\
++#mcmc.c mcmc.h\
++#times.c times.h\
++#tiporder.c tiporder.h\
++#mg.c mg.h\
++#m4.c m4.h\
++#xml.c xml.h
++
++#test_LDADD = -lm
++phyml_SOURCES = main.c \
++utilities.c  utilities.h\
++optimiz.c  optimiz.h\
++lk.c  lk.h\
++bionj.c bionj.h\
++models.c  models.h\
++free.c  free.h\
++help.c  help.h\
++simu.c  simu.h\
++eigen.c  eigen.h\
++pars.c  pars.h\
++alrt.c  alrt.h\
++interface.c  interface.h\
++cl.c  cl.h\
++spr.c spr.h\
++draw.c  draw.h\
++stats.c stats.h\
++rates.c rates.h\
++mcmc.c mcmc.h\
++times.c times.h\
++tiporder.c tiporder.h\
++mg.c mg.h\
++m4.c m4.h
++
++phyml_LDADD = -lm 
++all: all-am
++
++.SUFFIXES:
++.SUFFIXES: .c .o .obj
++$(srcdir)/Makefile.in:  $(srcdir)/Makefile.am  $(am__configure_deps)
++	@for dep in $?; do \
++	  case '$(am__configure_deps)' in \
++	    *$$dep*) \
++	      ( cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh ) \
++	        && { if test -f $@; then exit 0; else break; fi; }; \
++	      exit 1;; \
++	  esac; \
++	done; \
++	echo ' cd $(top_srcdir) && $(AUTOMAKE) --foreign src/Makefile'; \
++	$(am__cd) $(top_srcdir) && \
++	  $(AUTOMAKE) --foreign src/Makefile
++.PRECIOUS: Makefile
++Makefile: $(srcdir)/Makefile.in $(top_builddir)/config.status
++	@case '$?' in \
++	  *config.status*) \
++	    cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh;; \
++	  *) \
++	    echo ' cd $(top_builddir) && $(SHELL) ./config.status $(subdir)/$@ $(am__depfiles_maybe)'; \
++	    cd $(top_builddir) && $(SHELL) ./config.status $(subdir)/$@ $(am__depfiles_maybe);; \
++	esac;
++
++$(top_builddir)/config.status: $(top_srcdir)/configure $(CONFIG_STATUS_DEPENDENCIES)
++	cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh
++
++$(top_srcdir)/configure:  $(am__configure_deps)
++	cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh
++$(ACLOCAL_M4):  $(am__aclocal_m4_deps)
++	cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh
++$(am__aclocal_m4_deps):
++install-binPROGRAMS: $(bin_PROGRAMS)
++	@$(NORMAL_INSTALL)
++	test -z "$(bindir)" || $(MKDIR_P) "$(DESTDIR)$(bindir)"
++	@list='$(bin_PROGRAMS)'; test -n "$(bindir)" || list=; \
++	for p in $$list; do echo "$$p $$p"; done | \
++	sed 's/$(EXEEXT)$$//' | \
++	while read p p1; do if test -f $$p; \
++	  then echo "$$p"; echo "$$p"; else :; fi; \
++	done | \
++	sed -e 'p;s,.*/,,;n;h' -e 's|.*|.|' \
++	    -e 'p;x;s,.*/,,;s/$(EXEEXT)$$//;$(transform);s/$$/$(EXEEXT)/' | \
++	sed 'N;N;N;s,\n, ,g' | \
++	$(AWK) 'BEGIN { files["."] = ""; dirs["."] = 1 } \
++	  { d=$$3; if (dirs[d] != 1) { print "d", d; dirs[d] = 1 } \
++	    if ($$2 == $$4) files[d] = files[d] " " $$1; \
++	    else { print "f", $$3 "/" $$4, $$1; } } \
++	  END { for (d in files) print "f", d, files[d] }' | \
++	while read type dir files; do \
++	    if test "$$dir" = .; then dir=; else dir=/$$dir; fi; \
++	    test -z "$$files" || { \
++	      echo " $(INSTALL_PROGRAM_ENV) $(INSTALL_PROGRAM) $$files '$(DESTDIR)$(bindir)$$dir'"; \
++	      $(INSTALL_PROGRAM_ENV) $(INSTALL_PROGRAM) $$files "$(DESTDIR)$(bindir)$$dir" || exit $$?; \
++	    } \
++	; done
++
++uninstall-binPROGRAMS:
++	@$(NORMAL_UNINSTALL)
++	@list='$(bin_PROGRAMS)'; test -n "$(bindir)" || list=; \
++	files=`for p in $$list; do echo "$$p"; done | \
++	  sed -e 'h;s,^.*/,,;s/$(EXEEXT)$$//;$(transform)' \
++	      -e 's/$$/$(EXEEXT)/' `; \
++	test -n "$$list" || exit 0; \
++	echo " ( cd '$(DESTDIR)$(bindir)' && rm -f" $$files ")"; \
++	cd "$(DESTDIR)$(bindir)" && rm -f $$files
++
++clean-binPROGRAMS:
++	-test -z "$(bin_PROGRAMS)" || rm -f $(bin_PROGRAMS)
++m4$(EXEEXT): $(m4_OBJECTS) $(m4_DEPENDENCIES) 
++	@rm -f m4$(EXEEXT)
++	$(AM_V_CCLD)$(LINK) $(m4_OBJECTS) $(m4_LDADD) $(LIBS)
++part$(EXEEXT): $(part_OBJECTS) $(part_DEPENDENCIES) 
++	@rm -f part$(EXEEXT)
++	$(AM_V_CCLD)$(LINK) $(part_OBJECTS) $(part_LDADD) $(LIBS)
++phycont$(EXEEXT): $(phycont_OBJECTS) $(phycont_DEPENDENCIES) 
++	@rm -f phycont$(EXEEXT)
++	$(AM_V_CCLD)$(LINK) $(phycont_OBJECTS) $(phycont_LDADD) $(LIBS)
++phyml$(EXEEXT): $(phyml_OBJECTS) $(phyml_DEPENDENCIES) 
++	@rm -f phyml$(EXEEXT)
++	$(AM_V_CCLD)$(LINK) $(phyml_OBJECTS) $(phyml_LDADD) $(LIBS)
++phyml-mpi$(EXEEXT): $(phyml_mpi_OBJECTS) $(phyml_mpi_DEPENDENCIES) 
++	@rm -f phyml-mpi$(EXEEXT)
++	$(AM_V_CCLD)$(LINK) $(phyml_mpi_OBJECTS) $(phyml_mpi_LDADD) $(LIBS)
++phytime$(EXEEXT): $(phytime_OBJECTS) $(phytime_DEPENDENCIES) 
++	@rm -f phytime$(EXEEXT)
++	$(AM_V_CCLD)$(LINK) $(phytime_OBJECTS) $(phytime_LDADD) $(LIBS)
++rf$(EXEEXT): $(rf_OBJECTS) $(rf_DEPENDENCIES) 
++	@rm -f rf$(EXEEXT)
++	$(AM_V_CCLD)$(LINK) $(rf_OBJECTS) $(rf_LDADD) $(LIBS)
++test$(EXEEXT): $(test_OBJECTS) $(test_DEPENDENCIES) 
++	@rm -f test$(EXEEXT)
++	$(AM_V_CCLD)$(LINK) $(test_OBJECTS) $(test_LDADD) $(LIBS)
++tiporder$(EXEEXT): $(tiporder_OBJECTS) $(tiporder_DEPENDENCIES) 
++	@rm -f tiporder$(EXEEXT)
++	$(AM_V_CCLD)$(LINK) $(tiporder_OBJECTS) $(tiporder_LDADD) $(LIBS)
++
++mostlyclean-compile:
++	-rm -f *.$(OBJEXT)
++
++distclean-compile:
++	-rm -f *.tab.c
++
++include ./$(DEPDIR)/alrt.Po
++include ./$(DEPDIR)/bionj.Po
++include ./$(DEPDIR)/cl.Po
++include ./$(DEPDIR)/draw.Po
++include ./$(DEPDIR)/eigen.Po
++include ./$(DEPDIR)/free.Po
++include ./$(DEPDIR)/help.Po
++include ./$(DEPDIR)/interface.Po
++include ./$(DEPDIR)/lk.Po
++include ./$(DEPDIR)/m4.Po
++include ./$(DEPDIR)/main.Po
++include ./$(DEPDIR)/mcmc.Po
++include ./$(DEPDIR)/mg.Po
++include ./$(DEPDIR)/models.Po
++include ./$(DEPDIR)/mpi_boot.Po
++include ./$(DEPDIR)/optimiz.Po
++include ./$(DEPDIR)/pars.Po
++include ./$(DEPDIR)/rates.Po
++include ./$(DEPDIR)/simu.Po
++include ./$(DEPDIR)/spr.Po
++include ./$(DEPDIR)/stats.Po
++include ./$(DEPDIR)/times.Po
++include ./$(DEPDIR)/tiporder.Po
++include ./$(DEPDIR)/utilities.Po
++include ./$(DEPDIR)/xml.Po
++
++.c.o:
++	$(AM_V_CC)$(COMPILE) -MT $@ -MD -MP -MF $(DEPDIR)/$*.Tpo -c -o $@ $<
++	$(AM_V_at)$(am__mv) $(DEPDIR)/$*.Tpo $(DEPDIR)/$*.Po
++#	$(AM_V_CC) \
++#	source='$<' object='$@' libtool=no \
++#	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) \
++#	$(COMPILE) -c $<
++
++.c.obj:
++	$(AM_V_CC)$(COMPILE) -MT $@ -MD -MP -MF $(DEPDIR)/$*.Tpo -c -o $@ `$(CYGPATH_W) '$<'`
++	$(AM_V_at)$(am__mv) $(DEPDIR)/$*.Tpo $(DEPDIR)/$*.Po
++#	$(AM_V_CC) \
++#	source='$<' object='$@' libtool=no \
++#	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) \
++#	$(COMPILE) -c `$(CYGPATH_W) '$<'`
++
++ID: $(HEADERS) $(SOURCES) $(LISP) $(TAGS_FILES)
++	list='$(SOURCES) $(HEADERS) $(LISP) $(TAGS_FILES)'; \
++	unique=`for i in $$list; do \
++	    if test -f "$$i"; then echo $$i; else echo $(srcdir)/$$i; fi; \
++	  done | \
++	  $(AWK) '{ files[$$0] = 1; nonempty = 1; } \
++	      END { if (nonempty) { for (i in files) print i; }; }'`; \
++	mkid -fID $$unique
++tags: TAGS
++
++TAGS:  $(HEADERS) $(SOURCES)  $(TAGS_DEPENDENCIES) \
++		$(TAGS_FILES) $(LISP)
++	set x; \
++	here=`pwd`; \
++	list='$(SOURCES) $(HEADERS)  $(LISP) $(TAGS_FILES)'; \
++	unique=`for i in $$list; do \
++	    if test -f "$$i"; then echo $$i; else echo $(srcdir)/$$i; fi; \
++	  done | \
++	  $(AWK) '{ files[$$0] = 1; nonempty = 1; } \
++	      END { if (nonempty) { for (i in files) print i; }; }'`; \
++	shift; \
++	if test -z "$(ETAGS_ARGS)$$*$$unique"; then :; else \
++	  test -n "$$unique" || unique=$$empty_fix; \
++	  if test $$# -gt 0; then \
++	    $(ETAGS) $(ETAGSFLAGS) $(AM_ETAGSFLAGS) $(ETAGS_ARGS) \
++	      "$$@" $$unique; \
++	  else \
++	    $(ETAGS) $(ETAGSFLAGS) $(AM_ETAGSFLAGS) $(ETAGS_ARGS) \
++	      $$unique; \
++	  fi; \
++	fi
++ctags: CTAGS
++CTAGS:  $(HEADERS) $(SOURCES)  $(TAGS_DEPENDENCIES) \
++		$(TAGS_FILES) $(LISP)
++	list='$(SOURCES) $(HEADERS)  $(LISP) $(TAGS_FILES)'; \
++	unique=`for i in $$list; do \
++	    if test -f "$$i"; then echo $$i; else echo $(srcdir)/$$i; fi; \
++	  done | \
++	  $(AWK) '{ files[$$0] = 1; nonempty = 1; } \
++	      END { if (nonempty) { for (i in files) print i; }; }'`; \
++	test -z "$(CTAGS_ARGS)$$unique" \
++	  || $(CTAGS) $(CTAGSFLAGS) $(AM_CTAGSFLAGS) $(CTAGS_ARGS) \
++	     $$unique
++
++GTAGS:
++	here=`$(am__cd) $(top_builddir) && pwd` \
++	  && $(am__cd) $(top_srcdir) \
++	  && gtags -i $(GTAGS_ARGS) "$$here"
++
++distclean-tags:
++	-rm -f TAGS ID GTAGS GRTAGS GSYMS GPATH tags
++
++distdir: $(DISTFILES)
++	@srcdirstrip=`echo "$(srcdir)" | sed 's/[].[^$$\\*]/\\\\&/g'`; \
++	topsrcdirstrip=`echo "$(top_srcdir)" | sed 's/[].[^$$\\*]/\\\\&/g'`; \
++	list='$(DISTFILES)'; \
++	  dist_files=`for file in $$list; do echo $$file; done | \
++	  sed -e "s|^$$srcdirstrip/||;t" \
++	      -e "s|^$$topsrcdirstrip/|$(top_builddir)/|;t"`; \
++	case $$dist_files in \
++	  */*) $(MKDIR_P) `echo "$$dist_files" | \
++			   sed '/\//!d;s|^|$(distdir)/|;s,/[^/]*$$,,' | \
++			   sort -u` ;; \
++	esac; \
++	for file in $$dist_files; do \
++	  if test -f $$file || test -d $$file; then d=.; else d=$(srcdir); fi; \
++	  if test -d $$d/$$file; then \
++	    dir=`echo "/$$file" | sed -e 's,/[^/]*$$,,'`; \
++	    if test -d "$(distdir)/$$file"; then \
++	      find "$(distdir)/$$file" -type d ! -perm -700 -exec chmod u+rwx {} \;; \
++	    fi; \
++	    if test -d $(srcdir)/$$file && test $$d != $(srcdir); then \
++	      cp -fpR $(srcdir)/$$file "$(distdir)$$dir" || exit 1; \
++	      find "$(distdir)/$$file" -type d ! -perm -700 -exec chmod u+rwx {} \;; \
++	    fi; \
++	    cp -fpR $$d/$$file "$(distdir)$$dir" || exit 1; \
++	  else \
++	    test -f "$(distdir)/$$file" \
++	    || cp -p $$d/$$file "$(distdir)/$$file" \
++	    || exit 1; \
++	  fi; \
++	done
++check-am: all-am
++check: check-am
++all-am: Makefile $(PROGRAMS)
++installdirs:
++	for dir in "$(DESTDIR)$(bindir)"; do \
++	  test -z "$$dir" || $(MKDIR_P) "$$dir"; \
++	done
++install: install-am
++install-exec: install-exec-am
++install-data: install-data-am
++uninstall: uninstall-am
++
++install-am: all-am
++	@$(MAKE) $(AM_MAKEFLAGS) install-exec-am install-data-am
++
++installcheck: installcheck-am
++install-strip:
++	$(MAKE) $(AM_MAKEFLAGS) INSTALL_PROGRAM="$(INSTALL_STRIP_PROGRAM)" \
++	  install_sh_PROGRAM="$(INSTALL_STRIP_PROGRAM)" INSTALL_STRIP_FLAG=-s \
++	  `test -z '$(STRIP)' || \
++	    echo "INSTALL_PROGRAM_ENV=STRIPPROG='$(STRIP)'"` install
++mostlyclean-generic:
++
++clean-generic:
++
++distclean-generic:
++	-test -z "$(CONFIG_CLEAN_FILES)" || rm -f $(CONFIG_CLEAN_FILES)
++	-test . = "$(srcdir)" || test -z "$(CONFIG_CLEAN_VPATH_FILES)" || rm -f $(CONFIG_CLEAN_VPATH_FILES)
++
++maintainer-clean-generic:
++	@echo "This command is intended for maintainers to use"
++	@echo "it deletes files that may require special tools to rebuild."
++clean: clean-am
++
++clean-am: clean-binPROGRAMS clean-generic mostlyclean-am
++
++distclean: distclean-am
++	-rm -rf ./$(DEPDIR)
++	-rm -f Makefile
++distclean-am: clean-am distclean-compile distclean-generic \
++	distclean-tags
++
++dvi: dvi-am
++
++dvi-am:
++
++html: html-am
++
++html-am:
++
++info: info-am
++
++info-am:
++
++install-data-am:
++
++install-dvi: install-dvi-am
++
++install-dvi-am:
++
++install-exec-am: install-binPROGRAMS
++
++install-html: install-html-am
++
++install-html-am:
++
++install-info: install-info-am
++
++install-info-am:
++
++install-man:
++
++install-pdf: install-pdf-am
++
++install-pdf-am:
++
++install-ps: install-ps-am
++
++install-ps-am:
++
++installcheck-am:
++
++maintainer-clean: maintainer-clean-am
++	-rm -rf ./$(DEPDIR)
++	-rm -f Makefile
++maintainer-clean-am: distclean-am maintainer-clean-generic
++
++mostlyclean: mostlyclean-am
++
++mostlyclean-am: mostlyclean-compile mostlyclean-generic
++
++pdf: pdf-am
++
++pdf-am:
++
++ps: ps-am
++
++ps-am:
++
++uninstall-am: uninstall-binPROGRAMS
++
++.MAKE: install-am install-strip
++
++.PHONY: CTAGS GTAGS all all-am check check-am clean clean-binPROGRAMS \
++	clean-generic ctags distclean distclean-compile \
++	distclean-generic distclean-tags distdir dvi dvi-am html \
++	html-am info info-am install install-am install-binPROGRAMS \
++	install-data install-data-am install-dvi install-dvi-am \
++	install-exec install-exec-am install-html install-html-am \
++	install-info install-info-am install-man install-pdf \
++	install-pdf-am install-ps install-ps-am install-strip \
++	installcheck installcheck-am installdirs maintainer-clean \
++	maintainer-clean-generic mostlyclean mostlyclean-compile \
++	mostlyclean-generic pdf pdf-am ps ps-am tags uninstall \
++	uninstall-am uninstall-binPROGRAMS
++
++
++# all-am:	intro $(bin_PROGRAMS)$(EXEEXT)
++all-am:	intro $(bin_PROGRAMS)
++	@echo ""
++	@echo "Done."
++
++intro:	
++	@echo ""
++	@echo ""
++	@echo ":: Building [$(bin_PROGRAMS)]. Version $(VERSION) ::"
++	@echo ""
++	@echo ""
++
++# Tell versions [3.59,3.63) of GNU make to not export all variables.
++# Otherwise a system limit (for SysV at least) may be exceeded.
++.NOEXPORT:
+unchanged:
+--- ./phyml-20120412/src/mcmc.c	2012-03-20 14:10:49.000000000 +1300
++++ ./phyml-20120412-new-bak/src/mcmc.c	2013-02-18 21:51:51.459662430 +1300
+@@ -1977,7 +1977,7 @@
+       if(tree->mod->n_catg > 1)
+ 	{
+ 	  if(tree->mod->free_mixt_rates == NO)
+-	    PhyML_Fprintf(fp,"%G\t",tree->mod->alpha);
++	    PhyML_Fprintf(fp,"%G\t",tree->mod->alpha->v);
+ 	  else
+ 	    {
+ 	      For(i,tree->mod->n_catg) PhyML_Fprintf(fp,"%G\t",tree->mod->gamma_r_proba->v[i]);
+unchanged:
+--- ./phyml-20120412/src/optimiz.c	2012-03-13 19:04:03.000000000 +1300
++++ ./phyml-20120412-new-bak/src/optimiz.c	2013-02-18 21:51:51.462662427 +1300
+@@ -979,29 +979,18 @@
+ 
+   if((tree->mod->s_opt->opt_free_mixt_rates) && (tree->mod->free_mixt_rates == YES))
+     {
+-      /* int failed; */
++      int failed;
+       int i;
+         
+-      /* failed = 0; */
+-      /* tree->mod->update_eigen = 1; */
+-      /* BFGS(tree,tree->mod->gamma_r_proba_unscaled,tree->mod->n_catg,1.e-5,1.e-5, */
+-      /* 	   &Return_Abs_Lk, */
+-      /* 	   &Num_Derivative_Several_Param, */
+-      /* 	   &Lnsrch_Free_Mixt_Rates,&failed); */
+-
+       if(verbose) Print_Lk(tree,"[Rate class freqs.  ]");
+ 
++      failed = 0;
++      tree->mod->update_eigen = 1;
++      BFGS(tree,tree->mod->gamma_r_proba_unscaled,tree->mod->n_catg,1.e-5,1.e-5,
++      	   &Return_Abs_Lk,
++      	   &Num_Derivative_Several_Param,
++      	   &Lnsrch_Free_Mixt_Rates,&failed);
+ 
+-      /* For(i,tree->mod->n_catg-1) */
+-      /* 	{ */
+-      /* 	  Generic_Brent_Lk(&(tree->mod->gamma_r_proba_unscaled->v[i]), */
+-      /* 			   0., */
+-      /* 			   100, */
+-      /* 			   tree->mod->s_opt->min_diff_lk_global, */
+-      /* 			   tree->mod->s_opt->brent_it_max, */
+-      /* 			   tree->mod->s_opt->quickdirty, */
+-      /* 			   Wrap_Lk,NULL,tree,NULL); */
+-      /* 	} */
+ 
+       tree->mod->gamma_r_proba_unscaled->v[tree->mod->n_catg-1] = 100.;
+       For(i,tree->mod->n_catg-1)
+unchanged:
+--- ./phyml-20120412/src/tree	1970-01-01 12:00:00.000000000 +1200
++++ ./phyml-20120412-new-bak/src/tree	2013-02-19 10:09:36.962892889 +1300
+@@ -0,0 +1 @@
++((((((tax6:0.01,(tax5:0.01,tax1:0.026):0.008):0.003,(tax7:0.003,tax3:0.003):0.01271127):0.01,tax2:0.01):0.01,tax10:0.01):0.01,tax9:0.01):0.01,tax8:0.01,tax4:0.1);
+\ No newline at end of file
+unchanged:
+--- ./phyml-20120412/src/utilities.c	2012-04-11 11:48:30.000000000 +1200
++++ ./phyml-20120412-new-bak/src/utilities.c	2013-02-19 10:16:57.770281203 +1300
+@@ -1668,6 +1668,7 @@
+ 
+   Detect_Align_File_Format(io);
+ 
++
+   switch(io->data_file_format)
+     {
+     case PHYLIP: 
+@@ -1861,7 +1862,7 @@
+     {
+       if(!Get_Token(fp,token)) break;
+ 
+-/*       PhyML_Printf("\n+ Token: '%s' next_token=%d cur_token=%d",token,nxt_token_t,cur_token_t); */
++      PhyML_Printf("\n+ Token: '%s' next_token=%d cur_token=%d",token,nxt_token_t,cur_token_t);
+ 
+       if(token[0] == ';') 
+ 	{
+@@ -1940,7 +1941,7 @@
+ 	}
+       while(c != '"');
+       *token = c;
+-      c = fgetc(fp);
++      /* c = fgetc(fp); */
+       if(c == EOF) return 0;
+       *(token+1) = '\0';
+       return 1;
+@@ -1950,6 +1951,8 @@
+     {
+       Skip_Comment(fp);
+       c = fgetc(fp);
++      *token = c;
++      token++;
+       if(c == EOF) return 0;
+       return 1;
+     }
+@@ -1967,7 +1970,7 @@
+   else if(c == '-') { *token = c; token++; }
+   else
+     {
+-      while(isgraph(c) && c != ';' && c != '-' && c != ',')
++      while(isgraph(c) && c != ';' && c != '-' && c != ',' && c != '=')
+ 	{
+ 	  *(token++) = c;
+ 	  c = fgetc(fp);
+@@ -2308,8 +2311,6 @@
+ 
+ }
+ 
+-
+-
+ //////////////////////////////////////////////////////////////
+ //////////////////////////////////////////////////////////////
+ 
+@@ -7756,7 +7757,7 @@
+ 
+ #ifndef PHYTIME
+   mod->l_min = 1.E-8;
+-  mod->l_max = 100.0;
++  mod->l_max = 10.0;
+ #else
+   mod->l_min = 1.E-8;
+   mod->l_max = 2.0;
+@@ -8049,11 +8050,23 @@
+   t_node **bip1, **bip2;
+   int bip_size1, bip_size2, bip_size;
+   int different,identical;
++  int n_edges;
+ 
+   /* !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!*/
+   /* WARNING: call Match_Tip_Numbers and Get_Bip before using this function. */
+   /* !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!*/
+   
++  if(on_existing_edges_only == YES)
++    {
++      n_edges = 0;
++      For(i,2*tree1->n_otu-3) 
++	if(tree1->t_edges[i]->does_exist && tree2->t_edges[i]->does_exist) n_edges++;
++      n_edges -= tree1->n_otu;
++    }
++  else
++    {
++      n_edges = tree1->n_otu-3;
++    }
+   
+   identical = 0;
+   different = 0;
+@@ -8070,7 +8083,7 @@
+ 	      bip_size2 = MIN(b2->left->bip_size[b2->l_r],b2->rght->bip_size[b2->r_l]);
+ 
+ 	      if(bip_size2 > 1 && ((on_existing_edges_only == YES && b2->does_exist) || (on_existing_edges_only == NO)))
+-		{
++		{		  
+ 		  if(bip_size1 == bip_size2)
+ 		    {
+ 		      bip_size = bip_size1;
+@@ -8139,17 +8152,21 @@
+ 			  b1->bip_score++;
+ 			  b2->bip_score++;
+ 			  identical++;			  
+-			  break;
++			  goto out;
++			}
++		      else 
++			{
++			  different++; // Bipartitions have identical sizes but distinct elements
+ 			}
+-		      else different++; // Bipartitions has identical size but distinct elements
+ 		    }
+ 		  else different++; // Biparition have different sizes
+ 		}
+ 	    }
+ 	}
++    out: ;
+     }
+-  
+-  return different;
++  return n_edges - identical;
++  /* return different; */
+ }
+ 
+ //////////////////////////////////////////////////////////////
+@@ -12242,6 +12259,9 @@
+ 
+   aLRT(tree);
+   
++  Br_Len_Involving_Invar(tree);
++  Rescale_Br_Len_Multiplier_Tree(tree);
++
+   Free(s_tree);
+   s_tree = Write_Tree(tree,NO);
+ 
+@@ -12481,6 +12501,9 @@
+       Warn_And_Exit("");
+     }
+ 
++
++  if(node_list_size == 1) return node_list[0];
++
+   list = (t_node ***)mCalloc(node_list_size,sizeof(t_node **));
+   For(i,node_list_size) list[i] = (t_node **)mCalloc(2*tree->n_otu-1,sizeof(t_node *));
+   size = (int *)mCalloc(node_list_size,sizeof(int));
+unchanged:
+--- ./phyml-20120412/src/utilities.h	2012-04-11 11:39:02.000000000 +1200
++++ ./phyml-20120412-new-bak/src/utilities.h	2013-02-18 21:51:51.452662435 +1300
+@@ -272,6 +272,8 @@
+ /* #define P_LK_LIM_SUP 2.037035976e+90 /\* R: format(2^(+300),digits=10) *\/ */
+ #define  P_LK_LIM_INF   3.054936e-151 /* 2^-500 */
+ #define  P_LK_LIM_SUP   3.273391e+150 /* 2^500 */
++/* #define  P_LK_LIM_INF   1.499697e-241 /\* 2^-800 *\/ */
++/* #define  P_LK_LIM_SUP   6.668014e+240 /\* 2^800 *\/ */
+ 
+ 
+ /*!********************************************************/
+unchanged:
+--- utilities.c	2013-10-31 08:57:09.622430940 +1300
++++ utilities.new.c	2013-10-31 08:55:49.522433717 +1300
+@@ -4769,8 +4769,6 @@
+   t_node *n_e, *n_f;
+ 
+   n_otu                   = t_full->n_otu;
+-  t_empt->n_root          = t_full->n_root;
+-  t_empt->e_root          = t_full->e_root;
+   t_empt->c_lnL_sorted    = t_full->c_lnL_sorted;
+   t_empt->log_site_lk_cat = t_full->log_site_lk_cat;
+   t_empt->cur_site_lk     = t_full->cur_site_lk;
+@@ -6993,6 +6991,7 @@
+ 	  boot_tree->mat = boot_mat;
+ 	}
+ 
++
+       boot_tree->mod                = boot_mod;
+       boot_tree->io                 = tree->io;
+       boot_tree->data               = boot_data;
+@@ -7002,6 +7001,7 @@
+       boot_tree->io->print_site_lnl = 0;
+       boot_tree->io->print_trace    = 0;
+ 
++
+       if((boot_tree->mod->s_opt->random_input_tree) && (boot_tree->mod->s_opt->topo_search == SPR_MOVE)) Random_Tree(boot_tree);
+       Order_Tree_CSeq(boot_tree,boot_data);
+       Share_Lk_Struct(tree,boot_tree);
+@@ -7010,6 +7010,7 @@
+       Fill_Dir_Table(boot_tree);
+       Update_Dirs(boot_tree);
+ 
++
+       if(tree->mod->s_opt->greedy) Init_P_Lk_Tips_Double(boot_tree);
+       else                         Init_P_Lk_Tips_Int(boot_tree);
+       Init_Ui_Tips(boot_tree);
+@@ -13307,11 +13308,15 @@
+       Warn_And_Exit("");
+     }
+ 
++  has_outgrp = NO;
++
++  /* 'Root' node is t_nodes[0] */ 
+   if(strstr(tree->t_nodes[0]->name,"*"))
+     {
+       /* PhyML_Printf("\n. Found outgroup taxon: %s",tree->t_nodes[0]->name); */
+       tree->t_nodes[0]->s_ingrp[0]  = 0;
+       tree->t_nodes[0]->s_outgrp[0] = 1;
++      has_outgrp = YES;
+     }
+   else
+     {
+@@ -13319,7 +13324,6 @@
+       tree->t_nodes[0]->s_outgrp[0] = 0;
+     }
+ 
+-  has_outgrp = NO;
+   Get_Best_Root_Position_Post(tree->t_nodes[0],tree->t_nodes[0]->v[0],&has_outgrp,tree);  
+   Get_Best_Root_Position_Pre(tree->t_nodes[0],tree->t_nodes[0]->v[0],tree);  
+ 
+only in patch2:
+unchanged:
+--- src/utilities.c	2012-04-11 11:48:30.000000000 +1200
++++ src.new/utilities.c	2013-11-12 14:19:11.600689913 +1300
+@@ -13995,7 +13995,7 @@ int Check_Topo_Constraints(t_tree *big_t
+   For(i,2*big_tree_cpy->n_otu-3) big_tree_cpy->t_edges[i]->bip_score = 0;
+   For(i,2*small_tree->n_otu-3) small_tree->t_edges[i]->bip_score = 0;
+ 
+-  diffs = Compare_Bip(small_tree,big_tree_cpy,YES);
++  diffs = Compare_Bip(small_tree,big_tree_cpy,NO);
+ 
+   /* printf("\n"); */
+   /* printf("\n. %s",Write_Tree(big_tree_cpy,NO)); */