##### #test.show(quiver=False) ##### start = clock() #test.precondition(101,tol=1e-2,verbose=False) #test.precondition(81,tol=1e-1,verbose=False) test.precondition((60, 70), tol=1e-5, verbose=False) #test.precondition(251,tol=1e-3,verbose=False) precondtime = clock() - start times.append(precondtime) ##### #test.show(title='after preconditioning') ##### test.SOR(tol=tol, max_iter=max_iter, verbose=False) ##### time2 = clock() - start times.append(time2) #test.show(title='result with precon') ##### plt.figure() plt.imshow(np.abs(raw.potentials - test.potentials).T, origin='lower') plt.colorbar() plt.title('raw - test') plt.tight_layout() plt.show()