Hej, I have been using a very simple Scheduled Relaxation Jacobi (SRJ) method for tests purpose in one of my code, and would now like to implement into the big version that uses PETSc. So far, I have figured that I can use the weighted Jacobi method since it is nothing but a Jacobi-preconditioned Richardson. While in the weighted Jacobi, the relaxation weight omega is fixed, in the SRJ method, the value of the relaxation factor depends on the grid size and the iteration. Is there any simple way to implement such a iteration-varying relaxation weight given that I have text files with their appropriate values?
Thanks a lot. JC
