@amabdelrehim I think it is a very bad idea to double all the code here that is already in D_psi.c. There should be a solution similar to what we do in the even/odd case, where we use macros to build these operators. Here we could just re-use the functions local_H for instance, etc. But it must not be doubled!