def _calculateUAa(self, b, W): A = 1.0 + W*self._sign2 V = W/A Lk = self._calculateLk(self._G01, V) Gtb = dot(self._G01.T, b) GtV = ddot(self._G01.T, V, left=False) LtLGtV = stu(Lk.T, stl(Lk, GtV)) LtLGtVG = dot(LtLGtV, self._G01) bn = self._sign2*b a = b + dot(dot(GtV.T, LtLGtVG) - GtV.T, Gtb)\ + dot(GtV.T, dot(LtLGtV, bn)) - V*bn return a
def _calculateUAa(self, b, W): A = 1.0 + W * self._sign2 V = W / A Lk = self._calculateLk(self._G01, V) Gtb = dot(self._G01.T, b) GtV = ddot(self._G01.T, V, left=False) LtLGtV = stu(Lk.T, stl(Lk, GtV)) LtLGtVG = dot(LtLGtV, self._G01) bn = self._sign2 * b a = b + dot(dot(GtV.T, LtLGtVG) - GtV.T, Gtb)\ + dot(GtV.T, dot(LtLGtV, bn)) - V*bn return a
def _calculateUAa(self, b, W): Wsq = NP.sqrt(W) Ln = self._calculateLn(self._K, Wsq) a = b - Wsq * stu(Ln.T, stl(Ln, Wsq*dot(self._K,b))) return a
def _calculateUAa(self, b, W): Wsq = NP.sqrt(W) Ln = self._calculateLn(self._K, Wsq) a = b - Wsq * stu(Ln.T, stl(Ln, Wsq * dot(self._K, b))) return a