{"title":"Dual EP matrices and their properties and characterizations","authors":"Yubing Xie, Hongxing Wang","doi":"10.1016/j.cam.2025.116785","DOIUrl":"10.1016/j.cam.2025.116785","url":null,"abstract":"<div><div>In this paper, we introduce the notions of the perfect-EP, stronger-EP, X-EP, weak-EP and P-EP matrices, by applying DMPGI, MPDGI, WDGI, dual rank and range. We discuss the properties and characterizations of these five types of dual EP matrices, as well as their relationships.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116785"},"PeriodicalIF":2.1,"publicationDate":"2025-06-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144204474","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}
Zofia Długosz , Michał Długosz , Michał Rajewski , Piotr Bogacki , Aleksandra Śniatkowska , Janusz Winiecki , Tomasz Talaśka , Alicja Bartkowska-Śniatkowska , Rafał Długosz
{"title":"A computationally optimized procedure based on discrete wavelet transform supporting analysis of medical images","authors":"Zofia Długosz , Michał Długosz , Michał Rajewski , Piotr Bogacki , Aleksandra Śniatkowska , Janusz Winiecki , Tomasz Talaśka , Alicja Bartkowska-Śniatkowska , Rafał Długosz","doi":"10.1016/j.cam.2025.116774","DOIUrl":"10.1016/j.cam.2025.116774","url":null,"abstract":"<div><div>The paper presents a concept and implementation of an algorithm whose task is to automatically determine the outline of a spot representing cancerous tissue. We assume that the location of this tissue was previously determined based on data obtained during computed tomography (CT). In medical images obtained on the basis of CT, active areas are represented by spots of appropriate intensity. For the purposes of properly selected radiotherapy, it is necessary to determine the outline of these areas. At the same time, it is very important to maintain an appropriate margin around the main area where the cancer occurs, so as not to miss small cancer foci located near the larger area. The aim of the proposed method is to process the biomedical image in such a way that the mentioned small foci are not lost, but are located inside the outline designated with an appropriate margin. The algorithm is based on a modified wavelet transform and appropriately selected linear and nonlinear filters with the possibility of their returning. The algorithm was, as a prototype, implemented in the Matlab environment and verified on images selected so that they were representative for various situations.</div><div>The system for which the proposed algorithm is being developed operates in real time. For this reason, we paid special attention to the problem of computational complexity of its particular stages. This is also important because during a medical procedure we deal with a series of biomedical images, creating a 3-D structure. Savings in this area obtained for a single image therefore give measurable benefits for the entire set of images.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116774"},"PeriodicalIF":2.1,"publicationDate":"2025-06-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144195355","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}
Jia-Lin Zhang , Guo-Feng Zhang , Zhao-Zheng Liang , Li-Dan Liao , Rui-Xia Li , Ming-Guang Liu
{"title":"Anderson acceleration preconditioned tensor splitting approach for tensor absolute value equation","authors":"Jia-Lin Zhang , Guo-Feng Zhang , Zhao-Zheng Liang , Li-Dan Liao , Rui-Xia Li , Ming-Guang Liu","doi":"10.1016/j.cam.2025.116794","DOIUrl":"10.1016/j.cam.2025.116794","url":null,"abstract":"<div><div>In this paper, we introduce the Anderson acceleration-based preconditioned tensor splitting (AAPTS) method to solve a class of tensor absolute value equation (TAVE). This approach uniquely combines Anderson acceleration (AA) with tensor splitting (TS) and preconditioning techniques, enhancing both convergence rates and computational efficiency. By dynamically adjusting iteration steps through AA, the AAPTS method accelerates the fixed-point iteration process, reduces solution time, and improves stability, particularly in dealing with complex tensor problems. This integration preserves the flexibility and generalizability of the PTS approach while significantly enhancing its convergence properties. Theoretically, we present the convergence theory for both the TS approach and the preconditioned tensor splitting (PTS) approach for TAVE, which distinguishes our work from previous analyses. On this basis, by introducing the strong semi-smoothness of the absolute value function, the local convergence theory of the AAPTS algorithms are discussed. Furthermore, the numerical experiments are conducted to clarify the efficiency and practicability of the AAPTS methods, and a detailed analysis has been carried out to examine the effects of two parameters, both derived from the AA, on the performance of the AAPTS.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116794"},"PeriodicalIF":2.1,"publicationDate":"2025-06-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144204476","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 semi-adaptive finite difference method for simulating two-sided fractional convection–diffusion quenching problems","authors":"Rumin Dong , Lin Zhu , Qin Sheng , Bingxin Zhao","doi":"10.1016/j.cam.2025.116796","DOIUrl":"10.1016/j.cam.2025.116796","url":null,"abstract":"<div><div>This paper investigates quenching solutions of an one-dimensional, two-sided Riemann–Liouville fractional order convection–diffusion problem. Fractional order spatial derivatives are discretized using weighted averaging approximations in conjunction with standard and shifted Grünwald formulas. The advective term is handled utilizing a straightforward Euler formula, resulting in a semi-discretized system of nonlinear ordinary differential equations. The conservativeness of the proposed scheme is rigorously proved and validated through simulation experiments. The study is further advanced to a fully discretized, semi-adaptive finite difference method. Detailed analysis is implemented for the monotonicity, positivity and stability of the scheme. Investigations are carried out to assess the potential impacts of the fractional order on quenching location, quenching time, and critical length. The computational results are thoroughly discussed and analyzed, providing a more comprehensive understanding of the quenching phenomena modeled through two-sided fractional order convection–diffusion problems.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116796"},"PeriodicalIF":2.1,"publicationDate":"2025-06-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144195354","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}
Zelalem Arega Worku , Jason E. Hicken , David W. Zingg
{"title":"Very high-order symmetric positive-interior quadrature rules on triangles and tetrahedra","authors":"Zelalem Arega Worku , Jason E. Hicken , David W. Zingg","doi":"10.1016/j.cam.2025.116782","DOIUrl":"10.1016/j.cam.2025.116782","url":null,"abstract":"<div><div>We present novel fully-symmetric quadrature rules with positive weights and strictly interior nodes that are exact for polynomials of total degrees up to 84 on triangles and 40 on tetrahedra. Initial guesses for solving the nonlinear systems of equations needed to derive quadrature rules are generated by forming tensor-product structures on quadrilateral/hexahedral subdomains of the simplices using the Legendre–Gauss nodes on the first half of the line reference element. In combination with a methodology for node elimination, these initial guesses lead to the development of highly efficient quadrature rules, even for very high polynomial degrees. Using existing estimates of the minimum number of quadrature points for a given degree, we show that the derived quadrature rules on triangles and tetrahedra are more than 95% and 80% efficient, respectively, for almost all degrees. The accuracy of the quadrature rules is demonstrated through numerical examples.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116782"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144204477","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 local parallel fully mixed finite element method for superposed fluid and porous layers","authors":"Jian Li , Zhuoyu Gao , Luling Cao , Zhangxin Chen","doi":"10.1016/j.cam.2025.116798","DOIUrl":"10.1016/j.cam.2025.116798","url":null,"abstract":"<div><div>Numerical simulation of the geothermal energy extraction process, described by a system of the transient Navier–Stokes–Darcy–Boussinesq equations, is important for achieving the economic extraction and utilization of geothermal resources. A fluid velocity is a crucial parameter for improving geothermal energy recovery and ensuring efficient and stable operation of geothermal systems. A fully mixed finite element method not only solves for the fluid velocity <span><math><mi>u</mi></math></span> and pressure <span><math><mi>p</mi></math></span> simultaneously, but also captures fluid changes in a sustainable manner. In addition, practical geothermal systems often involve large spatial regions, which makes it necessary to use efficient parallel algorithms in space. The key idea of a local parallel mixed finite element method based on two grids is to first use a coarse grid globally to obtain low-frequency solutions decoupled across domains, and then, at each time step, solve a series of local residual equations on overlapping subdomains with a fine grid using a local parallel algorithm to obtain high-frequency solutions. This method significantly improves computational efficiency. Finally, its effectiveness and efficiency is validated through numerical experiments.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116798"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144204473","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":"An ADMM-based interior point method for solving nonnegative tensor least squares problems and its applications","authors":"Shenghao Feng , Yimin Wei , Renjie Xu","doi":"10.1016/j.cam.2025.116793","DOIUrl":"10.1016/j.cam.2025.116793","url":null,"abstract":"<div><div>This paper considers the nonnegative tensor least squares (NTLS) problem to find a nonnegative solution to tensor equations. To overcome the shortcomings of the traditional interior point method (IPM), we add additional equality constraints to transform the nonconvex NTLS problem into a multi-convex problem. Then, we apply the alternating direction method of multipliers (ADMM) to inexactly minimize the log-barrier penalty function at each iteration. Each subproblem is convex and has an explicit solution that can be easily obtained, and the resulting algorithm is the ADMM-based interior point method (ABIP). More importantly, we only need to compute third-order tensor products and matrix products by tensor-train (TT) decomposition and inexact computing instead of high-order tensor products, which greatly reduces the computational cost. We prove that each limit point of the sequences generated by the ABIP satisfies the corresponding Karush–Kuhn–Tucker (KKT) conditions under mild assumptions. Moreover, we generalize the ABIP to the complex case to solve the complex NTLS problem. In addition, we use the ABIP to solve the <span><math><msub><mrow><mi>l</mi></mrow><mrow><mn>1</mn></mrow></msub></math></span>-norm regularized NTLS problem and the tensor absolute value equation (TAVE). Finally, some numerical experiments are provided to illustrate the validity and efficiency of the ABIP algorithm.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116793"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144204470","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":"Structure-preserving, weighted implicit–explicit schemes for multi-phase incompressible Navier–Stokes/Darcy coupled nonlocal Allen–Cahn model","authors":"Meng Li, Ke Wang, Nan Wang","doi":"10.1016/j.cam.2025.116784","DOIUrl":"10.1016/j.cam.2025.116784","url":null,"abstract":"<div><div>A multitude of substances exist as mixtures comprising multiple chemical components in natural. These substances undergo morphological changes under external influences. In the phase-field model coupled with fluid flow, the dynamic movement and evolution of the phase interface intricately interact with the fluid motion. This article focuses on the N-component models that couple the conservative Allen–Cahn equation with two types of incompressible fluid flow systems: the Navier–Stokes equation and the Darcy equation. By utilizing the scalar auxiliary variable method and the projection method, we innovatively construct two types of structure-preserving weighted implicit–explicit schemes for the coupled models, resulting in fully decoupled linear systems and second-order accuracy in time. The schemes are proved to be mass-conservative. In addition, with the application of <span><math><mi>G</mi></math></span>-norm inspired by the idea of <span><math><mi>G</mi></math></span>-stability, we rigorously establish their unconditional energy stability. Finally, the performance of the proposed schemes is verified by some numerical simulations.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116784"},"PeriodicalIF":2.1,"publicationDate":"2025-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144204475","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":"Impact of inaccuracies of injection positioning for hyperthermia treatment for cancer","authors":"Gustavo Resende Fatigate , Gustavo Coelho Martins , Marcelo Lobosco , Ruy Freitas Reis","doi":"10.1016/j.cam.2025.116780","DOIUrl":"10.1016/j.cam.2025.116780","url":null,"abstract":"<div><div>Cancer remains a global health concern, necessitating novel treatment strategies. Magnetic nanoparticle-mediated hyperthermia, which induces localized heating of cancerous tissues, is a promising adjunct to conventional therapies. This study utilizes the Pennes bioheat model to evaluate tissue damage during hyperthermia treatments. So, once we obtain the optimal heating location based on the literature, this work uses forward uncertainty quantification (UQ) integrated into the bioheat model to identify a region of interest where minor deviations in injection positions have minimal impact on treatment outcomes. Given the high computational demands for this simulation, this study employs a parallel computation strategy utilizing CUDA architecture to accelerate execution. The study also compares the computational efficiency of CPU and GPU architectures, highlighting the advantages of parallelization in reducing computational time. Finally, the study’s main contribution is the integration of uncertainty quantification with the GPU-accelerated resolution of the Pennes bioheating model, which significantly improves simulation efficiency – achieving a <span><math><mrow><mn>68</mn><mo>.</mo><mn>18</mn><mo>×</mo></mrow></math></span> speedup compared to the serial CPU implementation – and allows robust analysis of hyperthermia treatments in cancer therapy.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116780"},"PeriodicalIF":2.1,"publicationDate":"2025-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144204471","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 boundary element formulation without singular integrals","authors":"M. Tadi","doi":"10.1016/j.cam.2025.116783","DOIUrl":"10.1016/j.cam.2025.116783","url":null,"abstract":"<div><div>This note is on a boundary element formulation for Laplace and Helmholtz equations. It introduces a practical approach to avoid the computations of singular integrals. This note considers 2-D problems, however the method can also be applied to 3-D problems. It introduces first-order and second order accurate formulation. The accuracy can also be improved. Numerical examples are presented for Laplace, Helmholtz, and modified Helmholtz operators.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"472 ","pages":"Article 116783"},"PeriodicalIF":2.1,"publicationDate":"2025-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144211984","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}