[naffrancois]

Hello,

I would like to write a very simple implicit scheme to solve Euler equations, backward Euler. Have one of you ever used operator overloading in Fortran to build the Jacobian matrix ? This is the line I am following but:

Second question, how do you usually debug the Jacobian ? I mean a good clue is when you reach bigger cfl than explicit. But if it does not, or if the time integration simply crashes, how do you check you are accurately computing the Jacobian ?

By the way, I am testing the algorithm on Sod shock tube problem

Cheers