Outer product

In linear algebra, the outer product of two coordinate vectors is a matrix. If the two vectors have dimensions ''n'' and ''m'', then their outer product is an ''n'' × ''m'' matrix. More generally, given two tensors (multidimensional arrays of numbers), their outer product is a tensor. The outer product of tensors is also referred to as their tensor product, and can be used to define the tensor algebra.

The outer product contrasts with:
* The dot product (a special case of "inner product"), which takes a pair of coordinate vectors as input and produces a scalar
* The Kronecker product, which takes a pair of matrices as input and produces a block matrix
* Standard matrix multiplication

Definition

Given two vectors of size $m \times 1$ and $n \times 1$ respectively
$\mathbf = \begin u_1 \\ u_2 \\ \vdots \\ u_m \end, \quad \mathbf = \begin v_1 \\ v_2 \\ \vdots \\ v_n \end$
their outer product, denoted $\mathbf \otimes \mathbf,$ is defined as the $m \times n$ matrix $\mathbf$ obtained by multiplying each element of $\mathbf$ by each element of $\mathbf$:
$\mathbf \otimes \mathbf = \mathbf = \begin u_1v_1 & u_1v_2 & \dots & u_1v_n \\ u_2v_1 & u_2v_2 & \dots & u_2v_n \\ \vdots & \vdots & \ddots & \vdots \\ u_mv_1 & u_mv_2 & \dots & u_mv_n \end$

Or in index notation:
$(\mathbf \otimes \mathbf)_ = u_i v_j$

Denoting the dot product by $\,\cdot,\,$ if given an $n \times 1$ vector $\mathbf,$ then $(\mathbf \otimes \mathbf) \mathbf = (\mathbf \cdot \mathbf) \mathbf.$ If given a $1 \times m$ vector $\mathbf,$ then $\mathbf (\mathbf \otimes \mathbf) = (\mathbf \cdot \mathbf) \mathbf^.$

If $\mathbf$ and $\mathbf$ are vectors of the same dimension bigger than 1, then $\det (\mathbf \otimes\mathbf) = 0$.

The outer product $\mathbf \otimes \mathbf$ is equivalent to a matrix multiplication $\mathbf \mathbf^,$ provided that $\mathbf$ is represented as a $m \times 1$ column vector and $\mathbf$ as a $n \times 1$ column vector (which makes $\mathbf^$ a row vector). For instance, if $m = 4$ and $n = 3,$ then
$\mathbf \otimes \mathbf = \mathbf\mathbf^\textsf = \beginu_1 \\ u_2 \\ u_3 \\ u_4\end \beginv_1 & v_2 & v_3\end = \begin u_1 v_1 & u_1 v_2 & u_1 v_3 \\ u_2 v_1 & u_2 v_2 & u_2 v_3 \\ u_3 v_1 & u_3 v_2 & u_3 v_3 \\ u_4 v_1 & u_4 v_2 & u_4 v_3 \end.$

For complex vectors, it is often useful to take the conjugate transpose of $\mathbf,$ denoted $\mathbf^\dagger$ or $\left(\mathbf^\textsf\right)^*$:
$\mathbf \otimes \mathbf = \mathbf \mathbf^\dagger = \mathbf \left(\mathbf^\textsf\right)^*.$

Contrast with Euclidean inner product

If $m = n,$ then one can take the matrix product the other way, yielding a scalar (or $1 \times 1$ matrix):
$\left\langle\mathbf, \mathbf\right\rangle = \mathbf^\textsf \mathbf$
which is the standard inner product for Euclidean vector spaces, better known as the dot product. The dot product is the trace of the outer product. Unlike the dot product, the outer product is not commutative.

Multiplication of a vector $\mathbf$ by the matrix $\mathbf \otimes \mathbf$ can be written in terms of the inner product, using the relation $\left(\mathbf \otimes \mathbf\right)\mathbf = \mathbf\left\langle\mathbf, \mathbf\right\rangle$.

The outer product of tensors

