I was not clear in previous mail the issue is that by rayleight C non-diagonal i get matrix-matrix multilication thatbreaks the sparsity. I am not sure what kind of approximation i should use, approximate matrix-matrix, or go use diagonal c matrix or do anyone of you a better suggestion to involve damping into the dynamical system.
Best Regards Johan _______________________________________________
