--- trunk/escript/py_src/pdetools.py 2008/04/15 08:47:57 1513 +++ trunk/escript/py_src/pdetools.py 2008/04/16 00:15:44 1514 @@ -579,7 +579,7 @@ vrot[i:i+2]=w1,w2 return vrot -def GMRES(b, Aprod, Msolve, bilinearform, stoppingcriterium, x=None, iter_max=100, iter_restart=10): +def GMRES(b, Aprod, Msolve, bilinearform, stoppingcriterium, x=None, iter_max=100, iter_restart=20): m=iter_restart iter=0 while True: @@ -589,7 +589,7 @@ if stopped: break return x -def GMRESm(b, Aprod, Msolve, bilinearform, stoppingcriterium, x=None, iter_max=100, iter_restart=10): +def GMRESm(b, Aprod, Msolve, bilinearform, stoppingcriterium, x=None, iter_max=100, iter_restart=20): iter=0 r=Msolve(b) r_dot_r = bilinearform(r, r)