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)
