예제 #1
0
void sample_hyperpriors(double res){
    sample_alpha(res);
    sample_U();
    sample_V();
    //if (tensor) 
    //  sample_T();
}
예제 #2
0
HHRESULT CGaussianMDP::sample()
{
   int j = 0;
   HHRESULT hr = HH_OK;

   numconfig = 0;

   for(j=0; j<cY; j++)
   {
      hr = sample_config(config, j, sigma2, cY, Y, phi, alpha);
      numconfig = imax2(numconfig, config[j]);
   }
   
   hr = sample_phi(config, Y, sigma2, s2, m, cY, phi, numconfig);
   hr = sample_m(s2, phi, numconfig, m);
   hr = sample_s2(s2, phi, m, numconfig, w, W);
   hr = sample_theta(config, phi, cY, theta);
   hr = sample_alpha(par1, par2, cY,numconfig, alpha);
      
Cleanup:
   return hr;
Error:
   goto Cleanup;
}