{"title":"Ellipsoidal density-equalizing map for genus-0 closed surfaces","authors":"Zhiyuan Lyu, Lok Ming Lui, Gary P. T. Choi","doi":"10.1007/s10444-026-10304-9","DOIUrl":"10.1007/s10444-026-10304-9","url":null,"abstract":"<div><p>Surface parameterization is a fundamental task in geometry processing and plays an important role in many science and engineering applications. In recent years, the density-equalizing map, a shape deformation technique based on the physical principle of density diffusion, has been utilized for the parameterization of simply connected and multiply connected open surfaces. More recently, a spherical density-equalizing mapping method has been developed for the parameterization of genus-0 closed surfaces. However, for genus-0 closed surfaces with extreme geometry, using a spherical domain for the parameterization may induce large geometric distortion. In this work, we develop a novel method for computing density-equalizing maps of genus-0 closed surfaces onto an ellipsoidal domain. This allows us to achieve ellipsoidal area-preserving parameterizations and ellipsoidal parameterizations with controlled area change. We further propose an energy minimization approach that combines density-equalizing maps and quasi-conformal maps, which allows us to produce ellipsoidal density-equalizing quasi-conformal maps for achieving a balance between density equalization and quasi-conformality. Using our proposed methods, we can significantly improve the performance of surface remeshing for genus-0 closed surfaces. Experimental results on a large variety of genus-0 closed surfaces are presented to demonstrate the effectiveness of our proposed methods.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"52 3","pages":""},"PeriodicalIF":2.1,"publicationDate":"2026-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10444-026-10304-9.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147686675","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}
Nicola Rares Franco, Andrea Manzoni, Paolo Zunino, Jan S. Hesthaven
{"title":"Deep orthogonal decomposition: a continuously adaptive neural network approach to model order reduction of parametrized partial differential equations","authors":"Nicola Rares Franco, Andrea Manzoni, Paolo Zunino, Jan S. Hesthaven","doi":"10.1007/s10444-026-10295-7","DOIUrl":"10.1007/s10444-026-10295-7","url":null,"abstract":"<div><p>We develop a novel deep learning technique, termed deep orthogonal decomposition (DOD), for dimensionality reduction and reduced order modeling of parameter-dependent partial differential equations. The approach involves constructing a deep neural network model that approximates the solution manifold using a continuously adaptive local basis. In contrast to global methods, such as proper orthogonal decomposition (POD), this adaptivity allows the DOD to mitigate the Kolmogorov barrier when dealing with space-interacting parameters, making the approach applicable to a wide spectrum of parametric problems. Leveraging this idea, we use the DOD to construct an adaptive alternative to the so-called POD-NN method, here termed DOD-NN. The approach is fully data-driven and non-intrusive but, at the same time, allows for a tight control on error propagation and remains highly interpretable thanks to the rich structure present in the latent space. For this reason, the proposed approach stands out as a valuable alternative to other nonlinear model order reduction techniques, such as those based on deep autoencoders. The methodology is discussed both theoretically and practically, evaluating its performances on problems involving nonlinear PDEs, parametrized geometries, and high-dimensional parameter spaces. Finally, we conclude with a brief discussion on potential applications of the DOD beyond DOD-NN, featuring, for instance, the integration of our approach within intrusive reduced order models such as the reduced basis method.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"52 3","pages":""},"PeriodicalIF":2.1,"publicationDate":"2026-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10444-026-10295-7.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147686690","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":"Convergence analysis of time-filtered backward differentiation methods up to fourth-order for the molecular beam epitaxial model without slope selection","authors":"Jiexin Wang, Hong-lin Liao","doi":"10.1007/s10444-026-10302-x","DOIUrl":"10.1007/s10444-026-10302-x","url":null,"abstract":"<div><p>This work presents the convergence of time-filtered backward differentiation formula (BDF) methods up to fourth-order time accuracy for the molecular beam epitaxial (MBE) equation without slope selection. With the help of the discrete orthogonal convolution kernels, some concise <span>(L^2)</span> norm error estimates are established via the discrete energy technique. To the best of our knowledge, this is the first time such type <span>(L^2)</span> norm error estimates of time-filtered BDF2 and time-filtered BDF3 schemes are obtained for a nonlinear parabolic equation. Numerical examples are included to show the effectiveness of the proposed methods.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"52 2","pages":""},"PeriodicalIF":2.1,"publicationDate":"2026-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147642115","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":"On randomized explicit block Kaczmarz method for solving large linear systems","authors":"Cun-Qiang Miao, Xue-Yuan Tan","doi":"10.1007/s10444-026-10301-y","DOIUrl":"10.1007/s10444-026-10301-y","url":null,"abstract":"<div><p>The randomized block Kaczmarz method (Linear Algebra Appl., 441: 199-221, 2014) proposed by Needell and Tropp is efficient for solving large consistent linear systems. However, each iteration of the randomized block Kaczmarz method carries a high cost, as it calls for the computation of a pseudoinverse, or equivalently, the solution of a least-squares problem. In this paper, we propose a randomized explicit block Kaczmarz method that avoids the direct computation of pseudoinverses by exploiting the structure of the block updates. This explicit formulation allows for a more flexible selection of the rows in the working block at each iteration, without relying on a predefined partition of the row indices. Based on a randomized block selection strategy, we further establish the convergence properties of the proposed method. It indicates that the randomized explicit block Kaczmarz method exhibits faster convergence compared to the multi-step standard randomized Kaczmarz method and the randomized block Kaczmarz method. Finally, numerical experiments are carried out to show great superiority and robustness over some state-of-the-art randomized block Kaczmarz methods.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"52 2","pages":""},"PeriodicalIF":2.1,"publicationDate":"2026-04-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147619931","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":"Spatially semi-discretely uniformly exponential stability approximation for one-dimensional thermoelastic system","authors":"Jiankang Liu, Bao-Zhu Guo","doi":"10.1007/s10444-026-10298-4","DOIUrl":"10.1007/s10444-026-10298-4","url":null,"abstract":"<div><p>This paper investigates the uniformly exponential stability of a spatially semi-discrete finite difference scheme for a one-dimensional thermoelastic system. The frequency domain approach is employed due to the difficulty of constructing a suitable time-domain multiplier for the continuous system. A novel semi-discrete finite difference scheme is developed using the order-reduction method. The uniform exponential stability of the resulting semi-discrete system is then established through a discrete counterpart of the frequency domain method.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"52 2","pages":""},"PeriodicalIF":2.1,"publicationDate":"2026-04-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147620161","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":"Analyzing the function (zuparrow uparrow !a) for fractional a","authors":"William Paulsen","doi":"10.1007/s10444-026-10300-z","DOIUrl":"10.1007/s10444-026-10300-z","url":null,"abstract":"<div><p>In this paper, we will extend Knuth’s up arrow notation <span>(b uparrow uparrow n)</span>, which for integer <i>n</i> is <span>(b^{b^{cdot ^{cdot ^{b}}}})</span> with <i>n</i> <i>b</i>’s, to allow <i>n</i> to be a fraction, or even complex. We do this by first considering the unique complex tetration <span>(kappa _b(z))</span>, which satisfies the equations <span>(kappa _b(z+1) = b^{kappa _b(z)})</span> with <span>(kappa _b(0)=1)</span>, and then fix the value of <i>z</i> as we vary the base <i>b</i>. In particular, we will consider the “half-iterate,” which is the case with <span>(z=1/2)</span>, and analytically continue the function into the complex <i>b</i>-plane. To consider “third-iterates” and “fourth-iterates,” we also use 11 other values of <i>z</i>, namely 1/12, 1/6, 1/4, 1/3, 5/12, 7/12, 2/3, 3/4, 5/6, 11/12, and finally <span>(z=i)</span>. We make several discoveries by analyzing these functions. First of all, we discover that there is a branch cut at <span>(b=e^{1/e})</span>, and find that as we approach this point, the functions approach one of the known super-exponentials for the base <span>(e^{1/e})</span>. Also, we find that <span>(b=1)</span> is also a branch point, and discover what happens as the base approaches 1.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"52 2","pages":""},"PeriodicalIF":2.1,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147586399","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":"GP-CMRH: an inner product free iterative method for block two-by-two nonsymmetric linear systems","authors":"Kui Du, Jia-Jun Fan","doi":"10.1007/s10444-026-10299-3","DOIUrl":"10.1007/s10444-026-10299-3","url":null,"abstract":"<div><p>We propose an inner product free iterative method called GP-CMRH for solving block two-by-two nonsymmetric linear systems. GP-CMRH relies on a new simultaneous Hessenberg process that reduces two rectangular matrices to upper Hessenberg form simultaneously, without employing inner products. Compared with GPMR [SIAM J. Matrix Anal. Appl., 44 (2023), pp. 293–311], GP-CMRH requires less computational cost per iteration and may be more suitable for high-performance computing and low or mixed precision arithmetic due to its inner product free property. Our numerical experiments demonstrate that GP-CMRH and GPMR exhibit comparable convergence behavior (with GP-CMRH requiring slightly more iterations), yet GP-CMRH consumes less computational time in most cases. GP-CMRH significantly outperforms GMRES and CMRH in terms of convergence rate and runtime efficiency.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"52 2","pages":""},"PeriodicalIF":2.1,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147586398","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":"Smoothed distance kernels for MMDs and applications in Wasserstein gradient flows","authors":"Nicolaj Rux, Michael Quellmalz, Gabriele Steidl","doi":"10.1007/s10444-026-10289-5","DOIUrl":"10.1007/s10444-026-10289-5","url":null,"abstract":"<div><p>Negative distance kernels <span>(K(x,y) := - Vert x-yVert )</span> were used in the definition of maximum mean discrepancies (MMDs) in statistics and lead to favorable numerical results in various applications. In particular, so-called slicing techniques for handling high-dimensional kernel summations profit from the simple parameter-free structure of the distance kernel. However, due to its non-smoothness in <span>(x=y)</span>, most of the classical theoretical results, e.g., on Wasserstein gradient flows of the corresponding MMD functional, no longer hold true. In this paper, we propose a new kernel which keeps the favorable properties of the negative distance kernel as being conditionally positive definite of order one with a nearly linear increase towards infinity and a simple slicing structure, but is Lipschitz differentiable now. Our construction is based on a simple 1D smoothing procedure of the absolute value function followed by a Riemann–Liouville fractional integral transform. Numerical results demonstrate that the new kernel performs similarly well as the negative distance kernel in gradient descent methods, but now with theoretical guarantees.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"52 2","pages":""},"PeriodicalIF":2.1,"publicationDate":"2026-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10444-026-10289-5.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147536003","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":"Constrained rational fitting for (mathcal {D})-stable model order reduction","authors":"Tommaso Bradde, Stefano Grivet-Talocia","doi":"10.1007/s10444-026-10297-5","DOIUrl":"10.1007/s10444-026-10297-5","url":null,"abstract":"<div><p>A classical problem in data-driven model order reduction (MOR) of linear time-invariant (LTI) systems is the preservation of structural properties of the underlying large-scale dynamics. When dealing with MOR based on transfer function measurements, one relevant problem is how to force the reduced-order model (ROM) to inherit the asymptotic stability of the reference system, i.e., to enforce the poles of the ROM transfer function to have strictly negative real part. In this work, we tackle the more general problem of placing such poles in arbitrary linear matrix inequality (LMI) regions of the complex plane, which include a rich class of convex sets symmetric with respect to the real axis. LTI systems with poles constrained to this kind of regions are called <span>(mathcal {D})</span>-stable. Combining well-established results from control theory and recent developments in asymptotically stable rational approximation algorithms, we show that the problem can be solved efficiently via standard convex optimization routines. Several numerical testbenches of engineering interest confirm the effectiveness of the proposed methodology in practical applications.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"52 2","pages":""},"PeriodicalIF":2.1,"publicationDate":"2026-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10444-026-10297-5.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147507849","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":"Anderson acceleration of a Picard solver for the Oldroyd-B model of viscoelastic fluids","authors":"Duygu Vargun, Igor O. Monteiro, Leo G. Rebholz","doi":"10.1007/s10444-025-10270-8","DOIUrl":"10.1007/s10444-025-10270-8","url":null,"abstract":"<div><p>We study an iterative nonlinear solver for the Oldroyd-B system describing incompressible viscoelastic fluid flow. We establish a range of attributes of the fixed-point-based solver, together with the conditions under which it becomes contractive, and examine the smoothness properties of its corresponding fixed-point function. Under these properties, we demonstrate that the solver meets the necessary conditions for the recent Anderson acceleration (AA) framework, thereby showing that AA enhances the solver’s linear convergence rate. Results from three benchmark tests illustrate how AA improves the solver’s ability to converge as the Weissenberg number is increased.</p></div>","PeriodicalId":50869,"journal":{"name":"Advances in Computational Mathematics","volume":"52 2","pages":""},"PeriodicalIF":2.1,"publicationDate":"2026-03-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147441121","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}