From 024867f3da16c7b51a768399401de1edcd417384 Mon Sep 17 00:00:00 2001 From: Peter Mikkelsen Date: Wed, 26 Jan 2022 13:51:28 +0000 Subject: Add a better implementation of inner product, and don't simplify nested arrays if the contents are also arrays --- operators.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'operators.c') diff --git a/operators.c b/operators.c index 2280e13..8fd9911 100644 --- a/operators.c +++ b/operators.c @@ -173,9 +173,9 @@ opInnerProduct(Datum *lefto, Datum *righto, Array *left, Array *right) throwerror(nil, EType); if(left->rank > 0 && right->rank > 0 && left->shape[left->rank-1] != right->shape[0]) - throwerror(L"Last dimension of A must match first dimension of B in A f.g B", EShape); + throwerror(nil, ELength); - return rundfn(L"⍶/¨(↓⍺)⍹⌾⍉↓⍉⍵", lefto, righto, left, right); + return rundfn(L"↑(↓⍺)(⍶/⍹¨)⌾⍉↓⍉⍵", lefto, righto, left, right); } Array * -- cgit v1.2.3