On Saturday, January 3, 2015 11:56:20 PM UTC-5, Jiahao Chen wrote: > > > On Sat, Jan 3, 2015 at 11:27 PM, Joshua Adelman <[email protected] > <javascript:>> wrote: > >> PS - maybe it's my experience with numpy, but I wonder if any thought has >> been given to following more of a numpy-style convention of always >> returning a view of an array when possible rather than a copy? > > > Search the Github issues for ArrayViews. > > All 3 versions of your Julia code have assignments of the form > > v2[:, k] = p > > which could be faster if you explicitly devectorize. > > Thanks, > > Jiahao Chen > Staff Research Scientist > MIT Computer Science and Artificial Intelligence Laboratory >
Hi Jiahao, Actually when I devectorized vander2 (explicitly looping over the first dimension as the inner-most loop) as suggested, while cases where N < 100 are indeed faster, for N > 100 the devectorized code is slower. Josh
