---
title: Rotating Gaussian Model
url: https://www.emergentmind.com/topics/rotating-gaussian-model
type: topic
---

# Rotating Gaussian Model

The term **Rotating Gaussian Model** does not denote a single universally standardized formalism. In the arXiv literature, it refers to several mathematically distinct constructions in which a Gaussian object is coupled to a notion of rotation: a Gaussian pressure anomaly in a rotating, stratified Boussinesq flow; a Gaussian surrogate obtained after an orthogonal data rotation for posterior approximation; iterative Gaussianization schemes that alternate marginal Gaussianization with orthonormal rotations; rotationally invariant planar Gaussian processes; rotating quantum Gaussian wave packets; and computer-vision models in which rotated boxes or dynamic 3D primitives are represented by Gaussian distributions [1605.06859] [1909.06753] [1602.00229].

## 1. Terminological scope and recurring structure

Across these usages, the common element is the Gaussian as the primary analytic or computational primitive, while “rotation” may refer to physical rotation, coordinate rotation, orthonormal mixing, or explicit orientation parameters.

| Research area | Gaussian object | Role of rotation |
|---|---|---|
| Rotating stratified flow | Gaussian pressure anomaly | Physical vortex equilibrium and stability |
| Bayesian inference | Gaussian approximation for nuisance projection | Orthogonal likelihood factorization |
| Gaussianization | Standard Gaussian target | Orthonormal mixing between marginal transforms |
| Stochastic processes | Planar Gaussian process | Rotational invariance and winding statistics |
| Quantum mechanics | Gaussian packet | Mean angular momentum and rotating dynamics |
| Computer vision | Gaussian box or 3D Gaussian primitive | Orientation, heading, or time-varying rotation |

This multiplicity matters because formally similar phrases encode different mathematical operations. In some cases rotation acts on the **state space** of the unknown density, as in RBIG; in others it acts on the **data space** to isolate nuisance structure, as in integrated rotated Gaussian approximation; and in others it is a **physical symmetry or kinematic variable**, as in vortices, wave packets, and dynamic Gaussian primitives [2206.03860] [2209.10839] [2412.04282].

## 2. Rotating Gaussian vortices in stratified Boussinesq flow

In geophysical and astrophysical fluid dynamics, the rotating Gaussian model is an exact steady vortex solution of the inviscid, non-hydrostatic Boussinesq equations on an \(f\)-plane with constant Coriolis parameter \(f\) and constant background Brunt–Väisälä frequency \(\bar N\). The total pressure and density are decomposed as
\[
p_{\rm tot}(r,z,t)=\bar p(z)+p(r,z,t),\qquad
\rho_{\rm tot}(r,z,t)=\bar\rho(z)+\rho(r,z,t),
\]
with
\[
\frac{d\bar p}{dz}=-\,\bar\rho\,g,\qquad
\bar N^2=-\,\frac{g}{\rho_o}\,\frac{d\bar\rho}{dz}=\text{const.}
\]
The vortex is specified by a Gaussian pressure anomaly
\[
p(r,z)=p_0\exp\!\Bigl[-(r/R)^2-(z/H)^2\Bigr],
\]
where \(p_0\) is the amplitude at the center, \(R\) is the horizontal length scale, and \(H\) is the vertical length scale [1605.06859].

Cyclo-geostrophic and hydrostatic balance yield the exact steady fields
\[
u_r=0,\quad u_z=0,
\]
\[
u_\phi(r,z)=\frac{f r}{2}\Biggl(-1+\sqrt{1-\frac{8p_0\,\exp[-(r/R)^2-(z/H)^2]}{\rho_o f^2 R^2}}\Biggr),
\]
and
\[
b(r,z)=-\,\frac{2p_0}{\rho_o H^2}\,z\,\exp\!\Bigl[-(r/R)^2-(z/H)^2\Bigr].
\]
At the center,
\[
\omega_c=\omega(r=0)=\frac{\partial (r u_\phi)}{r\,\partial r}\Big|_{r=0},\qquad
Ro=\frac{\omega_c}{2f},\qquad
N_c^2=\bar N^2+\Bigl(\partial b/\partial z\Bigr)_{r=0},
\]
with
\[
Ro=-\tfrac12+\sqrt{\tfrac14-\frac{2p_0}{\rho_o f^2 R^2}},\qquad
N_c^2=\bar N^2-\frac{2p_0}{\rho_o H^2}.
\]
The aspect ratio satisfies the universal scaling
\[
\Bigl(\tfrac{H}{R}\Bigr)^2
=\frac{-\,Ro(1+Ro)f^2}{\bar N^2\bigl[1-(N_c/\bar N)^2\bigr]}.
\]

