mosesdecoder/moses/TranslationModel/WordCoocTable.cpp

#include "moses/TranslationModel/WordCoocTable.h"
using namespace std;
namespace Moses {

  WordCoocTable::
  WordCoocTable()
  {
    m_cooc.reserve(1000000);
    m_marg1.reserve(1000000);
    m_marg2.reserve(1000000);
  }
  
  WordCoocTable::
  WordCoocTable(wordID_t const VocabSize1, wordID_t const VocabSize2)
    : m_cooc(VocabSize1), m_marg1(VocabSize1,0), m_marg2(VocabSize2, 0)
  {}
  
  void
  WordCoocTable::
  Count(size_t const a, size_t const b)
  {
    while (a >= m_marg1.size()) 
      {
	m_cooc.push_back(my_map_t());
	m_marg1.push_back(0);
      }
    while (b >= m_marg2.size()) 
      m_marg2.push_back(0);
    ++m_marg1[a];
    ++m_marg2[b];
    ++m_cooc[a][b];
  }
  
  uint32_t 
  WordCoocTable::
  GetJoint(size_t const a, size_t const b) const
  {
    if (a >= m_marg1.size() || b >= m_marg2.size()) return 0;
    my_map_t::const_iterator m = m_cooc.at(a).find(b);
    if (m == m_cooc[a].end()) return 0;
    return m->second;
  }

  uint32_t 
  WordCoocTable::
  GetMarg1(size_t const x) const
  {
    return x >= m_marg1.size() ? 0 : m_marg1[x];
  }

  uint32_t 
  WordCoocTable::
  GetMarg2(size_t const x) const
  {
    return x >= m_marg2.size() ? 0 : m_marg2[x];
  }

  float 
  WordCoocTable::
  pfwd(size_t const a, size_t const b) const
  {
    return float(GetJoint(a,b))/GetMarg1(a);
  }

  float 
  WordCoocTable::
  pbwd(size_t const a, size_t const b) const
  {
    // cerr << "at " << __FILE__ << ":" << __LINE__ << endl;
    return float(GetJoint(a,b))/GetMarg2(b);
  }
}