---
title: 'COS Procedure: Convex Oscillation Suppression'
url: https://www.emergentmind.com/topics/convex-oscillation-suppressing-cos-procedure
type: topic
---

# COS Procedure: Convex Oscillation Suppression

Searching arXiv for the cited COS-related papers and terminology.
arxiv_search(query="Convex-Oscillation-Suppressing OR COS Procedure oscillation suppression convex", max_results=10)
arxiv_search(query="1707.05342 OR 2604.00301 OR 2409.09620 OR 1810.10439 OR 2512.02745", max_results=10)
Convex-Oscillation-Suppressing (COS) Procedure denotes, across several research literatures, a class of constructions that combine convexity or convexification with mechanisms intended to suppress undesirable oscillatory behavior. In statistical learning, the label has been used retrospectively for an unrestricted procedure that suppresses empirical-process oscillations through robust blockwise comparisons and exploits convex averaging to recover convex-case rates even for non-convex classes. In high-order finite volume and discontinuous Galerkin discretizations, it appears as an explicit oscillation-suppressing convex constraint enforced by radial scaling. In non-convex optimization, it describes sequential convex programming based on inner-convex majorizers that prevent feasibility oscillations and cost increases. In computational probability, however, COS refers to the Fourier–cosine expansion method, where convergence is governed by tail cosine energy rather than convex admissibility. This suggests that COS is best understood as a cross-domain technical motif rather than a single universally standardized algorithm [1707.05342] [2604.00301] [1810.10439] [2512.02745].

## 1. Terminological scope and principal usages

The expression “COS” does not have a single fixed meaning across arXiv-indexed research. In some works it is an explicit part of the formal framework; in others it is a reconstruction or interpretive lens applied to a method that is not originally named that way. A common source of confusion is therefore terminological rather than mathematical.

| Domain | COS meaning | Representative source |
|---|---|---|
| Statistical learning | A reconstructed “Convex-Oscillation-Suppressing” view of an unrestricted learner | [1707.05342] |
| FV/DG stabilization | An explicit oscillation-suppressing convex constraint inside EPO | [2604.00301] |
| DG on triangular meshes | A reconstructed COS-like combination of RIOE damping and BP limiting | [2409.09620] |
| Non-convex optimization | A tailored COS interpretation of sequential convex programming with inner-convex surrogates | [1810.10439] |
| Computational probability | The Fourier–cosine expansion method | [2512.02745] |

In the learning paper, the procedure is “not explicitly named COS in the paper,” but is described through that lens because it “suppresses empirical-process oscillations” and “exploits convexity” via essential subsets, star-hulls, and convex averaging [1707.05342]. The triangular-mesh DG paper similarly “does not use the term ‘COS’,” and reconstructs a COS-like procedure by combining oscillation-eliminating methodology with optimal convex decomposition [2409.09620]. By contrast, the EPO framework explicitly defines an “oscillation-suppressing (COS) constraint” as one of three convex constraints enforced on a cellwise candidate state [2604.00301]. The computational-probability literature uses COS in a different sense altogether: a Fourier–cosine expansion method for densities and transforms [2512.02745].

## 2. Statistical learning: unrestricted learning through convex localization

For learning under squared loss, the COS interpretation in "An optimal unrestricted learning procedure" considers an arbitrary class \(F\) of square-integrable functions on \((\Omega,\mu)\), a distribution \(X \sim \mu\), and an unknown square-integrable target \(Y\). The loss is
\[
L(f(X),Y)=(f(X)-Y)^2,
\]
with risk
\[
R(f)=\mathbb{E}(f(X)-Y)^2,
\]
and oracle
\[
f^*=\arg\min_{f\in F}R(f).
\]
The analysis introduces \(\sigma^2=\mathbb{E}(f^*(X)-Y)^2\), a uniform integrability assumption on \(U-U\) with \(U=(F+F)/2\), and the target-class interaction parameter
\[
L_T=\sup_{f\in F}\left(\mathbb{E}\left[\left(\frac{(f-f^*)(X)}{\|f-f^*\|_{L_2}}\right)^2\left(\frac{f^*(X)-Y}{\sigma}\right)^2\right]\right)^{1/2},
\]
which ensures
\[
\mathbb{E}\big[(f-f^*)^2(X)(f^*(X)-Y)^2\big]\le L_T^2\sigma^2\|f-f^*\|_{L_2}^2.
\]
Localization is carried out through the star-shaped set
\[
F_{h,r}=\mathrm{star}(F-h,0)\cap rD,
\]
which enables convex localization even when \(F\) is non-convex [1707.05342].

