- private final PhylogenyNode getExternalPhylogenyNode( final int i ) {
- return _external_nodes[ _mappings[ i ] ];
+ private final void calculateDistancesFromNewNode( final int otu1, final int otu2, final double d ) {
+ final int m_otu1 = _mappings[ otu1 ];
+ final int m_otu2 = _mappings[ otu2 ];
+ for( int i = 0; i < _n; ++i ) {
+ if ( ( i == otu1 ) || ( i == otu2 ) ) {
+ continue;
+ }
+ final int m_i = _mappings[ i ];
+ if ( otu1 < i ) {
+ if ( otu2 > i ) {
+ _d_values[ m_otu1 ][ m_i ] = ( _d_values[ m_otu1 ][ m_i ] + _d_values[ m_i ][ m_otu2 ] - d ) / 2;
+ }
+ else {
+ _d_values[ m_otu1 ][ m_i ] = ( _d_values[ m_otu1 ][ m_i ] + _d_values[ m_otu2 ][ m_i ] - d ) / 2;
+ }
+ }
+ else {
+ if ( otu2 > i ) {
+ _d_values[ m_i ][ m_otu1 ] = ( _d_values[ m_i ][ m_otu1 ] + _d_values[ m_i ][ m_otu2 ] - d ) / 2;
+ }
+ else {
+ _d_values[ m_i ][ m_otu1 ] = ( _d_values[ m_i ][ m_otu1 ] + _d_values[ m_otu2 ][ m_i ] - d ) / 2;
+ }
+ }
+ }