Commit f1faedb1 authored by jansen31's avatar jansen31
Browse files

sign change in iterate? Ad hoc change of which I am not sure

parent 18538e9f
......@@ -46,6 +46,8 @@ class P_orthonorm(ProxOperator):
u_new = zeros_like(u_norm)
u_new[1] = u_norm[0] - (y / (y ** 2 - 1)) * (y * u_norm[0] - u_norm[1])
u_new[0] = (1 / (y ** 2 - 1)) * (y * u_norm[0] - u_norm[1])
if a < 0: # Not sure if or why this is needed...
u_new *= -1
return u_new
......
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