---
title: Operator-Adapted Wavelet Decomposition
url: https://www.emergentmind.com/topics/operator-adapted-wavelet-decomposition-algorithm
type: topic
---

# Operator-Adapted Wavelet Decomposition

Operator-adapted wavelet decomposition denotes a class of multiscale constructions in which the decomposition is tailored to a given operator rather than chosen independently of it. In the literature represented here, this adaptation appears in several mathematically distinct forms: the wavelet-vaguelette decomposition for linear inverse problems, the gamblet transform for symmetric positive operators, diagonal frame decompositions that generalize the SVD to redundant frames, generalized wavelets built from a Fourier multiplier operator \(L\), and fast product-convolution-to-wavelet algorithms that construct sparse wavelet matrices of spatially varying operators. A common objective is to obtain sparse, localized, or block-diagonal representations in which approximation, inversion, denoising, eigenpair computation, or multiscale finite-element correction become more efficient than in operator-agnostic coordinates [1310.7149], [1805.10736], [1210.1808], [2005.09870].

## 1. Principal meanings of operator adaptation

A first meaning of operator adaptation is the wavelet-vaguelette decomposition (WVD). For a linear operator \(K\), the decomposition assumes wavelets \((\phi_{j_0k},\psi_{jk})\), vaguelette systems \(\mathcal U=\{u_{jk}\}\) and \(\mathcal V=\{v_{jk}\}\), and scale-dependent pseudo-singular values \(\kappa_j\) such that
\[
K\psi_{jk}=\kappa_j v_{jk}, \qquad Ku_{jk}=\kappa_j \psi_{jk},
\]
with
\[
\langle u_{jk},v_{j'k'}\rangle=\delta_{j-j'}\delta_{k-k'}.
\]
This converts the white-noise inverse problem
\[
Y_\epsilon(dt)=(Kf)(t)\,dt+\epsilon W(dt)
\]
into the heteroscedastic Gaussian sequence model
\[
y_{jk}=\theta_{jk}+\epsilon_j z_{jk}, \qquad \epsilon_j=\epsilon \kappa_j^{-1},
\]
or, under the normalization \(\kappa_j^{-1}=2^{\beta j}\),
\[
y_{jk}=\theta_{jk}+\epsilon 2^{\beta j} z_{jk}.
\]
Here the decomposition is operator-adapted because the analysis functions are not raw wavelets of the data but vaguelettes tied to \(K\) [1310.7149].

A second meaning is the gamblet construction for symmetric positive operators \(L:H_0^s(\Omega)\to H^{-s}(\Omega)\). Starting from nested measurement functions \(\phi_i^{(k)}\), one defines operator-adapted pre-wavelets
\[
\psi_i^{(k)}=\sum_{j\in I^{(k)}}A_{i,j}^{(k)}L^{-1}\phi_j^{(k)},
\]
and operator-adapted wavelets
\[
\chi_i^{(k)}:=\sum_{j\in I^{(k)}}W_{i,j}^{(k)}\psi_j^{(k)}.
\]
The corresponding spaces satisfy
\[
\Psi^{(k)}=\Psi^{(k-1)}\oplus X^{(k)}
\]
orthogonally in the energy inner product
\[
\langle u,v\rangle=\int_\Omega uLv.
\]
Here adaptation is with respect to the operator-induced energy geometry rather than the \(L^2\) inner product [1805.10736], [1806.00565].

A third meaning is diagonalization by operator-adapted frames. A diagonal frame decomposition (DFD) for \(A:X\to Y\) is a triple \((u_\lambda,v_\lambda,\kappa_\lambda)\) such that \((u_\lambda)\) is a frame of \(X\), \((v_\lambda)\) is a frame of \(\operatorname{ran}(A)\), and
\[
A^*v_\lambda=\kappa_\lambda u_\lambda.
\]
This yields the reconstruction identity
\[
x=\bar U\,M_\kappa^+\,V^*(Ax),
\]
and a direct thresholding rule in operator-adapted coordinates. The construction generalizes the SVD to redundant systems such as wavelets, curvelets, and shearlets [1909.09364].

A fourth meaning is basis construction from a Fourier multiplier operator. For an admissible operator \(L\) with symbol \(\widehat L(\omega)\), the generalized scaling functions are
\[
\widehat\varphi_j(\omega)=\widehat L_{\mathrm d,j}(\omega)\widehat L(\omega)^{-1},
\]
the \(L^*L\)-interpolating spline is
\[
\widehat\phi_j(\omega)=\det(\mathbf D)^j
\frac{|\widehat\varphi_j(\omega)|^2}{
\sum_{k\in\mathbb Z^d}|\widehat\varphi_j(\omega+2\pi(\mathbf D^T)^{-j}k)|^2},
\]
and the wavelet is defined by
\[
\psi_{j+1}=L^*\phi_j.
\]
In this setting, the wavelet acts like the operator across scales [1210.1808].

A fifth meaning is representation-theoretic decomposition for wavelet representations associated with a finite-to-one endomorphism \(r:X\to X\). After lifting to the solenoid \(X_\infty\), one obtains a covariant representation \((H,U,\pi)\) and a direct-integral decomposition
\[
[H,U,\pi]=\int_F^\oplus [H_z,U_z,\pi_z]\,d\mu_\infty(z),
\]
where \(F\) is a fundamental domain and the fiber representations are irreducible for \(\mu_\infty\)-almost every \(z\). This is an operator-adapted decomposition in a representation-theoretic sense rather than a numerical sparse-matrix sense [1105.3442].

## 2. Fast construction of sparse wavelet matrices from product-convolution structure

A particularly explicit algorithmic realization is given for spatially varying operators through product-convolution expansions. Let \(\mathcal E\) be the space of discrete \(d\)-dimensional signals on
\[
\Omega=\{1,\ldots,N_1\}\times\cdots\times\{1,\ldots,N_d\},
\qquad N=\prod_{i=1}^d N_i,
\]
with periodic boundary conditions, and let
\[
G:f\mapsto \sum_{y\in\Omega}K(\cdot,y)f(y).
\]
The target object is the wavelet-domain matrix
\[
\Theta=\Psi^*H\Psi, \qquad \theta[\lambda,\mu]=\langle H\psi_\lambda,\psi_\mu\rangle,
\]
for a product-convolution approximation
\[
Hf=\sum_{k=1}^m u_k\star(v_k\odot f).
\]
The naive cost of generic wavelet decomposition is stated as \(O(N^3)\) in general. The product-convolution model replaces this by a structurally adapted factorization
\[
H=\sum_{k=1}^m U_kV_k, \qquad
\Theta=\sum_{k=1}^m A_kB_k,
\]
with
\[
A_k=\Psi^*U_k\Psi,\qquad B_k=\Psi^*V_k\Psi.
\]
The adaptation occurs because the local kernels \(u_k\) and coefficient maps \(v_k\) are derived from the operator’s spatially varying impulse responses, by direct interpolation, adaptive sampling, PCA on sampled impulse responses, or randomized SVD of the SVIR [2005.09870].

The algorithm computes an \(\eta\)-accurate sparse approximation \(\widetilde\Theta\) of \(\Theta\) such that
\[
\|\Theta-\widetilde\Theta\|_{2\to2}\le \eta.
\]
Its pipeline is explicit. One first builds or obtains the product-convolution expansion. For each term \(k\), one sets
\[
\epsilon_k=\frac{\eta}{m\|v_k\|_\infty}.
\]
The convolution part
\[
A_k=\Psi^*U_k\Psi
\]
is then computed by exploiting the fact that the wavelet subbands of a convolution operator are circulant; each subband of \(A_k\) is a rectangular circulant matrix. Consequently, \(A_k\) can be encoded using only \(O(N\log_2 N)\) coefficients and computed in
\[
O(N\log_2^2 N)
\]
operations. One then truncates \(A_k\) to \(\widetilde A_k\) so that
\[
\|A_k-\widetilde A_k\|_{2\to2}\le \epsilon_k.
\]
The multiplier part
\[
B_k=\Psi^*V_k\Psi
\]
is sparse because
\[
B_k[\lambda,\mu]=\langle v_k\odot \psi_\lambda,\psi_\mu\rangle
\]
vanishes unless \(\psi_\lambda\) and \(\psi_\mu\) overlap. If
\[
\mathcal I=\{(\lambda,\mu)\in\Lambda\times\Lambda\;|\;\psi_\lambda\cap\psi_\mu\neq\emptyset\},
\]
then
\[
\mathrm{supp}(B_k)\subseteq \mathcal I,
\qquad
\#\mathcal I \le \frac{c(d)}{d}\delta^d N\log_2 N.
\]
The matrix \(B_k\) is computed by a two-stage sparse cascade algorithm with cost
\[
O(\delta^{d+1}N\log_2 N).
\]
Finally, one forms
\[
C_k=\widetilde A_kB_k,
\qquad
\widetilde\Theta\leftarrow \widetilde\Theta+C_k.
\]

