anchor.diff:      INTEGER, DIMENSION(:), ALLOCATABLE :: values_per_column_l
       INTEGER, DIMENSION(:,:), ALLOCATABLE :: index_l
       REAL(NTREAL), DIMENSION(:,:), ALLOCATABLE :: values_l
  -    INTEGER, DIMENSION(:), ALLOCATABLE :: values_per_column_pi
  -    INTEGER, DIMENSION(:,:), ALLOCATABLE :: index_pi
  -    REAL(NTREAL), DIMENSION(:,:), ALLOCATABLE :: values_pi
       !! Broadcasted Column
       INTEGER, DIMENSION(:), ALLOCATABLE :: recv_index
       REAL(NTREAL), DIMENSION(:), ALLOCATABLE :: recv_values
  
positive.code:     ALLOCATE(index_l(sparse_a%rows, sparse_a%columns))
      ALLOCATE(values_l(sparse_a%rows, sparse_a%columns))
      values_per_column_l = 0
      ALLOCATE(values_per_column_pi(sparse_a%columns))
      ALLOCATE(index_pi(sparse_a%rows, sparse_a%columns))
      ALLOCATE(values_pi(sparse_a%rows, sparse_a%columns))
      values_per_column_pi = 0
      ALLOCATE(a_buf(sparse_a%columns))
  
      !! Allocate space for a received column
  
negative.code:        !! Loop over other columns
         fill_counter = 0
         DO II = 1, num_local_pivots
            local_pi_i = local_pivots(II)
            !! Insert Into L
            Aval = a_buf(II)
            insert_value = inverse_factor * (Aval - dot_values(II))
  