I have a little question concerning the implementation of gsl_siman_solve. The variable n_evals is used to count the number of evaluations of Ef, but n_evals is incremented only after the computation of new_E, not after after the computation of E in the first line of the while loop. Can you explain why you don't take it into account? Thank you, Fabian Bastin