anchor.diff:      DEALLOCATE(recv_index)
       DEALLOCATE(recv_values)
       DEALLOCATE(dot_values)
  +    DEALLOCATE(col_root_lookup)
       CALL DestructSparseMatrix(sparse_a)
     END SUBROUTINE CholeskyDecomposition
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  
positive.code:     INTEGER :: II, JJ, local_JJ, local_II, local_row
      REAL(NTREAL) :: Aval, insert_value, inverse_factor
      REAL(NTREAL), DIMENSION(:), ALLOCATABLE :: dot_values
      INTEGER :: root
  
      !! Optional Parameters
      IF (PRESENT(solver_parameters_in)) THEN
  
negative.code:           inverse_factor = 1.0/insert_value
            !! Insert
            IF (JJ .GE. AMat%start_row .AND. JJ .LT. AMat%end_row) THEN
               CALL AppendToVector(values_per_column_l(local_pi_j), &
                    & index_l(:,local_pi_j), values_l(:, local_pi_j), &
                    & local_JJ, insert_value)
            END IF
            !! Extract column J for sending later
            recv_num_values = values_per_column_l(local_pi_j)
            recv_index(:recv_num_values) = index_l(:recv_num_values,local_pi_j)
            recv_values(:recv_num_values) = values_l(:recv_num_values, local_pi_j)
         END IF
         col_root = col_root_lookup(pi_j)
  
  