newick.cc

Go to the documentation of this file.

#include "ball.h"
#include "node.h"
#include "nodeseq.h"
#include "genealogy.h"
#include "internal.h"
 
#include <R.h>
#include <Rmath.h>
#include <Rdefines.h>
 
#include <cstring>
#include <charconv>
#include <iostream>
 
std::string
double2string
(double value)
{
  char buffer[32];
  auto [ptr, ec] = std::to_chars(buffer,buffer+sizeof(buffer),value);
  if (ec == std::errc()) {
    return string_t(buffer,ptr);
  } else {
    err("error in %s",__func__); // #nocov
  }
}

string_t
ball_t::newick
(const slate_t &t, bool showdeme) const
{
  assert(color==black);
  string_t o = "[&&PhyloPOMP type=extant";
  if (showdeme && deme() != undeme)
    o += " deme=" + std::to_string(deme());
  o += "]"
    + std::to_string(uniq) +
    ":" + double2string(t);
  return o;
}

string_t
node_t::newick
(const slate_t& tnow, const slate_t& tpar,
 bool showdeme, bool extended) const
{
  string_t o1 = "", o2 = "", o3 = "";
  int n = nchildren();
  if (n > 0) {
    o1 = "("; o3 = ")";
  }
  if (extended) {
    o3 += "[&&PhyloPOMP ";
    if (holds(blue))
      o3 += "type=sample";
    else if (holds_own())
      o3 += "type=root";
    else
      o3 += "type=node";
    if (showdeme && deme() != undeme)
      o3 += " deme=" + std::to_string(deme());
    o3 += "]";
  }
  n = 0;
  for (ball_t *b : *this) {
    node_t *p = 0;
    switch (b->color) {
    case green:
      p = b->child();
      if (p != this) {
        if (n++ > 0) o2 += ",";
        o2 += p->newick(tnow,slate,showdeme,extended);
      }
      break;
    case black:
      assert(extended);
      if (n++ > 0) o2 += ",";
      o2 += b->newick(tnow-slate,showdeme);
      break;
    case blue:
      break;
    }
  }
  return o1 + o2 + o3
    + std::to_string(uniq)
    + ":" + double2string(slate - tpar);
}

string_t
nodeseq_t::newick
(slate_t t, slate_t te, bool showdeme, bool extended) const
{
  string_t o = "";
  for (node_t *p : *this) {
    if (p->is_root()) {
      o += p->newick(t,te,showdeme,extended) + ";";
    }
  }
  return o;
}

string_t
genealogy_t::newick
(bool extended) const
{
  return nodeseq_t::newick(time(),timezero(),(ndeme() > 0),extended);
}
 
extern "C" {

  SEXP newick (SEXP State, SEXP Extended) {
    PROTECT(Extended = AS_LOGICAL(Extended));
    bool extended = *LOGICAL(Extended);
    genealogy_t A(State);
    UNPROTECT(1);
    return mkString(A.newick(extended).c_str());
  }
 
}