On Tue, 19 Jul 2022, Thomas McGuire wrote:
Now the wierd thing is the timing for matmul3 is slow on par with emulating K
matrices in J
But if I use the fork directly it is as fast as Raul’s solution
You are measuring two completely different things.
10 timex'bmat1 ([: +/ *)"1 _ bmat2'
0.0482293
10 timex'bmat1 ([: +/ *"1 _) bmat2'
0.115583
(bmat1 ([: +/ *)"1 _ bmat2) -: (bmat1 ([: +/ *"1 _) bmat2)
0
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm