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

Reply via email to