---
title: Dual Tamed Scheme for Index-1 SDAEs
url: https://www.emergentmind.com/topics/dual-tamed-scheme
type: topic
---

# Dual Tamed Scheme for Index-1 SDAEs

A dual tamed scheme is an equivalent reformulation of a semi-implicit tamed discretization for index-1 stochastic differential algebraic equations (SDAEs) under non-global Lipschitz coefficients. In the usage documented in the recent SDAE literature, the term does not denote a second independent numerical algorithm; rather, it denotes the same semi-implicit tamed method rewritten in reduced variables so that the analysis aligns with the inherent stochastic differential equation obtained after eliminating the algebraic constraints. The central setting is the index-1 SDAE with singular matrix structure, superlinear drift growth, and strong convergence order \(1/2\) for the resulting approximation [2509.09032].

## 1. Terminology and scope

The phrase “dual tamed scheme” is used explicitly for SDAEs with singular matrix \(A(t)\), where the original numerical method is first posed in the full variable \(X_n^N\) and then rewritten in reduced coordinates \((u_n,\hat v^N)\). In that setting, “dual” means that the numerical scheme is expressed in the coordinates induced by projector-based elimination of the algebraic constraints, so that it can be compared directly with the reduced exact dynamics \(u(t)\) [2509.09032].

This terminology is narrower than the broader literature on tamed discretizations. In related work on McKean–Vlasov equations driven by Lévy noise, the core method is a tamed Euler scheme for an interacting particle system, and the paper states explicitly that it does **not** study a “Dual Tamed Scheme” in the sense of a separate duality-based numerical method; the approximation is instead two-stage, consisting of particle approximation and time discretization [2010.08585]. Likewise, the “semi-tamed Euler scheme” for jump–diffusion SDEs and the “tamed \(\theta\)-Euler-Maruyama scheme” for neutral SDDEs are distinct constructions, not dual schemes in this specific SDAE sense [1510.06067], [1612.02800].

A plausible implication is that the term should be reserved for contexts where a constrained stochastic system is first reduced to an inherent SDE and the numerical method is then rewritten in those reduced variables. Outside that setting, the literature more commonly distinguishes between tamed, semi-tamed, modified tamed, and \(\theta\)-tamed schemes rather than dual schemes.

## 2. SDAE setting and reduction to an inherent SDE

The underlying model is the SDAE
\[
A(t)dX(t)=\left[B(t)X(t)+f(t,X(t))\right]dt + g(t,X(t))dW(t), \quad t\in \left[0,T\right] ,\quad X(0)=X_0,
\]
where \(A(t)\in\mathbb{R}^{d\times d}\) is singular, \(B(t)\in\mathbb{R}^{d\times d}\), \(f\) is the drift, \(g\) the diffusion, and \(W(t)\) is an \(m\)-dimensional Wiener process [2509.09032].

Because \(A(t)\) is singular, the system is not a standard SDE. The analysis proceeds by using Itô’s formula on \(k(t,x)=A(t)x\), which yields
\[
A(t)X(t)=A(0)X_0+\int_0^tA'(s)X(s)+B(s)X(s)+f(s, X(s))ds
+\int_0^tg(s, X(s))dW(s).
\]
To define an index-1 SDAE, the paper requires that the noise does not enter the algebraic constraints,
\[
Img(t,X)\subseteq ImA(t),
\]
and that the constraint equations are globally uniquely solvable [2509.09032].

