def residuals(self, calib=None): if self.model_type in ("fg",'fga'): from dolo.algos.steady_state import residuals return residuals(self, calib) elif self.model_type in ('mfg','mfga'): from dolo.algos.mfg.steady_state import residuals return residuals(self, calib)
def residuals(self, calib=None): if self.model_type in ("fg", 'fga'): from dolo.algos.steady_state import residuals return residuals(self, calib) elif self.model_type in ('mfg', 'mfga'): from dolo.algos.mfg.steady_state import residuals return residuals(self, calib)
# create parameter range parm_vec = linspace(*range) # solve all calibrations (For the range of calibrations) for i in parm_vec: # Set calibrated parater value pp = {freeparam: i} model.set_calibration(**pp) ####################################### # Solve for the competitive equilibrium ######################################### model.set_calibration(po=0) calib_ce = find_steady_state(model) res = residuals(model, calib=calib_ce) if abs(sum(res['arbitrage'])) + abs(sum(res['transition'])) > 0.00001: print('model', case, 'ME', freeparam, '=', i, '--> steady state cant be found: non-zero residuals') dr_me_true = 0 else: dr_me_true = 1 dr_ce = perturb(model, order=2, steady_state=calib_ce) ########################################## # Solve for the social planner equilibrium ########################################## model.set_calibration(po=1) calib_po = find_steady_state(model)
def residuals(self, calib=None): from dolo.algos.steady_state import residuals return residuals(self, calib)
def residuals(self, calib=None): from dolo.algos.steady_state import residuals return residuals(self, calib)