ures the relative error desired in the sum of
          squares, for the initialization step of nonlinear part.
          Termination occurs when the actual relative reduction in
          the sum of squares is at most TOL1. In addition, if
          ALG = 'I', TOL1 also measures the relative residual of
          the solutions computed by the CG algorithm (for the
          initialization step). Termination of a CG process occurs
          when the relative residual is at most TOL1.
          If the user s