anchor.diff:         !! Loop over other columns
          DO II = 1, num_local_pivots
             !! Insert Into L
  -          Aval = a_buf(II)
  +          local_pi_i = local_pivots(II)
  +          Aval = a_buf(local_pi_i)
             insert_value = inverse_factor * (Aval - dot_values(II))
  -          local_pi_i = local_pivots(II)
             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:     INTEGER :: local_pi_i, local_pi_j
      REAL(NTREAL) :: Aval, insert_value, inverse_factor
      REAL(NTREAL), DIMENSION(:), ALLOCATABLE :: dot_values
      INTEGER :: row_root, col_root
  
      !! Optional Parameters
      IF (PRESENT(solver_parameters_in)) THEN
  
negative.code:        END IF
         col_root = col_root_lookup(pi_j)
  
         !! Broadcast column JJ, and Inverse Factor
         CALL BroadcastVector(recv_num_values, recv_index, recv_values, &
              & col_root, row_comm)
         CALL MPI_Bcast(inverse_factor, 1, MPINTREAL, col_root, row_comm, &
              & grid_error)
  
         !! Extract the row of A to a dense matrix for easy lookup
         DO II = MAX(acol%outer_index(pi_j),1), acol%outer_index(pi_j+1)
            a_buf(acol%inner_index(II)) = acol%values(II)
         END DO
  
         !! Compute Dot Products
         CALL DotAllPivoted(recv_num_values, recv_index, recv_values, &
  