Sundials KLU Linear solver

You have to add the diagonal to it since it’s I - gamma*J