gibbs.cpp

来自「The package includes 3 Matlab-interfaces」· C++ 代码 · 共 33 行

#include "Gibbs.h"
#include <stdlib.h>
#include <iostream>

void Gibbs::transition() {

  int i = chooseInteger(mc_Pi, ia_mrf->N);
  int qi = ia_mrf->V[i];
  for (int xi=0; xi<qi; xi++) {
    mc_Pxi[i][xi] = ia_mrf->localMat[i][xi];
  }
  
  for (int n=0; n<ia_mrf->neighbNum(i); n++) {
    int j = ia_mrf->adjMat[i][n];
    int xj = mc_currX[j];
    for (int xi=0; xi<qi; xi++) {
      mc_Pxi[i][xi] *= ia_mrf->pairPotential(i,n,xi,xj);
    }
  }

  // normalize Pxi
  double sum_Pxi = 0.0;
  for (int xi=0; xi<qi; xi++) {
    sum_Pxi += mc_Pxi[i][xi];
  }
  for (int xi=0; xi<qi; xi++) {
    mc_Pxi[i][xi] /= sum_Pxi;
  }
  
  mc_currX[i] = chooseInteger(mc_Pxi[i], qi);
  
}