The error control is spectral-norm explicit:
\[
\|\Theta-\widetilde\Theta\|_{2\to2}
\le
\sum_{k=1}^m \|B_k\|_{2\to2}\|A_k-\widetilde A_k\|_{2\to2}.
\]
Because \(\Psi\) is orthogonal and \(V_k\) is diagonal,
\[
\|B_k\|_{2\to2}=\|v_k\|_\infty,
\]
so the choice \(\epsilon_k=\eta/(m\|v_k\|_\infty)\) gives
\[
\|\Theta-\widetilde\Theta\|_{2\to2}\le \eta.
\]
Under the regularity assumptions, Algorithm 1 computes \(\widetilde\Theta\) in
\[
O\!\left(\delta^{d+1} m^{\max(1,d/\alpha)} N \log_2^2 N\, \eta^{-d/\alpha}\right)
\]
operations, and the number of nonzero coefficients in \(\widetilde\Theta\) is
\[
O\!\left(\delta^d N \log_2^2 N\, \eta^{-d/\alpha}\right).
\]
In many smooth applications \(\alpha\ge d\), so \(\max(1,d/\alpha)=1\), which is the quasi-linear regime emphasized by the authors [2005.09870].

## 3. Operator-adapted decompositions in inverse problems, denoising, eigenanalysis, and basis design

In inverse problems, the WVD provides an operator-adapted coordinate model in which estimation becomes levelwise penalized regression. The estimator at level \(j\) is
\[
\widehat \theta_j=\arg\min_{\mu\in\mathbb R^{2^j}}
\Bigl\{\|y_j-\mu\|^2+\epsilon_j^2\,\mathrm{pen}_j(N(\mu))\Bigr\},
\]
with scale-dependent penalty calibrated by the ill-posedness of \(K\). The main theorem partitions the parameter domain into dense, sparse, and critical regimes and states exact minimax rates over Besov sequence balls. In this setting, operator adaptation enters through the vaguelettes and pseudo-singular values, while adaptation to unknown regularity enters through levelwise complexity penalization [1310.7149].

For denoising with prior information on \(Lu\) rather than on \(u\), gamblet coefficients replace ordinary wavelet coefficients. The partial reconstruction at level \(l\) is
\[
\eta^{(l)}:=\sum_{k=1}^{l}[\phi^{(k),\chi},\eta]\cdot \chi^{(k)},
\]
and the main estimator is
\[
v^\dagger(\eta)=\eta^{(l^\dagger)},
\]
where \(l^\dagger\) minimizes
\[
\beta_l=
\begin{cases}
h^{2s}M^2, & l=0,\\[1mm]
\sigma^2 h^{-(2s+d)l}+h^{2s(l+1)}M^2, & 1\le l\le q-1,\\[1mm]
h^{-(2s+d)q}\sigma^2, & l=q.
\end{cases}
\]
The paper also gives hard and soft thresholding in gamblet coordinates. The principal guarantee is near minimax optimality, up to a multiplicative constant, over
\[
V_M^{(q)}=\{u\in \Psi^{(q)}:\|Lu\|_{L^2(\Omega)}\le M\},
\]
in the energy norm induced by the operator [1805.10736].

The same gamblet decomposition supports eigenpair computation. The multiresolution splitting
\[
V= W^{(1)}\oplus W^{(2)}\oplus\cdots\oplus W^{(q)}
\]
is energy-orthogonal, with level matrices
\[
B^{(k)}=W^{(k)}A^{(k)}W^{(k),T}.
\]
The algorithm first solves the eigenproblem on a coarse gamblet space and then performs hierarchical subspace correction. At each level \(k\), the correction solve is
\[
\langle \widetilde v^{(k,\ell+1)}, w\rangle
=
\lambda^{(k,\ell)} [v^{(k,\ell)},w], \qquad \forall w\in V^{(k)},
\]
followed by a Rayleigh-Ritz step in
\[
V^{(1,k)}=V^{(1)}+\operatorname{span}\{\widehat v^{(k,\ell+1)}\}.
\]
The method is designed for bijective positive symmetric operators and is shown to have near-linear complexity when \(\mathcal L\) is a local operator mapping \(\mathcal H_0^s(\Omega)\) to \(\mathcal H^{-s}(\Omega)\) [1806.00565].