After non-dimensionalization with
\[
Ro=\frac{\omega_c}{2f},\qquad
Bu=\bigl(NH/(fR)\bigr)^2,\qquad
\varepsilon=\frac{f}{\bar N},
\]
the linearized stability problem is formulated by normal modes
\[
(\tilde v_r,\tilde v_\phi,\tilde v_z,\tilde b,\tilde p)
=\boldsymbol g(r,z)\,\exp\bigl(\lambda t+i m\phi\bigr),\qquad
\lambda=\sigma-i m c,
\]
leading to a two-dimensional eigenvalue problem in \((r,z)\) for each azimuthal wavenumber \(m\) and each symmetry in \(z\). Numerically, the paper studies \(-0.5<Ro<0.5\) and \(0.02<Bu<2.3\), finding neutrally stable vortices only in a small region of parameter space: cyclones with \(Ro\sim0.02\text{–}0.05\) and \(Bu\sim0.85\text{–}0.95\). Anticyclones generally have slower growth rates than cyclones, and for \(Ro<0\) and \(0.5\lesssim Bu\lesssim1.3\) the most unstable eigenmode is slower than \(50\) turn-around times of the vortex. The results are numerically insensitive to reducing \(f/\bar N\) from \(0.1\) to \(0.01\), with eigenvalues and eigenvectors varying by \(<4\%\), because the fastest-growing modes satisfy to very good accuracy the local hydrostatic relation \(\partial_z\tilde p\approx\tilde b\) [1605.06859].

## 3. Rotated Gaussian approximations in Bayesian inference

In Bayesian regression with low-dimensional parameters of interest and high-dimensional nuisance structure, the rotated Gaussian model appears as the **Integrated Rotated Gaussian Approximation**. The starting point is
\[
y\sim N(X\theta+\eta,\sigma^2 I_n),
\]
where \(\theta\in\mathbb R^p\) is the parameter of interest and \(\eta\in\mathbb R^n\) is a high-dimensional nuisance term. Let \(X=Q R\) be a QR decomposition, with \(Q=[M\mid S]\), \(M^\top M=I_p\), \(S^\top S=I_{n-p}\), and \(M^\top S=0\). Rotating the data,
\[
\tilde y=Q^\top y=(M^\top y;\,S^\top y),
\]
gives
\[
M^\top y\sim N(M^\top X\theta+M^\top\eta,\sigma^2 I_p),\qquad
S^\top y\sim N(S^\top\eta,\sigma^2 I_{n-p}),
\]
so that
\[
p(y\mid \theta,\eta)=p(M^\top y\mid \theta,M^\top\eta)\cdot p(S^\top y\mid \eta).
\]
If \(u=M^\top\eta\), the key approximation is
\[
p(u\mid S^\top y)\approx N(u\mid m,\Sigma).
\]
Integrating out \(u\) yields
\[
\tilde\pi(\theta\mid y)\propto \pi(\theta)\,N(M^\top y\mid M^\top X\theta+m,\sigma^2 I_p+\Sigma),
\]
reducing the problem to a \(p\)-dimensional Gaussian linear model [1909.06753].

The paper establishes an approximation-error bound in KL divergence: the loss in approximating \(\theta\)'s posterior is controlled by how well the rotated nuisance law \(p(u\mid S^\top y)\) is approximated after convolution with Gaussian noise. Under the special case \(\eta=Z\alpha\), with possibly highly non-Gaussian prior on \(\alpha\), the projected law \(u=M^\top Z\alpha\) is shown to be approximately Gaussian when \(p\ll q\) and the conditional covariance of \(\alpha\) is not too degenerate. With a \(g\)-prior on each submodel and \(g_n\to\infty\) slowly, the method preserves variable-selection consistency under standard design assumptions, provided the Gaussian approximation for \(u\) concentrates around the truth at rate \(o(\log g_n)\). Empirically, the paper reports a nonparametric Gaussian-process example with \(n=100,p=3\) in which full MCMC takes \(\approx 6\) min while IRGA completes in \(2\) sec; a diabetes variable-selection example with \(n=442,r=64\) in which IRGA with \(p=4\) and VAMP for the \(q=60\) nuisance runs in \(\approx 4\) sec versus \(11\) min for Gibbs; and a gene-expression example with \(n=462,p=5,q=10\,000\) in which IRGA with VAMP runs in \(\sim 10\) sec [1909.06753].

