anchor.diff:              last_print_time = cur_time
           end if
           if (cur_time - last_save_time >= save_time) then
  -            call save_point_file(points, iteration_count)
  +            call save_point_file(points, energy, force, iteration_count)
               last_save_time = cur_time
           end if
       end subroutine finish_iteration
  
positive.code:     end subroutine print_table_header
  
  
      subroutine save_point_file(points, idx)
          real(rk), intent(in) :: points(d + 1, num_points)
          integer, intent(in) :: idx
  
          integer :: i, j, u
          character(len=20) :: fname
  
          write(fname,"(a,i10.10,a)") "pcreo_", idx, ".csv"
          open(file=fname, newunit=u)
          do j = 1, num_points
              do i = 1, d + 1
                  if (i > 1) write(u,'(A)',advance="no") ", "
                  write(u,'(SP,'//rf//')',advance="no") points(i,j)
              end do
              write(u,*)
          end do
  
negative.code: 
      subroutine terminate_iteration
          call print_optimization_status
          call save_point_file(points, energy, force, iteration_count)
          write(*,*) "Convergence has been achieved (up to numerical&
                  & round-off error). Exiting."
          stop
  