nojhan/paradiseo
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Pack and unpack functions

Browse source 
git-svn-id: svn://scm.gforge.inria.fr/svnroot/paradiseo@1394 331e1502-861f-0410-8da2-ba01fb791d7f
This commit is contained in:
wcancino 2009-02-03 14:33:04 +00:00
commit 75675fb66b
2 changed files with 90 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

69
contribution/branches/PhyloMOEA/PhyloMOEA/PhyloMOEO_packunpack.cpp Normal file
View file

@ -0,0 +1,69 @@
//
// C++ Implementation: PhyloMOEO_packunpack
//
// Description:
//
//
// Author: <>, (C) 2009
//
// Copyright: See COPYING file that comes with this distribution
//
//
#include <PhyloMOEO_packunpack.h>
void pack( PhyloMOEO & ind)
{
cout << "packing individual" << endl;
phylotreeIND & tree = ind.get_tree();
string s = tree.newick_traverse2( false, false);
::pack(s);
edge edgeaux = *( tree.taxon_number(0).inout_edges_begin() );
node root_traverse = edgeaux.opposite( tree.taxon_number(0) );
postorder_Iterator it = tree.postorder_begin( root_traverse );
postorder_Iterator it2 = tree.postorder_end( root_traverse );
double *blens = new double[(2*tree.number_of_taxons() -3 )];
int top = tree.TREE.number_of_edges() -1;
// reverse order ..... for GTL internals
for( ; it!=it2; top--)
{
if( *it != root_traverse )
blens[top] = tree.get_branch_length( it.branch() );
++it;
}
for(int i=0; i< (2*tree.number_of_taxons() -3 ) ; i++) ::pack(blens[i]);
delete [] blens;
cout << "packing finished..." << endl;
//::pack(blens);
}
void unpack( PhyloMOEO &ind )
{
cout << "unpacking individual" << endl;
phylotreeIND tree = ind.get_tree();
string newickstring;
::unpack(newickstring);
tree.read_newick2( newickstring );
edge edgeaux = *( tree.taxon_number(0).inout_edges_begin() );
node root_traverse = edgeaux.opposite( tree.taxon_number(0) );
postorder_Iterator it = tree.postorder_begin( root_traverse );
postorder_Iterator it2 = tree.postorder_end( root_traverse );
//double blens[2*Tind->number_of_taxons() -3];
//MPI_Unpack(buffer, bufsize, pos, blens, 2*Tind->number_of_taxons() -3, MPI_DOUBLE, com);
// reverse order ..... for GTL internals
while(it!=it2)
{
double blen;
::unpack(blen);
if( *it != root_traverse )
tree.set_branch_length( it.branch(), blen );
++it;
}
cout << "unpacking finished..." << endl;
}

21
contribution/branches/PhyloMOEA/PhyloMOEA/PhyloMOEO_packunpack.h Normal file
View file

@ -0,0 +1,21 @@
//
// C++ Interface: PhyloMOEO_packunpack
//
// Description:
//
//
// Author: <>, (C) 2009
//
// Copyright: See COPYING file that comes with this distribution
//
//
#ifndef PHYLOMOEO_PACKUNPACK_H_
#define PHYLOMOEO_PACKUNPACK_H_
#include <PhyloMOEO.h>
#include <peo>
void pack( PhyloMOEO &);
void unpack( PhyloMOEO & );
#endif