A related, but distinct, variational formulation performs mean-field variational inference in a rotated coordinate system chosen by PCA of a cross-covariance matrix involving the target’s score. For an unnormalized density \(p(x)\propto e^{-U(x)}\), with \(\gamma=\mathcal N(0,I_d)\), define
\[
h(x)=\nabla\log p(x)+x,
\]
and
\[
C=\mathbb E_{x\sim\gamma}\bigl[x\,h(x)^T\bigr]
=\mathbb E_{x\sim\gamma}\bigl[x(\nabla\log p(x)+x)^T\bigr].
\]
If \(C=V D V^\top\), then \(R=V^\top\) maximizes
\[
J(R)=\sum_{i=1}^d \bigl[(R C R^\top)_{ii}\bigr]^2.
\]
In the rotated coordinates \(z=R^\top x\), one applies MFVI with \(q(z)=\prod_i q_i(z_i)\), where the coordinatewise optimum satisfies
\[
\nabla_{z_i}\log q_i^*(z_i)
=\mathbb E_{z_{-i}\sim q_{-i}}\bigl[\nabla_{z_i}\log p_R(z)\bigr].
\]
Iterating rotation and coordinatewise Gaussianization produces a compositional transport map with monotonic KL reduction,
\[
KL(\gamma\|p^{(k)})\le KL(\gamma\|p^{(k-1)}),
\]
and, for Gaussian targets, geometric contraction under random rotations [2510.07732].

## 4. Rotation-based iterative Gaussianization and density transport

In signal processing and generative modeling, the rotating Gaussian model is exemplified by **Rotation-based Iterative Gaussianization** (RBIG). Starting from \(x^{(0)}\in\mathbb R^d\) with unknown density \(p_x(x)\), RBIG constructs
\[
f_k(x^{(k-1)})=R_k\circ G_k(x^{(k-1)}),
\]
where \(G_k\) is a marginal Gaussianization acting componentwise and \(R_k\) is an orthonormal rotation. Writing
\[
u^{(k)}=G_k(x^{(k-1)}),\qquad x^{(k)}=R_k u^{(k)},
\]
the full map is
\[
f_K(x)=x^{(K)}=(R_K\circ G_K)\circ\cdots\circ(R_1\circ G_1)(x).
\]
For each coordinate,
\[
G_k(x^{(k-1)})_i=\Phi^{-1}(F_{i,k}(x_i^{(k-1)})),
\]
where \(F_{i,k}\) is the one-dimensional marginal CDF and \(\Phi^{-1}\) is the probit transform [1602.00229].

The rotation stage mixes coordinates whose marginals are already \(N(0,1)\) but remain statistically dependent. The literature considers PCA, ICA, and random orthonormal rotations. Because each layer is differentiable and invertible, density recovery follows by change of variables. The Jacobian of \(G_k\) is diagonal, with
\[
[J_{G_k}(x)]_{ii}=g(z_i)^{-1}p_{i,k}(x_i),
\]
and the original density is recovered from the terminal Gaussian density and the product of Jacobian determinants. Convergence is expressed through negentropy
\[
J(x)=D_{KL}[p_x(x)\|\mathcal N(0,I)].
\]
RBIG satisfies
\[
J(x)-J(R\Psi(x))=J_m(x)\ge 0,\qquad \forall R,
\]
where \(J_m(x)\) is the sum of one-dimensional negentropies, and also
\[
I(\Psi(x))-I(R\Psi(x))=J_m(R\Psi(x))\ge 0,\qquad \forall R,
\]
so total negentropy strictly decreases at every iteration unless the variable is already a spherical Gaussian [1602.00229].

