This line would really benefit from parameterizing q in NN.
Right now this.q * Yd is unstable and the append is slow.
Would also help to pre-allocate Ydata to it's full size. One would just need to first count the non-zero entries in this.outTimes with something like:
length(find(x -> x==1, this.outTimes))
https://github.com/XtractOpen/Meganet.jl/blob/master/src/integrators/NN.jl#L82
This line would really benefit from parameterizing
q
inNN
.Right now
this.q * Yd
is unstable and the append is slow.Would also help to pre-allocate Ydata to it's full size. One would just need to first count the non-zero entries in this.outTimes with something like:
length(find(x -> x==1, this.outTimes))