Hi Florian,

    Thank you for the note. I agree this need to be fixed.

Best regards,
Abdou


On Mon, Apr 7, 2014 at 5:43 PM, Florian <notifications@github.com> wrote:

> In solver/eigcg.c:
>
> > +        _vector_null(x[i].s2);
> > +        _vector_null(x[i].s3);
> > +      }
> > +
> > +      *flag = 0;		
> > +      *reshist = 0.0;
> > +      if( g_proc_id == g_stdio_proc)
> > +        displayInfo(eps_sq,maxit,*flag,*iter,*reshist);
> > +      return;
> > +     }
> > +
> > +  }
> > +
> > +  /* Set up for the method */
> > +  *flag = 1;
> > +  tolb = eps_sq * (*normb)*(*normb);	/* Relative to b tolerance */
>
> Since tolb is multiplied with normb^2 the stopping criterion is
> automatically thought of as relative. More important however: when called
> from incr_eigcg normb depends on the history of the algorithm.
>
> --
> Reply to this email directly or view it on GitHub<https://github.com/etmc/tmLQCD/pull/246/files#r11346274>
> .
>