For image-sized data, dense rotations do not scale. **Convolutional RBIG** replaces each dense orthonormal rotation \(R_i\) with a convolutional operator \(C_i\), using
\[
x_i=C_i * \Psi_i(x_{i-1}),
\qquad
\Psi_i(x_{i-1})\approx C_i^T\!\,\overset{\small\smile}{*}\,x_i,
\]
and learns filters by minimizing
\[
L_c(c)=\|x-C^T(Cx)\|_2+\lambda_{\rm act}\|Cx\|_1+\lambda_c\|c\|_1.
\]
This preserves the layerwise structure of Gaussianization while making image Gaussianization tractable at high dimension; it also retains information-theoretic bookkeeping through the drop in total correlation
\[
\Delta I_i=\sum_d h(x^{(i-1)}_d)-\sum_d h(x^{(i)}_d)
\]
when \(C_i\) is approximately orthonormal [2206.03860].

A later convergence analysis studies Gaussianization with random rotations in dimension \(D\). If \(x_\ell=\phi(R_\ell x_{\ell-1})\), then for Gaussian input \(p(x)=\mathcal N(0,\Sigma)\) with \(\mathrm{tr}\,\Sigma=D\), exact representation requires almost surely
\[
L\ge \frac{D+1}{2}=\Omega(D),
\]
and the iterative KL loss satisfies
\[
\mathbb E[\mathrm{KL}_L]\le \Bigl(1-\frac{2}{D+2}\Bigr)^L \mathrm{KL}_0.
\]
Hence reducing \(\mathrm{KL}_0\) to \(\varepsilon\) requires
\[
L=\Theta\bigl(D\ln(\mathrm{KL}_0/\varepsilon)\bigr)=\Omega(D).
\]
This suggests that random-rotation Gaussianization is theoretically simple and invertible, but its layer complexity grows linearly with dimension unless additional structure is exploited [2306.13520].

## 5. Stochastic-process and Gaussian-process formulations

