Longze Tan , Xueping Guo , Mingyu Deng , Jingrun Chen
{"title":"On the adaptive deterministic block Kaczmarz method with momentum for solving large-scale consistent linear systems","authors":"Longze Tan , Xueping Guo , Mingyu Deng , Jingrun Chen","doi":"10.1016/j.cam.2024.116328","DOIUrl":"10.1016/j.cam.2024.116328","url":null,"abstract":"<div><div>The Kaczmarz method is a traditional and widely used iterative algorithm for solving large-scale consistent linear systems, while its improved block Kaczmarz-type methods have received much attention and research in recent years due to their excellent numerical performance. Hence, in this paper, we present a deterministic block Kaczmarz method with momentum, which is based on Polyak’s heavy ball method and a row selection criterion for a set of block-controlled indices defined by the Euclidean norm of the residual vector. The proposed method does not need to compute the pseudo-inverses of a row submatrix at each iteration and it adaptively selects and updates the set of block control indices, thus this is different from the block Kaczmarz-type methods that are based on projection and pre-partitioning of row indices. The theoretical analysis of the proposed method shows that it converges linearly to the unique least-norm solutions of the consistent linear systems. Numerical experiments demonstrate that the deterministic block Kaczmarz method with momentum method is more efficient than the existing block Kaczmarz-type methods.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116328"},"PeriodicalIF":2.1,"publicationDate":"2024-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142554470","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":"Detective generalized multiscale hybridizable discontinuous Galerkin(GMsHDG) method for porous media","authors":"Do Yang Park, Minam Moon","doi":"10.1016/j.cam.2024.116320","DOIUrl":"10.1016/j.cam.2024.116320","url":null,"abstract":"<div><div>The Detective Generalized Multiscale Hybridizable Discontinuous Galerkin (Detective GMsHDG) method aims to further reduce the computational cost of the GMsHDG method. The GMsHDG method itself reduces the computational cost of the HDG method by employing an upscaled structure on a two-grid mesh. Given a PDE within a specified domain, we subdivide the domain into polygonal subdomains and transforms a HDG problem into globular and local problems. Globular problem concerns whether the solutions on smaller domains glue well to form a globular solution. The process involves generation of multiscale spaces, which is a vector space of functions defined on edges of the polygonal regions. A naive approximation by polynomials fails, especially in porous media, necessitating the generation of problem-specific spaces. The Detective GMsHDG method improves this process by replacing the generation of the multiscale space with the detective algorithm. The Detective GMsHDG method has two stages. First is called an offline stage. During the offline stage, we construct a detective function which, given a permeability distribution, it gives a multiscale space. Later stage is called the offline stage where, given the multiscale space, we use GMsHDG method to solve a given PDE numerically. We show numerical results to argue the liability of the solution using the detective GMsHDG method.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116320"},"PeriodicalIF":2.1,"publicationDate":"2024-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142532047","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 novel post-processed finite element method and its convergence for partial differential equations","authors":"Wenming He , Jiming Wu , Zhimin Zhang","doi":"10.1016/j.cam.2024.116319","DOIUrl":"10.1016/j.cam.2024.116319","url":null,"abstract":"<div><div>In this article, by combining high-order interpolation on coarse meshes and low-order finite element solutions on fine meshes, we propose a novel approach to improve the accuracy of the finite element method. The new method is in general suitable for most partial differential equations. For simplicity, we use the second-order elliptic problem as an example to show how the novel approach improves the accuracy of the finite element method. Numerical tests are also conducted to validate the main theoretical results.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116319"},"PeriodicalIF":2.1,"publicationDate":"2024-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142532048","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 efficient Newton-type matrix splitting algorithm for solving generalized absolute value equations with application to ridge regression problems","authors":"Xuehua Li, Cairong Chen","doi":"10.1016/j.cam.2024.116329","DOIUrl":"10.1016/j.cam.2024.116329","url":null,"abstract":"<div><div>A generalized Newton-based matrix splitting (GNMS) method is proposed for solving the generalized absolute value equations (GAVEs). Under mild conditions, the GNMS method converges to the unique solution of GAVEs. Moreover, we can obtain a few weaker convergence conditions for some existing methods. Numerical results verify the effectiveness of the proposed method.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116329"},"PeriodicalIF":2.1,"publicationDate":"2024-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142554554","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":"Dimension reduction based on time-limited cross Gramians for bilinear systems","authors":"Zhi-Hua Xiao , Yao-Lin Jiang , Zhen-Zhong Qi","doi":"10.1016/j.cam.2024.116302","DOIUrl":"10.1016/j.cam.2024.116302","url":null,"abstract":"<div><div>The cross Gramian is a useful tool in model order reduction but only applicable to square dynamical systems. Throughout this paper, time-limited cross Gramians is firstly extended to square bilinear systems that satisfies a generalized Sylvester equation, and then concepts from decentralized control are used to approximate a cross Gramian for non-square bilinear systems. In order to illustrate these cross Gramians, they are calculated efficiently based on shifted Legendre polynomials and applied to dimension reduction, which leads to a lower dimensional model by truncating the states that are associated with smaller approximate generalized Hankel singular values. In combination of the dominant subspace projection method, our reduction procedure is modified to produce a bounded-input bounded-output stable-preserved reduced model under some certain conditions. At last, the performance of numerical experiments indicates the validity of our reduction methods.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116302"},"PeriodicalIF":2.1,"publicationDate":"2024-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142554551","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 analysis of a thermoelastic problem of Moore–Gibson–Thompson type with history dependence in the thermal displacement","authors":"N. Bazarra , J.R. Fernández , R. Quintanilla","doi":"10.1016/j.cam.2024.116317","DOIUrl":"10.1016/j.cam.2024.116317","url":null,"abstract":"<div><div>In this work, we study, from the numerical point of view, a heat conduction model which is described by the history dependent version of the Moore–Gibson–Thompson equation. First, we consider the thermal problem, introducing a fully discrete approximation by means of the finite element method and the implicit Euler scheme. The discrete stability of its solution is proved, and an a priori error analysis is provided, which leads to the linear convergence imposing suitable regularity conditions. Secondly, we deal with the natural extension to the thermoelastic case. Following the analysis of the thermal problem, similar results are shown. Finally, we present some one-dimensional numerical simulations for both problems which demonstrate the accuracy of the approximations and the behavior of the discrete energies and the solutions.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116317"},"PeriodicalIF":2.1,"publicationDate":"2024-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142532044","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":"Convergence analysis of a weak Galerkin finite element method on a Shishkin mesh for a singularly perturbed fourth-order problem in 2D","authors":"Shicheng Liu , Xiangyun Meng , Qilong Zhai","doi":"10.1016/j.cam.2024.116324","DOIUrl":"10.1016/j.cam.2024.116324","url":null,"abstract":"<div><div>In this paper, we apply the weak Galerkin (WG) finite element method to solve the singularly perturbed fourth-order boundary value problem in a 2D domain. A Shishkin mesh is used to ensure that the method exhibits uniform convergence, regardless of the singular perturbation parameter. Asymptotically optimal order error estimate in a <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> discrete norm is established for the corresponding WG solutions. Numerical tests are provided to verify the convergence theory.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116324"},"PeriodicalIF":2.1,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142531790","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":"Non-Intrusive Reduced Basis two-grid method for flow and transport problems in heterogeneous porous media","authors":"Wansheng Gao , Ludovic Chamoin , Insa Neuweiler","doi":"10.1016/j.cam.2024.116321","DOIUrl":"10.1016/j.cam.2024.116321","url":null,"abstract":"<div><div>Due to its non-intrusive nature and ease of implementation, the Non-Intrusive Reduced Basis (NIRB) two-grid method has gained significant popularity in numerical computational fluid dynamics simulations. The efficiency of the NIRB method hinges on separating the procedure into offline and online stages. In the offline stage, a set of high-fidelity computations is performed to construct the reduced basis functions, which is time-consuming but is only executed once. In contrast, the online stage adapts a coarse-grid model to retrieve the expansion coefficients of the reduced basis functions. Thus it is much less costly than directly solving a high-fidelity model. However, coarse grids in heterogeneous porous media of flow models are often accompanied by upscaled hydraulic parameters (e.g. hydraulic conductivity), thus introducing upscaling errors. In this work, we introduce the two-scale idea to the existing NIRB two-grid method: when dealing with coarse-grid models, we also employ upscaled model parameters. Both the discretization and upscaling errors are compensated by the rectification post-processing. The numerical examples involve flow and heat transport problems in heterogeneous hydraulic conductivity fields, which are generated by self-affine random fields. Our research findings indicate that the modified NIRB method can effectively capture the large-scale features of numerical solutions, including pressure, velocity, and temperature. However, accurately retrieving velocity fields with small-scale features remains highly challenging.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116321"},"PeriodicalIF":2.1,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142532046","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Analysis of dependent complementary competing risks data from a generalized inverted family of lifetime distributions under a maximum ranked set sampling procedure with unequal samples","authors":"Liang Wang , Chunfang Zhang , Yogesh Mani Tripathi , Yuhlong Lio","doi":"10.1016/j.cam.2024.116309","DOIUrl":"10.1016/j.cam.2024.116309","url":null,"abstract":"<div><div>This paper explores analysis of a dependent complementary competing risks model when the failure causes are distributed by the proposed generalized inverted family of lifetime distributions. Under maximum ranked set sampling with unequal samples (MRSSU), statistical inference of model parameters and reliability indices is discussed under classical frequentist and Bayesian approaches, respectively. Maximum likelihood estimators along with their existence and uniqueness are obtained for model parameters, and associated approximate confidence intervals are constructed in consequence. Bayesian estimation is also performed with respect to general flexible priors, and the Markov Chain Monte Carlo (MCMC) algorithm is proposed for complex posterior computation. The study further examines classical and Bayesian estimations with order restriction of parameters when additional historical information is available in the MRSSU scenario. Finally, the performance of different results is evaluated through numerical simulations and a real data example is presented for demonstrating the application of our methods.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116309"},"PeriodicalIF":2.1,"publicationDate":"2024-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142441065","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 mathematical model for the role of vaccination and treatment in measles transmission in Turkey","authors":"Osman Rasit Isik , Necibe Tuncer , Maia Martcheva","doi":"10.1016/j.cam.2024.116308","DOIUrl":"10.1016/j.cam.2024.116308","url":null,"abstract":"<div><div>A previously developed and analyzed deterministic model for the transmission dynamics of measles, which takes into account the possibility of vaccinated people also contracting the disease, has been developed for Turkey. The model consists of nine compartments. The structural identifiability of the model was examined using software and detailed tables are given assuming that the incidence is given for structural identifiability. As a result of this analysis, the model is found to be structurally identifiable if at least two parameters are given along with the incidence. The parameters in this non-autonomous model are determined by considering the 1970–2021 measles case numbers in Turkey. We realize that the changes in immigration rates in Turkey, especially since the early 2000s, the changes in vaccination rates from 1970 to the present, and the changes in the vaccination rates of susceptible individuals, are significant changes in terms of time, and so we assume that these three parameters are time dependent. The practical identifiability of the model with the determined parameters is examined and it is found that if two parameters are given, all parameters except five parameters are practical identifiable. Unidentified parameters are fixed to a value by taking reference sources into account, and a model with all parameters practically identifiable is achieved. With the obtained values, the associated reproduction number of the model was obtained as 1.07 which means that the disease will persist in Turkey.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116308"},"PeriodicalIF":2.1,"publicationDate":"2024-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142532043","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}