How can we create a Mathematica function that can accept 0 or more arguments in the same way the Times function does?

Specially, I'd like to replace the Times function to work with Reals, Complexes, Quaternions or 4x4 matrices.

It looks like we can do the following:

MySpecialFunction[x1_, x2_, ..., xn_] := ... (* I don't know how to create this function yet *)

Times = MySpecialFunction;

True == (Times[1, 2, 3, 4, 5] == MySpecialFunction[1, 2, 3, 4, 5]);
True == (Times[2, 3] == MySpecialFunction[2, 3]);
True == (Times[] == MySpecialFunction[]);

If MySpecialFunction is, for example, Plus, it works.


(* It looks like ## can help us to get all arguments, but I don't know how to use it to create the function. *)
In[1]:= f[##] &[a, b, c, d]
Out[1]= f[a, b, c, d]
  • $\begingroup$ "Quaternions or 4x4 matrices" - Times[] is commutative, so you probably want to be careful there. $\endgroup$ – J. M.'s discontentment May 6 at 14:44
  • $\begingroup$ Looks like a duplicate of this. $\endgroup$ – Rohit Namjoshi May 6 at 14:48
  • $\begingroup$ @J.M., I agree. It is just that I'd like to unify the * (Times), ** (NonCommutativeMultiply) and . (Dot) operations at one. But I agree that I need to be careful. I'd like to apply the same concept to other functions as well. Like Sin[matrix] be actually MatrixFunction[Sin, matrix] instead, to simplify my analysis. $\endgroup$ – GarouDan May 6 at 14:54
  • $\begingroup$ For another approach, look up BlankSequence and BlankNullSequence for pattern matching the many arguments in the function's definition. $\endgroup$ – eyorble May 6 at 15:13


I assume in the following that we have a multiplication operation mytimes which is non-commutative, an addition operation Plus which is abelian and scalar multiplication Times.

We will want our operation to be distributive and linear this is the tricky thing in my eyes.

Let's code it

My favorite Mathematica book is Power Programming with Mathematica - The Kernel by David B. Wagner. You can find a nice example how one could code a linear operator there. This case is very similar.

If we have something like mytimes[c*somehead[z]] and if ourc is a scalar (any element of $\mathbb{C}$), then we want to move it out mytimes (homogeneity). We also want to preserve the order in case we have multiple arguments.

pattern1 = c1_*c2_ /; Element[c1, Complexes]; 
mytimes[y___] /; Length[Cases[{y}, pattern1]] != 0 := {y} /. {x1___, x2_*x3_, x4___} /; Element[x2, Complexes] :> x2*mytimes[x1, x3, x4];

The left hand side of the upper code ensures that we only apply the definition if we have at least one element where at least one scalar is multiplied with something else.

This will already simplify mytimes[3*somehead[y], 4*somehead[z]] to

12 mytimes[somehead[y],somehead[z]]

Next is distributivity:

pattern2 = c1_ + c2_;
mytimes[y___] /; Length[Cases[{y}, pattern2]] != 0 := {y} /. {x1___, x2_ + x3_, x4___} :> mytimes[x1, x2, x4] + mytimes[x1, x3, x4];

Analogous to the previous case we only want to match if at least one argument consists of a sum. If we have such an argument we let the pattern matching engine do the complicated stuff on the right hand side.

Lets see some stuff

With the upper definitions

mytimes[3*somehead[y] + 4*somehead[z], 5*somehead[x], 6*somehead[w]]

simplifies to:

30 (3 mytimes[somehead[y],somehead[x],somehead[w]]+4 mytimes[somehead[z],somehead[x],somehead[w]])

With just one argument to mytimes it gives

 mytimes[3*somehead[y] + 4*somehead[z]]
3 mytimes[somehead[y]]+4 mytimes[somehead[z]]


You could now add more rules to mytimes to define multiplication of quaternions, 4x4 matrices or other non-commmutative operations. You can define mytimes for all of them at once as stated in the question by introducing rules for specific heads.

| improve this answer | |

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.