Hi Sebastien
There appears to be a slight problem with the current implementation of k_order_pert.m.
It is called from the top of dr1.m and dr1.m then returns after running it but it then does not appear to provide for the log-linear adjustment of ghx and ghu usually performed further down the line in dr1.m when if the model is log-linear - though I may be missing something.
There is also a C++ log-linear adjustment built in the C++ version of the dr1 for C++ DsgeLikelihood and that code may be integrated into k_order_perturbation.cpp too to run adjustment within C++ when if needed.
Best regards
George