Commit 2395b49a authored by Marco Selig's avatar Marco Selig

bug in conjugate_gradient._without fixed.

parent 9421e010
......@@ -717,7 +717,8 @@ class conjugate_gradient(object):
if(limii is None):
limii = 10*self.b.domain.dim(split=False)
d = r = self.b-self.A(self.x)
r = self.b-self.A(self.x)
d = field(self.b.domain,val=np.copy(r.val),target=self.b.target)
gamma = r.dot(d)
delta_ = np.absolute(gamma)**(-0.5)
......
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