anchor.diff:    integer                        :: shift, shift2, itermax
     include 'constants.include.F'
     
  -  !DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: U, W, R, S, y, h, lambda
  +  !DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: U, W, S, y, h, lambda
     if (N_st_diag*3 > sze) then
        print *,  'error in Davidson :'
        print *,  'Increase n_det_max_jacobi to ', N_st_diag*3
  
positive.code:   U = 0.d0
    W = 0.d0
    S = 0.d0
    R = 0.d0
    y = 0.d0
  
  
  
negative.code:   deallocate (                                                       &
        W, residual_norm,                                              &
        U,                                                             &
        c, S,                                                          &
        h,                                                             &
        y, s_, s_tmp,                                                  &
        lambda                                                         &
  