{"title":"Efficient algorithms for Tucker decomposition via approximate matrix multiplication","authors":"Maolin Che, Yimin Wei, Hong Yan","doi":"10.1007/s10444-025-10232-0","DOIUrl":"10.1007/s10444-025-10232-0","url":null,"abstract":"<div><p>This paper develops fast and efficient algorithms for computing Tucker decomposition with a given multilinear rank. By combining random projection and the power scheme, we propose two efficient randomized versions for the truncated high-order singular value decomposition (T-HOSVD) and the sequentially T-HOSVD (ST-HOSVD), which are two common algorithms for approximating Tucker decomposition. To reduce the complexities of these two algorithms, fast and efficient algorithms are designed by combining two algorithms and approximate matrix multiplication. The theoretical results are also achieved based on the bounds of singular values of standard Gaussian matrices and the theoretical results for approximate matrix multiplication. Finally, the efficiency of these algorithms is illustrated via some test tensors from synthetic and real datasets.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"51 3","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143856476","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Computing the action of the matrix generating function of Bernoulli polynomials on a vector with an application to non-local boundary value problems","authors":"Lidia Aceto, Luca Gemignani","doi":"10.1007/s10444-025-10231-1","DOIUrl":"10.1007/s10444-025-10231-1","url":null,"abstract":"<div><p>This paper deals with efficient numerical methods for computing the action of the matrix generating function of Bernoulli polynomials, say <span>(q(tau ,A))</span>, on a vector when <i>A</i> is a large and sparse matrix. This problem occurs when solving some non-local boundary value problems. Methods based on the Fourier expansion of <span>(q(tau ,w))</span> have already been addressed in the scientific literature. The contribution of this paper is twofold. First, we place these methods in the classical framework of Krylov-Lanczos (polynomial-rational) techniques for accelerating Fourier series. This allows us to apply the convergence results developed in this context to our function. Second, we design a new acceleration scheme. Some numerical results are presented to show the effectiveness of the proposed algorithms.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"51 2","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10444-025-10231-1.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143809364","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A discontinuous plane wave neural network method for Helmholtz equation and time-harmonic Maxwell’s equations","authors":"Long Yuan, Qiya Hu","doi":"10.1007/s10444-025-10229-9","DOIUrl":"10.1007/s10444-025-10229-9","url":null,"abstract":"<div><p>In this paper, we propose a <i>discontinuous</i> plane wave neural network (DPWNN) method with <span>(hp-)</span>refinement for approximately solving Helmholtz equation and time-harmonic Maxwell equations. In this method, we define a quadratic functional as in the plane wave least square (PWLS) method with <span>(h-)</span>refinement and introduce new discretization sets spanned by element-wise neural network functions with a single hidden layer, where the activation function on each element is chosen as a complex-valued exponential function like the plane wave function. The desired approximate solution is recursively generated by iteratively solving a quasi-minimization problem associated with the functional and the sets described above, which is defined by a sequence of approximate minimizers of the underlying residual functionals, where plane wave direction angles and activation coefficients are alternatively computed by iterative algorithms. For the proposed DPWNN method, the plane wave directions are adaptively determined in the iterative process, which is different from that in the standard PWLS method (where the plane wave directions are preliminarily given). Numerical experiments will confirm that this DPWNN method can generate approximate solutions with higher accuracy than the PWLS method.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"51 2","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-04-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143786662","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Low-rank exponential integrators for stiff differential Riccati equations","authors":"Hao Chen, Alfio Borzì","doi":"10.1007/s10444-025-10228-w","DOIUrl":"10.1007/s10444-025-10228-w","url":null,"abstract":"<div><p>Exponential integrators are an efficient alternative to implicit schemes for the time integration of stiff system of differential equations. In this paper, low-rank exponential integrators of orders one and two for stiff differential Riccati equations are proposed and investigated. The error estimates of the proposed schemes are established. The proposed approach allows to overcome the main difficulties that lay in the interplay of time integration and low-rank approximation in the numerical schemes, which is uncommon in standard discretization of differential equations. Results of numerical experiments demonstrate the validity of the convergence analysis and show the performance of the proposed low-rank approximations with different settings.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"51 2","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143749225","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A quasi-boundary-value method for solving a nonlinear space-fractional backward diffusion problem","authors":"Xiaoli Feng, Xiaoyu Yuan, Yun Zhang","doi":"10.1007/s10444-025-10230-2","DOIUrl":"10.1007/s10444-025-10230-2","url":null,"abstract":"<div><p>In this paper, we adopt a quasi-boundary-value method to solve the nonlinear space-fractional backward problem with perturbed both final value and variable diffusion coefficient in general dimensional space, which is a severely ill-posed problem. The existence, uniqueness and stability of the solution for the quasi-boundary-value problem are proved. Convergence estimates are presented under an <i>a-priori</i> bound assumption of the exact solution. Finally, several numerical examples are given by the finite difference scheme and the fixed-point iteration method to show the effectiveness of the theoretical results.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"51 2","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-03-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143736958","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Strong convergence of a fully discrete scheme for stochastic Burgers equation with fractional-type noise","authors":"Yibo Wang, Wanrong Cao","doi":"10.1007/s10444-025-10227-x","DOIUrl":"10.1007/s10444-025-10227-x","url":null,"abstract":"<div><p>We investigate numerical approximations for the stochastic Burgers equation driven by an additive cylindrical fractional Brownian motion with Hurst parameter <span>(H in (frac{1}{2}, 1))</span>. To discretize the continuous problem in space, a spectral Galerkin method is employed, followed by the presentation of a nonlinear-tamed accelerated exponential Euler method to yield a fully discrete scheme. By showing the exponential integrability of the stochastic convolution of the fractional Brownian motion, we present the boundedness of moments of semidiscrete and full-discrete approximations. Building upon these results and the convergence of the fully discrete scheme in probability proved by a stopping time technique, we derive the strong convergence of the proposed scheme.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"51 2","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143676326","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Product kernels are efficient and flexible tools for high-dimensional scattered data interpolation","authors":"Kristof Albrecht, Juliane Entzian, Armin Iske","doi":"10.1007/s10444-025-10226-y","DOIUrl":"10.1007/s10444-025-10226-y","url":null,"abstract":"<div><p>This work concerns the construction and characterization of product kernels for multivariate approximation from a finite set of discrete samples. To this end, we consider composing different component kernels, each acting on a low-dimensional Euclidean space. Due to Aronszajn (Trans. Am. Math. Soc. <b>68</b>, 337–404 1950), the product of positive <i>semi-</i>definite kernel functions is again positive <i>semi-</i>definite, where, moreover, the corresponding native space is a particular instance of a tensor product, referred to as Hilbert tensor product. We first analyze the general problem of multivariate interpolation by product kernels. Then, we further investigate the tensor product structure, in particular for <i>grid-like</i> samples. We use this case to show that the product of positive definite kernel functions is again positive definite. Moreover, we develop an efficient computation scheme for the well-known Newton basis. Supporting numerical examples show the good performance of product kernels, especially for their flexibility.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"51 2","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10444-025-10226-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143655244","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The Kolmogorov N-width for linear transport: exact representation and the influence of the data","authors":"Florian Arbes, Constantin Greif, Karsten Urban","doi":"10.1007/s10444-025-10224-0","DOIUrl":"10.1007/s10444-025-10224-0","url":null,"abstract":"<div><p>The Kolmogorov <i>N</i>-width describes the best possible error one can achieve by elements of an <i>N</i>-dimensional linear space. Its decay has extensively been studied in approximation theory and for the solution of partial differential equations (PDEs). Particular interest has occurred within model order reduction (MOR) of parameterized PDEs, e.g., by the reduced basis method (RBM). While it is known that the <i>N</i>-width decays exponentially fast (and thus admits efficient MOR) for certain problems, there are examples of the linear transport and the wave equation, where the decay rate deteriorates to <span>(N^{-1/2})</span>. On the other hand, it is widely accepted that a smooth parameter dependence admits a fast decay of the <i>N</i>-width. However, a detailed analysis of the influence of properties of the data (such as regularity or slope) on the rate of the <i>N</i>-width seems to be lacking. In this paper, we state that the optimal linear space is a direct sum of shift-isometric eigenspaces corresponding to the largest eigenvalues, yielding an exact representation of the <i>N</i>-width as their sum. For the linear transport problem, which is modeled by half-wave symmetric initial and boundary conditions <i>g</i>, we obtain such an optimal decomposition by sorted trigonometric functions with eigenvalues that match the Fourier coefficients of <i>g</i>. Further, for normalized <i>g</i> in the Sobolev space <span>(H^r)</span> of broken order <span>(r>0)</span>, the sorted eigenfunctions give the sharp upper bound of the <i>N</i>-width, which is a reciprocal of a certain power sum. Yet, for ease, we also provide the decay <span>((pi N)^{-r})</span>, obtained by the non-optimal space of ordering the trigonometric functions by frequency rather than by eigenvalue. Our theoretical investigations are complemented by numerical experiments which confirm the sharpness of our bounds and give additional quantitative insight.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"51 2","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-03-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10444-025-10224-0.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143546172","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"On the recovery of two function-valued coefficients in the Helmholtz equation for inverse scattering problems via neural networks","authors":"Zehui Zhou","doi":"10.1007/s10444-025-10225-z","DOIUrl":"10.1007/s10444-025-10225-z","url":null,"abstract":"<div><p>Recently, deep neural networks (DNNs) have become powerful tools for solving inverse scattering problems. However, the approximation and generalization rates of DNNs for solving these problems remain largely under-explored. In this work, we introduce two types of combined DNNs (uncompressed and compressed) to reconstruct two function-valued coefficients in the Helmholtz equation for inverse scattering problems from the scattering data at two different frequencies. An analysis of the approximation and generalization capabilities of the proposed neural networks for simulating the regularized pseudo-inverses of the linearized forward operators in direct scattering problems is provided. The results show that, with sufficient training data and parameters, the proposed neural networks can effectively approximate the inverse process with desirable generalization. Preliminary numerical results show the feasibility of the proposed neural networks for recovering two types of isotropic inhomogeneous media. Furthermore, the trained neural network is capable of reconstructing the isotropic representation of certain types of anisotropic media.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"51 1","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-02-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10444-025-10225-z.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143379751","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"On a non-uniform (alpha )-robust IMEX-L1 mixed FEM for time-fractional PIDEs","authors":"Lok Pati Tripathi, Aditi Tomar, Amiya K. Pani","doi":"10.1007/s10444-025-10221-3","DOIUrl":"10.1007/s10444-025-10221-3","url":null,"abstract":"<div><p>A non-uniform implicit-explicit L1 mixed finite element method (IMEX-L1-MFEM) is investigated for a class of time-fractional partial integro-differential equations (PIDEs) with space-time-dependent coefficients and non-self-adjoint elliptic part. The proposed fully discrete method combines an IMEX-L1 method on a graded mesh in the temporal variable with a mixed finite element method in spatial variables. The focus of the study is to analyze stability results and to establish optimal error estimates, up to a logarithmic factor, for both the solution and the flux in <span>(L^2)</span>-norm when the initial data <span>(u_0in H_0^1(Omega )cap H^2(Omega ))</span>. Additionally, an error estimate in <span>(L^infty )</span>-norm is derived for 2D problems. All the derived estimates and bounds in this article remain valid as <span>(alpha rightarrow 1^{-})</span>, where <span>(alpha )</span> is the order of the Caputo fractional derivative. Finally, the results of several numerical experiments conducted at the end of this paper are confirming our theoretical findings.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"51 1","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-02-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143373297","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}