Diffusion wavelets

Diffusion wavelets are a fast multiscale framework for the analysis of functions on discrete (or discretized continuous) structures like graphs, manifolds, and point clouds in Euclidean space. Diffusion wavelets are an extension of classical wavelet theory from harmonic analysis. Unlike classical wavelets whose basis functions are predetermined, diffusion wavelets are adapted to the geometry of a given diffusion operator [math]\displaystyle{ T }[/math] (e.g., a heat kernel or a random walk). Moreover, the diffusion wavelet basis functions are constructed by dilation using the dyadic powers (powers of two) of [math]\displaystyle{ T }[/math]. These dyadic powers of [math]\displaystyle{ T }[/math] diffusion over the space and propagate local relationships in the function throughout the space until they become global. And if the rank of higher powers of [math]\displaystyle{ T }[/math] decrease (i.e., its spectrum decays), then these higher powers become compressible. From these decaying dyadic powers of [math]\displaystyle{ T }[/math] comes a chain of decreasing subspaces. These subspaces are the scaling function approximation subspaces, and the differences in the subspace chain are the wavelet subspaces. Diffusion wavelets were first introduced in 2004 by Ronald Coifman and Mauro Maggioni at Yale University.

Algorithm

This algorithm constructs the scaling basis functions and the wavelet basis functions along with the representations of the diffusion operator [math]\displaystyle{ T }[/math] at these scales.

In the algorithm below, the subscript notation [math]\displaystyle{ \Phi_a }[/math] and [math]\displaystyle{ \Psi_b }[/math] represents the scaling basis functions at scale [math]\displaystyle{ a }[/math] and the wavelet basis functions at scale [math]\displaystyle{ b }[/math] respectively. The notation [math]\displaystyle{ [\Phi_b]_{\Phi_a} }[/math] denotes the matrix representation of the scaling basis [math]\displaystyle{ \Phi_b }[/math] represented with respect to the basis [math]\displaystyle{ \Phi_a }[/math]. Lastly, the notation [math]\displaystyle{ [T]_{\Phi_a}^{\Phi_b} }[/math] denotes the matrix represents of the operator [math]\displaystyle{ T }[/math], where the row space of [math]\displaystyle{ T }[/math] is represented with respect to the basis [math]\displaystyle{ \Phi_a }[/math], and the column space of [math]\displaystyle{ T }[/math] is represented with respect to the basis [math]\displaystyle{ \Phi_b }[/math]. Otherwise put, the domain of operator [math]\displaystyle{ T }[/math] is represented with respect to the basis [math]\displaystyle{ \Phi_a }[/math] and the range is represented with respect to the basis [math]\displaystyle{ \Phi_b }[/math]. The function [math]\displaystyle{ QR }[/math] is a sparse QR decomposition with [math]\displaystyle{ \epsilon }[/math] precision.

// Input:
//    [math]\displaystyle{ T }[/math] is the matrix representation of the diffusion operator.
//    [math]\displaystyle{ \epsilon }[/math] is the precision of the QR decomposition, e.g., 1e-6.
//    [math]\displaystyle{ J }[/math] is the maximum number of scale levels (note: this is an optional upper bound, it may converge sooner.)
// Output:
//    [math]\displaystyle{ \lbrace\Phi_j\rbrace }[/math] is the set of scaling basis functions indexed by scale [math]\displaystyle{ j }[/math].
//    [math]\displaystyle{ \lbrace\Psi_j\rbrace }[/math] is the set of wavelet basis functions indexed by scale [math]\displaystyle{ j }[/math].
  
[math]\displaystyle{ \lbrace\Phi_j\rbrace, \lbrace\Psi_j\rbrace \leftarrow \text{function DiffusionWaveletTree} ( T , \epsilon , J ): }[/math]
    [math]\displaystyle{ \textbf{for } j\leftarrow 0 \text{ to } J-1 }[/math]:
        [math]\displaystyle{ [\Phi_{j+1}]_{\Phi_j}, [T^{2^j}]_{\Phi_j}^{\Phi_{j+1}} \leftarrow QR\left([T^{2^j}]_{\Phi_j}^{\Phi_{j}}, \epsilon\right)   }[/math]
        [math]\displaystyle{ [T^{2^{j+1}}]_{\Phi_{j+1}}^{\Phi_{j+1}} \leftarrow  \left([T^{2^j}]_{\Phi_j}^{\Phi_{j+1}} [\Phi_{j+1}]_{\Phi_j}\right)^2 }[/math]
        [math]\displaystyle{ [\Psi_j]_{\Phi_j} \leftarrow QR\left(I_{\langle\Phi_j\rangle}-[\Phi_{j+1}]_{\Phi_j}\left([\Phi_{j+1}]_{\Phi_j}\right)^*, \epsilon\right) }[/math] 
    [math]\displaystyle{ \textbf{end for} }[/math]

Applications

Mathematics

Diffusion wavelets are of general interest in mathematics. Specifically, they allow for the direct calculation of the Green′s function and the inverse graph Laplacian.

Computer science

Diffusion wavelets have been used extensively in computer science, especially in machine learning. They have been applied to the following fields:

• solving Markov decision processes and Markov chains for machine learning,
• transfer learning,
• value function approximation in reinforcement learning,
• dimensionality reduction,
• mesh compression for 3D graphics,
• topic model analysis of document corpora.
• relation extraction.

See also

• Wavelets

External links

• Mauro Maggioni's MATLAB code implementation
• Chang Wang's diffusion wavelet page