public Lattice(model m, inference inf, dataSeq x) { _w = x.Count; _h = m.NTag; _logBel = new belief(_w, _h); List <dMatrix> YYlist = new List <dMatrix>(); List <List <double> > Ylist = new List <List <double> >(); inf.getYYandY(m, x, YYlist, Ylist); for (int i = 0; i < _w; i++) { _logBel.belState[i] = new List <double>(Ylist[i]); if (i > 0) { _logBel.belEdge[i] = new dMatrix(YYlist[i]); } } _heuListList = new List <List <double> >(); for (int i = 0; i < _w; i++) { _heuListList.Add(new List <double>(new double[_h])); } Viterbi _bwdViterbi = new Viterbi(_w, _h); for (int i = 0; i < _w; i++) { _bwdViterbi.setScores(i, Ylist[i], YYlist[i]); } List <int> tags = new List <int>(); _bwdViterbi.runViterbi(ref tags); //update the viterbiHeuristicMap for (int i = 0; i < _w; i++) { for (int j = 0; j < _h; j++) { double h = _bwdViterbi.getPathScore(i, j); setHeuMap(i, j, h); } } //get zGold ZGold = 0; for (int i = 0; i < x.Count; i++) { int s = x.getTags(i); ZGold += Ylist[i][s]; if (i > 0) { int sPre = x.getTags(i - 1); ZGold += YYlist[i][sPre, s]; } } }
virtual public double getGradCRF(List <double> vecGrad, model m, dataSeq x, baseHashSet <int> idSet) { if (idSet != null) { idSet.Clear(); } int nTag = m.NTag; //compute beliefs belief bel = new belief(x.Count, nTag); belief belMasked = new belief(x.Count, nTag); //store the YY and Y List <dMatrix> YYlist = new List <dMatrix>(), maskYYlist = new List <dMatrix>(); List <List <double> > Ylist = new List <List <double> >(), maskYlist = new List <List <double> >(); _inf.getYYandY(m, x, YYlist, Ylist, maskYYlist, maskYlist); _inf.getBeliefs(bel, m, x, YYlist, Ylist); _inf.getBeliefs(belMasked, m, x, maskYYlist, maskYlist); double ZGold = belMasked.Z; double Z = bel.Z; List <featureTemp> fList; //Loop over nodes to compute features and update the gradient for (int i = 0; i < x.Count; i++) { fList = _fGene.getFeatureTemp(x, i); foreach (featureTemp im in fList) { for (int s = 0; s < nTag; s++) { int f = _fGene.getNodeFeatID(im.id, s); if (idSet != null) { idSet.Add(f); } vecGrad[f] += bel.belState[i][s] * im.val; vecGrad[f] -= belMasked.belState[i][s] * im.val; } } } //Loop over edges to compute features and update the gradient for (int i = 1; i < x.Count; i++) { for (int s = 0; s < nTag; s++) { for (int sPre = 0; sPre < nTag; sPre++) { int f = _fGene.getEdgeFeatID(sPre, s); if (idSet != null) { idSet.Add(f); } vecGrad[f] += bel.belEdge[i][sPre, s]; vecGrad[f] -= belMasked.belEdge[i][sPre, s]; } } } return(Z - ZGold);//-log{P(y*|x,w)} }