DOC:

**Issue [Documentation] `pylops.lsqr` documentation lists 9 results even though there is 11 #637**
- added `iiter` and `xnorm` to the documentation of the `pylops.lsqr` solver and related methods

Author: IruNikZe

pylops/optimization/basic.py

@@ -242,6 +242,8 @@ def lsqr(
 
         ``7`` means the iteration limit has been reached
 
+    iiter : :obj:`int`
+        Iteration number upon termination
     r1norm : :obj:`float`
         :math:`||\mathbf{r}||_2^2`, where
         :math:`\mathbf{r} = \mathbf{y} - \mathbf{Op}\,\mathbf{x}`
@@ -257,6 +259,8 @@ def lsqr(
     arnorm : :obj:`float`
         Estimate of norm of :math:`\cond(\mathbf{Op}^H\mathbf{r}-
         \epsilon^2\mathbf{x})`
+    xnorm : :obj:`float`
+        :math:`||\mathbf{x}||_2`
     var : :obj:`float`
         Diagonals of :math:`(\mathbf{Op}^H\mathbf{Op})^{-1}` (if ``damp=0``)
         or more generally :math:`(\mathbf{Op}^H\mathbf{Op} +

pylops/optimization/cls_basic.py

@@ -1043,6 +1043,8 @@ def solve(
 
             ``7`` means the iteration limit has been reached
 
+        iiter : :obj:`int`
+            Iteration number upon termination
         r1norm : :obj:`float`
             :math:`||\mathbf{r}||_2^2`, where
             :math:`\mathbf{r} = \mathbf{y} - \mathbf{Op}\,\mathbf{x}`
@@ -1058,6 +1060,8 @@ def solve(
         arnorm : :obj:`float`
             Estimate of norm of :math:`\cond(\mathbf{Op}^H\mathbf{r}-
             \epsilon^2\mathbf{x})`
+        xnorm : :obj:`float`
+            :math:`||\mathbf{x}||_2`
         var : :obj:`float`
             Diagonals of :math:`(\mathbf{Op}^H\mathbf{Op})^{-1}` (if ``damp=0``)
             or more generally :math:`(\mathbf{Op}^H\mathbf{Op} +