Non-symmetric IRFs in recursive utility

Dear friends,
I try to investigate the economic dynamics under recursive utility set in the paper Caldara et al. (2012). I use the code shared by prof.Pfeifer from his homepage.
I omit the volatility shock and focus on level of productivity shock. I set calculate the IRFs to one positive standard error shocks and one negative standard error shocks under the order=3.
In theory, as I take 3 order perturbation, the irfs to positive and negative shocks will be different, but the code present that they are same. My question is that, should they be different or same? If it is different, how can I get the result? If it is same, how to explain?
Caldara_et_al_2012.mod (5.8 KB)

They should not be the same. See

for code how to generate state dependent IRFs.

Thank you very much!