Maria Robert , Suresh Kumar Nadupuri , Nagaiah Chamakuri
{"title":"Moreau-Yosida regularization to optimal control of the monodomain model with pointwise control and state constraints in cardiac electrophysiology","authors":"Maria Robert , Suresh Kumar Nadupuri , Nagaiah Chamakuri","doi":"10.1016/j.cam.2024.116306","DOIUrl":"10.1016/j.cam.2024.116306","url":null,"abstract":"<div><div>In this work, we study the optimal control problem of a coupled reaction-diffusion system, which is a monodomain model in cardiac electrophysiology with pointwise bilateral control and state constraints. We adopt the Moreau-Yosida regularization as a penalization technique to deal with the state constraints. The regularized problem’s first-order optimality condition is derived. In addition, sufficient second-order optimality condition is derived for the regularized problem using the virtual control concept by proving equivalence between Moreau-Yosida regularization and the virtual control concept. The convergence of optimal controls of the regularized problems to the optimal control of the original problem is proved. Moreover, the semi-smooth Newton method for numerically finding the optimal solution to the regularization problem is presented. Finally, numerical experiments are conducted, and the results allow us to understand the extinction of the wave excitation in cardiac defibrillation in the presence of both control and state constraints.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116306"},"PeriodicalIF":2.1,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142441064","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}
Algirdas Lančinskas , Julius Žilinskas , Pascual Fernández , Blas Pelegrín
{"title":"Population-based algorithm for discrete facility location with ranking of candidate locations","authors":"Algirdas Lančinskas , Julius Žilinskas , Pascual Fernández , Blas Pelegrín","doi":"10.1016/j.cam.2024.116304","DOIUrl":"10.1016/j.cam.2024.116304","url":null,"abstract":"<div><div>Facility location problems are mathematical optimization problems that involve finding the best locations for facilities (e.g., factories, warehouses, stores) to serve customers within a given geographic area. The goal is typically to minimize costs, maximize efficiency, or optimize other objectives. Facility location problems can vary in several ways, including customer behavior rules, the type of search space, and constraints on locations for new facilities being located. These variations directly impact the complexity of the problem and the appropriate solution methods that can be used to tackle the problem. This research is focused on the discrete competitive facility location problem for an entering firm, which is a crucial scenario for new firms entering the existing market. The goal is to strategically locate new facilities to maximize their profit, while considering existing competitors. A new random search heuristic algorithm to approximate the optimal solution for discrete competitive facility location problems for firm expansion has been developed. The algorithm extends its precursor that ranks potential locations for the new facilities depending on their usefulness and uselessness in creating new solutions in the past. The new algorithm uses a population to handle and reuse the best solutions found so far and new strategies for ranking potential locations, considering features of the solutions in the population. The designed algorithm has been investigated by solving competitive facility location problems actual for an entering firms using real geographical data.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116304"},"PeriodicalIF":2.1,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142426745","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 two-step relaxed-inertial derivative-free projection based algorithm for solving standard nonlinear pseudo-monotone equations and logistic regression problems","authors":"Wenli Liu , Jinbao Jian , Jianghua Yin","doi":"10.1016/j.cam.2024.116327","DOIUrl":"10.1016/j.cam.2024.116327","url":null,"abstract":"<div><div>This paper explores a two-step inertial derivative-free projection method with a relaxation factor <span><math><mrow><mi>γ</mi><mo>∈</mo><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mn>2</mn><mo>)</mo></mrow></mrow></math></span> for solving nonlinear pseudo-monotone equations. Unlike existing inertial algorithms for the system of nonlinear pseudo-monotone equations, the inertial step of our method involves the current iteration point and the previous two iteration points. In particular, one of the inertial parameters is nonpositive. In the proposed algorithm, the search direction possesses not only the sufficient descent property but also the trust region property, independent of the line search technique. Moreover, we also establish the global convergence and the convergence rate of the algorithm without the Lipschitz continuity of the underlying mapping. Finally, our method provides competitive results on standard nonlinear monotone and pseudo-monotone equations and logistic regression problems compared with two inertial algorithms existing in the literature.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116327"},"PeriodicalIF":2.1,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142531791","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":"Domain decomposition with local time discretization for the nonlinear Stokes–Biot system","authors":"Hemanta Kunwar , Hyesuk Lee","doi":"10.1016/j.cam.2024.116311","DOIUrl":"10.1016/j.cam.2024.116311","url":null,"abstract":"<div><div>This work presents a domain decomposition method for the fluid-poroelastic structure interaction (FPSI) system, which utilizes local time integration for subproblems. To derive the domain decomposition scheme, we introduce a Lagrange multiplier and define time-dependent Steklov–Poincaré-type operators based on the interface conditions. These operators are employed to transform the coupled system into an evolutionary nonlinear interface problem, which is then solved using an iterative algorithm. This approach provides the flexibility to use different time discretization schemes and step sizes in subdomains, making it an efficient method for simulating multiphysics systems. We present numerical tests for both non-physical and physical problems to demonstrate the accuracy and efficiency of this method.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116311"},"PeriodicalIF":2.1,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142445770","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}
Yves-Cédric Bauwelinckx , Jan Dhaene , Milan van den Heuvel , Tim Verdonck
{"title":"On the causality-preservation capabilities of generative modelling","authors":"Yves-Cédric Bauwelinckx , Jan Dhaene , Milan van den Heuvel , Tim Verdonck","doi":"10.1016/j.cam.2024.116312","DOIUrl":"10.1016/j.cam.2024.116312","url":null,"abstract":"<div><div>Modelling is essential in both the financial and insurance industries. The emergence of machine learning and deep learning models offers new tools for this, but they often require large datasets that are typically unavailable in business fields due to privacy and ethical concerns. This lack of data is currently one of the main hurdles in developing better models. Generative modelling, such as Generative Adversarial Networks (GANs), can address this issue by creating synthetic data that can be freely shared. While GANs are widely studied in fields like computer vision, their use in business is limited, primarily because business questions often focus on identifying causal effects, whereas GANs and neural networks typically emphasise high-dimensional correlations. This paper explores whether GANs can produce synthetic data that reliably answers causal questions by performing causal analyses on GAN-generated data under varying assumptions. The study includes cross-sectional, time series, and complete structural model scenarios. Findings show that while basic GANs replicate causal relationships in simple cross-sectional data, they struggle with more complex structural models. In contrast, CausalGAN effectively replicates the original causal model, and TimeGAN modifies the causal representation in time series data.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116312"},"PeriodicalIF":2.1,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142445744","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":"Self representation based methods for tensor completion problem","authors":"Faezeh Aghamohammadi, Fatemeh Shakeri","doi":"10.1016/j.cam.2024.116297","DOIUrl":"10.1016/j.cam.2024.116297","url":null,"abstract":"<div><div>Tensor, the higher-order data array, naturally arises in many fields, such as information sciences, seismic data reconstruction, physics, video inpainting and so on. In this paper, we intend to provide a new model to recover a tensor, based on self-representation, for the all-mode unfoldings of the desired tensor, regardless of the tensor rank. The suggested idea generalizes self-representation to tensor and recovers an incomplete tensor by reconstructing one fiber by others in such a way that they all belong to the same subspace. We design least-square and low-rank self-representation algorithms based on the Linearized Alternating Direction Method utilizing this concept. We show that the proposed algorithms converge to the rank-minimization of the incomplete tensor.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116297"},"PeriodicalIF":2.1,"publicationDate":"2024-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142445769","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":"On semi-continuity and continuity of solution maps of parametric generalized multiobjective games in fuzzy environments","authors":"Nguyen Van Hung","doi":"10.1016/j.cam.2024.116303","DOIUrl":"10.1016/j.cam.2024.116303","url":null,"abstract":"<div><div>The purpose of this article is to study new results on the continuity of solution maps for parametric generalized multiobjective games in fuzzy environments. Firstly, we revisit parametric fuzzy generalized multiobjective games. Secondly, we establish some sufficient conditions for upper semi-continuity, Hausdorff upper semi-continuity, closedness and compactness of solution maps for such problem under suitable conditions. Thirdly, we introduce the auxiliary problem of parametric fuzzy generalized multiobjective games and the concept of strong <span><math><msub><mrow><mi>K</mi></mrow><mrow><mi>i</mi></mrow></msub></math></span>-convexity of objective functions for these problems. Finally, results on the lower semi-continuity, Hausdorff lower semi-continuity, continuity and Hausdorff continuity of solution maps to parametric generalized multiobjective games in fuzzy environments are established and studied. Many examples are given to illustrate our results.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116303"},"PeriodicalIF":2.1,"publicationDate":"2024-10-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142426743","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 computational technique for solving Volterra integro-differential equations of the third kind using meshless collocation method","authors":"E. Aourir , H. Laeli Dastjerdi","doi":"10.1016/j.cam.2024.116294","DOIUrl":"10.1016/j.cam.2024.116294","url":null,"abstract":"<div><div>The primary goal of this study is to give an approximate algorithm for solving Volterra integro-differential equations (VIDEs) of the third kind using meshless collocation techniques. The basic framework of the novel approach is based on a collocation scheme and radial basis functions (RBFs) created on scattered points. This technique requires no background approximation cells, and the algorithm is powerful, has greater stability, and does not require much computer memory. This approach represents the solution of VIDEs of the third kind by interpolating the RBFs based on the Gauss–Legendre quadrature formula. The problem is reduced to a system of algebraic equations that can be easily solved. A description of the technique for the proposed equations is provided. Furthermore, the error analysis of this scheme is examined. A few numerical experiments are presented to prove the reliability and precision of the suggested approach for solving VIDEs of the third kind. Certain problems were compared with analytical solutions, the moving least squares method, and other methods to prove the effectiveness and applicability of the approach described.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116294"},"PeriodicalIF":2.1,"publicationDate":"2024-10-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142426742","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":"Relaxation RKN-type integrators that preserve two invariants for second-order (oscillatory) systems","authors":"Kai Liu , Bin Wang , Ting Fu","doi":"10.1016/j.cam.2024.116300","DOIUrl":"10.1016/j.cam.2024.116300","url":null,"abstract":"<div><div>Recently, the relaxation technique has been widely used to impose conservation of invariants while retaining the full accuracy of the original method. So far, only a single invariant of a system has been considered. In this work, by a mild generalization of the relaxation technique, the Runge–Kutta–Nyström (RKN) integrators are modified to preserve two invariants for second-order system of Ordinary Differential Equations (ODEs). The proposed integrators can be explicit and of arbitrarily high order. The accuracy of the relaxation RKN integrators and the existence of valid relaxation parameters have been proved. The construction of the new integrators is under the framework of adapted RKN (ARKN) integrators which are specially designed for numerical solving second-order oscillatory systems. Therefore, the proposed integrators could be oscillation-preserving in the sense that they exactly integrate homogeneous oscillatory system <span><math><mrow><msup><mrow><mi>q</mi></mrow><mrow><mo>′</mo><mo>′</mo></mrow></msup><mo>+</mo><mi>K</mi><mi>q</mi><mo>=</mo><mn>0</mn></mrow></math></span>. Some numerical experiments are conducted to show the advantage and efficiency of the proposed integrators in comparison with the standard (A)RKN integrators.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116300"},"PeriodicalIF":2.1,"publicationDate":"2024-10-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142426744","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 primal–dual optimization-by-continuation algorithm","authors":"Ignace Loris , Simone Rebegoldi","doi":"10.1016/j.cam.2024.116299","DOIUrl":"10.1016/j.cam.2024.116299","url":null,"abstract":"<div><div>We present a numerical iterative optimization algorithm for the minimization of a cost function consisting of a linear combination of three convex terms, one of which is differentiable, a second one is prox-simple and the third one is the composition of a linear map and a prox-simple function. The algorithm’s special feature lies in its ability to approximate, in a single iteration run, the minimizers of the cost function for many different values of the parameters determining the relative weight of the three terms in the cost function. A proof of convergence of the algorithm, based on an inexact variable metric approach, is also provided. As a special case, one recovers a generalization of the primal–dual algorithm of Chambolle and Pock, and also of the proximal-gradient algorithm. Finally, we show how it is related to a primal–dual iterative algorithm based on inexact proximal evaluations of the non-smooth terms of the cost function.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116299"},"PeriodicalIF":2.1,"publicationDate":"2024-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142426768","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}