anchor.diff:          if (ex) then
               write(*,*) "Loading initial point configuration from file..."
               open(newunit=u, file="pcreo_input.csv")
  +            read(u,*)
  +            read(u,*)
               read(u,*) points
               close(u)
           else
  
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 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
  