Mark,
interesting idea. Given the fact that in 2-d euclidean metric, the
Einstein summation conventions are only a way to write out
conventional matrix multiplications, do you consider at some point to
include a non-euclidean metric in this thing? (As you have in special
relativity, for example)
Something along the lines:
eta = np.diag(-1,1,1,1)
a = np.array(1,2,3,4)
b = np.array(1,1,1,1)
such that
einsum('i,i', a,b, metric=eta) = -1 + 2 + 3 + 4
I don't know how useful it would be, just a thought,
Hanno
Am 26.01.2011 um 21:27 schrieb Mark Wiebe:
> I wrote a new function, einsum, which implements Einstein summation
> notation, and I'd like comments/thoughts from people who might be
> interested in this kind of thing.
_______________________________________________
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion
