I was wondering if that would get quicker if it was first run on a smaller matrix? JIT and all that...
On Thursday, January 30, 2014 6:15:40 PM UTC+1, Stefan Karpinski wrote: > > Yes, that's the plan. It's a tricky change. > > > On Thu, Jan 30, 2014 at 12:11 PM, Johan Sigfrids > <[email protected]<javascript:> > > wrote: > >> Couldn't Julia provide a view by default when indexing? You could then >> copy if you want a new one. >> >> >> On Thursday, January 30, 2014 6:55:34 PM UTC+2, Tim Holy wrote: >>> >>> On Thursday, January 30, 2014 07:32:18 AM John Myles White wrote: >>> > This is pretty standard fare for Julia. Things like sum are really >>> wasteful >>> > with memory, whereas the nuclear option is very conservative when >>> > implemented right. >>> >>> To be fair, it's not sum() that's to blame, the problem is allocating a >>> new >>> array with A[m-f:m+f, indx]. >>> >>> --Tim >>> >>> >