Given two tensors $\mathbf, \mathbf$ with dimensions $(k_1, k_2, \dots, k_m)$ and $(l_1, l_2, \dots, l_n)$, their outer product $\mathbf \otimes \mathbf$ is a tensor with dimensions $(k_1, k_2, \dots, k_m, l_1, l_2, \dots, l_n)$ and entries
$(\mathbf \otimes \mathbf)_ = u_ v_$

For example, if $\mathbf$ is of order 3 with dimensions $(3, 5, 7)$ and $\mathbf$ is of order 2 with dimensions $(10, 100),$ then their outer product $\mathbf$ is of order 5 with dimensions $(3, 5, 7, 10, 100).$ If $\mathbf$ has a component and $\mathbf$ has a component , then the component of $\mathbf$ formed by the outer product is .

Connection with the Kronecker product

The outer product and Kronecker product are closely related; in fact the same symbol is commonly used to denote both operations.

If $\mathbf = \begin1 & 2 & 3\end^\textsf$ and $\mathbf = \begin4 & 5\end^\textsf$, we have:
$\begin \mathbf \otimes_\text \mathbf &= \begin 4 \\ 5 \\ 8 \\ 10 \\ 12 \\ 15\end, & \mathbf \otimes_\text \mathbf &= \begin 4 & 5 \\ 8 & 10 \\ 12 & 15\end \end$

In the case of column vectors, the Kronecker product can be viewed as a form of vectorization (or flattening) of the outer product. In particular, for two column vectors $\mathbf$ and $\mathbf$, we can write:
$\mathbf \otimes_ \mathbf = \operatorname(\mathbf \otimes_\text \mathbf)$

Note that the order of the vectors is reversed in the right side of the equation.

Another similar identity that further highlights the similarity between the operations is
$\mathbf \otimes_ \mathbf^\textsf = \mathbf u \mathbf^\textsf = \mathbf \otimes_ \mathbf$

where the order of vectors needs not be flipped. The middle expression uses matrix multiplication, where the vectors are considered as column/row matrices.

Connection with the matrix product

Given a pair of matrices $\mathbf$ of size $m\times p$ and $\mathbf$ of size $p\times n$, consider the matrix product $\mathbf = \mathbf\,\mathbf$ defined as usual as a matrix of size $m\times n$.

Now let $\mathbf a^\text_k$ be the $k$-th column vector of $\mathbf A$ and let $\mathbf b^\text_k$ be the $k$-th row vector of $\mathbf B$. Then $\mathbf$ can be expressed as a sum of column-by-row outer products:

$\mathbf = \mathbf\, \mathbf = \left( \sum_^p _\, _ \right)_ = \begin & & \\ \mathbf a^\text_ & \cdots & \mathbf a^\text_ \\ & & \end \begin & \mathbf b^\text_ & \\ & \vdots & \\ & \mathbf b^\text_ & \end = \sum_^p \mathbf a^\text_k \otimes \mathbf b^\text_k$
Note the duality of this expression with the more common one as a matrix built with row-by-column inner product entries (or dot product): $C_ = \langle\rangle$

This relation is relevant in the application of the Singular Value Decomposition (SVD) (and Spectral Decomposition as a special case). In particular, the decomposition can be interpreted as the sum of outer products of each left ($\mathbf_k$) and right ($\mathbf_k$) singular vectors, scaled by the corresponding nonzero singular value $\sigma_k$:
$\mathbf = \mathbf = \sum_^(\mathbf_k \otimes \mathbf_k) \, \sigma_k$

This result implies that $\mathbf$ can be expressed as a sum of rank-1 matrices with spectral norm $\sigma_k$ in decreasing order. This explains the fact why, in general, the last terms contribute less, which motivates the use of the Truncated SVD as an approximation. The first term is the least squares fit of a matrix to an outer product of vectors.

Properties

