Jianghui Du

Where do I get the latest? Is it already registered?

Still doesn't work... Again, CVODE+LapackBand takes 3 mins to get to t=13000. FBDF + BandedMatrices It seems to be working. But very slow. I waited for half an hour and...

> > I used parametric type in BeelerReuterCpu to make the type consistent between u and the fields in BeelerReuterCpu . I assume that's necessary for the Dual Number to...

So basically all the arrays in the `struct` are caches? I thought it refers to temporary caches allocated during calculations. ```julia mutable struct BeelerReuterCpu{T}

I created all the caches and pass them as parameters, and now it works. However, I had trouble with the `log` functions in the code ( inside `calc_iCa` and `update_C_cpu`)....

Decreasing `reltol` and `abstol` down to `1e-15` doesn't make a difference. However, decreasing `dtmax` does make the negative number in `ERROR: DomainError with (a negative number）` smaller , like with...

Also, it seems the `DiffEqBase.dualcache` can only be passed as parameters `p` to the ODE function `(f::BeelerReuterCpu)(du, u, p, t)`, and cannot go inside the definition of `mutable struct BeelerReuterCpu`...

I can't get it work to put the cache inside the struct. If I only add a cache inside like ```Julia mutable struct BeelerReuterCpu{T}

I realized that in `XI=DiffEqBase.get_tmp(f.XI_cache, u)` the left hand side doesn't need to be allocated. I was trying to allocate for that and that was the problem. Now this works...

I cann't get it work on non-uniform grid ```julia # grid set up L = 50.0 # domain [0,L] N = 1000 # number of grid points # a uniform...