diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index d97fcce462a74bc9d56e27feecb16b5321ad5890..04e780cbe357afb938884c4bae20ccd67c19a971 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -538,11 +538,12 @@ int main(int argc, char *argv[]) { SingletonVectorType computed_state; stateUpdater->extractState(computed_state); double const correction = stateEnergyNorm.diff(computed_state, alpha); - if (state_fpi <= 2 // Let the first two steps pass through unchanged - || correction < requiredReduction * lastCorrection) { + + if (state_fpi <= 2 or correction < + requiredReduction * lastCorrection) { alpha = computed_state; relaxationWriter << "N "; - } else { // alpha is still the old time step here + } else { alpha *= relaxation; Arithmetic::addProduct(alpha, 1.0 - relaxation, computed_state); relaxationWriter << "Y ";