In stochastic-process theory, a rotating Gaussian model may refer to a smooth, rotationally invariant, centered Gaussian process in the plane,
\[
X_t=(\xi^x_t,\xi^y_t),\qquad \xi_t=\xi^x_t+i\xi^y_t=r_t e^{i\phi_t},
\]
with covariance
\[
\langle \xi^i_t\xi^j_{t'}\rangle=\delta_{ij}C_{tt'}.
\]
The central observable is the winding angle \(\phi_t\). The angular-velocity correlator is
\[
{\cal C}_v(t,t')
=\bigl\langle \dot\phi_t\,\dot\phi_{t'}\bigr\rangle
=-\tfrac12\Bigl(\partial_t\partial_{t'}\ln|c_{tt'}|\Bigr)\ln(1-c_{tt'}^2),
\]
where \(c_{tt'}=C_{tt'}/\sqrt{C_{tt}C_{t't'}}\). The winding-angle variance is
\[
\Phi_{tt'}=\bigl\langle(\phi_t-\phi_{t'})^2\bigr\rangle
=\int_{t'}^t ds_1\int_{t'}^t ds_2\,{\cal C}_v(s_1,s_2).
\]
For most stationary processes \(C_{tt'}=C(t-t')\), the large-time behavior is diffusive,
\[
\Phi(\tau)\simeq 2D\tau,\qquad
D=\int_0^\infty ds\,\frac{[C'(s)]^2}{C(0)^2-C(s)^2}.
\]
For smooth processes with stationary increments, the variance grows as \(\tfrac12(\ln t)^2\), and the paper also analyzes fractional Brownian motion, correlators \(\langle e^{in(\phi_t-\phi_{t'})}\rangle\), the distribution of \(\dot\phi_t\), and the variance of algebraic area [0904.0582].

A different Gaussian-process construction arises in stellar variability modeling. The **starry-process** models the light curve of a rotating, evolving stellar surface by marginalizing over spot configurations. In the spherical-harmonic representation,
\[
f=1+A(I,P,u)\,y,
\]
where \(A\) is the design matrix depending on inclination \(I\), rotation period \(P\), limb-darkening coefficients \(u\), and observation times. After marginalizing over spot realizations parameterized by \(\theta_\bullet=(n,c,\mu_\phi,\sigma_\phi,r,\Delta r)\), the light curve is approximated by
\[
f\mid P,u,\theta_\bullet\sim \mathcal N(\mu,\Sigma),
\]
with
\[
\mu=1+A(I,P,u)\mu_y,\qquad
\Sigma=A(I,P,u)\Sigma_y A^\top(I,P,u).
\]
The moments \(\mu_y\) and \(\Sigma_y\) are obtained from nested closed-form integrals over radius, latitude, longitude, and contrast, using Wigner-rotation matrices \(R_x\) and \(R_y\). Optional inclination marginalization and a normalization correction are also derived in closed form. The implementation reports \(\sim 20\) ms for GP mean, covariance, and log-likelihood at \(K\sim 10^3\) and \(l_{\max}=15\) on a modern laptop [2102.01697].

These two stochastic-process uses share Gaussian-process structure but differ sharply in emphasis. In the planar winding problem, rotational invariance is a property of the **state law**; in the stellar model, rotation is part of the **forward map** from a random surface to a light curve.

## 6. Quantum and few-body rotating Gaussian states

In quantum mechanics, the rotating Gaussian model appears as a two-dimensional Gaussian packet with fixed mean angular momentum,
\[
\psi(x,y)=N\exp[-\mu(a x^2+b x y+c y^2)+F x+G y],
\]
where \(a,b,c,F,G\) are complex parameters subject to positivity conditions on the real parts. The expectation of angular momentum decomposes into
\[
\langle L_z\rangle=\hbar\,\mathcal L=\hbar(\mathcal L_c+\mathcal L_i),
\]
with \(\mathcal L_c\) the external contribution from motion of the packet center and \(\mathcal L_i\) the intrinsic contribution due to quantum fluctuations. For the isotropic oscillator, minimizing the mean energy at fixed \(\mathcal L_c\) and \(\mathcal L_i\) yields
\[
E^{\min}=\hbar\omega[1+|\mathcal L_c|+|\mathcal L_i|].
\]
In the co-rotating case, where \(\mathcal L=\mathcal L_c+\mathcal L_i\), this becomes
\[
E^{\min}=\hbar\omega(1+|\mathcal L|).
\]
The minimizing packets have nonzero coordinate-momentum correlations and moderate quadrature squeezing, with
\[
S_x=S_y=\frac{1}{1+\eta}<1.
\]
The same framework is developed for free evolution and for a charged particle in a homogeneous magnetic field, where co-rotating packets in the same sense as the Larmor precession are static in time [1509.07992].

A semiclassical extension studies rotating Gaussian wave packets in weak external potentials in two and three dimensions. A minimal packet is parameterized by its center \(\mathbf q\), mean velocity \(\mathbf v\), and a complex symmetric matrix
\[
\Omega=\Omega_{\rm Re}+i\Omega_{\rm Im},\qquad \Omega_{\rm Im}>0.
\]
The internal angular momentum is
\[
\mathcal L_{\rm i}
=\frac{\hbar}{2}\,[\Omega_{\rm Re},(\Omega_{\rm Im})^{-1}]^{\rm ax},
\]
so it is nonzero only when \([\Omega_{\rm Re},\Omega_{\rm Im}]\neq 0\). Using an eikonal approximation, the paper derives an explicit first-order change in internal angular momentum induced by a weak potential. For a two-dimensional particle crossing a tilted ridge potential,
\[
[\Delta\mathcal L_{\rm i}]_3
=\hbar\,\frac{qV_0}{\mu v^3}(\omega_1-\omega_2),
\]
showing that anisotropy of the initial packet fixes the sense of rotation [1703.06413].

In few-body physics, the Gaussian enters through the interaction rather than the wave packet. Two identical spinless bosons in a rotating 2D harmonic trap interact through the finite-range Gaussian potential
\[
V_\sigma(r)=\frac{1}{2\pi\sigma^2}\exp\!\Bigl(-\frac{r^2}{2\sigma^2}\Bigr),
\]
which tends to a contact \(\delta\)-function as \(\sigma\to 0\). After separation into center-of-mass and relative motion, the relative Hamiltonian yields a transcendental equation for each angular-momentum sector \(m\). The paper reports that, for a given relative angular momentum and interaction strength \(g_2>0\), the ground-state energy increases with interaction range; below \(g_2V(r)\le -1\), the \(m=0\) ground-state energy diverges to physically unacceptable negative infinity; and for \(|m|=1\), the ground-state energy becomes independent of the interaction strength. In the \(\delta\)-limit, convergence of the ground-state energy requires a considerably large critical Hilbert space, whereas for Gaussian interaction potential with \(\sigma\to 1\) convergence occurs for a considerably small critical Hilbert space [2205.11958].

## 7. Geometric and computer-vision formulations

In rotated-object detection, a rotated rectangle with center \((x,y)\), width \(w\), height \(h\), and angle \(\theta\) is modeled as a Gaussian distribution
\[
p(\mathbf u)=\mathcal N(\mathbf u;\mu,\Sigma),\qquad \mu=(x,y)^\top,
\]
with
\[
\Sigma^{1/2}=R(\theta)\Lambda,\qquad
\Lambda=\mathrm{diag}(w/2,h/2),
\]
equivalently
\[
\Sigma=R(\theta)\Lambda^2R(\theta)^\top.
\]
This representation removes several pathologies of angle-parameterized box regression. Because \(\Sigma(\theta)=\Sigma(\theta-\pi)\), the model is blind to \(180^\circ\) flips; swapping \(w\leftrightarrow h\) and \(\theta\to\theta-\pi/2\) leaves \(\Sigma\) unchanged; and when \(w\approx h\), \(\Sigma\) becomes nearly isotropic and almost independent of \(\theta\). Regression is performed with Gaussian-to-Gaussian distances, especially KLD,
\[
D_{KL}(p_1\|p_2)
=\tfrac12\bigl[(\mu_2-\mu_1)^\top\Sigma_2^{-1}(\mu_2-\mu_1)
+\mathrm{tr}(\Sigma_2^{-1}\Sigma_1)
-\ln\det(\Sigma_2^{-1}\Sigma_1)-2\bigr].
\]
The paper further proposes Gaussian-metric label assignment with
\[
G_{ij}=1/(\tau+D_{KL}(p_j\|p_i)),\qquad t_i=m_i+v_i,
\]
using ATSS-style thresholding. It extends the construction to 3D by embedding a box \((x,y,z,w,h,l,\theta)\) in a \(3\times 3\) Gaussian with
\[
\Sigma_3^{1/2}=R_y(\theta)\,\mathrm{diag}(w/2,h/2,l/2).
\]
Reported gains include \(+23.97\) AP\(_{75}\) points on HRSC2016, a rise from \(65.73\%\) to \(71.28\%\) mAP\(_{50}\) on DOTA-v1.0 with RetinaNet, \(+1.9\) 3D mAP on KITTI val for PointPillars+KLD, and \(+1.8\) points in overall mAP on Waymo L2 [2209.10839].

In dynamic novel-view rendering, each static Gaussian primitive has mean \(\mu\in\mathbb R^3\), covariance \(\Sigma\in\mathbb R^{3\times 3}\), color \(c\in\mathbb R^3\), and opacity \(\alpha\), with
\[
\Sigma=R S R^\top,
\]
where \(S=\mathrm{diag}(s_1^2,s_2^2,s_3^2)\) and \(R\in SO(3)\). In the dynamic setting,
\[
G_i(t)=[\mu_i(t),\Sigma_i(t),\alpha_i(t)],\qquad
\Sigma_i(t)=R_i(t)S_i(t)R_i(t)^\top.
\]
The temporal evolution of position, quaternion, and scale is modeled by infinite-order Taylor expansions around a reference time \(t_r\), truncated in practice to third order and supplemented by a learnable Peano remainder,
\[
T_i(t)=f^{\rm polynomial}_i(t)+H_i(t).
\]
The paper uses a sparse set of “Global Primitives” and linear blend skinning for “Local Primitives”, with all Taylor coefficients and remainder weights trained end-to-end by photometric reconstruction loss. An ablation removing the time-rotation module drops PSNR on the Cut Roast Beef scene from \(36.96\) dB to \(32.94\) dB and SSIM from \(0.984\) to \(0.972\), quantifying the effect of explicit rotation modeling [2412.04282].

A common source of confusion is terminological rather than mathematical. In object detection and dynamic view rendering, rotation parameterizes a covariance or quaternion field attached to a Gaussian primitive. In Gaussianization and rotated posterior approximation, rotation is a coordinate transform used to simplify dependence structure. In rotating vortices and quantum packets, it is part of the physical dynamics. This suggests that “Rotating Gaussian Model” is best understood not as one model class, but as a family of Gaussian-based constructions linked by the repeated appearance of rotation as either symmetry, transport, factorization, or orientation.

Source: https://www.emergentmind.com/topics/rotating-gaussian-model