acroy

@jiahao : Thanks for the comments. I will fix the typos in the next iteration. The indentation issue is odd, because locally I have four spaces. Is there any git(hub)...

@kmsquire : Thanks! Indentation fixed in last commit.

@tlycken : that would be great. I should add that from reading a bit about adaptive symplectic methods, I got the impression that methods based on rescaling time (like [here](http://www.math.ku.edu/~huang/research/paper/HL97.pdf))...

Good question ;-) I guess there might be two reasons: 1) Although the method works it is very crude IMO and 2) We never settled the API - for example,...

I agree that this might be useful to have, although I am not sure about the relations between tolerances, minstep and maxIterations (actually, refinements would be a better name). Could...

@ahwillia : Returning `Vector{typeof(y0)}` was actually a conscious choice for the API, discussed in #20. As Tomas explained recently in #78, the reason is that we also want to support...

> Anyway, I also find this a bit annoying. Could the 2D array be the standard for 1D array IC input? And only the vector of vectors for scalar-like things?...

I am not saying that matrices are bad, it's just that `Vector{typeofy0)}` is a reasonable choice and gives a consistent API. I would also argue that the performance of the...

@yuyichao : Don't get me wrong, I don't insist on the current behavior. However, we went through some discussions (#20) to define the current API, so we shouldn't make fast...

The plan was to use an iterator interface for the low level functions. Then one could do `for t,y in OdeXX(F,y0,tspan, args) ... end` and have and option for "event...