Pub Date : 2025-11-01DOI: 10.1016/j.rinam.2025.100666
Qiang Niu , Mianmian Chen , Jinheng Wu
The Lanczos algorithm is a well-known three-term recurrence that can be used to generate an orthogonal basis for a Krylov subspace derived by a symmetric matrix. In the paper, we present a statistical interpretation of the entries of the tridiagonal matrix generated by the Lanczos process with a diagonal matrix and an initial vector . We show that the entries on the main diagonal line can be interpreted as weighted mean and the entries on the super-diagonal line can be understood as weighted sum of variance. Besides, a recurrence for producing the entries on the off-diagonal entries of the tridiagonal matrix is discovered, which leads to a new implementation of the Lanczos process. Finally, numerical examples are provided to investigate the preservation of orthogonality and efficiency in data fitting.
{"title":"Lanczos algorithm explained in statistics","authors":"Qiang Niu , Mianmian Chen , Jinheng Wu","doi":"10.1016/j.rinam.2025.100666","DOIUrl":"10.1016/j.rinam.2025.100666","url":null,"abstract":"<div><div>The Lanczos algorithm is a well-known three-term recurrence that can be used to generate an orthogonal basis for a Krylov subspace derived by a symmetric matrix. In the paper, we present a statistical interpretation of the entries of the tridiagonal matrix generated by the Lanczos process with a diagonal matrix <span><math><mi>X</mi></math></span> and an initial vector <span><math><mi>e</mi></math></span>. We show that the entries on the main diagonal line can be interpreted as <em>weighted mean</em> and the entries on the super-diagonal line can be understood as <em>weighted sum of variance</em>. Besides, a recurrence for producing the entries on the off-diagonal entries of the tridiagonal matrix is discovered, which leads to a new implementation of the Lanczos process. Finally, numerical examples are provided to investigate the preservation of orthogonality and efficiency in data fitting.</div></div>","PeriodicalId":36918,"journal":{"name":"Results in Applied Mathematics","volume":"28 ","pages":"Article 100666"},"PeriodicalIF":1.3,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145519627","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-01DOI: 10.1016/j.rinam.2025.100663
Zhikun Tian , Jianyun Wang , Zixin Zhong
In this paper, we investigate the two-step backward differentiation formula (BDF2) finite element method for a two-dimensional time-dependent Schrödinger equation. By applying the finite element method for space discretization and the BDF2 for time discretization, we derive a fully discrete scheme for the Schrödinger equation. The errors of the exact solution with the finite element solution are divided into temporal and spatial errors for separate analysis. We obtain the optimal error estimate in both space and time for the fully discrete scheme. Finally, a numerical experiment is performed to demonstrate the accuracy and efficiency of the proposed numerical scheme.
{"title":"Optimal error estimates of BDF2 finite element method for the two-dimensional time-dependent Schrödinger equation","authors":"Zhikun Tian , Jianyun Wang , Zixin Zhong","doi":"10.1016/j.rinam.2025.100663","DOIUrl":"10.1016/j.rinam.2025.100663","url":null,"abstract":"<div><div>In this paper, we investigate the two-step backward differentiation formula (BDF2) finite element method for a two-dimensional time-dependent Schrödinger equation. By applying the finite element method for space discretization and the BDF2 for time discretization, we derive a fully discrete scheme for the Schrödinger equation. The errors of the exact solution with the finite element solution are divided into temporal and spatial errors for separate analysis. We obtain the optimal error estimate in both space and time for the fully discrete scheme. Finally, a numerical experiment is performed to demonstrate the accuracy and efficiency of the proposed numerical scheme.</div></div>","PeriodicalId":36918,"journal":{"name":"Results in Applied Mathematics","volume":"28 ","pages":"Article 100663"},"PeriodicalIF":1.3,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145466098","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-01DOI: 10.1016/j.rinam.2025.100655
Ikha Magdalena , Anwar Efendi Nasution
Wave resonance in coastal basins can cause significant environmental and structural damage, highlighting the need for effective mitigation strategies in coastal engineering. This study explores resonance phenomena in a rectangular semi-closed basin protected by submerged triangular breakwaters. The wave dynamics are modeled utilizing modified Linear Shallow Water Equations (LSWEs), which incorporate a breakwater-induced friction factor. Analytical and numerical approaches are employed in determining the basin’s natural period—a key parameter governing the onset of resonance. Numerical simulations, formulated utilizing the finite volume method, are conducted to identify the conditions that trigger resonance. The findings reveal that submerged triangular breakwaters substantially affect the natural period and help attenuate resonance effects, providing valuable insights for the design of resilient coastal infrastructure.
{"title":"Analytical and numerical investigation of wave resonance over rectangular basin with submerged triangular breakwaters","authors":"Ikha Magdalena , Anwar Efendi Nasution","doi":"10.1016/j.rinam.2025.100655","DOIUrl":"10.1016/j.rinam.2025.100655","url":null,"abstract":"<div><div>Wave resonance in coastal basins can cause significant environmental and structural damage, highlighting the need for effective mitigation strategies in coastal engineering. This study explores resonance phenomena in a rectangular semi-closed basin protected by submerged triangular breakwaters. The wave dynamics are modeled utilizing modified Linear Shallow Water Equations (LSWEs), which incorporate a breakwater-induced friction factor. Analytical and numerical approaches are employed in determining the basin’s natural period—a key parameter governing the onset of resonance. Numerical simulations, formulated utilizing the finite volume method, are conducted to identify the conditions that trigger resonance. The findings reveal that submerged triangular breakwaters substantially affect the natural period and help attenuate resonance effects, providing valuable insights for the design of resilient coastal infrastructure.</div></div>","PeriodicalId":36918,"journal":{"name":"Results in Applied Mathematics","volume":"28 ","pages":"Article 100655"},"PeriodicalIF":1.3,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145417120","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-01DOI: 10.1016/j.rinam.2025.100672
Yanyan Song , Yanqiu Wang , Qiao Xin
In this paper, we construct two finite element discretizations on convex polygonal meshes for the Stokes equations, using the generalized barycentric coordinates (GBCs). Both constructions use a bubble-enhanced second-order Floater-Lai GBC to discretize the velocity field. The pressure field is discretized by discontinuous piecewise constants and discontinuous piecewise linears, respectively. They can be viewed as the generalization of the - and the - elements to convex polygonal meshes. We prove the inf–sup stability and the optimal convergence for both discretizations. Supporting numerical results are presented.
本文利用广义质心坐标在凸多边形网格上构造了Stokes方程的两个有限元离散化。两种结构都使用气泡增强的二阶float - lai GBC来离散速度场。压力场分别采用不连续分段常数和不连续分段线性进行离散。它们可以看作是P2-P0和Q2-P1元素对凸多边形网格的推广。证明了这两种离散化方法的稳定性和最优收敛性。给出了相应的数值结果。
{"title":"Second-order finite element discretization of Stokes equations on convex polygonal meshes","authors":"Yanyan Song , Yanqiu Wang , Qiao Xin","doi":"10.1016/j.rinam.2025.100672","DOIUrl":"10.1016/j.rinam.2025.100672","url":null,"abstract":"<div><div>In this paper, we construct two finite element discretizations on convex polygonal meshes for the Stokes equations, using the generalized barycentric coordinates (GBCs). Both constructions use a bubble-enhanced second-order Floater-Lai GBC to discretize the velocity field. The pressure field is discretized by discontinuous piecewise constants and discontinuous piecewise linears, respectively. They can be viewed as the generalization of the <span><math><msub><mrow><mi>P</mi></mrow><mrow><mn>2</mn></mrow></msub></math></span>-<span><math><msub><mrow><mi>P</mi></mrow><mrow><mn>0</mn></mrow></msub></math></span> and the <span><math><msub><mrow><mi>Q</mi></mrow><mrow><mn>2</mn></mrow></msub></math></span>-<span><math><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub></math></span> elements to convex polygonal meshes. We prove the inf–sup stability and the optimal convergence for both discretizations. Supporting numerical results are presented.</div></div>","PeriodicalId":36918,"journal":{"name":"Results in Applied Mathematics","volume":"28 ","pages":"Article 100672"},"PeriodicalIF":1.3,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145614876","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-01DOI: 10.1016/j.rinam.2025.100676
Li Wang , Jiawei Wang , Ren-Cang Li
This paper is concerned with sparse PCA via the matrix (2,1)-norm regularization (). It can produce a row-sparse projection, a useful tool in machine learning when it comes to, for example, feature selection, that aims to choose most relevant features. Mathematically, is a non-smooth optimization problem on the Stiefel manifold. For a suitably chosen regularization parameter, the optimal projection matrix has many negligible rows. A practical NEPv approach (nonlinear eigenvalue problem with eigenvector dependency) is proposed to iteratively compute the optimal projection matrix. It is shown that the approach is globally convergent in the sense that the objective is monotonically increasing during the iterative process and any accumulation point of the iterates is a stationary point to the optimization problem. Extensive numerical experiments, with an application to feature selection, have been conducted to demonstrate the performance of the practical NEPv approach, with comparison against existing feature selection methods in terms of classification accuracy. The numerical results demonstrate that is highly effective and often produces superior classification results to existing feature selection methods that are in use today.
{"title":"Sparse PCA via matrix (2,1)-norm regularization with an application to feature selection","authors":"Li Wang , Jiawei Wang , Ren-Cang Li","doi":"10.1016/j.rinam.2025.100676","DOIUrl":"10.1016/j.rinam.2025.100676","url":null,"abstract":"<div><div>This paper is concerned with sparse PCA via the matrix (2,1)-norm regularization (<span><math><msub><mrow><mo>PCA</mo></mrow><mrow><mn>2</mn><mo>,</mo><mn>1</mn></mrow></msub></math></span>). It can produce a row-sparse projection, a useful tool in machine learning when it comes to, for example, feature selection, that aims to choose most relevant features. Mathematically, <span><math><msub><mrow><mo>PCA</mo></mrow><mrow><mn>2</mn><mo>,</mo><mn>1</mn></mrow></msub></math></span> is a non-smooth optimization problem on the Stiefel manifold. For a suitably chosen regularization parameter, the optimal projection matrix has many negligible rows. A practical NEPv approach (nonlinear eigenvalue problem with eigenvector dependency) is proposed to iteratively compute the optimal projection matrix. It is shown that the approach is globally convergent in the sense that the objective is monotonically increasing during the iterative process and any accumulation point of the iterates is a stationary point to the optimization problem. Extensive numerical experiments, with an application to feature selection, have been conducted to demonstrate the performance of the practical NEPv approach, with comparison against existing feature selection methods in terms of classification accuracy. The numerical results demonstrate that <span><math><msub><mrow><mo>PCA</mo></mrow><mrow><mn>2</mn><mo>,</mo><mn>1</mn></mrow></msub></math></span> is highly effective and often produces superior classification results to existing feature selection methods that are in use today.</div></div>","PeriodicalId":36918,"journal":{"name":"Results in Applied Mathematics","volume":"28 ","pages":"Article 100676"},"PeriodicalIF":1.3,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145680919","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-01DOI: 10.1016/j.rinam.2025.100673
Elyas Shivanian , Ahmad Jafarabadi , Mousa J. Huntul
This study focuses on retrieving a time-dependent source term in the heat equation governed by two distinct nonlocal boundary conditions. The inverse problem is structured with an interior energy over-specification constraint. The proposed computational framework combines the partition of unity approach for spatial discretization with the finite difference scheme for temporal advancement. Through energy analysis, the semi-discrete time-stepping formulation is proven to be unconditionally stable and convergent at a rate of . Despite being linear and uniquely solvable, the problem is inherently ill-posed, as slight disturbances in input data can induce significant errors in the reconstructed solution. To counteract this instability, Tikhonov regularization is implemented, yielding a stable approximation even under noisy data conditions. Moreover, a novel parameter selection strategy for the regularization is introduced, which surpasses standard methods by delivering substantially improved results. Numerical simulations corroborate the scheme’s robustness, demonstrating its accuracy with noise-free inputs and its resilience when handling contaminated measurements.
{"title":"A local meshless technique for recovering dual forms of time-varying sources in the nonlocal inverse heat equation","authors":"Elyas Shivanian , Ahmad Jafarabadi , Mousa J. Huntul","doi":"10.1016/j.rinam.2025.100673","DOIUrl":"10.1016/j.rinam.2025.100673","url":null,"abstract":"<div><div>This study focuses on retrieving a time-dependent source term in the heat equation governed by two distinct nonlocal boundary conditions. The inverse problem is structured with an interior energy over-specification constraint. The proposed computational framework combines the partition of unity approach for spatial discretization with the finite difference scheme for temporal advancement. Through energy analysis, the semi-discrete time-stepping formulation is proven to be unconditionally stable and convergent at a rate of <span><math><mrow><mi>O</mi><mrow><mo>(</mo><mi>δ</mi><mi>t</mi><mo>)</mo></mrow></mrow></math></span>. Despite being linear and uniquely solvable, the problem is inherently ill-posed, as slight disturbances in input data can induce significant errors in the reconstructed solution. To counteract this instability, Tikhonov regularization is implemented, yielding a stable approximation even under noisy data conditions. Moreover, a novel parameter selection strategy for the regularization is introduced, which surpasses standard methods by delivering substantially improved results. Numerical simulations corroborate the scheme’s robustness, demonstrating its accuracy with noise-free inputs and its resilience when handling contaminated measurements.</div></div>","PeriodicalId":36918,"journal":{"name":"Results in Applied Mathematics","volume":"28 ","pages":"Article 100673"},"PeriodicalIF":1.3,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145614874","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-01DOI: 10.1016/j.rinam.2025.100665
Jingwei Chen
This paper develops an -parametrized framework for analyzing the strong convergence of the stochastic theta (ST) method for stochastic differential equations driven by time-changed Lévy noise (TCSDEwLNs). The analysis accommodates time–space-dependent coefficients satisfying local Lipschitz conditions. Properties of the inverse subordinator are investigated and explicit moment bounds for the exact solution are derived with jump rate incorporated. The analysis demonstrates that the ST method converges strongly with order of , establishing a precise relationship between numerical accuracy and the time-change mechanism. This theoretical advancement extends existing results and facilitates applications in finance, physics and biology where time-changed Lévy models are prevalent.
{"title":"α-scaled strong convergence of stochastic theta method for stochastic differential equations driven by time-changed Lévy noise beyond Lipschitz continuity","authors":"Jingwei Chen","doi":"10.1016/j.rinam.2025.100665","DOIUrl":"10.1016/j.rinam.2025.100665","url":null,"abstract":"<div><div>This paper develops an <span><math><mi>α</mi></math></span>-parametrized framework for analyzing the strong convergence of the stochastic theta (ST) method for stochastic differential equations driven by time-changed Lévy noise (TCSDEwLNs). The analysis accommodates time–space-dependent coefficients satisfying local Lipschitz conditions. Properties of the inverse subordinator <span><math><mi>E</mi></math></span> are investigated and explicit moment bounds for the exact solution are derived with jump rate incorporated. The analysis demonstrates that the ST method converges strongly with order of <span><math><mrow><mo>min</mo><mrow><mo>{</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>F</mi></mrow></msub><mo>,</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>G</mi></mrow></msub><mo>,</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>H</mi></mrow></msub><mo>,</mo><mi>α</mi><mo>/</mo><mn>2</mn><mo>}</mo></mrow></mrow></math></span>, establishing a precise relationship between numerical accuracy and the time-change mechanism. This theoretical advancement extends existing results and facilitates applications in finance, physics and biology where time-changed Lévy models are prevalent.</div></div>","PeriodicalId":36918,"journal":{"name":"Results in Applied Mathematics","volume":"28 ","pages":"Article 100665"},"PeriodicalIF":1.3,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145519625","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We analyze pattern formation in a two-component system within an isotropically growing or shrinking domain. By studying the evolution of a Lyapunov-like function, we derive time-dependent Turing bifurcation conditions through a stability analysis of linear perturbations across all Fourier modes. This general framework enables explicit characterization of pattern formation dynamics. Numerically, we consider two cases: a steady base state (exponential growth) and a time-dependent state (linear growth). First, we validate our approach by recovering the well-known conditions for fixed domains. Then, we simulate the Brusselator reaction system in dynamic domains, obtaining excellent agreement with our model’s predictions. These simulations highlight key pattern features, including evolution, amplitude growth, and wavenumber inertia. Our findings provide a novel energetic and geometrical perspective on the Turing bifurcation.
{"title":"Turing conditions for a two-component isotropic growing system from a potential function","authors":"Aldo Ledesma-Durán, Consuelo García-Alcántara, Iván Santamaría-Holek","doi":"10.1016/j.rinam.2025.100664","DOIUrl":"10.1016/j.rinam.2025.100664","url":null,"abstract":"<div><div>We analyze pattern formation in a two-component system within an isotropically growing or shrinking domain. By studying the evolution of a Lyapunov-like function, we derive time-dependent Turing bifurcation conditions through a stability analysis of linear perturbations across all Fourier modes. This general framework enables explicit characterization of pattern formation dynamics. Numerically, we consider two cases: a steady base state (exponential growth) and a time-dependent state (linear growth). First, we validate our approach by recovering the well-known conditions for fixed domains. Then, we simulate the Brusselator reaction system in dynamic domains, obtaining excellent agreement with our model’s predictions. These simulations highlight key pattern features, including evolution, amplitude growth, and wavenumber inertia. Our findings provide a novel energetic and geometrical perspective on the Turing bifurcation.</div></div>","PeriodicalId":36918,"journal":{"name":"Results in Applied Mathematics","volume":"28 ","pages":"Article 100664"},"PeriodicalIF":1.3,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145466100","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
An analytical approach to solve a time-fractional Cauchy problem of order based on the Ornstein–Uhlenbeck (OU), Cox–Ingersoll–Ross (CIR) and Jacobi processes with time-dependent parameters by transforming it into a system of linear fractional differential equations is established. We consider the process as an inhomogeneous Pearson diffusion and derive the analytical formulas for conditional expectations via the Volterra fractional integral equation. We also provide the -conditional moments of the OU, CIR and Jacobi processes where . Finally, we illustrate with examples of the first and second moments of the extended OU and extended CIR processes by obtaining solutions with different values and comparing to .
{"title":"Analytical solutions for time-fractional Cauchy problem based on OU, CIR and Jacobi processes with time-dependent parameters","authors":"Muntiranee Mongkolsin , Khamron Mekchay , Phiraphat Sutthimat","doi":"10.1016/j.rinam.2025.100657","DOIUrl":"10.1016/j.rinam.2025.100657","url":null,"abstract":"<div><div>An analytical approach to solve a time-fractional Cauchy problem of order <span><math><mrow><mn>0</mn><mo><</mo><mi>α</mi><mo>≤</mo><mn>1</mn></mrow></math></span> based on the Ornstein–Uhlenbeck (OU), Cox–Ingersoll–Ross (CIR) and Jacobi processes with time-dependent parameters by transforming it into a system of linear fractional differential equations is established. We consider the process as an inhomogeneous Pearson diffusion and derive the analytical formulas for conditional expectations via the Volterra fractional integral equation. We also provide the <span><math><mi>β</mi></math></span>-conditional moments of the OU, CIR and Jacobi processes where <span><math><mrow><mi>β</mi><mo>∈</mo><mi>R</mi></mrow></math></span>. Finally, we illustrate with examples of the first and second moments of the extended OU and extended CIR processes by obtaining solutions with different <span><math><mi>α</mi></math></span> values and comparing to <span><math><mrow><mi>α</mi><mo>=</mo><mn>1</mn></mrow></math></span>.</div></div>","PeriodicalId":36918,"journal":{"name":"Results in Applied Mathematics","volume":"28 ","pages":"Article 100657"},"PeriodicalIF":1.3,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145466217","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-01DOI: 10.1016/j.rinam.2025.100671
Hassan Ranjbar, Afshin Babaei
This paper develops the balanced Euler–Maruyama integrator for stochastic Volterra integral equations. First, an upper bound for the designed integrator is rigorously established in the mean square sense. Next, the scheme is proved to give a strong convergence rate of 1/2 for general diffusion matrices. Furthermore, for a special case of diffusion matrices, we theoretically detect that the established integrator super-converges with strong order 1.0. Numerical experiments are provided to confirm the theoretical findings.
{"title":"Enhancing the Euler–Maruyama integrator via a balancing strategy for stochastic Volterra integral equations","authors":"Hassan Ranjbar, Afshin Babaei","doi":"10.1016/j.rinam.2025.100671","DOIUrl":"10.1016/j.rinam.2025.100671","url":null,"abstract":"<div><div>This paper develops the balanced Euler–Maruyama integrator for stochastic Volterra integral equations. First, an upper bound for the designed integrator is rigorously established in the mean square sense. Next, the scheme is proved to give a strong convergence rate of 1/2 for general diffusion matrices. Furthermore, for a special case of diffusion matrices, we theoretically detect that the established integrator super-converges with strong order 1.0. Numerical experiments are provided to confirm the theoretical findings.</div></div>","PeriodicalId":36918,"journal":{"name":"Results in Applied Mathematics","volume":"28 ","pages":"Article 100671"},"PeriodicalIF":1.3,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145614875","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}