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
