anchor.diff:               END IF
                CALL MPI_Allreduce(MPI_IN_PLACE, Aval, 1, MPINTREAL, MPI_SUM, &
                     & column_comm, grid_error)
  -             CALL DotOneHelper(recv_num_values,recv_index,recv_values, &
  -                  & values_per_column_l(local_pi_i), index_l(:,local_pi_i), &
  -                  & values_l(:,local_pi_i),insert_value,column_comm)
                !! Insert Into L
  -             insert_value = inverse_factor * (Aval - insert_value)
  +             fill_counter = fill_counter + 1
  +             insert_value = inverse_factor * (Aval - dot_values(fill_counter))
                IF (JJ .GE. AMat%start_row .AND. JJ .LT. AMat%end_row) THEN
                   CALL AppendToVector(values_per_column_l(local_pi_i), &
                        & index_l(:,local_pi_i), values_l(:, local_pi_i), &
  
positive.code:        CALL MPI_Allreduce(MPI_IN_PLACE, inverse_factor, 1, MPINTREAL, MPI_SUM, &
              & row_comm, grid_error)
  
         !! Loop over other columns
         DO II = JJ + 1, AMat%actual_matrix_dimension
            pi_i = pivot_vector(II)
            local_pi_i = pi_i - AMat%start_column + 1
  
negative.code:     DEALLOCATE(values_per_column_l)
      DEALLOCATE(index_l)
      DEALLOCATE(values_l)
      DEALLOCATE(values_per_column_pi)
      DEALLOCATE(index_pi)
      DEALLOCATE(values_pi)
      DEALLOCATE(recv_index)
      DEALLOCATE(recv_values)
      DEALLOCATE(dot_values)
  