For direct sparse regularization, the DFD yields the explicit thresholding formula
\[
B_\alpha^{\mathrm{DFD}}(y)
=
(\bar U\circ M_\kappa^+\circ S_{\alpha d/\kappa}\circ V^*)(y).
\]
This is a non-iterative operator-adapted frame thresholding method. If the frame is a basis, then synthesis and analysis \(\ell^1\)-regularization as well as DFD thresholding are equivalent; in the redundant case, those three approaches are pairwise different [1909.09364].

At the level of basis construction, operator-like wavelets are generated from a Fourier multiplier operator by first defining generalized B-splines and then setting
\[
\psi_{j+1}=L^*\phi_j.
\]
In the stochastic model
\[
Ls=w,
\]
the coefficient identity
\[
\langle s,\psi\rangle=\langle Ls,\phi\rangle=\langle w,\phi\rangle
\]
shows why sparsity of the innovation can be inherited by the wavelet coefficients [1210.1808].

## 4. Hierarchical finite-element realizations on unstructured meshes

A recent finite-element line of work constructs operator-adapted wavelet decompositions on unstructured triangular meshes with Whitney edge elements for multiscale electromagnetic problems. The continuous equation is
\[
\mathcal L\mathbf u=\mathbf g,
\]
and in the principal experiments the governing PDE is the vector Helmholtz equation
\[
\left(\nabla\times \left[\frac{1}{\alpha}\nabla\times\right]-k_0^2\theta\right)\mathbf{u}(\mathbf r)=\mathbf g(\mathbf r)\quad \text{in }\Omega,
\]
with boundary conditions
\[
\mathcal{B}_i\{\mathbf{u}(\mathbf r)\}=\mathbf b_i(\mathbf r)\quad \text{on }\Gamma_i.
\]
Starting from nested FEM spaces \(V^1\subset V^2\subset\cdots\subset V^q\), the operator-adapted hierarchy imposes
\[
\mathbb V^{j+1}=\mathbb V^j \oplus_{\mathcal L} \mathbb W^j,
\]
so that
\[
\mathbb V^q=\mathbb V^1\oplus_{\mathcal L}\mathbb W^1\oplus_{\mathcal L}\cdots\oplus_{\mathcal L}\mathbb W^{q-1}.
\]
The transformed stiffness matrix becomes block diagonal,
\[
\mathbb L=
\begin{bmatrix}
\mathbb A^1 & 0 & \cdots & 0\\
0 & \mathbb B^1 & \cdots & 0\\
\vdots & \vdots & \ddots & \vdots\\
0 & 0 & \cdots & \mathbb B^{q-1}
\end{bmatrix},
\]
which decouples the scales. The hierarchy is built from finest to coarser levels using sparse refinement matrices \(\tilde{\mathbf C}^j\), null-space matrices \(\tilde{\mathbf W}^j\), and recursively defined operator-adapted transfers \(\mathbb C^j\). The coarse and detail right-hand sides are
\[
\mathbbm g^1=\mathbb C^1\mathbb C^2\cdots \mathbb C^{q-1}\mathbbm g^q,
\qquad
\mathbbm b^j=\tilde{\mathbf W}^j \mathbb C^{j+1}\cdots \mathbb C^{q-1} \mathbbm g^q,
\]
and the independent systems are
\[
\mathbb A^1\mathbbm v^1=\mathbbm g^1,
\qquad
\mathbb B^j\mathbbm w^j=\mathbbm b^j.
\]
The full solution is reconstructed as
\[
\mathbbm u^\zeta
=
\mathbbm s^{\text{coarse}}
+\sum_{j=1}^{\zeta-1}\mathbbm s^{\text{detail},j}.
\]
The principal numerical claim is that the overall computational complexity is nearly linear, with sparse matrix-vector multiplications as the dominant operations [2507.17989].

A related polygonal-element formulation begins from a finest triangular mesh and generates coarser meshes by agglomerating triangles into convex polygons. The operator-adapted hierarchy is again
\[
\mathbb{V}^{j+1}=\mathbb{V}^{j}\oplus_{\mathcal L}\mathbb{W}^{j},
\]
with operator orthogonality
\[
\int_{\Omega}\mathoutline{\phi}_{i}^{\,j}\,\mathcal{L}\,\mathoutline{\psi}_{\ell}^{\,j,T}\,d\Omega = 0.
\]
On a convex polygon \(P\subset\mathbb R^2\), generalized barycentric coordinates \(\{\lambda_i(x)\}\) are used to define Whitney one-form edge bases
\[
\mathbf{W}_{\ell,j}(x)=\lambda_{\ell}(x)\,\nabla \lambda_{j}(x)-\lambda_{j}(x)\,\nabla \lambda_{\ell}(x).
\]
The operator-agnostic refinement matrices \(\tilde{\mathbf C}^j\) and null-space matrices \(\tilde{\mathbf W}^j\) are assembled first, after which the operator-adapted matrices \(\mathbb C^j\) and compressed operators \(\mathbb A^j\) are defined recursively. As in the triangular-mesh formulation, the detail levels can be added without recomputing the coarser levels [2512.16004].

## 5. Computational properties and demonstrated applications

The product-convolution-to-wavelet algorithm was validated on two space-varying image blurs on images up to \(4096\times 4096\): a vertically varying isotropic Gaussian blur and a rotated/skewed Gaussian blur varying with location. The product-convolution ranks were chosen as \(m=5\) and \(m=25\), and the decomposition precision was \(\eta=5\times 10^{-4}\). For \(N=16\) million pixels, the reported speed-up over naive decomposition using product-convolution was \(4680\times\) for blur 1 and \(2330\times\) for blur 2. In wavelet-\(\ell_1\) deblurring on \(1024\times 1024\) images, wavelet-domain FISTA using the sparse wavelet operator was much faster than spatial-domain or product-convolution FISTA. For blur 1, FISTA-W gave \(16.2\times\) speed-up over exact FISTA, preconditioned FISTA-WP gave \(30.7\times\), and GPU \(670\times\); for blur 2, the corresponding figures were \(62\times\), \(120\times\), and \(4089\times\). The loss in image quality was less than \(0.2\) dB relative to the exact model in that setup [2005.09870].

In gamblet denoising, the central guarantee is not an empirical benchmark but a statistical one:
\[
\sup_{u\in V_M^{(q)}} E\big[\|u-v^\dagger(\eta)\|^2\big]
<
C\inf_{v(\eta)} \sup_{u\in V_M^{(q)}} E\big[\|u-v(\eta)\|^2\big].
\]
The paper also proves a high-probability energy-norm bound stating that, if \(l^\dagger\neq 0\), then with probability at least \(1-\varepsilon\),
\[
\|\eta^{(l^\dagger)}\|
\le
\|u\|
+
C\sqrt{ \log\frac1\varepsilon\; \sigma^{\frac{2s+d}{4s+d}} M^{\frac{2s+2d}{4s+d}} }.
\]
The transform itself is computable to precision \(\epsilon\) in
\[
O\!\left(N\log^{3d}\frac{N}{\epsilon}\right)
\]
complexity, with an incomplete Cholesky version stated as
\[
O\!\left(N\log^2N\,\log^{2d}\frac{N}{\epsilon}\right).
\]
This places gamblet thresholding among near-linear operator-adapted decomposition methods [1805.10736].

For low eigenpairs of rough operators, gamblet multilevel correction was tested on the SPE10 elliptic benchmark, random checkerboard coefficients, and Anderson localization or disordered Schrödinger-type problems. The reported observations were that gamblet-based multilevel correction converges much faster than geometric multigrid-based correction, that the first 12 eigenvalues improve dramatically level by level, and that in the SPE10 case the gamblet multilevel correction method showed about a ten-fold speedup in online CPU time versus ARPACK in the reported implementation [1806.00565].

In multiscale electromagnetics on unstructured triangular meshes, the L-shaped and U-shaped waveguide problems were used as validation. For the L-shaped waveguide, the relative \(L^2\) error between the finest-level FEM solution and the hierarchical operator-adapted solution with coarse plus three detail levels was about
\[
2.5\times 10^{-6},
\]
and the relative \(L^2\) error versus numerical mode matching was
\[
1.21\times 10^{-4}.
\]
For the U-shaped waveguide, the corresponding errors were approximately
\[
8\times 10^{-6}
\quad\text{and}\quad
7.6\times 10^{-4}.
\]
For finest-level DoFs ranging from \(12800\) to \(1638400\), the reported cost was about \(\mathcal O(N^{0.91})\) per iteration in a two-level scenario, about \(\mathcal O(N^{0.92})\) excluding precomputation in a six-level scenario, and about \(\mathcal O(N^{1.07})\) including precomputation [2507.17989].

In the polygonal-element formulation, the wedge singularity problem and a microporous Si slab were the main examples. The reported fitted slopes from the wedge example were approximately \(\mathcal O(N^{0.94})\) for the coarsest-level calculation, \(\mathcal O(N^{0.59})\) for the first detail, \(\mathcal O(N^{0.37})\) for the second detail, \(\mathcal O(N^{0.23})\) for the third detail, and \(\mathcal O(N^{0.11})\) for the fourth detail; including precomputation, the global fitted slope was about \(1.04\). With four detail levels, the relative \(L^2\)-error in the near-tip region dropped below \(2\%\) relative to finest-level FEM, while the error versus the analytic series solution was about \(3\%\) or less in the same region. The paper also reports peak memory savings up to about \(40\%\) [2512.16004].

## 6. Adjacent methods, common misconceptions, and structural limitations

A recurrent misconception is that any use of wavelets inside an operator-learning or PDE pipeline constitutes an operator-adapted wavelet decomposition. The Wavelet Neural Operator is explicitly described as **not a classical operator-adapted wavelet decomposition algorithm in the numerical-analysis sense**. It uses a fixed multilevel discrete wavelet transform, with the operator-specific part learned through trainable kernels acting on wavelet coefficients. The basis is fixed in advance, with wavelet families such as \(\mathrm{db}4\) and \(\mathrm{db}6\), and there is no algorithm that derives wavelets from the PDE operator. What is adapted is the learned mapping in wavelet space, not the decomposition itself [2205.02191].

A different nearby construction is the monogenic synchrosqueezed wavelet transform for AM-FM image decomposition. It is adapted to the monogenic/Riesz operator and to a structured class of oscillatory image models, with local wave-vector estimators
\[
\Lambda_1(a,b)=\partial_{b_1} c_F(a,b)\times(c_F(a,b))^{-1},
\qquad
\Lambda_2(a,b)=\partial_{b_2} c_F(a,b)\times(c_F(a,b))^{-1},
\]
and a reassigned representation
\[
S_{F,\varepsilon}^{\delta,\nu}(b,k,n).
\]
The principal theorem states that significant wavelet coefficients occur only near the scale matching one mode and that each separated mode can be reconstructed by integrating the synchrosqueezed representation over a small neighborhood of its estimated local wave vector. This is an operator-adapted decomposition in the sense of analysis geometry and reassignment, but not in the sense of building wavelets from an arbitrary elliptic, convolution, or inverse operator [1211.5082].

Across the numerical-analysis literature, the main limitations are structural. The product-convolution route assumes that the operator either is already in product-convolution form or can be approximated well by one; if the SVIR is not low-rank or local translation invariance is weak, \(m\) may need to be large. The guarantees in that setting concern decomposition of \(H\), not the original operator \(G\), unless \(H\approx G\). The WVD approach applies only to operators admitting a wavelet-vaguelette decomposition, with scale-only pseudo-singular values and frame bounds. Gamblet constructions assume a symmetric positive bijection and rely on hierarchical measurements, operator locality, and exponentially localized bases for near-linear computation. DFD thresholding requires a diagonal frame decomposition and yields a direct thresholding method, but in the redundant case it is distinct from both analysis and synthesis \(\ell^1\)-regularization. The finite-element electromagnetic formulations depend on preserving sparsity under recursive operator-adapted transformations and presently focus on 2D mesh hierarchies with Whitney edge elements or convex polygonal elements [2005.09870], [1310.7149], [1805.10736], [1909.09364], [2507.17989], [2512.16004].

Taken together, these works show that “operator-adapted wavelet decomposition algorithm” does not denote a single canonical procedure. It denotes a family of constructions in which the multiscale representation is tied to an operator through vaguelettes, energy-orthogonal wavelets, diagonal frame relations, Fourier-multiplier design, product-convolution structure, or operator-orthogonal finite-element detail spaces. This suggests that the most precise use of the term depends on which object is being adapted: the coefficient-extraction map, the basis functions, the transformed stiffness matrix, or the representation itself.

Source: https://www.emergentmind.com/topics/operator-adapted-wavelet-decomposition-algorithm