public Factor duplicate () { return new LogTableFactor (this); }
public AbstractTableFactor recenter () { // return (AbstractTableFactor) normalize (); int loc = argmax (); double lval = probs.valueAtLocation(loc); // should be refactored timesEqualsLog (-lval); return this; } }
private void comparePotentialValues (LogTableFactor ptl, LogTableFactor ptl2) { AssignmentIterator it1 = ptl.assignmentIterator (); AssignmentIterator it2 = ptl2.assignmentIterator (); while (it1.hasNext ()) { assertTrue (ptl.value (it1) == ptl.value (it2)); it1.advance (); it2.advance (); } }
protected double rawValue (Assignment assn) { int numVars = getNumVars (); int[] indices = new int[numVars]; for (int i = 0; i < numVars; i++) { Variable var = getVariable (i); indices[i] = assn.get (var); } return rawValue (indices); }
public void testPlusEquals () { Variable var = new Variable (4); // log 0, log 1, log 2, log 3 double[] vals = new double[] { Double.NEGATIVE_INFINITY, 0, 0.6931471805599453, 1.0986122886681098 }; LogTableFactor factor = LogTableFactor.makeFromLogValues (var, vals); factor.plusEquals (0.1); // log 0.1, log 1.1, log 2.1, log 3.1 double[] expected = new double[] { -2.3025850929940455, 0.09531017980432493, 0.7419373447293773, 1.1314021114911006 }; LogTableFactor ans = LogTableFactor.makeFromLogValues (var, expected); assertTrue ("Error: expected "+ans.dumpToString ()+" but was "+factor.dumpToString (), factor.almostEquals (ans)); }
public void testRecenter2 () { Variable var = new Variable (4); double[] vals = new double[]{ 0, 1.4, 1.4, 0 }; LogTableFactor ltbl1 = LogTableFactor.makeFromLogValues (var, vals); ltbl1.recenter (); double[] expected = new double[]{ -1.4, 0, 0, -1.4 }; LogTableFactor ans = LogTableFactor.makeFromLogValues (var, expected); assertTrue (!ltbl1.isNaN ()); assertTrue ("Error: expected "+ans.dumpToString ()+"but was "+ltbl1.dumpToString (), ans.almostEquals (ltbl1)); }
private DiscreteFactor constructConditionalCpt (FactorGraph mdl, Variable var, Assignment fullAssn) { List ptlList = mdl.allFactorsContaining (var); LogTableFactor ptl = new LogTableFactor (var); for (AssignmentIterator it = ptl.assignmentIterator (); it.hasNext(); it.advance ()) { Assignment varAssn = it.assignment (); fullAssn.setValue (var, varAssn.get (var)); ptl.setRawValue (varAssn, sumValues (ptlList, fullAssn)); } ptl.normalize (); return ptl; }
public void testLogValue () { Variable[] vars = new Variable[] { new Variable (2), new Variable (2) }; LogTableFactor ptl = LogTableFactor.makeFromValues (vars, new double[] { 1, 2, 3, 4 }); Assignment assn = new Assignment (vars, new int [vars.length]); assertEquals (0, ptl.logValue (assn), 1e-5); assertEquals (0, ptl.logValue (ptl.assignmentIterator()), 1e-5); assertEquals (0, ptl.logValue (0), 1e-5); assertEquals (1, ptl.value (assn), 1e-5); assertEquals (1, ptl.value (ptl.assignmentIterator()), 1e-5); assertEquals (1, ptl.value (0), 1e-5); LogTableFactor ptl2 = LogTableFactor.makeFromLogValues (vars, new double[] { 0, Math.log (2), Math.log (3), Math.log (4) }); Assignment assn2 = new Assignment (vars, new int [vars.length]); assertEquals (0, ptl2.logValue (assn2), 1e-5); assertEquals (0, ptl2.logValue (ptl2.assignmentIterator()), 1e-5); assertEquals (0, ptl2.logValue (0), 1e-5); assertEquals (1, ptl2.value (assn2), 1e-5); assertEquals (1, ptl2.value (ptl2.assignmentIterator()), 1e-5); assertEquals (1, ptl2.value (0), 1e-5); }
public void testRecenter () { Variable var = new Variable (4); double[] vals = new double[]{ 2.0, 4.0, 6.0, 8.0 }; LogTableFactor ltbl1 = LogTableFactor.makeFromValues (var, vals); ltbl1.recenter (); double[] expected = new double[] { Math.log (0.25), Math.log(0.5), Math.log (0.75), 0 }; LogTableFactor ans = LogTableFactor.makeFromLogValues (var, expected); assertTrue ("Error: expected "+ans.dumpToString ()+"but was "+ltbl1.dumpToString (), ans.almostEquals (ltbl1)); }
protected Factor slice_onevar (Variable var, Assignment observed) { Assignment assn = (Assignment) observed.duplicate (); double[] vals = new double [var.getNumOutcomes ()]; for (int i = 0; i < var.getNumOutcomes (); i++) { assn.setValue (var, i); vals[i] = logValue (assn); } return LogTableFactor.makeFromLogValues (var, vals); }
public void testSparseValueAndLogValue () { Variable[] vars = new Variable[] { new Variable (2), new Variable (2) }; int[] szs = { 2, 2 }; int[] idxs1 = new int[] { 1, 3 }; double[] vals1 = new double[]{ 4.0, 8.0 }; LogTableFactor ptl1 = LogTableFactor.makeFromMatrix (vars, new SparseMatrixn (szs, idxs1, vals1)); AssignmentIterator it = ptl1.assignmentIterator (); assertEquals (1, it.indexOfCurrentAssn ()); assertEquals (Math.log (4), ptl1.logValue (it), 1e-5); assertEquals (Math.log (4), ptl1.logValue (it.assignment ()), 1e-5); assertEquals (4, ptl1.value (it), 1e-5); assertEquals (4, ptl1.value (it.assignment ()), 1e-5); it = ptl1.varSet ().assignmentIterator (); assertEquals (0, it.indexOfCurrentAssn ()); assertEquals (Double.NEGATIVE_INFINITY, ptl1.logValue (it), 1e-5); assertEquals (Double.NEGATIVE_INFINITY, ptl1.logValue (it.assignment ()), 1e-5); assertEquals (0, ptl1.value (it), 1e-5); assertEquals (0, ptl1.value (it.assignment ()), 1e-5); }
public void testLogMultiVarSlice () { double[] vals = { 0.0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 }; Variable v1 = new Variable (2); Variable v2 = new Variable (2); Variable v3 = new Variable (2); Variable v4 = new Variable (2); Variable[] vars = new Variable[]{v1, v2, v3, v4}; Factor ptl = LogTableFactor.makeFromValues (vars, vals); System.out.println (ptl.dumpToString ()); Assignment assn = new Assignment (v4, 0); LogTableFactor sliced = (LogTableFactor) ptl.slice (assn); LogTableFactor expected = LogTableFactor.makeFromValues (new Variable[] { v1,v2,v3 }, new double[] { 0, 2, 4, 6, 8, 10, 12, 14 }); assertTrue ("Test failed. Expected: "+expected.dumpToString ()+"\nActual: "+sliced.dumpToString (), expected.almostEquals (sliced)); }
public double value (Assignment assn) { Flops.exp (); if (getNumVars () == 0) return 1.0; return Math.exp (rawValue (assn)); }
protected void plusEqualsAtLocation (int loc, double v) { Flops.log (); Flops.sumLogProb (1); double oldVal = logValue (loc); setRawValue (loc, Maths.sumLogProb (oldVal, Math.log (v))); }
public Factor normalize () { double sum = logspaceOneNorm (); if (sum < -500) System.err.println ("Attempt to normalize all-0 factor "+this.dumpToString ()); for (int i = 0; i < probs.numLocations (); i++) { double val = probs.valueAtLocation (i); probs.setValueAtLocation (i, val - sum); } return this; }
public static LogTableFactor multiplyAll (Collection phis) { /* Get all the variables */ VarSet vs = new HashVarSet (); for (Iterator it = phis.iterator (); it.hasNext ();) { Factor phi = (Factor) it.next (); vs.addAll (phi.varSet ()); } /* define a new potential over the neighbors of NODE */ LogTableFactor newCPF = new LogTableFactor (vs); for (Iterator it = phis.iterator (); it.hasNext ();) { Factor phi = (Factor) it.next (); newCPF.multiplyBy (phi); } return newCPF; }
protected void divideByInternal (DiscreteFactor ptl) { int[] projection = largeIdxToSmall (ptl); int numLocs = probs.numLocations (); for (int singleLoc = 0; singleLoc < numLocs; singleLoc++) { int smallIdx = projection[singleLoc]; double prev = this.probs.valueAtLocation (singleLoc); double newVal = ptl.logValue (smallIdx); double product = prev - newVal; /* by convention, let -Inf + Inf (corresponds to 0/0) be -Inf */ if (Double.isInfinite (newVal)) { product = Double.NEGATIVE_INFINITY; } this.probs.setValueAtLocation (singleLoc, product); } Flops.increment (numLocs); // handle the pluses }
public double logsum () { return logspaceOneNorm (); }
private Factor sliceForAlpha (Assignment assn) { double alph = assn.getDouble (alpha); int[] sizes = sizesFromVarSet (xs); Matrix diag = Matrices.diag (sizes, alph); Matrix matrix = Matrices.constant (sizes, -alph); matrix.plusEquals (diag); return LogTableFactor.makeFromLogMatrix (xs.toVariableArray (), (SparseMatrixn) matrix); }
private DiscreteFactor constructConditionalCpt (FactorGraph mdl, Variable var, Assignment fullAssn) { List ptlList = mdl.allFactorsContaining (var); LogTableFactor ptl = new LogTableFactor (var); for (AssignmentIterator it = ptl.assignmentIterator (); it.hasNext(); it.advance ()) { Assignment varAssn = it.assignment (); fullAssn.setValue (var, varAssn.get (var)); ptl.setRawValue (varAssn, sumValues (ptlList, fullAssn)); } ptl.normalize (); return ptl; }