@Override
public void incrementSufficientStatisticsFromMarginal(
SufficientStatistics gradient,
SufficientStatistics currentParameters,
Factor marginal,
Assignment conditionalAssignment,
double count,
double partitionFunction) {
if (conditionalAssignment.containsAll(getVars().getVariableNumsArray())) {
// Short-circuit the slow computation below if possible.
double multiplier = marginal.getTotalUnnormalizedProbability() * count / partitionFunction;
incrementSufficientStatisticsFromAssignment(
gradient, currentParameters, conditionalAssignment, multiplier);
} else {
VariableNumMap conditionedVars =
initialWeights.getVars().intersection(conditionalAssignment.getVariableNumsArray());
TableFactor productFactor =
(TableFactor)
initialWeights
.product(
TableFactor.pointDistribution(
conditionedVars, conditionalAssignment.intersection(conditionedVars)))
.product(marginal);
Tensor productFactorWeights = productFactor.getWeights();
double[] productFactorValues = productFactorWeights.getValues();
int tensorSize = productFactorWeights.size();
double multiplier = count / partitionFunction;
TensorSufficientStatistics tensorGradient = (TensorSufficientStatistics) gradient;
for (int i = 0; i < tensorSize; i++) {
int builderIndex = (int) productFactorWeights.indexToKeyNum(i);
tensorGradient.incrementFeatureByIndex(productFactorValues[i] * multiplier, builderIndex);
}
}
}
