Commit efd86dfb authored by s.gretchko's avatar s.gretchko
Browse files

Fixed bug in DyRePr algorithm

parent 9abb1308
......@@ -67,7 +67,7 @@ class DyRePr(AvP):
if self.accelerator is not None:
u_intermed = self.accelerator.evaluate(u_intermed, self)
for j in range(self.n_product_Prox):
weight = norm(u_intermed[j], ord=2)**2
weight = norm(u - u_intermed[j], ord=2)**2
u_intermed[j] /= weight + 1
u_new = self.p_avg.eval(u_intermed)
else:
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment