On Thu, 15 Mar 2018, oleg_komashko(a)ukr.net wrote:
Dear all,
Today's git help entry for max
does not mention matrix argument
But
set seed 13
matrix a_matrix = mnormal(10^4,10^4)
set stopwatch
ma1 = max(vec(a_matrix))
eval $stopwatch
set stopwatch
ma2 = maxc(maxr(a_matrix))
eval $stopwatch
I have 1.6487076 sec for max(vec..
and 7.9388368 sec for the second version
And the first version seems to be doomed
for deprecation
Not necessarily. We have in mind (well, not fully discussed yet) that
max(a_matrix) as revised would return the largest element in a_matrix
(regardless of rows and columns). It would therefore be equivalent to
but faster than max(vec(a_matrix)).
Allin