AdaptiveMap.cpp

#include "AdaptiveMap.hpp"

AdaptiveMap::AdaptiveMap(const Shared* const sh, const int n) : shared(sh), t(n) {
  dt = DivisionTable::getDT();
}

void AdaptiveMap::update(uint32_t *const p, const int limit) {
  assert(limit > 0 && limit < 1024); 
  uint32_t p0 = p[0];
  const int n = p0 & 1023; //count
  const int pr = p0 >> 10; //prediction (22-bit fractional part)
  if( n < limit ) {
    ++p0;
  } else {
    p0 = (p0 & 0xfffffc00) | limit;
  }
  INJECT_SHARED_y
  const int target = y << 22; //(22-bit fractional part)
  const int delta = ((target - pr) >> 3) * dt[n]; //the larger the count (n) the less it should adapt pr+=(target-pr)/(n+1.5)
  p0 += delta & 0xfffffc00;
  p[0] = p0;
}