you just need to end the line with an opertor +, * etc so that Julia knows it must parse the next line
Thank you for the tips, I’ve learned something.
So, I decided to benchmark properly with the other code from Petsc.jl (I should have done that from the beginning…) and… The solutions are the same (slight difference, but I guess it is because my method is more accurate in time):
So I was looking for a bug that doesn’t exist.
Sorry for that, I didn’t expect the solution of my equations to behave exactly that way at these conditions. I’ve learned by experience that generally when it works on the first try, there is something fishy but it wasn’t the case here 
Thank you for your time!