Changes in version 2025-01-17                      

  - NULL is replaced by a vector of 0 when estimating b norm
  - fixed returned norm of residual for convergence on x0

                     Changes in version 2024-10-17                      

  - initial R implemetation