The procedure \(P\) has two components, \(P_1\) and \(P_2\), and is run twice. \(P_1\) is a distance oracle: for \(h,f\in H\), it computes \(P_1(h,f)=v_\ell^*\), the \(\ell\)-th largest entry of the nonincreasing rearrangement of \((|f-h|(X_i))_{i=1}^N\). \(P_2\) is a tournament on a second independent sample partitioned into \(n\) blocks of size \(m=N/n\). For each block \(j\),
\[
Q_{h,f}(j)=\frac{1}{m}\sum_{i\in I_j}(h-f)^2(X_i),\qquad
M_{h,f}(j)=\frac{2}{m}\sum_{i\in I_j}(h-f)(X_i)(f(X_i)-Y_i),
\]
and
\[
B_{h,f}(j)=Q_{h,f}(j)+M_{h,f}(j).
\]
The “home-and-away” decision rule declares \(f\succ h\) if, on more than \(n/2\) blocks, either \(B_{h,f}(j)\ge -\theta_2 r^2\) when \(P_1(h,f)\le \theta_4 r\), or \(B_{h,f}(j)\ge -\theta_3 P_1(h,f)^2\) when \(P_1(h,f)>\theta_4 r\). The winners form
\[
P_2(H)=\{f\in H:\text{for all }h\in H,\ f\succ h\}.
\]
The algorithm first applies \(P_1,P_2\) on \(F\) to obtain \(F_1\), then forms \(\bar F_1=(F_1+F_1)/2\), repeats \(P_1,P_2\) on \(\bar F_1\), obtains \(F_2\), and outputs any \(\tilde f\in F_2\) [1707.05342].

