123 |
return True |
return True |
124 |
else: |
else: |
125 |
return False |
return False |
126 |
def stoppingcriterium_GMRES(self,rho,r): |
def stoppingcriterium_GMRES(self,norm_r,norm_b): |
127 |
if self.verbose: print "GMRES step %s: L2(rho) = %s, L2(b)*TOL=%s"%(self.iter,rho,r*self.getTolerance()) |
if self.verbose: print "GMRES step %s: L2(r) = %s, L2(b)*TOL=%s"%(self.iter,norm_r,norm_b*self.getTolerance()) |
128 |
self.iter+=1 |
self.iter+=1 |
129 |
if rho <= r*self.getTolerance(): |
if norm_r <= norm_b*self.getTolerance(): |
130 |
if self.verbose: print "GMRES terminated after %s steps."%self.iter |
if self.verbose: print "GMRES terminated after %s steps."%self.iter |
131 |
return True |
return True |
132 |
else: |
else: |