added preallocation of productivity ot singlefirm to avoid error at runtime
Merging in work from last night at home
Merge branch 'master' of https://github.com/veryshuai/EEJKT_est
increased possible matrix size for a fairer calibration