anchor.diff:  !------------------------------------------------------------------------------
   
       ! Construct auxiliary objects
  -    num_obs = SIZE(data_array, 1)
  +    num_obs = SIZE(data_evaluate, 1)
   
       shocks_cov = MATMUL(shocks_cholesky, TRANSPOSE(shocks_cholesky))
       is_deterministic = ALL(shocks_cov .EQ. zero_dble)
  
positive.code: 
      REAL(our_dble), INTENT(IN)      :: periods_payoffs_systematic(num_periods, max_states_period, 4)
      REAL(our_dble), INTENT(IN)      :: periods_draws_prob(num_periods, num_draws_prob, 4)
      REAL(our_dble), INTENT(IN)      :: data_array(:, :)
      REAL(our_dble), INTENT(IN)      :: shocks_cholesky(4, 4)
      REAL(our_dble), INTENT(IN)      :: periods_emax(num_periods, max_states_period)
  
  
negative.code:  CONTAINS
  !******************************************************************************
  !******************************************************************************
  SUBROUTINE fort_evaluate(rslt, periods_payoffs_systematic, mapping_state_idx, periods_emax, states_all, shocks_cholesky, data_evaluate, periods_draws_prob)
  
      !   DEVELOPMENT NOTES
      !
  