Discussion: new operators for numerical computation

[Gregory Lielens]

I was probably wrong in my first attempt to define kronecker product,
but something like 

 def __specialproduct__(self,b):
   rank4tensor = self.i(2,4) (*) b.i(1,3)
   return(rank4tensor[RemoveAxis,:,RemoveAxis,:])

may be ok... if for a(...,n,m,...), b=a[...,RemoveAxis,:,...] produce a
tensor b(...,n*m,...)

Enough of this, my head is spinning :-)
Anyway, element and inner product should be sufficient as infix operator
in a first step imho, in a second we should
be able to specify generalized inner product where the indices are
specified, and be able to create/destroy dimensions in arrays, a la
NumPy. 

Greg.