Commit aeef62a7 by Martin Reinecke

### fix

parent 6a35d71f
Pipeline #14523 passed with stage
in 6 minutes and 15 seconds
 ... @@ -204,7 +204,7 @@ class InformationStore(object): ... @@ -204,7 +204,7 @@ class InformationStore(object): kmj = (k-m+j)%m kmj = (k-m+j)%m result[i, j] = self._ss_store[kmi, kmj] result[i, j] = self._ss_store[kmi, kmj] result[i, m+j] = result[m+j, i] = self._sy_store[kmi, kmj] result[i, m+j] = result[m+j, i] = self._sy_store[kmi, kmj] result[m+i, m+j] = self._yy_store[(k-m+i)%m, (k-m+j)%m] result[m+i, m+j] = self._yy_store[kmi, kmj] sgrad_i = self.s[k-m+i].vdot(self.last_gradient) sgrad_i = self.s[k-m+i].vdot(self.last_gradient) result[2*m, i] = result[i, 2*m] = sgrad_i result[2*m, i] = result[i, 2*m] = sgrad_i ... @@ -242,7 +242,7 @@ class InformationStore(object): ... @@ -242,7 +242,7 @@ class InformationStore(object): for i in xrange(2*m+1): for i in xrange(2*m+1): delta[i] *= b_dot_b[m-1, 2*m-1]/b_dot_b[2*m-1, 2*m-1] delta[i] *= b_dot_b[m-1, 2*m-1]/b_dot_b[2*m-1, 2*m-1] for j in xrange(m-1, -1, -1): for j in xrange(m): delta_b_b = sum([delta[l]*b_dot_b[m+j, l] for l in xrange(2*m+1)]) delta_b_b = sum([delta[l]*b_dot_b[m+j, l] for l in xrange(2*m+1)]) beta = delta_b_b/b_dot_b[j, m+j] beta = delta_b_b/b_dot_b[j, m+j] delta[j] += (alpha[j] - beta) delta[j] += (alpha[j] - beta) ... ...
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!