Invert the forward matrix using conjugate gradient scheme on the normal equation. Regularization is provided by truncated the operation after a fixed number of iterations.
| Syntax: | [X, resid] = tcgls(A, Y, X0, nIter[, W]); | |
|---|---|---|
| Inputs: | A | The forward matrix |
| Y | The residue appropriate to |
|
| X0 | Initial guess or zeros if passed as the empty list |
|
| nIter | How many iterations to compute | |
| W | Relaxation parameter. Optional: 1.0 if not specified | |
| Outputs: | X | The reconstructed image |
| resid | The final residual vector | |