giving more strict convergence criterias for both linear and non-linear problems.
mainly adding the DratesDxi terms to the Jacobian formulation.
changing the residual calculation so that the 1/dt * M + K form is taken.
turn off the comments. Attention the local problem is not correct yet.
create memory space to store drates_dxi.
set looser criteria for local problem.