{"title":"Discontinuous Galerkin Methods for the Vlasov–Stokes System","authors":"Harsha Hutridurga, Krishan Kumar, Amiya K. Pani","doi":"10.1515/cmam-2023-0243","DOIUrl":"https://doi.org/10.1515/cmam-2023-0243","url":null,"abstract":"This paper develops and analyses a semi-discrete numerical method for the two-dimensional Vlasov–Stokes system with periodic boundary condition. The method is based on the coupling of the semi-discrete discontinuous Galerkin method for the Vlasov equation with discontinuous Galerkin scheme for the stationary incompressible Stokes equation. The proposed method is both mass and momentum conservative. Since it is difficult to establish non-negativity of the discrete local density, the generalized discrete Stokes operator become non-coercive and indefinite, and under the smallness condition on the discretization parameter, optimal error estimates are established with help of a modified the Stokes projection to deal with the Stokes part and, with the help of a special projection, to tackle the Vlasov part. Finally, numerical experiments based on the dG method combined with a splitting algorithm are performed.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"5 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139646424","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Robust Multigrid Methods for Discontinuous Galerkin Discretizations of an Elliptic Optimal Control Problem","authors":"Sijing Liu","doi":"10.1515/cmam-2023-0132","DOIUrl":"https://doi.org/10.1515/cmam-2023-0132","url":null,"abstract":"We consider discontinuous Galerkin methods for an elliptic distributed optimal control problem, and we propose multigrid methods to solve the discretized system. We prove that the 𝑊-cycle algorithm is uniformly convergent in the energy norm and is robust with respect to a regularization parameter on convex domains. Numerical results are shown for both 𝑊-cycle and 𝑉-cycle algorithms.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"32 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139646768","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Machine Learning Estimators: Implementation and Comparison in Python","authors":"Fabian Merle","doi":"10.1515/cmam-2023-0198","DOIUrl":"https://doi.org/10.1515/cmam-2023-0198","url":null,"abstract":"We compare different <jats:italic>machine learning</jats:italic> estimators and present details about their implementation in Python. The computational studies are conducted for <jats:italic>classification</jats:italic> as well as <jats:italic>regression</jats:italic> problems. Moreover, as one of the founding problems of machine learning, we present the specific <jats:italic>classification</jats:italic> task of <jats:italic>handwritten digit recognition</jats:italic>. In this connection, we discuss the mathematical formulation and of course the implementation details of this problem. All corresponding Python code is fully provided on request and can be downloaded from the author’s GitHub page <jats:ext-link xmlns:xlink=\"http://www.w3.org/1999/xlink\" ext-link-type=\"uri\" xlink:href=\"https://github.com/Fab1Fatal\">https://github.com/Fab1Fatal</jats:ext-link>.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"7 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139557710","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Efficient P1-FEM for Any Space Dimension in Matlab","authors":"Stefanie Beuter, Stefan A. Funken","doi":"10.1515/cmam-2022-0239","DOIUrl":"https://doi.org/10.1515/cmam-2022-0239","url":null,"abstract":"This paper deals with the efficient implementation of the finite element method with continuous piecewise linear functions (P1-FEM) in <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msup> <m:mi mathvariant=\"double-struck\">R</m:mi> <m:mi>d</m:mi> </m:msup> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2022-0239_ineq_0001.png\" /> <jats:tex-math>mathbb{R}^{d}</jats:tex-math> </jats:alternatives> </jats:inline-formula> (<jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:mrow> <m:mi>d</m:mi> <m:mo>∈</m:mo> <m:mi mathvariant=\"double-struck\">N</m:mi> </m:mrow> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2022-0239_ineq_0002.png\" /> <jats:tex-math>dinmathbb{N}</jats:tex-math> </jats:alternatives> </jats:inline-formula>). Although at present there does not seem to be a very high practical demand for finite element methods that use higher-dimensional simplicial partitions, there are some advantages in studying the efficient implementation of the method independent of the dimension. For instance, it provides additional insights into necessary data structures and the complexity of implementations. Throughout, the focus is on an efficient realization using <jats:sc>Matlab</jats:sc> built-in functions and vectorization. The fast and vectorized <jats:sc>Matlab</jats:sc> function can be easily implemented in many other vector languages and is provided in Julia, too. The complete implementation of the adaptive FEM is given, including assembling stiffness matrix, building load vector, error estimation, and adaptive mesh-refinement. Numerical experiments underline the efficiency of our freely available code which is observed to be of a slightly more than linear complexity with respect to the number of elements when memory limits are not exceeded.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"7 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139557522","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Simultaneous Reconstruction of Speed of Sound and Nonlinearity Parameter in a Paraxial Model of Vibro-Acoustography in Frequency Domain","authors":"Barbara Kaltenbacher, Teresa Rauscher","doi":"10.1515/cmam-2023-0076","DOIUrl":"https://doi.org/10.1515/cmam-2023-0076","url":null,"abstract":"In this paper, we consider the inverse problem of vibro-acoustography, a technique for enhancing ultrasound imaging by making use of nonlinear effects. It amounts to determining two spatially variable coefficients in a system of PDEs describing propagation of two directed sound beams and the wave resulting from their nonlinear interaction. To justify the use of Newton’s method for solving this inverse problem, on one hand, we verify well-definedness and differentiability of the forward operator corresponding to two versions of the PDE model; on the other hand, we consider an all-at-once formulation of the inverse problem and prove convergence of Newton’s method for its solution.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"1 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139557530","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Robust PRESB Preconditioning of a 3-Dimensional Space-Time Finite Element Method for Parabolic Problems","authors":"Ladislav Foltyn, Dalibor Lukáš, Marco Zank","doi":"10.1515/cmam-2023-0085","DOIUrl":"https://doi.org/10.1515/cmam-2023-0085","url":null,"abstract":"We present a recently developed preconditioning of square block matrices (PRESB) to be used within a parallel method to solve linear systems of equations arising from tensor-product discretizations of initial boundary-value problems for parabolic partial differential equations. We consider weak formulations in Bochner–Sobolev spaces and tensor-product finite element approximations for the heat and eddy current equations. The fast diagonalization method is employed to decouple the arising linear system of equations into auxiliary spatial complex-valued linear systems that can be solved concurrently. We prove that the real part of the system matrix is positive definite, which allows us to accelerate the flexible generalized minimal residual method (FGMRES) by the PRESB preconditioner. The action of PRESB on a vector includes two solutions with positive definite matrices. The spectrum of the preconditioned system lies between 1/2 and 1. Finally, we combine the PRESB-FGMRES method with multigrid-CG iterations and illustrate the numerical efficiency and the robustness for spatial discretizations up to 12 millions degrees of freedom.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"77 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139557707","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Quasi-Optimality of an AFEM for General Second Order Elliptic PDE","authors":"Arnab Pal, Thirupathi Gudi","doi":"10.1515/cmam-2023-0238","DOIUrl":"https://doi.org/10.1515/cmam-2023-0238","url":null,"abstract":"In this article, convergence and quasi-optimal rate of convergence of an adaptive finite element method (in short, AFEM) is shown for a general second-order non-selfadjoint elliptic PDE with convection term <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:mrow> <m:mi>b</m:mi> <m:mo>∈</m:mo> <m:msup> <m:mrow> <m:mo stretchy=\"false\">[</m:mo> <m:mrow> <m:msup> <m:mi>L</m:mi> <m:mi mathvariant=\"normal\">∞</m:mi> </m:msup> <m:mo></m:mo> <m:mrow> <m:mo stretchy=\"false\">(</m:mo> <m:mi mathvariant=\"normal\">Ω</m:mi> <m:mo stretchy=\"false\">)</m:mo> </m:mrow> </m:mrow> <m:mo stretchy=\"false\">]</m:mo> </m:mrow> <m:mi>d</m:mi> </m:msup> </m:mrow> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2023-0238_eq_0388.png\" /> <jats:tex-math>{bin[L^{infty}(Omega)]^{d}}</jats:tex-math> </jats:alternatives> </jats:inline-formula> and using minimal regularity of the dual problem, i.e., the solution of the dual problem has only <jats:inline-formula> <jats:alternatives> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msup> <m:mi>H</m:mi> <m:mn>1</m:mn> </m:msup> </m:math> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_cmam-2023-0238_eq_0235.png\" /> <jats:tex-math>{H^{1}}</jats:tex-math> </jats:alternatives> </jats:inline-formula>-regularity, which extends the result [J. M. Cascon, C. Kreuzer, R. H. Nochetto and K. G. Siebert, Quasi-optimal convergence rate for an adaptive finite element method, SIAM J. Numer. Anal. 46 2008, 5, 2524–2550]. The theoretical results are illustrated by numerical experiments.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"36 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-01-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139373452","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Numerical Approximation of Gaussian Random Fields on Closed Surfaces","authors":"Andrea Bonito, Diane Guignard, Wenyu Lei","doi":"10.1515/cmam-2022-0237","DOIUrl":"https://doi.org/10.1515/cmam-2022-0237","url":null,"abstract":"We consider the numerical approximation of Gaussian random fields on closed surfaces defined as the solution to a fractional stochastic partial differential equation (SPDE) with additive white noise. The SPDE involves two parameters controlling the smoothness and the correlation length of the Gaussian random field. The proposed numerical method relies on the Balakrishnan integral representation of the solution and does not require the approximation of eigenpairs. Rather, it consists of a sinc quadrature coupled with a standard surface finite element method. We provide a complete error analysis of the method and illustrate its performances in several numerical experiments.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"23 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-01-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139104590","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An Optimal Method for High-Order Mixed Derivatives of Bivariate Functions","authors":"Evgeniya V. Semenova, Sergiy G. Solodky","doi":"10.1515/cmam-2023-0137","DOIUrl":"https://doi.org/10.1515/cmam-2023-0137","url":null,"abstract":"The problem of optimal recovering high-order mixed derivatives of bivariate functions with finite smoothness is studied. Based on the truncation method, an algorithm for numerical differentiation is constructed, which is order-optimal both in the sense of accuracy and in terms of the amount of involved Galerkin information. Numerical examples are provided to illustrate the fact that our approach can be implemented successfully.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"1 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139079091","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Wave Propagation in High-Contrast Media: Periodic and Beyond","authors":"Élise Fressart, Barbara Verfürth","doi":"10.1515/cmam-2023-0066","DOIUrl":"https://doi.org/10.1515/cmam-2023-0066","url":null,"abstract":"This work is concerned with the classical wave equation with a high-contrast coefficient in the spatial derivative operator. We first treat the periodic case, where we derive a new limit in the one-dimensional case. The behavior is illustrated numerically and contrasted to the higher-dimensional case. For general unstructured high-contrast coefficients, we present the Localized Orthogonal Decomposition and show a priori error estimates in suitably weighted norms. Numerical experiments illustrate the convergence rates in various settings.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":"72 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139079450","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}