anchor.diff:  
       subroutine terminate_iteration
           call print_optimization_status
  -        call save_point_file(points, iteration_count)
  +        call save_point_file(points, energy, force, iteration_count)
           write(*,*) "Convergence has been achieved (up to numerical&
                   & round-off error). Exiting."
           stop
  
positive.code: 
  
      subroutine print_optimization_status
          write(*,'(I10,A)',advance="no") iteration_count, " |"
          write(*,'('//rf//',A)',advance="no") energy, " |"
          write(*,'('//rf//',A)',advance="no") &
                  & norm2(force) / sqrt(real(num_points, rk)), " |"
  #ifdef PCREO_TRACK_ANGLE
          write(*,'('//rf//',A)',advance="no") step_size, " |"
          write(*,'('//rf//')') step_angle
  #else
          write(*,'('//rf//')') step_size
  
negative.code:             last_print_time = cur_time
          end if
          if (cur_time - last_save_time >= save_time) then
              call save_point_file(points, energy, force, iteration_count)
              last_save_time = cur_time
          end if
      end subroutine finish_iteration
  