Commit 0522f7b6 authored by Carsten Kemena's avatar Carsten Kemena

fixed nj 'root'

parent 92ddcc8e
Pipeline #767 passed with stages
in 1 minute and 47 seconds
......@@ -240,10 +240,18 @@ PhylogeneticTree<DataType>::nj(Matrix<float> &distMat, const std::vector<std::st
}
}
this->root_.reset(nodes[minJ]);
//root_->parent_ = nullptr;
nodes[minJ]->addChild(nodes[minI]);
nodes[minI]->edgeLength=distMat[minI][minJ];
if (nodes[minJ]->nChildren() > nodes[minI]->nChildren())
{
this->root_.reset(nodes[minJ]);
nodes[minJ]->addChild(nodes[minI]);
nodes[minI]->edgeLength=distMat[minI][minJ];
}
else
{
this->root_.reset(nodes[minI]);
nodes[minI]->addChild(nodes[minJ]);
nodes[minJ]->edgeLength=distMat[minI][minJ];
}
isRooted_=false;
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment