{"title":"A fast Bregman projection method for linearly constrained optimization problems","authors":"Yu-Xin Ye , Jun-Feng Yin , Yu-Rui Jiang , Ze Wang","doi":"10.1016/j.cam.2025.116792","DOIUrl":"10.1016/j.cam.2025.116792","url":null,"abstract":"<div><div>A fast Bregman projection method is proposed for the solution of linearly constrained optimization problems by greedily making use of the residual to generate a weighted hyperplane. The convergence theories of the proposed method are established and studied under both noise-free and noisy cases. The linear convergence rate and its upper bound are derived in details, which is better than that of the sketched Bregman projection method. Numerical experiments verify the efficiency of the proposed method in terms of the number of iterations and CPU time.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116792"},"PeriodicalIF":2.1,"publicationDate":"2025-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144211983","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Further characterizations of W-weighted core-EP matrices","authors":"Ehsan Kheirandish , Abbas Salemi , Qing-Wen Wang","doi":"10.1016/j.cam.2025.116788","DOIUrl":"10.1016/j.cam.2025.116788","url":null,"abstract":"<div><div>Suppose that <span><math><mrow><mi>A</mi><mo>∈</mo><msub><mrow><mi>M</mi></mrow><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>ℂ</mi><mo>)</mo></mrow></mrow></math></span> and <span><math><mrow><mn>0</mn><mo>≠</mo><mi>W</mi><mo>∈</mo><msub><mrow><mi>M</mi></mrow><mrow><mi>n</mi><mo>,</mo><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>ℂ</mi><mo>)</mo></mrow></mrow></math></span>. The matrix <span><math><mi>A</mi></math></span> is called W-weighted core-EP matrix, if <span><math><mrow><msup><mrow><mi>A</mi></mrow><mrow><mi>†</mi></mrow></msup><msup><mrow><mi>A</mi></mrow><mrow><mi>c</mi><mo>,</mo><mi>W</mi></mrow></msup><mi>W</mi><mo>=</mo><mi>W</mi><msup><mrow><mi>A</mi></mrow><mrow><mi>c</mi><mo>,</mo><mi>W</mi></mrow></msup><msup><mrow><mi>A</mi></mrow><mrow><mi>†</mi></mrow></msup></mrow></math></span> where <span><math><msup><mrow><mi>A</mi></mrow><mrow><mi>c</mi><mo>,</mo><mi>W</mi></mrow></msup></math></span> is the <span><math><mi>W</mi></math></span>-weighted core part of <span><math><mi>A</mi></math></span>, and <span><math><msup><mrow><mi>A</mi></mrow><mrow><mi>†</mi></mrow></msup></math></span> represents the Moore–Penrose inverse of <span><math><mi>A</mi></math></span>. In this paper, some equivalent conditions of <span><math><mi>W</mi></math></span>-weighted core-EP matrices are investigated. Moreover, the W-weighted core part of <span><math><mi>A</mi></math></span> is introduced, and characterizations of <span><math><mi>W</mi></math></span>-weighted <em>DMP</em>, <span><math><mi>W</mi></math></span>-weighted <em>MPD</em>, and <span><math><mi>W</mi></math></span>-weighted <em>CMP</em> inverses are given. Finally, applications of the W-weighted core part of <span><math><mi>A</mi></math></span> in solving singular linear systems are studied.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116788"},"PeriodicalIF":2.1,"publicationDate":"2025-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144221733","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Mostafa Kadiri , Mohammed Louaked , Saber Trabelsi
{"title":"Tumor growth model with chemotaxis and active transport: Control and parameters recovery","authors":"Mostafa Kadiri , Mohammed Louaked , Saber Trabelsi","doi":"10.1016/j.cam.2025.116769","DOIUrl":"10.1016/j.cam.2025.116769","url":null,"abstract":"<div><div>In this paper, we formulate and analyze an optimal control problem for a system of Cahn–Hilliard equations modeling tumor growth, accounting for chemotaxis and active transport. The dynamical system was introduced in Hawkins-Daarud et al. (2012), and mathematical results of existence and uniqueness of weak solutions were obtained in Garcke and Yayla (2020). In this contribution, we prove the continuous dependence of the solutions on the physical parameters in addition to the initial data. In addition, we introduce an optimal control problem where the cost functional depends on a target function, but most importantly, on physical parameters targets. We establish the existence of a unique minimizer and provide optimality conditions. Eventually, we present simple numerical illustrations in full agreement with our theoretical results.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116769"},"PeriodicalIF":2.1,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144178826","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Numerical simulation of fluid–structure–acoustic interaction motivated by healthy human phonation","authors":"Jan Valášek , Petr Sváček","doi":"10.1016/j.cam.2025.116771","DOIUrl":"10.1016/j.cam.2025.116771","url":null,"abstract":"<div><div>This paper presents the mathematical model for fluid–structure–acoustic interaction (FSAI) problems with a particular focus on numerical simulations of healthy human phonation. The mathematical model consists of elastic body deformation, fluid flow, acoustics and their coupling conditions. A special attention is paid to challenges posed by the vocal folds contact, which is addressed in the model by the inlet penalization boundary conditions and an introduction of fictitious porous media.</div><div>The numerical model is based on the stabilized finite element method, implemented in an in-house solver. Numerical results are presented for the considered symmetric vocal fold vibrations, for which particularly the influence of penalization parameter and gap threshold is analyzed. Finally, the aeroacoustic simulations in the vocal tract model representing the vowel [u:] are performed with the help of the Lighthill analogy and the aeroacoustic wave equation approach.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116771"},"PeriodicalIF":2.1,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144185766","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Limit theorems for conditional U-statistics analysis on hyperspheres for missing at random data in the presence of measurement error","authors":"Salim Bouzebda, Nourelhouda Taachouche","doi":"10.1016/j.cam.2025.116811","DOIUrl":"10.1016/j.cam.2025.116811","url":null,"abstract":"<div><div>Missing data and measurement errors are prevalent challenges in modern statistical analyses, mainly when observations lie on complex structures like unit hyperspheres. To address these issues, we introduce a comprehensive framework for conditional <em>U</em>-statistics of general order, tailored explicitly for data missing at random and contaminated by measurement errors in such settings. We propose a novel deconvolution method for these conditional <em>U</em>-statistics and, for the first time, investigate its convergence rate and asymptotic distribution. Our unified approach establishes general asymptotic properties under broad model conditions, enabling us to derive asymptotic confidence intervals based on the estimator’s distribution. To demonstrate the practical significance of our framework, we provide new insights into the Kendall rank correlation coefficient and address discrimination problems.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116811"},"PeriodicalIF":2.1,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144168354","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The Tustin Method for approximating eigenvalues of delay systems","authors":"Markel Irastorza-Zabalegi , Felipe Ponce-Vanegas","doi":"10.1016/j.cam.2025.116772","DOIUrl":"10.1016/j.cam.2025.116772","url":null,"abstract":"<div><div>Numerical approximation of eigenvalues of Delay Differential Equations (DDEs) is an active field of research due to its impact in the modeling of many processes of industrial interest. In this work we introduce the <em>Tustin Method</em>, a new numerical method for Linear Time-Invariant (LTI) systems based on the Cayley transform. We introduce two alternatives, one using a trigonometric basis, and other using splines with maximum smoothness. For splines, we exploit the so called <span><math><mi>h</mi></math></span>-<span><math><mi>p</mi></math></span>-<span><math><mi>k</mi></math></span> refinement, achieving very high accuracy. We prove that the Tustin method recovers all the eigenvalues, and we give estimates for the rate of convergence. We make numerical experiments for generic systems, for systems with large delay, and for systems close to the stability boundary.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116772"},"PeriodicalIF":2.1,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144195356","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A combined Method Of Lines and Orthogonal Collocation with Second kind Chebyshev nodes for convection–diffusion–reaction equation with Danckwerts Conditions","authors":"L. Agud-Albesa , M. Boix-García , J.R. Torregrosa","doi":"10.1016/j.cam.2025.116731","DOIUrl":"10.1016/j.cam.2025.116731","url":null,"abstract":"<div><div>We analyze and develop two numerical methods to solve the convection–diffusion–reaction equation with Danckwerts boundary conditions. One of the methods is an approach based on the method of lines using spatial discretization by orthogonal collocation. While this method has been applied to other equations, it has not been previously studied for this particular case. Furthermore, the convergence of the method is demonstrated for various values of the Péclet and Damköhler numbers. We also describe the implementation of a weighted residual method by orthogonal collocation method, using Lagrange polynomials and Chebyshev nodes of the second kind to solve the same problem. Both methods are presented in matrix form to facilitate its implementation in <span>Matlab</span>. Finally, we compare the results with both the analytical solution and those from a previous conventional method of lines discretization based on finite differences, implemented by the authors. The computed errors demonstrate that the adapted method of lines with orthogonal collocation yields a more accurate overall approximation than the alternative approaches.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116731"},"PeriodicalIF":2.1,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144204472","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Adaptive finite element method for phase field fracture models based on recovery error estimates","authors":"Tian Tian , Chunyu Chen , Liang He , Huayi Wei","doi":"10.1016/j.cam.2025.116732","DOIUrl":"10.1016/j.cam.2025.116732","url":null,"abstract":"<div><div>The phase-field model is a widely used mathematical approach for describing crack propagation in continuum damage fractures. In the context of phase field fracture simulations, adaptive finite element methods (AFEM) are often employed to address the mesh size dependency of the model. However, existing AFEM approaches for this application frequently rely on heuristic adjustments and empirical parameters for mesh refinement. In this paper, we introduce an adaptive finite element method based on a recovery type posteriori error estimates approach grounded in theoretical analysis. This method transforms the gradient of the numerical solution into a smoother function space, using the difference between the recovered gradient and the original numerical gradient as an error indicator for adaptive mesh refinement. This enables the automatic capture of crack propagation directions without the need for empirical parameters. We have implemented this adaptive method for the Hybrid formulation of the phase-field model using the open-source software package FEALPy. The accuracy and efficiency of the proposed approach are demonstrated through simulations of classical 2D and 3D brittle fracture examples, validating the robustness and effectiveness of our implementation.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116732"},"PeriodicalIF":2.1,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144185787","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Truncated Euler–Maruyama method for hybrid stochastic functional differential equations with infinite time delay","authors":"Jingchao Zhou, Henglei Xu, Xuerong Mao","doi":"10.1016/j.cam.2025.116773","DOIUrl":"10.1016/j.cam.2025.116773","url":null,"abstract":"<div><div>Li et al. (2023) developed a new theory to approximate the solution of hybrid stochastic functional differential equations (SFDEs) with infinite time delay via the numerical solution of the corresponding hybrid SFDEs with finite time delay. But hybrid SFDEs were required to be globally Lipschitz continuous. In this paper, we will lift this restriction. Under the local Lipschitz condition and the Khasminskii-type condition, numerical solutions of hybrid SFDEs with infinite time delay will be designed by using the truncated Euler–Maruyama method. The strong convergence and convergence rate of the numerical solutions in <span><math><mrow><msup><mrow><mi>L</mi></mrow><mrow><mi>q</mi></mrow></msup><mrow><mo>(</mo><mi>q</mi><mo>≥</mo><mn>2</mn><mo>)</mo></mrow></mrow></math></span> will be obtained. Finally, an example to stochastic functional volatility model is given to demonstrate the effectiveness of our new theory.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116773"},"PeriodicalIF":2.1,"publicationDate":"2025-05-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144168355","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Error analysis of piecewise collocation method for delay integro-differential–algebraic equations","authors":"P. Teimoori, S. Pishbin","doi":"10.1016/j.cam.2025.116750","DOIUrl":"10.1016/j.cam.2025.116750","url":null,"abstract":"<div><div>Piecewise collocation method for integro-differential–algebraic equations (IDAEs) with non-vanishing delay is applied. The theory needed to understand the numerical approach and analyze the numerical treatment by collocation methods is developed. Based on the structure of solutions of delay IDAEs and identification of initial discontinuity points, numerical analysis and properties of the convergence are investigated. Here, collocation method which depends on the numerical solution in a fixed number of previous time steps is described by the constructive technique and dividing the definition domain into several subintervals according to the initial discontinuous points associated with the delay function. Finally, the numerical experiments are given to demonstrate the conclusions.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116750"},"PeriodicalIF":2.1,"publicationDate":"2025-05-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144147250","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}