Following the ``trial solution'' approach described in (), I change the variables 397#397 and solve equation () for the preconditioned variable 398#398:
399#399 | (158) |
Operator 400#400 in equation () is a cascade of two filters described in the previous section. After finding a solution for 398#398,I evaluate 397#397 to obtain the solution of the original problem. To avoid high frequency noise in the model, I introduce a regularization term into the problem and solve the system of equations:
401#401 |
In this paper I used the laplacian as the regularization operator 8#8 in equation ().
Figure shows the solution to the least-squares problem [equation ()] after 25 iterations of the conjugate-gradient method, the velocity stack after the filtering, and the solution to the preconditioned least-squares problem [equation ()] after 25 iterations of the conjugate-gradient.
Figure shows the residual for the preconditioned problem [equation ()] and the problem without preconditioning [equation ()].
res_ann
Figure 6 CG convergence with and without preconditioning |
As Figures and show, although the solution for the preconditioned problem does not have artifacts, it converges much slower than the solution to the problem without preconditioning, which probably makes this method of preconditioning impractical.