@asmeurer : Could you tell me how the prototype of the function should be, because I tested a few partial differential equations with the function, and it seemed to give the correct output. I should have missed something.