wehub-resource-sync
9 lines
348 B
Matlab
Executable File
9 lines
348 B
Matlab
Executable File
function lnZ = gibbsEnergy(nodePot, edgePot, nodeBel, edgeBel)
|
|
% Compute Gibbs free energy
|
|
% TBD: deal with log(0) for entropy
|
|
edgePot = reshape(edgePot,[],size(edgePot,3));
|
|
edgeBel = reshape(edgeBel,[],size(edgeBel,3));
|
|
Ex = dot(nodeBel,nodePot,1);
|
|
Exy = dot(edgeBel,edgePot,1);
|
|
Hx = dot(nodeBel,log(nodeBel),1);
|
|
lnZ = -(sum(Ex)+sum(Exy)+sum(Hx)); |