2025-03-06 | | Total: 13
Many equilibrated flux recovery methods for finite element solutions rely on ad hoc or method-specific techniques, limiting their generalizability and efficiency. In this work, we introduce the Equilibrated Averaging Residual Method (EARM), a unified framework for flux recovery that not only reproduces state-of-the-art locally conservative fluxes but also enables the derivation of new equilibrated fluxes with improved properties. In this paper, EARM is applied to conforming, nonconforming, and discontinuous Galerkin methods, ensuring local conservation and robust a posteriori error estimation. Despite the unified nature of the variational problem, the framework retains the flexibility to fully leverage the inherent properties of finite element spaces. Moreover, EARM offers explicit and computationally efficient flux reconstructions for all methods in two dimensions. In three dimensions, only simple local problems need to be solved for the conforming finite element methods.
In this paper we study the harmonic map heat flow problem for a radially symmetric case. The corresponding partial dfferential equation plays a key role in many analyses of harmonic map heat flow problems. We consider a basic discretization method for this problem, namely a second order finite difference discretization in space combined with a semi-implicit Euler method in time. The semi-implicit Euler method results in a linear problem in each time step. We restrict to the regime of smooth solutions of the continuous problem and present an error analysis of this discretization method. This results in optimal order discretization error bounds (apart from a logarithmic term). We also present discrete energy estimates that mimic the decrease of the energy of the continuous solution.
We consider the solution of the general Sylvester equation $AX+XB=C$ in mixed precision. First, we investigate the use of GMRES-based iterative refinement (GMRES-IR) to solve the equation using implicitly its Kronecker product form: we propose an efficient scheme to use the Schur factors of the coefficient matrices as preconditioners, but we demonstrate that this approach is not suitable in the case of the Sylvester equation. By revisiting a stationary iteration for linear systems, we therefore derive a new iterative refinement scheme for the quasi-triangular Sylvester equation, and our rounding error analysis provides sufficient conditions for convergence and a bound on the attainable relative residual. We leverage this iterative scheme to solve the general Sylvester equation in mixed precision. The new algorithms compute the Schur decomposition of the matrix coefficients in low precision, use the low-precision Schur factors to obtain an approximate solution to the quasi-triangular equation, and iteratively refine it to obtain a working-precision solution to the quasi-triangular equation. However, being only orthonormal to low precision, the unitary Schur factors of $A$ and $B$ cannot be used to recover the solution to the original equation. We propose two effective approaches to address this issue: one is based on re-orthonormalization in the working precision, and the other on explicit inversion of the almost-unitary factors. We test these mixed-precision algorithms on various Sylvester and Lyapunov equations from the literature. Our numerical experiments show that, for both classes of equations, the new algorithms are at least as accurate as existing ones. Our cost analysis, on the other hand, suggests that they would typically be faster than mono-precision alternatives if implemented on hardware that natively supports low precision.
It is demonstrated that one of the equations from the Lie classification list of second-order ODEs is a first integral of the Schwarz equation. As symmetry-preserving finite-difference schemes have been previously constructed for both equations, the preservation of a similar connection between these schemes is studied. It is shown that the schemes for the Schwarz equation and the second-order ODE can be related through a Bäcklund-type difference transformation. In addition, previously unexamined aspects of the difference scheme for the second-order ODE are discussed, including its singular solution and the complete set of difference first integrals.
We provide a finite element discretization of $\ell$-form-valued $k$-forms on triangulation in $\mathbb{R}^{n}$ for general $k$, $\ell$ and $n$ and any polynomial degree. The construction generalizes finite element Whitney forms for the de~Rham complex and their higher-order and distributional versions, the Regge finite elements and the Christiansen--Regge elasticity complex, the TDNNS element for symmetric stress tensors, the MCS element for traceless matrix fields, the Hellan--Herrmann--Johnson (HHJ) elements for biharmonic equations, and discrete divdiv and Hessian complexes in [Hu, Lin, and Zhang, 2025]. The construction discretizes the Bernstein--Gelfand--Gelfand (BGG) diagrams. Applications of the construction include discretization of strain and stress tensors in continuum mechanics and metric and curvature tensors in differential geometry in any dimension.
Tomography is a central tool in medical applications, allowing doctors to investigate patients' interior features. The Radon transform is commonly used to model the corresponding measurement process. Since in any practical application one can only obtain and process finitely many measurements, suitable discretization of the Radon transform and its adjoint (called the backprojection) is crucial. The most commonly used discretization approach combines the ray-driven Radon transform with the pixel-driven backprojection, as anecdotal reports describe these as showing the best approximation performance. However, there is little rigorous understanding of induced approximation errors. These methods involve three discretization parameters: the spatial-, detector-, and angular resolutions. Most commonly, balanced resolutions are used, i.e., the same (or similar) spatial- and detector resolutions. We present a novel interpretation of ray- and pixel-driven discretizations as `convolutional methods'. This allows for a structured analysis that can explain observed behavior. In particular, we prove convergence in the strong operator topology of the ray-driven Radon transform and the pixel-driven backprojection under suitable choice of discretization parameters. Notably, said suitable parameter choices include the combination of the ray-driven Radon transform with the pixel-driven backprojection under balanced resolutions, thus theoretically justifying this approach. Our theoretical investigations are complemented by numerical experiments that show behavior in simulations that the theory predicted or suggested.
In this paper, we investigate the stability and time-step constraints for solving advection-diffusion equations using exponential time differencing (ETD) Runge-Kutta (RK) methods in time and discontinuous Galerkin (DG) methods in space. We demonstrate that the resulting fully discrete scheme is stable when the time-step size is upper bounded by a constant. More specifically, when central fluxes are used for the advection term, the schemes are stable under the time-step constraint tau <= tau_0 * d / a^2, while when upwind fluxes are used, the schemes are stable if tau <= max{tau_0 * d / a^2, c_0 * h / a}. Here, tau is the time-step size, h is the spatial mesh size, and a and d are constants for the advection and diffusion coefficients, respectively. The constant c_0 is the CFL constant for the explicit RK method for the purely advection equation, and tau_0 is a constant that depends on the order of the ETD-RK method. These stability conditions are consistent with those of the implicit-explicit RKDG method. The time-step constraints are rigorously proved for the lowest-order case and are validated through Fourier analysis for higher-order cases. Notably, the constant tau_0 in the fully discrete ETD-RKDG schemes appears to be determined by the stability condition of their semi-discrete (continuous in space, discrete in time) ETD-RK counterparts and is insensitive to the polynomial degree and the specific choice of the DG method. Numerical examples, including problems with nonlinear convection in one and two dimensions, are provided to validate our findings.
We compute the rate of decay of the persistence probabilities of spherical fractional Brownian motion, which was defined by Lévy (1965) and Istas (2005). The rate resembles the Euclidean case treated in Molchan (1999). As a by-product we consider the coefficients of series representations of functions with algebraic endpoint singularities in terms of re-scaled Gegenbauer polynomials, which partly generalises Sidi (2009).
Drawing parallels with the way biological networks are studied, we adapt the treatment--control paradigm to explainable artificial intelligence research and enrich it through multi-parametric input alterations. In this study, we propose a framework for investigating the internal inference impacted by input data augmentations. The internal changes in network operation are reflected in activation changes measured by variance, which can be decomposed into components related to each augmentation, employing Sobol indices and Shapley values. These quantities enable one to visualize sensitivity to different variables and use them for guided masking of activations. In addition, we introduce a way of single-class sensitivity analysis where the candidates are filtered according to their matching to prediction bias generated by targeted damaging of the activations. Relying on the observed parallels, we assume that the developed framework can potentially be transferred to studying biological neural networks in complex environments.
Kernel approximation with exponentials is useful in many problems with convolution quadrature and particle interactions such as integral-differential equations, molecular dynamics and machine learning. This paper proposes a weighted balanced truncation to construct an optimal model reduction method for compressing the number of exponentials in the sum-of-exponentials approximation of kernel functions. This method shows great promise in approximating long-range kernels, achieving over 4 digits of accuracy improvement for the Ewald-splitting and inverse power kernels in comparison with the classical balanced truncation. Numerical results demonstrate its excellent performance and attractive features for practical applications.
The high cost of sequential time integration is one major constraint that limits the speedup of a time-parallel algorithm like the Parareal algorithm due to the difficulty of coarsening time steps in a stiff numerical problem. To address this challenge, we develop a parallel-in-time approach based on the Parareal algorithm, in which we construct a novel coarse solver using a data-driven method based on Dynamic Mode Decomposition in place of a classic time marching scheme. The proposed solver computes an approximation of the solution using two numerical schemes of different accuracies in parallel, and apply High-Order Dynamic Mode Decomposition (HODMD) to reduce the cost of sequential computations. Compared to the original Parareal algorithm, the proposed approach allows for the construction of low-cost coarse solvers for many complicated stiff problems. We demonstrate through several numerical examples in fluid dynamics that the proposed method can effectively reduce the serial computation cost and improve the parallel speedup of long-time simulations which are hard to accelerate using the original Parareal algorithm.
Objective: ascending aortic aneurysm growth prediction is still challenging in clinics. In this study, we evaluate and compare the ability of local and global shape features to predict ascending aortic aneurysm growth. Material and methods: 70 patients with aneurysm, for which two 3D acquisitions were available, are included. Following segmentation, three local shape features are computed: (1) the ratio between maximum diameter and length of the ascending aorta centerline, (2) the ratio between the length of external and internal lines on the ascending aorta and (3) the tortuosity of the ascending tract. By exploiting longitudinal data, the aneurysm growth rate is derived. Using radial basis function mesh morphing, iso-topological surface meshes are created. Statistical shape analysis is performed through unsupervised principal component analysis (PCA) and supervised partial least squares (PLS). Two types of global shape features are identified: three PCA-derived and three PLS-based shape modes. Three regression models are set for growth prediction: two based on gaussian support vector machine using local and PCA-derived global shape features; the third is a PLS linear regression model based on the related global shape features. The prediction results are assessed and the aortic shapes most prone to growth are identified. Results: the prediction root mean square error from leave-one-out cross-validation is: 0.112 mm/month, 0.083 mm/month and 0.066 mm/month for local, PCA-based and PLS-derived shape features, respectively. Aneurysms close to the root with a large initial diameter report faster growth. Conclusion: global shape features might provide an important contribution for predicting the aneurysm growth.
Hybrid quantum systems are fundamental in quantum materials and quantum information science. In this work, we demonstrate that Hamiltonian learning in hybrid spin-boson systems can achieve the Heisenberg limit. Specifically, we establish a rigorous theoretical framework proving that, given access to an unknown hybrid Hamiltonian system, our algorithm can estimate the Hamiltonian coupling parameters up to root mean square error (RMSE) $\epsilon$ with a total evolution time scaling as $T \sim \mathcal{O}(\epsilon^{-1})$ using only $\mathcal{O}({\rm polylog}(\epsilon^{-1}))$ measurements. Furthermore, it remains robust against small state preparation and measurement (SPAM) errors. In addition, we also provide an alternative algorithm based on distributed quantum sensing, which significantly reduces the maximum evolution time per measurement.To validate our method, we apply it to the generalized Dicke model for Hamiltonian learning and the spin-boson model for spectrum learning, demonstrating its efficiency in practical quantum systems. These results provide a scalable and robust framework for precision quantum sensing and Hamiltonian characterization in hybrid quantum platforms.