https://papers.cool/arxiv/math.NANumerical Analysis2024-08-09T00:00:00+00:00python-feedgenCool Papers - Immersive Paper Discoveryhttps://papers.cool/arxiv/2408.04044Asymptotically optimal $t$-design curves on $S^3$2024-08-09T00:00:00+00:00Ayodeji LindbladA $\textit{spherical $t$-design curve}$ was defined by Ehler and Gr\"{o}chenig to be a continuous, piecewise smooth, closed curve on the sphere with finitely many self-intersections whose associated line integral applied to any polynomial of degree $t$ or less evaluates to the average of this polynomial on the sphere. These authors posed the problem of proving that there exist sequences $\{\gamma_t\}_{t\in\Bbb N}$ of $t$-design curves on $S^d$ of asymptotically optimal length $\ell(\gamma_t)=\Theta(t^{d-1})$ as $t\to\infty$ and solved this problem for $d=2$. This work solves the problem for $d=3$ by proving existence of a constant $\mathcal C>0$ such that for any $C>\mathcal C$ and $t\in\Bbb N_+$, there exists a $t$-design curve with no self-intersections on $S^3$ of length $Ct^2$.https://papers.cool/arxiv/2408.04208Visibility Analysis of the Sun as Viewed from Multiple Spacecraft at the Sun-Earth Lagrange Points2024-08-09T00:00:00+00:00Jinsung LeeSung-Hong ParkArik PosnerKyung-Suk ChoJaemyung AhnBeyond the Sun-Earth line, spacecraft equipped with various solar telescopes are intended to be deployed at several different vantage points in the heliosphere to carry out coordinated, multi-view observations of the Sun and its dynamic activities. In this context, we investigate solar visibility by imaging instruments onboard the spacecraft orbiting the Sun-Earth Lagrange points L1, L4 and L5, respectively. An optimal arrival time for vertical periodic orbits stationed at L4 and L5 is determined based on geometric considerations that ensure maximum visibility of solar poles or higher latitudes per year. For a different set of orbits around the three Lagrange points (L1, L4 and L5), we calculate the visibility of the solar surface (i.e., observation days per year) as a function of the solar latitude. We also analyze where the solar limb viewed from one of the three Sun-Earth Lagrange points under consideration is projected onto the solar surface visible to the other two. This analysis particularly aims at determining the feasibility of studying solar eruptions, such as flares and coronal mass ejections, with coordinated observations of off-limb erupting coronal structures and their on-disk magnetic footpoints. In addition, visibility analysis of a feature (such as sunspots) on the solar surface is made for multiple spacecraft in various types of orbits with different inclinations to quantify the improvement in continuous tracking of the target feature for studying its long-term evolution from emergence, growth and to decay. A comprehensive comparison of observations from single (L1), double (L1 and L4) and multi-space missions (L1, L4 and L5) is carried out through our solar visibility analysis, and this may help us to design future space missions of constructing multiple solar observatories at the Sun-Earth Lagrange points.https://papers.cool/arxiv/2408.04009Error Bounds for Open Quantum Systems with Harmonic Bosonic Bath2024-08-09T00:00:00+00:00Kaizhao LiuJianfeng LuWe investigate the dependence of physical observable of open quantum systems with Bosonic bath on the bath correlation function. We provide an error estimate of the difference of physical observable induced by the variation of bath correlation function, based on diagrammatic and combinatorial arguments. This gives a mathematically rigorous justification of the result in [Mascherpa et al, Phys Rev Lett 2017].https://papers.cool/arxiv/2408.04080A hybrid interpolation ACA accelerated method for parabolic boundary integral operators2024-08-09T00:00:00+00:00Sivaram AmbikasaranRitesh KhanJohannes TauschSihao WangWe consider piecewise polynomial discontinuous Galerkin discretizations of boundary integral reformulations of the heat equation. The resulting linear systems are dense and block-lower triangular and hence can be solved by block forward elimination. For the fast evaluation of the history part, the matrix is subdivided into a family of sub-matrices according to the temporal separation. Separated blocks are approximated by Chebyshev interpolation of the heat kernel in time. For the spatial variable, we propose an adaptive cross approximation (ACA) framework to obtain a data-sparse approximation of the entire matrix. We analyse how the ACA tolerance must be adjusted to the temporal separation and present numerical results for a benchmark problem to confirm the theoretical estimates.https://papers.cool/arxiv/2408.04128Splitting Methods for Computing Matrix Functions (Elements) Based on Non-Zero Diagonals Positions2024-08-09T00:00:00+00:00Majed HamadiNezam Mahdavi-AmiriIn applications, we often need to compute functions of matrices, such as banded matrices, the Kronecker sum of banded matrices, Toeplitz matrices, and many other types, which all share the common feature that their non-zero elements are concentrated around certain diagonals. We approximate matrix functions by considering the positions of non-zero diagonals in the original matrix. Focusing on non-zero diagonals provides us with simple algorithms to be used as tools to reduce complexity of other algorithms for computing matrix functions. Here, we first establish a decay bound for elements of matrix functions using the non-zero diagonals. Then, we develop methods that involve dividing the problem of computing matrix functions into functions of some submatrices of the original matrix. The size of these submatrices depends on the positions and number of non-zero diagonals in a monomial of the original matrix, ranging from degree zero to a given degree. For Toeplitz matrices, we demonstrate that our method turns to a simpler algorithm and works more efficiently. The convergence analysis of our proposed methods is conducted by establishing connections to the best polynomial approximation. When only specific elements or the trace of matrix functions are required, we derive submatrices from the original matrix based solely on the indices of the elements of interest. Additionally, for the special case of banded-symmetric Toeplitz matrices, we derive an approximation for elements of matrix functions with geometrically reducing error, using closed-form formulas that depend solely on the indices of the elements.https://papers.cool/arxiv/2408.04157R-adaptive DeepONet: Learning Solution Operators for PDEs with Discontinuous Solutions Using an R-adaptive Strategy2024-08-09T00:00:00+00:00Yameng ZhuJingrun ChenWeibing DengDeepONet has recently been proposed as a representative framework for learning nonlinear mappings between function spaces. However, when it comes to approximating solution operators of partial differential equations (PDEs) with discontinuous solutions, DeepONet poses a foundational approximation lower bound due to its linear reconstruction property. Inspired by the moving mesh (R-adaptive) method, we propose an R-adaptive DeepONet method, which contains the following components: (1) the output data representation is transformed from the physical domain to the computational domain using the equidistribution principle; (2) the maps from input parameters to the solution and the coordinate transformation function over the computational domain are learned using DeepONets separately; (3) the solution over the physical domain is obtained via post-processing methods such as the (linear) interpolation method. Additionally, we introduce a solution-dependent weighting strategy in the training process to reduce the final error. We establish an upper bound for the reconstruction error based on piecewise linear interpolation and show that the introduced R-adaptive DeepONet can reduce this bound. Moreover, for two prototypical PDEs with sharp gradients or discontinuities, we prove that the approximation error decays at a superlinear rate with respect to the trunk basis size, unlike the linear decay observed in vanilla DeepONets. Therefore, the R-adaptive DeepONet overcomes the limitations of DeepONet, and can reduce the approximation error for problems with discontinuous solutions. Numerical experiments on PDEs with discontinuous solutions, including the linear advection equation, the Burgers' equation with low viscosity, and the compressible Euler equations of gas dynamics, are conducted to verify the advantages of the R-adaptive DeepONet over available variants of DeepONet.https://papers.cool/arxiv/2408.04209Solving high-dimensional Hamilton-Jacobi-Bellman equation with functional hierarchical tensor2024-08-09T00:00:00+00:00Xun TangNan ShengLexing YingThis work proposes a novel numerical scheme for solving the high-dimensional Hamilton-Jacobi-Bellman equation with a functional hierarchical tensor ansatz. We consider the setting of stochastic control, whereby one applies control to a particle under Brownian motion. In particular, the existence of diffusion presents a new challenge to conventional tensor network methods for deterministic optimal control. To overcome the difficulty, we use a general regression-based formulation where the loss term is the Bellman consistency error combined with a Sobolev-type penalization term. We propose two novel sketching-based subroutines for obtaining the tensor-network approximation to the action-value functions and the value functions, which greatly accelerate the convergence for the subsequent regression phase. We apply the proposed approach successfully to two challenging control problems with Ginzburg-Landau potential in 1D and 2D with 64 variables.https://papers.cool/arxiv/2408.04316Explicit expression and fast algorithms for the inverse of some matrices arising from implicit time integration2024-08-09T00:00:00+00:00Li ShishunWei HuileIn this paper, we first present an explicit expression for the inverse\emph{} of a type of matrices. As special applications, the inverse of some matrices arising from implicit time integration techniques, such as the well-known implicit Runge-Kutta schemes and block implicit methods, can also be explicitly determined. Adiitionally, we introduce three fast algorithms for computing the elements of the inverse of these matrices in $O(n^2)$ arithmetic operations, i.e., the first one is based on Traub algorithm for fast inversion of Vandermonde matrices, while the other two utilize the special structure of the matrices. Finally, some symbolic and numerical results are presented to show that our algorithms are both highly efficient and accurate.https://papers.cool/arxiv/2408.04351Sensitivity analysis of fractional linear systems based on random walks with negligible memory usage2024-08-09T00:00:00+00:00Andrés CentenoJuan A. AcebrónJosé MonteiroA random walk-based method is proposed to efficiently compute the solution of a large class of fractional in time linear systems of differential equations (linear F-ODE systems), along with the derivatives with respect to the system parameters. Such a method is unbiased and unconditionally stable, and can therefore be used to provide an unbiased estimation of individual entries of the solution, or the full solution. By using stochastic differentiation techniques, it can be used as well to provide unbiased estimators of the sensitivities of the solution with respect to the problem parameters without any additional computational cost. The time complexity of the algorithm is discussed here, along with suitable variance bounds, which prove in practice the convergence of the algorithm. Finally, several test cases were run to assess the validity of the algorithm.https://papers.cool/arxiv/2408.04371Efficient finite element schemes for a phase field model of two-phase incompressible flows with different densities2024-08-09T00:00:00+00:00Jiancheng WangMaojun LiCheng WangIn this paper, we present two multiple scalar auxiliary variable (MSAV)-based, finite element numerical schemes for the Abels-Garcke-Gr{\"u}n (AGG) model, which is a thermodynamically consistent phase field model of two-phase incompressible flows with different densities. Both schemes are decoupled, linear, second-order in time, and the numerical implementation turns out to be straightforward. The first scheme solves the Navier-Stokes equations in a saddle point formulation, while the second one employs the artificial compressibility method, leading to a fully decoupled structure with a time-independent pressure update equation. In terms of computational cost, only a sequence of independent elliptic or saddle point systems needs to be solved at each time step. At a theoretical level, the unique solvability and unconditional energy stability (with respect to a modified energy functional) of the proposed schemes are established. In addition, comprehensive numerical simulations are performed to verify the effectiveness and robustness of the proposed schemes.https://papers.cool/arxiv/2408.04372A Space-Time Multigrid Method for Space-Time Finite Element Discretizations of Parabolic and Hyperbolic PDEs2024-08-09T00:00:00+00:00Nils MargenbergPeter MunchWe present a space-time multigrid method based on tensor-product space-time finite element discretizations. The method is facilitated by the matrix-free capabilities of the {\ttfamily deal.II} library. It addresses both high-order continuous and discontinuous variational time discretizations with spatial finite element discretizations. The effectiveness of multigrid methods in large-scale stationary problems is well established. However, their application in the space-time context poses significant challenges, mainly due to the construction of suitable smoothers. To address these challenges, we develop a space-time cell-wise additive Schwarz smoother and demonstrate its effectiveness on the heat and acoustic wave equations. The matrix-free framework of the {\ttfamily deal.II} library supports various multigrid strategies, including $h$-, $p$-, and $hp$-refinement across spatial and temporal dimensions. Extensive empirical evidence, provided through scaling and convergence tests on high-performance computing platforms, demonstrate high performance on perturbed meshes and problems with heterogeneous and discontinuous coefficients. Throughputs of over a billion degrees of freedom per second are achieved on problems with more than a trillion global degrees of freedom. The results prove that the space-time multigrid method can effectively solve complex problems in high-fidelity simulations and show great potential for use in coupled problems.https://papers.cool/arxiv/2408.04503Row-aware Randomized SVD with applications2024-08-09T00:00:00+00:00Davide PalittaSascha PortaroThe randomized singular value decomposition proposed in [12] has certainly become one of the most well-established randomization-based algorithms in numerical linear algebra. The key ingredient of the entire procedure is the computation of a subspace which is close to the column space of the target matrix $\mathbf{A}$ up to a certain probabilistic confidence. In this paper we propose a modification to the standard randomized SVD procedure which leads, in general, to better approximations to $\text{Range}(\mathbf{A})$ at the same computational cost. To this end, we explicitly construct information from the row space of $\mathbf{A}$ enhancing the quality of our approximation. We also observe that very few pieces of information from $\text{Range}(\mathbf{A}^T)$ are indeed necessary. We thus design a variant of our algorithm equipped with a subsampling step which largely increases the efficiency of our procedure while attaining competitive accuracy records. Our findings are supported by both theoretical analysis and numerical results.https://papers.cool/arxiv/2408.04507Sharp error bounds for edge-element discretisations of the high-frequency Maxwell equations2024-08-09T00:00:00+00:00Théophile Chaumont-FreletJeffrey GalkowskiEuan A. SpenceWe prove sharp wavenumber-explicit error bounds for first- or second-type-N\'ed\'elec-element (a.k.a. edge-element) conforming discretisations, of arbitrary (fixed) order, of the variable-coefficient time-harmonic Maxwell equations posed in a bounded domain with perfect electric conductor (PEC) boundary conditions. The PDE coefficients are allowed to be piecewise regular and complex-valued; this set-up therefore includes scattering from a PEC obstacle and/or variable real-valued coefficients, with the radiation condition approximated by a perfectly matched layer (PML). In the analysis of the $h$-version of the finite-element method, with fixed polynomial degree $p$, applied to the time-harmonic Maxwell equations, the $\textit{asymptotic regime}$ is when the meshwidth, $h$, is small enough (in a wavenumber-dependent way) that the Galerkin solution is quasioptimal independently of the wavenumber, while the $\textit{preasymptotic regime}$ is the complement of the asymptotic regime. The results of this paper are the first preasymptotic error bounds for the time-harmonic Maxwell equations using first-type N\'ed\'elec elements or higher-than-lowest-order second-type N\'ed\'elec elements. Furthermore, they are the first wavenumber-explicit results, even in the asymptotic regime, for Maxwell scattering problems with a non-empty scatterer.