The reduction uses projector objects associated with \(A(t)\): the pseudo-inverse \(A^-(t)\), the projector \(P(t)=A^-(t)A(t)=I-Q(t)\), and \(R(t)\) with \(R(t)A(t)=0\). Under these projectors, the SDAE is reduced to an equivalent inherent SDE in a reduced variable \(u(t)\), with the algebraic component recovered through \(\hat v(t,u(t))\). The reduced form is
\[
\left \{ \begin{array}{c c c}
u(t)-u_0 &= \int_{0}^{t}P'(s)\left[ u(s)+\hat{v}(s,u(s))\right] +A^-(s) B(s)\left[u(s)+\hat{v}(s,u(s))\right]ds\\
&+\int_{0}^{t}A^-(s)f(s,u(s)+\hat{v}(s,u(s))ds\\
&+ \int_{0}^{t}A^-(s)g(s,u(s)+\hat{v}(s,u(s)))dW(s),\\
\hat{v}(t,u(t))&=-(A(t)+R(t)B(t))^{-1}\left[R(t)B(t)u(t)+f_1(t,u(t))\right],\\
u_0&=P(0)X_0,\\
\hat{v}(0,u_0)&=Q(0)X_0,
\end{array}\right.
\]
where
\[
f_1(t,u(t))=R(t)f(t,u(t)+\hat{v}(t,u(t))).
\]
This decomposition is the exact structural template that the dual tamed scheme mirrors [2509.09032].

## 3. Semi-implicit tamed discretization and its dual form

The original numerical method is the semi-implicit tamed scheme. Its design principle is stated explicitly: the linear component of the drift term is approximated implicitly, whereas its nonlinear component is tamed and approximated explicitly. The scheme is
\[
A(t_n)X_{n+1}^N-A(t_n)X_n^N=hB(t_n)X_{n+1}^N+\frac{f(X_n^N) h}{1+h\left\|f(X_n^N) \right\| }+g(X_n^N)\Delta W_n^N,
\]
or equivalently
\[
(A(t_n)-h B(t_n))X_{n+1}^N=A(t_n)X_{n}^N+\frac{f(X_n^N)h}{1+h\left\|f(X_n^N) \right\| }+g(X_n^N)\Delta W_n^N,
\]
with \(h=\frac{T}{N}\) and \(\Delta W_n^N=W(t_{n+1})-W(t_n)\) [2509.09032].

The explicit decomposition of the update is central. The linear drift term \(B(t_n)X_{n+1}^N\) is handled implicitly, the nonlinear drift \(f(X_n^N)\) is tamed through
\[
\frac{f(X_n^N)h}{1+h\|f(X_n^N)\|},
\]
and the diffusion remains explicit [2509.09032].

The dual tamed scheme is then stated in reduced coordinates. In Lemma 7, it is written as
\[
\left \{ \begin{array}{c c c}
u_{n+1}-u_n &= P'(t_n)\left[ u_{n+1}+\hat{v}^N(t_{n+1},u_{n+1})\right] h+A^-(t_n) B(t_n)\left[u_{n+1}+\hat{v}^N(t_{n+1},u_{n+1})\right]h\\
&+\cfrac{A^-(t_n)f(t_n,u_{n}+\hat{v}^N(t_{n},u_{n}))h}{1+h\|f(t_n,u_{n}+\hat{v}^N(t_{n},u_{n}))\|} +A^-(t_n)g(t_n,u_{n}+\hat{v}^N(t_{n},u_{n}))\Delta W_n,\\
\hat{v}^N(t_{n+1},u_{n+1})&=-(A(t_n)+R(t_n)B(t_n))^{-1}\left(R(t_n)B(t_n)u_{n+1}+\dfrac{R(t_n)f(t_n,u_{n}+\hat{v}^N(t_{n},u_{n}))}{1+h\|f(t_n,u_{n}+\hat{v}^N(t_{n},u_{n}))\|} \right),\\
u_0&=P(0)\zeta,\\
\hat{v}^N(0,u_0)&=Q(0)\zeta,\\
X_{n+1}&=u_{n+1} +\hat{v}^N(t_{n+1},u_{n+1}),
\end{array}\right.
\]
where \(\hat v^N\) is globally Lipschitz in \(u\) [2509.09032].

This dual form is the numerical counterpart of the exact decomposition \(X(t)=u(t)+\hat v(t,u(t))\). The paper states that the original semi-implicit scheme in \(X_n^N\) is equivalent to the dual scheme in \((u_n,\hat v^N)\), provided \(A(t_n)+R(t_n)B(t_n)\) is nonsingular. This is the precise sense in which the method is dual rather than separate [2509.09032].

## 4. Assumptions, stability structure, and strong convergence

The coefficient assumptions for the SDAE formulation are:
\[
\text{(A.1)} \quad f \text{ is locally Lipschitz, and } g \text{ is globally Lipschitz;}
\]
\[
\text{(A.2)} \quad f \text{ is one-sided Lipschitz: }  
\langle X-Y, f(X)-f(Y)\rangle \leq C\|X-Y\|^2;
\]
\[
\text{(A.3)} \quad f \text{ has superlinear growth control: } \|f(X)-f(Y)\| \leq C(1+\|X\|^c+\|Y\|^c)\|X-Y\|;
\]
\[
\text{(A.4)} \quad \mathbb{E}\|X_0\|^2\le M.
\]
These are the non-global Lipschitz hypotheses under which the method is analyzed [2509.09032].

The matrix and projector assumptions are equally important. Assumption 2 requires that \(A, A^-, P, P', B\) are bounded and Lipschitz; that \((A(t_n)+R(t_n)B(t_n))\) and \(S_h=A(t_n)-hB(t_n)\) are nonsingular with bounded inverses; and that
\[
I_d+hM_1(t_n)
\]
is nonsingular with
\[
|(I_{d}+hM_1(t_n))^{-1}|_1\leq \exp{(Kh)},\,\,n \in \{ 0,1,2,..., N-1 \},
\]
where
\[
\begin{align*}
M_1(t_n)&=-P'(t_n)+P'(t_n)[A(t_n)+R(t_n)B(t_n)]^{-1}R(t_n)B(t_n)\\
&-A^-(t_n) B(t_n)+A^-(t_n) B(t_n)[A(t_n)+R(t_n)B(t_n)]^{-1}R(t_n)B(t_n).
\end{align*}
\]
This algebraic structure is what makes the reduced formulation mathematically tractable [2509.09032].

The main convergence theorem states that, for all \(p \in [1, \infty)\), there exists a constant \(C >0\) such that
\[
\left(\mathbb{E}\sup_{t\in [0,T]}\|X_t-\bar{X}_t^N\|^p\right)^{\frac{1}{p}\leq Ch^{\frac{1}{2},~~~h=\frac{T}{N}.
\]
Thus the scheme converges strongly with order \(1/2\) [2509.09032].

The proof structure combines several intermediate results. Theorem 2 gives uniform moment bounds for the numerical solution,
\[
\sup_{N\in \mathbb{N}\sup_{n}\mathbb{E}[\|X_n^N\|^p]<\infty,
\]
Lemma 10 gives interpolation estimates for the reduced variable,
\[
\sup_{t\in [0,T]}\left(\mathbb{E}\left[\|\bar{u}^N_t-\bar{u}^N_{\lfloor t\rfloor}\|^p\right]\right)^{\frac{1}{p}\leq C_ph^{\frac{1}{2},
\]
and Lemma 9 controls the algebraic component through
\[
\left( \mathbb{E}\sup_{t\in [0,T]} \|\hat{v}(t,u(t))- \hat{v}^N(t,\bar{u}^N(t))\|^p\right )^{\frac{1}{p}
\leq Ch+C \left( \mathbb{E}\sup_{t\in [0,T]} \|u(t)-\bar{u}^N(t)\|^p\right )^{\frac{1}{p}.
\]
The final estimate for \(u_t-\bar u_t^N\) is then transferred to \(X_t-\bar X_t^N\) [2509.09032].

## 5. Relation to other tamed constructions

The dual tamed scheme belongs to a broader family of explicit or semi-implicit stabilizations for stochastic systems with superlinear drift, but its defining feature is the projector-based reduced-coordinate reformulation for SDAEs.

For jump–diffusion SDEs, one related construction is the semi-tamed Euler scheme, where the drift is split as \(f=u+v\) and only the nonglobally Lipschitz part is tamed:
\[
Y_{n+1} = Y_n + u(Y_n)\Delta t + \frac{\Delta t\, v(Y_n)}{1+\Delta t\|v(Y_n)\|} + g(Y_n)\Delta W_n + h(Y_n)\Delta N_n.
\]
The same paper also studies the fully tamed drift scheme
\[
X_{n+1} = X_n + \frac{\Delta t\, f(X_n)}{1+\Delta t\|f(X_n)\|} + g(X_n)\Delta W_n + h(X_n)\Delta N_n,
\]
and compares their linear and nonlinear mean-square stability under non-global Lipschitz conditions [1510.06067]. This is a duality only in the informal sense of comparing two taming placements, not in the reduced-coordinate SDAE sense.

For neutral stochastic delay differential equations, the tamed \(\theta\)-Euler-Maruyama scheme is semi-implicit in the drift through a \(\theta\)-weight:
\[
y_{k+1}-D(y_{k+1-m}) = y_k-D(y_{k-m}) +\theta\, b_\Delta(y_{k+1},y_{k+1-m})\Delta +(1-\theta)\,b_\Delta(y_k,y_{k-m})\Delta +\sigma_\Delta(y_k,y_{k-m})\Delta W_k.
\]
Its tamed coefficients satisfy bounds such as
\[
|b_\Delta(x,y)| \le \min\!\big(K_5 \Delta^{-a}(1+|x|+|y|),\,|b(x,y)|\big),
\]
and the strong convergence estimate
\[
\mathbb E\Big[\sup_{0\le t\le T}|X(t)-Y_\Delta(t)|^p\Big] \le C\,\Delta^{ap}
\]
under global monotonicity assumptions [1612.02800]. Again, this is not presented as a dual scheme.

For McKean–Vlasov SDEs driven by Lévy noise, the numerical method is a tamed Euler scheme applied to an interacting particle system,
\[
x_t^{i,N,n} = x_0^i + \int_0^t b_{\kappa(n,s)}\!\big(x_{\kappa(n,s)}^{i,N,n},\mu_{\kappa(n,s)}^{N,n}\big)\,ds + \cdots,
\]
with taming exemplified by
\[
b_n(x,\mu)=\frac{f(x,\cdot)}{1+n^{-1/2}|x|},\qquad \sigma_n(x,\mu)=\frac{g(x,\cdot)}{1+n^{-1/2}|x|},\qquad \gamma_n(x,\mu,z)=\frac{h(x,\cdot,z)}{1+n^{-1/2}|x|}.
\]
The paper states that there is no dual scheme there; the relevant viewpoint is only the two-level approximation by particle approximation and time discretization [2010.08585].

A further related variant is the modified tamed scheme
\[
b^h(t,x)=\frac{b(t,x)}{1+\psi(\gamma h^\alpha |b(t,x)|)},
\]
where the drift is unchanged whenever \(\gamma h^\alpha |b|\le 1\), and the paper proves strong order \(1/2\) and weak order \(1\) for Euler under the stated assumptions [2507.09475]. This suggests that contemporary taming research distinguishes methods mainly by where and when the damping acts: on the full drift, on a decomposed nonlinear component, on an implicit-explicit split, or on a reduced coordinate system.

## 6. Example, interpretation, and limitations

The representative numerical example for the dual tamed scheme is a three-dimensional SDAE,
\[
\begin{pmatrix} 1 & 0 &1  \\ 0 & 0& 0\\ 1&0 &0 \end{pmatrix}dX=[BX(t)+f(t,X(t))]dt+g(t,X(t))dW(t),~~ X(0)=(10^{-2};0;10^{-2}), ~~t\in\left[0,T \right],
\]
with
\[
B=\begin{pmatrix} 0 & 0 &0  \\ 0 & -1& 0\\ 0&0 &1 \end{pmatrix},
\]
\[
f(x_1,x_2,x_3) =\begin{pmatrix} x_1  \\ x_2-x_2^5\\ x_3-x_3^3 \end{pmatrix},
\qquad
g(x_1,x_2,x_3)=\begin{pmatrix} x_1-x_3 & x_2&x_3  \\ 0 & 0& 0\\ x_2-x_3&	x_1&x_2-x_1 \end{pmatrix}.
\]
The corresponding projectors satisfy
\[
A^-=\begin{pmatrix} 0 & 0 &1  \\ 0 & 0& 0\\ 1&0&-1 \end{pmatrix}, \qquad R=I_{3\times 3}-AA^-=\begin{pmatrix} 0 & 0 &0  \\ 0 & 1& 0\\ 0&0&0 \end{pmatrix},
\]
and the constraint equation yields
\[
v_1(t)=0,\qquad v_3(t)=0,\qquad v_2(t)=-u_2(t),
\]
showing that the SDAE is index-1 [2509.09032].

The paper further verifies the one-sided Lipschitz property
\[
\left\langle (X-Y)^T,f(t,X)-f(t,Y)\right\rangle\le \|X-Y\|^2
\]
and the matrix estimate
\[
(I_{3} +hM_1(t_n))^{-1}=\begin{pmatrix} 1 & 0 &\frac{h}{h+1}  \\ 0 & 1& 0\\ 0&0&\frac{1}{1+h} \end{pmatrix},
\qquad
|(I_{3} +hM_1(t_n))^{-1} |_1=1\leq \exp{(kh)}.
\]
The simulation uses \(N=2^{18}\), \(T=1\), \(X_0=(10^{-2},0,10^{-2})\), and 300 Monte Carlo samples; the resulting log-log plot shows strong convergence with slope \(1/2\), matching the theorem [2509.09032].

The principal limitation is terminological rather than mathematical. In the available literature, “dual tamed scheme” is tied specifically to the projector-based reduction of index-1 SDAEs. By contrast, in other settings the literature either uses a tamed Euler scheme without any dual formulation, or introduces different modifiers such as semi-taming, \(\theta\)-taming, or selective cut-off taming [2010.08585], [1510.06067], [1612.02800], [2507.09475]. This suggests that the term names an analytical representation of a scheme adapted to constrained stochastic dynamics, rather than a universal class of tamed integrators.

Source: https://www.emergentmind.com/topics/dual-tamed-scheme