Now the result of HC_decay_1Du exmaple is correct. Do not need full output. Just every 10 steps.
only need every 10 time steps output.
set the nonlinear iteration convergence criteria to 1e-12.
use the same output time span and nonlinear iteration convergence criteria.
leave the mass lumping terms on, but correct the localF factorization.