anchor.diff:  
           use calc, only: calc_type, ccmc_calc
           use system, only: set_spin_polarisation
  +        use report, only: wrapper_end_report
   
           integer(c_int) :: nresult
           type(c_ptr), value :: L
  
positive.code:         type(qmc_state_t), pointer :: qmc_state_out, qmc_state_restart
  
          logical :: have_restart_state
          integer :: opts, ierr
          real(p), allocatable :: sampling_probs(:)
          character(10), parameter :: keys(9) = [character(10) :: 'sys', 'qmc', 'dmqmc', 'ipdmqmc', 'operators', 'rdm', &
                                                                  'restart', 'reference', 'qmc_state']
  
          lua_state = flu_copyptr(L)
          call get_sys_t(lua_state, sys)
  
negative.code: 
          use calc, only: calc_type, mc_canonical_estimates
          use canonical_estimates, only: estimate_canonical
          use report, only: wrapper_end_report
  
          integer(c_int) :: nresult
          type(c_ptr), value :: l
  