Its oscillation-suppressing content is encoded by localized intrinsic and extrinsic complexity measures,
\[
N_{\mathrm{int}}(T,r,\kappa),\qquad N_{\mathrm{ext}}(T,r,\kappa),
\]
defined through localized Rademacher and multiplier averages over \(F_{f^*,r}\). The fixed-point prescription is to choose \(r^2\asymp \varepsilon\) so that these localized oscillations are at most proportional to \(r\) and \(r^2\), respectively. On high-probability events, \(P_2(H)\) becomes an \((\rho,r')\)-essential subset with \(\rho\le 1/20\) and \(r'\asymp r\), and the two-stage convexification yields
\[
\|\tilde f(X)-Y\|_{L_2}^2\le \inf_{f\in F}\|f(X)-Y\|_{L_2}^2+\frac{3r^2}{2}.
\]
The main theorem states that if
\[
N\ge 2(N_1+N_2),
\]
with \(N_1,N_2\) defined from \(N_{\rm int}\), \(N_{\rm ext}\), and the term
\[
c_2\Big(\frac{L_T^2\sigma^2}{\varepsilon}+1\Big)\log\Big(\frac{64}{\delta}\Big),
\]
then with probability at least \(1-\delta\),
\[
R(\tilde f)-\inf_{f\in F}R(f)\le \frac{3}{2}\varepsilon.
\]
The sample complexity matches, up to constants, the lower bounds required to overcome intrinsic and extrinsic obstructions. If \(F\) is convex, then \((F+F)/2=F\), so the procedure is proper; if \(F\) is non-convex, the output lies in a convex superset and the procedure is unrestricted. The heavy-tail robustness comes from order-statistic distance estimation in \(P_1\) and majority-of-block control in \(P_2\) [1707.05342].

## 3. Finite volume and discontinuous Galerkin stabilization: the EPO ray and COS radius

In "EPO: A Unified Framework for Entropy Stability, Positivity, and Oscillation Suppression," COS is an explicit local constraint acting on a candidate finite volume or discontinuous Galerkin update for a hyperbolic system \(\partial_t U+\partial_xF(U)=0\). On each cell \(I_i\), the candidate nodal array \(U_i^{\text{cand}}\) has conservative average
\[
\bar U_i^{\text{cand}}=\sum_k\omega_k U_{i,k}^{\text{cand}},
\]
and all corrections are restricted to the mean-preserving ray
\[
\mathcal{S}_{\bar U_i^{\text{cand}}}(\theta;U_i^{\text{cand}})
=\bar U_i^{\text{cand}}+\theta\big(U_i^{\text{cand}}-\bar U_i^{\text{cand}}\big),\qquad 0\le \theta\le 1.
\]
At quadrature nodes,
\[
U_{i,k}(\theta)=\bar U_i^{\text{cand}}+\theta\big(U_{i,k}^{\text{cand}}-\bar U_i^{\text{cand}}\big).
\]
The final limited state is
\[
U_i^{\text{lim}}=\bar U_i^{\text{cand}}+\theta_i\big(U_i^{\text{cand}}-\bar U_i^{\text{cand}}\big),
\]
where \(\theta_i\) is determined by convex constraints for admissible state, entropy, and oscillation suppression [2604.00301].

The admissible-state radius \(\theta_i^{\text{pos}}\) is the largest \(\theta\in[0,1]\) such that every node belongs to a closed convex invariant set \(G\). For the compressible Euler equations, a common choice is
\[
G_\varepsilon=\Bigl\{(\rho,m,E)^\top:\ \rho\ge\varepsilon,\quad
p(U)=(\gamma-1)\Big(E-\frac{m^2}{2\rho}\Big)\ge\varepsilon\Bigr\}.
\]
The entropy radius \(\theta_i^\eta\) is defined by a convex nodal inequality
\[
\sum_{k=1}^L\omega_k\,\eta(U_{i,k}(\theta))\le B_i,
\]
where \(B_i\) is a local budget coming from weak cell-average entropy stability. Along the ray,
\[
\Psi_i(\theta)=\sum_{k=1}^L\omega_k\,\eta(U_{i,k}(\theta))
\]
is convex and nondecreasing, so
\[
\theta_i^\eta=\sup\{\theta\in[0,1]:\Psi_i(\theta)\le B_i\}.
\]
For quadratic entropy \(\eta(U)=|U|^2/2\), the paper gives a closed form for \(\theta_i^\eta\). The construction also extends to any prescribed finite family of convex entropy pairs by computing one radius per pair and taking the minimum [2604.00301].

The COS constraint itself is posed as a bound-preserving convex local bound at quadrature or nodal points. For each component \(\ell\) and node \(k\),
\[
B_{i,k,\ell}^{\min}\le U_{i,k,\ell}\le B_{i,k,\ell}^{\max},
\]
with practical bounds derived from neighbor averages, a monotone low-order update, or two-point Riemann/Lax–Friedrichs averages. Writing
\[
U_{i,k,\ell}(\theta)=\bar U_{i,\ell}^{\text{cand}}+\theta\,\delta_{i,k,\ell}^{\text{cand}},
\]
the componentwise restriction yields closed-form upper bounds \(T_{i,k,\ell}\), and the COS radius is
\[
\theta_i^{\text{COS}}=\min\Bigl\{1,\ \min_{k,\ell}T_{i,k,\ell}\Bigr\}.
\]
The final limiter is then
\[
\theta_i=\min\Big\{\theta_i^{\text{pos}},\ \min_m\theta_i^{\eta_m},\ \theta_i^{\text{COS}}\Big\},
\]
or, in positivity-first form,
\[
\theta_i=\min\Big\{\theta_i^{\text{PE,all}},\ \theta_i^{\text{COS}}\Big\}.
\]
Because the admissible-state set, entropy sublevel set, and COS set are closed convex sets along the ray, the limiter is the maximal radial projection into their intersection. The framework proves preservation of cell averages, invariant-set preservation, local and global strong entropy inequalities, stagewise budgets for SSP Runge–Kutta methods, an SSP multistep variant, and extensions on rectangular and unstructured triangular meshes [2604.00301].

## 4. Unstructured triangular DG: RIOE damping and optimal convex decomposition

On unstructured triangular meshes, the COS-like construction in "Robust DG Schemes on Unstructured Triangular Meshes: Oscillation Elimination and Bound Preservation via Optimal Convex Decomposition" combines oscillation-eliminating damping with convex bound preservation. The scalar or system conservation law is discretized in the DG space
\[
V_h^k=\{v_h\in L^2(\Omega):v_h|_K\in P^k(K)\ \forall K\in\mathcal{T}_h\},
\]
with cell average
\[
\bar u_K(t)=\frac{1}{|K|}\int_Ku_h(x,y,t)\,dx\,dy=u_K^{(0)}(t).
\]
The semi-discrete method uses a Lax–Friedrichs flux
\[
\widehat F(u^-,u^+,n)=\frac12\big(F(u^-)\cdot n+F(u^+)\cdot n-\alpha^{LF}(u^+-u^-)\big),
\]
and SSPRK time stepping [2409.09620].

The oscillation-eliminating operator \(\mathfrak{F}_\tau\) is defined by a pseudo-time damping equation whose exact modal solution preserves the average and exponentially damps higher modes:
\[
u_\sigma^{(0)}(\Delta t)=u_K^{(0)}(t^n),\qquad
u_\sigma^{(\ell)}(\Delta t)=u_K^{(\ell)}(t^n)\exp\Big(-\Delta t\sum_{j=0}^m\sigma_K^j(u_h^n)\Big),
\]
for \(L_{m-1}+1\le \ell\le L_m\). The coefficients
\[
\sigma_K^j(u_h)=\sum_{i=1}^3\beta_K^{(i)}h_K^{(i)}\delta_K^{(i),j}(u_h)
\]
depend only on the cell and its edge neighbors, and the paper establishes conservation, compactness, scale invariance, and evolution invariance. For systems such as Euler, direct component-wise OE is not rotation-invariant. The rotation-invariant OE (RIOE) procedure therefore redefines momentum damping via normal and tangential edge components,
\[
m_{n,h}^{(i)}=m_h\cdot n_K^{(i)},\qquad
m_{t,h}^{(i)}=m_h\cdot n_{t,K}^{(i)},
\]
and replaces momentum entries by
\[
\hat\delta_K^{(i),j}(m_h)=\max\{\hat\delta_K^{(i),j}(m_{n,h}^{(i)}),\hat\delta_K^{(i),j}(m_{t,h}^{(i)})\}.
\]
With this change, \(\mathfrak{F}_\tau^{RI}\) is rotation-invariant whenever the PDE and numerical flux are rotation-invariant [2409.09620].

Bound preservation is obtained from feasible convex decompositions of cell averages into edge integrals and internal nodes. For \(P^1\) and \(P^2\) triangles, the paper constructs optimal decompositions and the associated maximal BP CFL numbers.

| Degree | Optimal BP CFL number | Reported gain over classic decomposition |
|---|---|---|
| \(P^1\) | \(C_{K,1}^{DCW}=2/[3(l_K^{(1)}+l_K^{(2)})]\) | \(100\%\)–\(200\%\) |
| \(P^2\) | \(C_{K,2}^{DCW}=2/(9\bar l_K+3\hat l_K)\) | \(280.38\%\)–\(350\%\) |

For \(P^1\), the optimal feasible decomposition uses
\[
w_i=\frac{2l_K^{(i)}}{3l_K^{(1)}+3l_K^{(2)}},\qquad
\omega_1=\frac{l_K^{(1)}+l_K^{(2)}-2l_K^{(3)}}{3l_K^{(1)}+3l_K^{(2)}},
\]
with an explicit internal node \((\xi_1,\eta_1)\). For \(P^2\), it uses
\[
w_i=\frac{2l_K^{(i)}}{9\bar l_K+3\hat l_K},
\]
and two internal nodes with equal weights
\[
\omega_s=\frac{\bar l_K+\hat l_K}{6\bar l_K+2\hat l_K}.
\]
A simplified two-step limiter then rescales the polynomial about the cell average, first to enforce density positivity and then to enforce positivity of specific internal energy [2409.09620].

The resulting stagewise scheme is
\[
\text{DG update} \rightarrow \text{RIOE damping} \rightarrow \text{BP limiter}.
\]
All coefficients use only edge-neighboring cell information, so compactness and parallel efficiency are retained. Numerical evidence includes optimal \((k+1)\)-th order accuracy for smooth scalar tests, non-oscillatory Burgers solutions, rotation-invariant Euler implosion results, and substantial efficiency gains in forward-facing step, double Mach reflection, and shock diffraction. Average time-step ratios satisfy \(\overline{\Delta t}^{DCW}/\overline{\Delta t}^{ZXS}\approx 2.85\) for \(P^1\) and \(\approx 4.49\) for \(P^2\), with corresponding CPU-time reductions such as \(95.46\) h versus \(250.07\) h for \(P^1\) and \(319.97\) h versus \(1274.51\) h for \(P^2\) in the forward-facing-step case [2409.09620].

## 5. Non-convex optimization: inner-convex SCP as an oscillation-suppressing convexification

In "A recursively feasible and convergent Sequential Convex Programming procedure to solve non-convex problems with linear equality constraints," the COS interpretation refers to a sequential convex programming procedure for
\[
\min_x f_0(x)\quad\text{s.t.}\quad f_i(x)\le 0,\ i=1,\dots,p,\qquad Ax=b,
\]
where \(f_0\) and \(f_i\) are non-convex and continuously differentiable, while the equalities are affine. At iteration \(k\), the method builds convex majorizers
\[
\tilde f_0(x;x_k),\qquad \tilde f_i(x;x_k),
\]
satisfying
\[
\tilde f(x;x_k)\ge f(x)\ \text{for all }x,\qquad
\tilde f(x_k;x_k)=f(x_k),\qquad
\nabla\tilde f(x_k;x_k)=\nabla f(x_k),
\]
and solves the convex subproblem
\[
\min_x \tilde f_0(x;x_k)\quad\text{s.t.}\quad Ax=b,\qquad \tilde f_i(x;x_k)\le 0.
\]
These are inner-convex approximations: feasibility for the convexified constraints implies feasibility for the original problem, and minimizing the surrogate guarantees monotone descent of the true objective [1810.10439].

For difference-of-convex functions \(f=c_1-c_2\), the surrogate reduces to linearization of the concave part,
\[
\tilde f(x;x_k)=c_1(x)-c_2(x_k)-\nabla c_2(x_k)^T(x-x_k).
\]
More generally, the paper develops Taylor-based inner-convex approximations. If
\[
H_f(x_k)=W\Lambda W^T,
\]
the quadratic term keeps only the positive semidefinite part,
\[
H_f^+(x_k)=W\Lambda^+W^T,
\]
and higher-order Taylor terms are convexly overestimated through diagonal terms and absolute coefficient sums. The resulting surrogate is
\[
\tilde f(x;x_k)=f(x_k)+\nabla f(x_k)^T(x-x_k)+\frac12(x-x_k)^TH_f^+(x_k)(x-x_k)+\text{convex higher-order terms}.
\]
If the Taylor series is truncated at order \(|\alpha|\), the method adds
\[
\frac{M}{(|\alpha|+1)!}\|x-x_k\|^{|\alpha|+1}
\]
to recover a global upper bound; if majorization fails at the tentative solution, \(M\) is increased and the convex subproblem is re-solved [1810.10439].

The guarantees are recursive feasibility, monotone descent, and convergence to a KKT point, provided the first convex subproblem is feasible. Under local strong convexity of the Lagrangian Hessian at the limit point and local Lipschitz continuity of the Hessian, the procedure has a quadratic local rate:
\[
\|x_{k+1}-x_k\|\le \gamma_k\|x_k-x_{k-1}\|^2.
\]
The “oscillation-suppressing” role of the convexification is therefore not nodal or modal but iterative: the method eliminates cost increases and constraint-violation oscillations that are typical of pure successive linearization. The paper also reports real-time performance in aerial-vehicle trajectory optimization. With \(N=25\) nodes and non-convex thrust and keep-out-zone constraints, the implementation was approximately \(100\times\) faster than GPOPS-II on a laptop; a penalty phase found an admissible solution in approximately \(98.1\%\) of cases, and once feasibility was obtained the SCP achieved \(100\%\) convergence [1810.10439].

## 6. Fourier–cosine COS method: tail-energy admissibility and convergence

In computational probability and option pricing, COS refers to the Fourier–cosine expansion method. Here the target density \(f\) is approximated on a truncation interval \([a,b]\) by
\[
f(x)\approx {\sum\nolimits\!'}_{k=0}^\infty A_k\cos\!\left(k\pi\frac{x-a}{b-a}\right),
\]
with coefficients
\[
A_k=\frac{2}{b-a}\int_a^b f(x)\cos\!\left(k\pi\frac{x-a}{b-a}\right)\,dx.
\]
If
\[
\phi(\omega)=\int_{\mathbb{R}}e^{i\omega x}f(x)\,dx
\]
is the characteristic function, practical coefficients are computed as
\[
F_k=\frac{2}{b-a}\,\Re\!\left[\phi\!\left(\frac{k\pi}{b-a}\right)\exp\!\left(-i\frac{ka\pi}{b-a}\right)\right],
\]
and the \(N\)-term approximant is
\[
f_1^N(x)={\sum\nolimits\!'}_{k=0}^{N-1}F_k\cos\!\left(k\pi\frac{x-a}{b-a}\right).
\]
The 2025 convergence note studies the symmetric case \([-L,L]\) and defines the tail cosine energy
\[
B(L)=\sum_{k=0}^{\infty}\frac{1}{L}\Bigg|\int_{\mathbb{R}\setminus[-L,L]}f(x)\cos\!\Big(\frac{k\pi(x+L)}{2L}\Big)\,dx\Bigg|^2.
\]
A density \(f\in L^1\) is COS-admissible if \(B(L)\to 0\) as \(L\to\infty\) [2512.02745].

The main sufficient condition is moment-based. In one dimension, if
\[
f\in L^1(\mathbb{R})\cap L^2(\mathbb{R})
\]
and
\[
\int_{\mathbb{R}}|x|^p|f(x)|^2\,dx<\infty
\]
for some \(p>1\), then \(f\) is COS-admissible and
\[
B(L)\le 2\zeta(p)\Bigg(L^{-p}\!\!\int_{|x|>L}|x|^p|f(x)|^2\,dx+\int_{|x|>L}|f(x)|^2\,dx\Bigg),
\]
hence
\[
B(L)\le 4\zeta(p)L^{-p}\!\!\int_{\mathbb{R}}|x|^p|f(x)|^2\,dx,
\qquad B(L)=O(L^{-p}).
\]
In \(d\) dimensions, if
\[
\int_{\mathbb{R}^d}|x|^p|f(x)|^2\,dx<\infty
\]
for some \(p>d\), then the \(d\)-dimensional tail energy \(B_d(L)\) satisfies
\[
B_d(L)\le C_{d,p}L^{-p}\!\!\int_{\mathbb{R}^d\setminus[-L,L]^d}|x|^p|f(x)|^2\,dx,
\]
where
\[
C_{d,p}=2^{d-1}(1+d^{p/2})\sum_{m\in\mathbb{Z}^d\setminus\{0\}}|m|^{-p}.
\]
The paper’s main convergence theorem is
\[
\lim_{L\to\infty}\limsup_{N\to\infty}\|f-f_1^N\|_2=0
\]
for COS-admissible \(f\in L^1\cap L^2\) [2512.02745].

This use of COS differs conceptually from the convex-admissibility frameworks in learning, PDE stabilization, and SCP. The relevant “oscillation suppression” here is the control of truncation-induced oscillatory tail contributions through weighted \(L^2\) moment bounds and blockwise orthogonality arguments. The note enlarges previous admissibility conditions by weakening the one-dimensional moment requirement from \(p=2\) to any \(p>1\), and from the multidimensional requirement to any \(p>d\). As a result, it includes heavy-tailed distributions such as Student \(t\) with small degrees of freedom [2512.02745].

## 7. Cross-cutting principles, distinctions, and common misconceptions

Across these literatures, the recurrent technical pattern is the use of convex structure to regularize unstable behavior. In unrestricted learning, convexity enters through star-hulls, essential subsets, and midpoint averaging, while oscillation suppression is achieved by majority-of-block control of quadratic and multiplier empirical processes [1707.05342]. In EPO, convexity is literal geometry on a mean-preserving ray, and oscillation suppression is a nodal bound-preserving constraint intersected with positivity and entropy sublevel sets [2604.00301]. In triangular-mesh DG, oscillation elimination is modal damping, while convexity appears in the feasible decomposition of cell averages and in the convex rescaling limiter [2409.09620]. In SCP, convexity arises from inner-convex surrogates, and the suppressed oscillations are iterative cost increases and feasibility losses [1810.10439].

A central misconception is to treat COS as a single canonical procedure. The sources do not support that interpretation. Several papers explicitly state that the relevant method is being presented “through the lens” of COS or reconstructed as “COS-like,” whereas the EPO framework contains an explicit COS constraint, and the Fourier–cosine literature uses COS as an acronym for a different numerical expansion method. This suggests that “COS procedure” is an umbrella label covering at least four distinct mathematical mechanisms: robust aggregation, radial projection into convex admissible sets, convex decomposition with rescaling, and inner-convex majorization; the Fourier–cosine COS method forms a separate acronymic lineage rather than a variant of those convex-geometric procedures [1707.05342] [2604.00301] [1810.10439] [2512.02745].

A second distinction concerns what is being suppressed. In learning, the target is empirical-process oscillation under heavy tails. In hyperbolic discretization, it is spurious oscillation at nodal or modal level, often tied to positivity and entropy stability. In optimization, it is oscillatory or non-monotone iterate behavior. In Fourier–cosine approximation, it is truncation error encoded by tail cosine energy. The shared terminology is therefore organizational rather than ontological: the underlying objects, guarantees, and proof techniques differ substantially even when convexity and oscillation suppression are both present.

Source: https://www.emergentmind.com/topics/convex-oscillation-suppressing-cos-procedure