anchor.diff:  
   
       subroutine print_optimization_status
  -        write(*,'(I10,A)',advance="no") iteration_count, " |"
  -        write(*,'('//rf//',A)',advance="no") energy, " |"
  -        write(*,'('//rf//',A)',advance="no") &
  +        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//',A)',advance='no') step_size, " |"
           write(*,'('//rf//')') step_angle
   #else
           write(*,'('//rf//')') step_size
  
positive.code:         if (ex) then
              write(*,*) "Loading initial point configuration from file..."
              open(newunit=u, file="pcreo_input.csv")
              read(u,*) points
              close(u)
          else
  
negative.code:     end subroutine print_table_header
  
  
      subroutine save_point_file(points, energy, force, idx)
          real(rk), intent(in) :: points(d + 1, num_points)
          real(rk), intent(in) :: energy, force(d + 1, num_points)
          integer, intent(in) :: idx
  
          integer :: i, j, u
          character(len=20) :: fname
          character(len=80) :: line
  
          write(fname,"(a,i10.10,a)") "pcreo_", idx, ".csv"
          open(file=fname, newunit=u)
  
          write(line,*) d
          write(u,'(A,A)',advance='no') trim(adjustl(line)), ", "
          write(line,*) s
          write(u,'(A,A)',advance='no') trim(adjustl(line)), ", "
  #ifdef PCREO_SYMMETRY
          write(line,*) num_external_points + symmetry_group_order * num_points
  #else
          write(line,*) num_points
  #endif
          write(u,'(A)') trim(adjustl(line))
          write(u,'(SP,'//rf//',A)',advance='no') energy, ", "
          write(u,'(SP,'//rf//')') norm2(force) / sqrt(real(num_points, rk))
  
          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
  