void sample_hyperpriors(double res){ sample_alpha(res); sample_U(); sample_V(); //if (tensor) // sample_T(); }
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; }