www.neuron.yale.edu

I guess I understand why sparse needs converge and the stop condition is |rhs| = 0 (L1 norm of rhs equals to zero). Sparse method uses "backward eular" to solve equations to get new state at t+dt. What's more, instead of setting the nonlinear system like "AX = b"(both A and b are...

When using method sparse to solve kinetic schame, I found that the code in scopemath/sparse.c solved the coefficient matrix several times to get converge when linflag=0. I am really confused with that. If you call the (*fun) which is the state function in the c file compiled from mod file N times, d...