The outer product of vectors satisfies the following properties:
$\begin (\mathbf \otimes \mathbf)^\textsf &= (\mathbf \otimes \mathbf) \\ (\mathbf + \mathbf) \otimes \mathbf &= \mathbf \otimes \mathbf + \mathbf \otimes \mathbf \\ \mathbf \otimes (\mathbf + \mathbf) &= \mathbf \otimes \mathbf + \mathbf \otimes \mathbf \\ c (\mathbf \otimes \mathbf) &= (c\mathbf) \otimes \mathbf = \mathbf \otimes (c\mathbf) \end$

The outer product of tensors satisfies the additional associativity property:
$(\mathbf \otimes \mathbf) \otimes \mathbf = \mathbf \otimes (\mathbf \otimes \mathbf)$

Rank of an outer product

If u and v are both nonzero, then the outer product matrix uv^T always has matrix rank 1. Indeed, the columns of the outer product are all proportional to the first column. Thus they are all linearly dependent on that one column, hence the matrix is of rank one.

("Matrix rank" should not be confused with "tensor order", or "tensor degree", which is sometimes referred to as "rank".)

Definition (abstract)

Let and be two vector spaces. The outer product of $\mathbf v \in V$ and $\mathbf w \in W$ is the element $\mathbf v \otimes \mathbf w \in V \otimes W$.

If is an inner product space, then it is possible to define the outer product as a linear map . In which case, the linear map $\mathbf x \mapsto \langle \mathbf v, \mathbf x\rangle$ is an element of the dual space of . The outer product is then given by
$(\mathbf v \otimes \mathbf w) (\mathbf x) = \left\langle \mathbf v, \mathbf x \right\rangle \mathbf w$

This shows why a conjugate transpose of is commonly taken in the complex case.

In programming languages

In some programming languages, given a two-argument function f (or a binary operator), the outer product of f and two one-dimensional arrays A and B is a two-dimensional array C such that C, j= f(A B . This is syntactically represented in various ways: in APL, as the infix binary operator ∘.f; in J, as the postfix adverb f/; in R, as the function outer(A, B, f) or the special %o%; in Mathematica, as Outer, A, B/syntaxhighlight>. In MATLAB, the function kron(A, B) is used for this product. These often generalize to multi-dimensional arguments, and more than two arguments.

In the Python library NumPy, the outer product can be computed with function np.outer(). In contrast, np.kron results in a flat array. The outer product of multidimensional arrays can be computed using np.multiply.outer.

Applications

As the outer product is closely related to the Kronecker product, some of the applications of the Kronecker product use outer products. These applications are found in quantum theory, signal processing, and image compression.

Spinors

Suppose so that and are in . Then the outer product of these complex 2-vectors is an element of , the 2 × 2 complex matrices:
$\begin sw & tw \\ sz & tz \end.$
The determinant of this matrix is because of the commutative property of C.

In the theory of spinors in three dimensions, these matrices are associated with isotropic vectors due to this null property. Élie Cartan described this construction in 1937, but it was introduced by Wolfgang Pauli in 1927 so that M(2, C) has come to be called Pauli algebra.

Concepts

The block form of outer products is useful in classification. Concept analysis is a study that depends on certain outer products:

When a vector has only zeros and ones as entries, it is called a ''logical vector'', a special case of a logical matrix. The logical operation and takes the place of multiplication. The outer product of two logical vectors and is given by the logical matrix $\left(a_\right) = \left(u_i \land v_j\right)$. This type of matrix is used in the study of binary relations, and is called a rectangular relation or a cross-vector.Ki Hang Kim (1982) ''Boolean Matrix Theory and Applications'', page 37, Marcel Dekker

See also

* Dyadics
* Householder transformation
* Norm (mathematics)
* Scatter matrix
* Ricci calculus

Products

* Cartesian product
* Cross product
* Exterior product
* Hadamard product

Duality

* Complex conjugate
* Conjugate transpose
* Transpose
* Bra–ket notation for outer product

References

Further reading

*

{{DEFAULTSORT:Outer Product
Bilinear maps
Operations on vectors
Higher-order functions
Articles with example Python (programming language) code