GeneralOuterProduct(X,f,g,Y) accepts two Arrays, X and Y, with N and M dimensions respectively. The result is an Array of N+M dimensions.

For each combination of indices of X and indices of Y, the function f is applied to the elements from X and Y.

For example, when X and Y are vectors of integers from 1 to 10, and f is the multiplication function (`*`), GeneralOuterProduct produces a matrix representing a traditional multiplication table.

If X or Y is a list, it is treated as one-dimensional Array.

If either X or Y is not an Array or list, it is interpreted as an Array with zero dimensions (and the value X or Y).

By default, the datatype of the resulting Array is determined by SuggestedDatatype based on the two inputs. Depending on f, this may not be suitable. For example, if X and Y have datatype integer[1] and f is `*`, the computed values might be too large to store in that datatype.

The option ':-datatype'=dt, where dt is a valid Array datatype, can be used to override the default datatype for the result.

If option evalhf or evalhf=true is passed, the arguments passed to each call to f are first evaluated to hardware floating-point using evalhf, as is the result returned by f. Furthermore, unless the datatype option was specified, the resulting Array has datatype float[8].

The subtype of the result depends on the subtypes of X and Y, as described by SuggestedSubtype. The option ':-subtype'=st', where st is a valid rtable subtype, overrides the default subtype.

The storage order of the result is based on the storage orders of X and Y, as described by SuggestedOrder. This can be overridden using the option :-order=ord, where ord is one of C_order or Fortran_order.

This function is part of the ArrayTools package, so it can be used in the short form GeneralOuterProduct(..) only after executing the command with(ArrayTools). However, it can always be accessed through the long form of the command by using ArrayTools[GeneralOuterProduct](..).

$\mathrm{with}\left(\mathrm{ArrayTools}\right)\:$

$V≔⟨\mathrm{`\$`}\left(1..10\right)⟩\:$

$\mathrm{GeneralOuterProduct}\left(V\,\mathrm{`*`}\,V\right)$

$R≔⟨1\,2\,3\,4\,5\,6⟩\:$

$\mathrm{R2}≔\mathrm{GeneralOuterProduct}\left(R\,\mathrm{`+`}\,R\right)$

$\mathrm{R2}≔\left[\right]$

$\mathrm{whattype}\left(\mathrm{R2}\right)$

$\mathrm{Matrix}$

$\mathrm{R3}≔\mathrm{GeneralOuterProduct}\left(\mathrm{R2}\,\mathrm{`+`}\,R\right)$

$\mathrm{whattype}\left(\mathrm{R3}\right)$

$\mathrm{Array}$

$\mathrm{R3}\left[1\right]$

$\mathrm{R3}\left[6\right]$

$\mathrm{map2}\left(\mathrm{numboccur}\,\mathrm{R3}\,\left[\mathrm{`\$`}\left(1..18\right)\right]\right)$

$\left[0\,0\,1\,3\,6\,10\,15\,21\,25\,27\,27\,25\,21\,15\,10\,6\,3\,1\right]$

$A≔\mathrm{Vector}\left[\mathrm{row}\right]\left(\left[4\,8\,12\right]\,\mathrm{datatype}=\mathrm{integer}\left[1\right]\right)$

$A≔\left[\right]$

$\mathrm{GeneralOuterProduct}\left(A\,\mathrm{`*`}\,A\right)$

Error, (in ArrayTools:-GeneralOuterProduct) unable to store '144' when datatype=integer[1]

$\mathrm{GeneralOuterProduct}\left(A\,\mathrm{`*`}\,A\,\mathrm{datatype}=\mathrm{integer}\left[2\right]\right)$

$V≔\mathrm{Array}\left(\left[\mathrm{seq}\left(i\,i=0..1\,\mathrm{numelems}=101\right)\right]\,\mathrm{datatype}=\mathrm{float}\left[8\right]\right)\:$

$V\left[1..4\right]$

$V\left[-4..-1\right]$

$\mathrm{img}≔\mathrm{GeneralOuterProduct}\left(V\,\max \,V\right)$

$\mathrm{img}\left[1..4,1..4\right]$

$\mathrm{img}\left[-4..-1,-4..-1\right]$

$\mathrm{with}\left(\mathrm{ImageTools}\right)\:$

$\mathrm{Embed}\left(\mathrm{img}\right)$

The ArrayTools[GeneralOuterProduct] command was introduced in Maple 2021.

For more information on Maple 2021 changes, see Updates in Maple 2021.