anchor.diff:      REAL(our_dble), ALLOCATABLE     :: W(:)     ! Working space    
       INTEGER(our_int), ALLOCATABLE   :: JW(:)    ! Working space 
       INTEGER(our_int)                :: L_JW     ! Working space 
  -    INTEGER(our_int)                :: L_JW     ! Working space 
  +    INTEGER(our_int)                :: L_W      ! Working space 
   
       INTEGER(our_int)                :: MINEQ    ! Locals
       INTEGER(our_int)                :: N1       ! Locals
  
positive.code:                         mapping_state_idx, states_all, num_periods, & 
                          periods_emax, delta, shocks_cholesky)
  
              A(1,:2) = divergence_derivative(X, shocks_cov, &
                              level, tiny)
  
          END IF
  
          !Call to SLSQP code
          CALL slsqp(M, MEQ, LA, N, X, XL, XU, F, C, G, A, ACC, &
                  ITER, MODE, W, L_W, JW, L_JW)
  
          ! Check if SLSQP has completed
  
negative.code:     LOGICAL                         :: is_finished
      LOGICAL                         :: is_success
  
      !/* SLSQP interface          */
  
      INTEGER(our_int)                :: M        ! Total number of constraints
      INTEGER(our_int)                :: MEQ      ! Total number of equality constraints
  