Pub Date : 2023-06-15DOI: 10.1007/s10915-023-02252-z
Xiangyun Bai, Jiebao Sun, Jie Shen, Wenjuan Yao, Zhichang Guo
{"title":"A Ginzburg-Landau-H-1 Model and Its SAV Algorithm for Image Inpainting","authors":"Xiangyun Bai, Jiebao Sun, Jie Shen, Wenjuan Yao, Zhichang Guo","doi":"10.1007/s10915-023-02252-z","DOIUrl":"https://doi.org/10.1007/s10915-023-02252-z","url":null,"abstract":"","PeriodicalId":21812,"journal":{"name":"SIAM J. Sci. Comput.","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-06-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76947690","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Generalized Singleton Bound and List-Decoding Reed-Solomon Codes Beyond the Johnson Radius","authors":"Chong Shangguan, Itzhak Tamo","doi":"10.1137/20m138795x","DOIUrl":"https://doi.org/10.1137/20m138795x","url":null,"abstract":"","PeriodicalId":21812,"journal":{"name":"SIAM J. Sci. Comput.","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90286649","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Uniform Restricted Chase Termination","authors":"Tomasz Gogacz, J. Marcinkowski, Andreas Pieris","doi":"10.1137/20m1377035","DOIUrl":"https://doi.org/10.1137/20m1377035","url":null,"abstract":"","PeriodicalId":21812,"journal":{"name":"SIAM J. Sci. Comput.","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91541223","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Higher Order Extended Dynamic Mode Decomposition Based on the Structured Total Least Squares","authors":"Weiyang Ding, J. Li","doi":"10.1137/21m1463665","DOIUrl":"https://doi.org/10.1137/21m1463665","url":null,"abstract":"","PeriodicalId":21812,"journal":{"name":"SIAM J. Sci. Comput.","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90710856","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Dynamic Finite-Budget Allocation of Stratified Sampling with Adaptive Variance Reduction by Strata","authors":"Chenxiao Song, Ray Kawai","doi":"10.1137/21m1430996","DOIUrl":"https://doi.org/10.1137/21m1430996","url":null,"abstract":"","PeriodicalId":21812,"journal":{"name":"SIAM J. Sci. Comput.","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"83826567","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Shape Optimization of the Stokes Eigenvalue Problem","authors":"Jiajie Li, Shengfeng Zhu","doi":"10.1137/21m1451543","DOIUrl":"https://doi.org/10.1137/21m1451543","url":null,"abstract":"","PeriodicalId":21812,"journal":{"name":"SIAM J. Sci. Comput.","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"89311227","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An Incomplete Cholesky Preconditioner Based on Orthogonal Approximations","authors":"Artem Napov","doi":"10.1137/21m1468334","DOIUrl":"https://doi.org/10.1137/21m1468334","url":null,"abstract":"","PeriodicalId":21812,"journal":{"name":"SIAM J. Sci. Comput.","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"79214859","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
. Whether due to the sheer size of a computational domain, the fine resolution required, 5 or the multiples scales and stochasticity of the dynamics, the dimensionality of a system must often 6 be reduced so that problems of interest become computationally tractable. In this paper, we develop 7 retractions for time-integration schemes that efficiently and accurately evolve the dynamics of a 8 system’s low-rank approximation. Through differential geometry, we analyze the error incurred 9 at each time step due to the high-order curvature of the manifold of fixed-rank matrices. We first 10 obtain a novel, explicit, computationally inexpensive set of algorithms that we refer to as perturbative 11 retractions and show that the set converges to an ideal retraction that projects optimally and exactly 12 to the manifold of fixed-rank matrices by reducing what we define as the projection-retraction error. 13 Furthermore, each perturbative retraction itself exhibits high-order convergence to the best low- 14 rank approximation of the full-rank solution. Using perturbative retractions, we then develop a new 15 class of integration techniques that we refer to as dynamically orthogonal Runge-Kutta (DORK) 16 schemes. DORK schemes integrate along the nonlinear manifold, updating the subspace upon which 17 we project the system’s dynamics as it is integrated. Through numerical test cases, we demonstrate 18 our schemes for matrix addition, real-time data compression, and deterministic and stochastic partial 19 differential equations. We find that DORK schemes are highly accurate by incorporating knowledge 20 of the dynamic, nonlinear manifold’s high-order curvature and computationally efficient by limiting 21 the growing rank needed to represent the evolving dynamics. 22 schemes that account for the 774 evolution of the reduced-order integration space during the time-step. We show that 775 DORK schemes are (i) highly accurate by incorporating knowledge of the dynamic, 776 nonlinear manifold’s high-order curvature from stage to stage and (ii) computationally 777 efficient by limiting the growing rank needed to represent the evolving dynamics L . In 778 our real-time data compression example, we show that these retractions may be used 779 when the dynamics are data-driven (rather than model-driven), and the perturbative 780 retractions with high-order corrections drift from the best low-rank approximation at 781 a much slower rate than retractions with low-order corrections. Note that adjusting 782 how the dynamics are calculated to correct the drift in the first place is important 783 to obtain the most accurate scheme. Lastly, we show that the retractions may be 784 used in deterministic and stochastic differential equations: the dynamical low-rank 785 approximation and the retractions are agnostic to the nature of the mathematical 786 spaces we choose to compress, whether deterministic or stochastic. In both cases, the 787 low-rank approximation converges quickly not just to th
{"title":"Dynamically Orthogonal Runge-Kutta Schemes with Perturbative Retractions for the Dynamical Low-Rank Approximation","authors":"A. Charous, Pierre FJ Lermusiaux","doi":"10.1137/21m1431229","DOIUrl":"https://doi.org/10.1137/21m1431229","url":null,"abstract":". Whether due to the sheer size of a computational domain, the fine resolution required, 5 or the multiples scales and stochasticity of the dynamics, the dimensionality of a system must often 6 be reduced so that problems of interest become computationally tractable. In this paper, we develop 7 retractions for time-integration schemes that efficiently and accurately evolve the dynamics of a 8 system’s low-rank approximation. Through differential geometry, we analyze the error incurred 9 at each time step due to the high-order curvature of the manifold of fixed-rank matrices. We first 10 obtain a novel, explicit, computationally inexpensive set of algorithms that we refer to as perturbative 11 retractions and show that the set converges to an ideal retraction that projects optimally and exactly 12 to the manifold of fixed-rank matrices by reducing what we define as the projection-retraction error. 13 Furthermore, each perturbative retraction itself exhibits high-order convergence to the best low- 14 rank approximation of the full-rank solution. Using perturbative retractions, we then develop a new 15 class of integration techniques that we refer to as dynamically orthogonal Runge-Kutta (DORK) 16 schemes. DORK schemes integrate along the nonlinear manifold, updating the subspace upon which 17 we project the system’s dynamics as it is integrated. Through numerical test cases, we demonstrate 18 our schemes for matrix addition, real-time data compression, and deterministic and stochastic partial 19 differential equations. We find that DORK schemes are highly accurate by incorporating knowledge 20 of the dynamic, nonlinear manifold’s high-order curvature and computationally efficient by limiting 21 the growing rank needed to represent the evolving dynamics. 22 schemes that account for the 774 evolution of the reduced-order integration space during the time-step. We show that 775 DORK schemes are (i) highly accurate by incorporating knowledge of the dynamic, 776 nonlinear manifold’s high-order curvature from stage to stage and (ii) computationally 777 efficient by limiting the growing rank needed to represent the evolving dynamics L . In 778 our real-time data compression example, we show that these retractions may be used 779 when the dynamics are data-driven (rather than model-driven), and the perturbative 780 retractions with high-order corrections drift from the best low-rank approximation at 781 a much slower rate than retractions with low-order corrections. Note that adjusting 782 how the dynamics are calculated to correct the drift in the first place is important 783 to obtain the most accurate scheme. Lastly, we show that the retractions may be 784 used in deterministic and stochastic differential equations: the dynamical low-rank 785 approximation and the retractions are agnostic to the nature of the mathematical 786 spaces we choose to compress, whether deterministic or stochastic. In both cases, the 787 low-rank approximation converges quickly not just to th","PeriodicalId":21812,"journal":{"name":"SIAM J. Sci. Comput.","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88368051","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"(text{PIN}^{mathcal L}) : Preconditioned Inexact Newton with Learning Capability for Nonlinear System of Equations","authors":"L. Luo, X. Cai","doi":"10.1137/22m1507942","DOIUrl":"https://doi.org/10.1137/22m1507942","url":null,"abstract":"","PeriodicalId":21812,"journal":{"name":"SIAM J. Sci. Comput.","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78016035","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Recursive Algorithms to Update a Numerical Basis Matrix of the Null Space of the Block Row, (Banded) Block Toeplitz, and Block Macaulay Matrix","authors":"C. Vermeersch, B. Moor","doi":"10.1137/22m1495858","DOIUrl":"https://doi.org/10.1137/22m1495858","url":null,"abstract":"","PeriodicalId":21812,"journal":{"name":"SIAM J. Sci. Comput.","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78804554","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：