anchor.diff:    U = 0.d0
     W = 0.d0
     S = 0.d0
  -  R = 0.d0
     y = 0.d0
   
   
  
positive.code:   deallocate (                                                       &
        W, residual_norm,                                              &
        U,                                                             &
        R, c, S,                                                       &
        h,                                                             &
        y, s_, s_tmp,                                                  &
        lambda                                                         &
  
negative.code:   integer                        :: k_pairs, kl
    
    integer                        :: iter2
    double precision, allocatable  :: W(:,:),  U(:,:), S(:,:)
    double precision, allocatable  :: y(:,:), h(:,:), lambda(:), s2(:)
    double precision, allocatable  :: c(:), s_(:,:), s_tmp(:,:)
    double precision               :: diag_h_mat_elem
  