sum.cc

Go to the documentation of this file.

#include "node.h"
#include "genealogy.h"
#include "generics.h"
#include "internal.h"
#include <Rinternals.h>

void
node_t::reuniqify
(name_t shift)
{
  this->uniq += shift;
  for (ball_t *b : *this) {
    b->uniq += shift;
  }
}

void
genealogy_t::reuniqify
(name_t shift)
{
  this->_unique += shift;
  for (node_t *p : *this) p->reuniqify(shift);
}

genealogy_t&
genealogy_t::operator+=
(const genealogy_t& other)
{
  genealogy_t G = other;
  slate_t t = time();
  slate_t t0 = timezero();
  t0 = (t0 > G.timezero()) ? t0 : G.timezero();
  t = (t < G.time()) ? t : G.time();
  G.reuniqify(_unique);
  merge(G,compare);
  timezero() = (timezero() < G.timezero()) ? timezero() : G.timezero();
  time() = (time() > G.time()) ? time() : G.time();
  ndeme() = (ndeme() > G.ndeme()) ? ndeme() : G.ndeme();
  _unique = G._unique;
  curtail(t,t0);
  return *this;
}
 
extern "C" {

  SEXP genealSum (SEXP args) {
    args = CDR(args);
    genealogy_t A(R_NegInf); // a "null" genealogy on [-inf,inf]
    A.time() = R_PosInf;
    while (args != R_NilValue) {
      A += CAR(args);
      args = CDR(args);
    }
    SEXP S;
    PROTECT(S = serial(A));
    SET_ATTR(S,install("class"),mkString("gpgen"));
    UNPROTECT(1);
    return S;
  }
 
}