首页 > 最新文献

Electronic Transactions on Numerical Analysis最新文献

英文 中文
Optimized surface parameterizations with applications to Chinese virtual broadcasting 优化表面参数化及其在中文虚拟广播中的应用
IF 1.3 4区 数学 Q2 Mathematics Pub Date : 2020-01-01 DOI: 10.1553/etna_vol53s383
M. Yueh, H. Huang, Tie-xiang Li, Wen-Wei Lin, S. Yau
Surface parameterizations have been widely applied in the computer-aided design for the geometric processing tasks of surface registration, remeshing, texture mapping, and so on. In this paper, we present an efficient balanced energy minimization (BEM) algorithm for the computation of simply connected open surface parameterizations with balanced angle and area distortions. The existence of a nontrivial accumulation function of the BEM algorithm is guaranteed under some mild conditions and the limiting function is shown to be one-to-one. Comparisons of the BEM algorithm with the angleand the area-preserving parameterizations show that the angular distortion is close to that of the angle-preserving parameterization while the area distortion is significantly improved. An application of the BEM on the Chinese virtual broadcasting technique is demonstrated thereafter, which is consisted of surface remeshing, registration, and morphing.
曲面参数化在计算机辅助设计中广泛应用于曲面配准、网格划分、纹理映射等几何处理任务。本文提出了一种有效的平衡能量最小化(BEM)算法,用于计算具有平衡角度和面积畸变的单连通开曲面参数化。在一些温和的条件下,保证了边界元算法的非平凡累积函数的存在性,并证明了极限函数是一对一的。将边界元算法与保角参数化和保面积参数化进行比较,结果表明,边界元算法的角畸变与保角参数化的角畸变接近,而面积畸变明显改善。然后演示了边界元在中文虚拟广播技术中的应用,包括曲面重划分、配准和变形。
{"title":"Optimized surface parameterizations with applications to Chinese virtual broadcasting","authors":"M. Yueh, H. Huang, Tie-xiang Li, Wen-Wei Lin, S. Yau","doi":"10.1553/etna_vol53s383","DOIUrl":"https://doi.org/10.1553/etna_vol53s383","url":null,"abstract":"Surface parameterizations have been widely applied in the computer-aided design for the geometric processing tasks of surface registration, remeshing, texture mapping, and so on. In this paper, we present an efficient balanced energy minimization (BEM) algorithm for the computation of simply connected open surface parameterizations with balanced angle and area distortions. The existence of a nontrivial accumulation function of the BEM algorithm is guaranteed under some mild conditions and the limiting function is shown to be one-to-one. Comparisons of the BEM algorithm with the angleand the area-preserving parameterizations show that the angular distortion is close to that of the angle-preserving parameterization while the area distortion is significantly improved. An application of the BEM on the Chinese virtual broadcasting technique is demonstrated thereafter, which is consisted of surface remeshing, registration, and morphing.","PeriodicalId":50536,"journal":{"name":"Electronic Transactions on Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2020-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"67352258","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}
引用次数: 4
Finite element discretization of semilinear acoustic wave equations with kinetic boundary conditions 具有动力学边界条件的半线性声波方程的有限元离散化
IF 1.3 4区 数学 Q2 Mathematics Pub Date : 2020-01-01 DOI: 10.1553/ETNA_VOL53S522
M. Hochbruck, Jan Leibold
We consider isoparametric finite element discretizations of semilinear acoustic wave equations with kinetic boundary conditions and derive a corresponding error bound as our main result. The difficulty is that such problems are stated on domains with curved boundaries and this renders the discretizations nonconforming. Our approach is to provide a unified error analysis for nonconforming space discretizations for semilinear wave equations. In particular, we introduce a general, abstract framework for nonconforming space discretizations in which we derive a-priori error bounds in terms of interpolation, data and conformity errors. The theory applies to a large class of problems and discretizations that fit into the abstract framework. The error bound for wave equations with kinetic boundary conditions is obtained from the general theory by inserting known interpolation and geometric error bounds into the abstract error result of the unified error analysis.
我们考虑具有动力学边界条件的半线性声波方程的等参有限元离散化,并推导出相应的误差界作为我们的主要结果。困难之处在于这些问题是在具有弯曲边界的域上表述的,这使得离散化不一致。我们的方法是为半线性波动方程的非一致性空间离散提供统一的误差分析。特别地,我们引入了非一致性空间离散化的一般抽象框架,在该框架中,我们推导了插值误差、数据误差和一致性误差的先验误差界。该理论适用于符合抽象框架的大量问题和离散化。在统一误差分析的抽象误差结果中插入已知的插值和几何误差边界,从一般理论推导出具有运动边界条件的波动方程的误差边界。
{"title":"Finite element discretization of semilinear acoustic wave equations with kinetic boundary conditions","authors":"M. Hochbruck, Jan Leibold","doi":"10.1553/ETNA_VOL53S522","DOIUrl":"https://doi.org/10.1553/ETNA_VOL53S522","url":null,"abstract":"We consider isoparametric finite element discretizations of semilinear acoustic wave equations with kinetic boundary conditions and derive a corresponding error bound as our main result. The difficulty is that such problems are stated on domains with curved boundaries and this renders the discretizations nonconforming. Our approach is to provide a unified error analysis for nonconforming space discretizations for semilinear wave equations. In particular, we introduce a general, abstract framework for nonconforming space discretizations in which we derive a-priori error bounds in terms of interpolation, data and conformity errors. The theory applies to a large class of problems and discretizations that fit into the abstract framework. The error bound for wave equations with kinetic boundary conditions is obtained from the general theory by inserting known interpolation and geometric error bounds into the abstract error result of the unified error analysis.","PeriodicalId":50536,"journal":{"name":"Electronic Transactions on Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2020-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"67352310","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}
引用次数: 9
A posteriori stopping criteria for space-time domain decomposition for the heat equation in mixed formulations 混合公式中热方程空时域分解的后验停止准则
IF 1.3 4区 数学 Q2 Mathematics Pub Date : 2018-06-18 DOI: 10.1553/ETNA_VOL49S151
S. Hassan, C. Japhet, M. Vohralík
We propose and analyse a posteriori estimates for global-in-time, nonoverlapping domain decomposition methods for heterogeneous and anisotropic porous media diffusion problems. We consider mixed formulations, with a lowest-order Raviart–Thomas–Nedelec discretization, often used for such problems. Optimized Robin transmission conditions are employed on the space-time interface between subdomains, and different time grids are used to adapt to different time scales in the subdomains. Our estimators allow to distinguish the spatial discretization, the temporal discretization, and the domain decomposition error components. We design an adaptive space-time domain decomposition algorithm, wherein the iterations are stopped when the domain decomposition error does not affect significantly the global error. Thus, a guaranteed bound on the overall error is obtained on each iteration of the space-time domain decomposition algorithm, and simultaneously important savings in terms of the number of domain decomposition iterations can be achieved. Numerical results for two-dimensional problems with strong heterogeneities and local time stepping are presented to illustrate the performance of our adaptive domain decomposition algorithm.
我们提出并分析了非均匀和各向异性多孔介质扩散问题的全局时间非重叠域分解方法的后验估计。我们考虑通常用于此类问题的具有最低阶Raviart–Thomas–Nedelec离散化的混合公式。子域之间的时空接口采用了优化的Robin传输条件,并使用不同的时间网格来适应子域中的不同时间尺度。我们的估计器允许区分空间离散化、时间离散化和域分解误差分量。我们设计了一种自适应时空域分解算法,其中当域分解误差对全局误差没有显著影响时,停止迭代。因此,在时空域分解算法的每次迭代中都获得了总体误差的保证界,并且同时可以实现在域分解迭代次数方面的重要节省。给出了具有强非均匀性和局部时间步长的二维问题的数值结果,以说明我们的自适应域分解算法的性能。
{"title":"A posteriori stopping criteria for space-time domain decomposition for the heat equation in mixed formulations","authors":"S. Hassan, C. Japhet, M. Vohralík","doi":"10.1553/ETNA_VOL49S151","DOIUrl":"https://doi.org/10.1553/ETNA_VOL49S151","url":null,"abstract":"We propose and analyse a posteriori estimates for global-in-time, nonoverlapping domain decomposition methods for heterogeneous and anisotropic porous media diffusion problems. We consider mixed formulations, with a lowest-order Raviart–Thomas–Nedelec discretization, often used for such problems. Optimized Robin transmission conditions are employed on the space-time interface between subdomains, and different time grids are used to adapt to different time scales in the subdomains. Our estimators allow to distinguish the spatial discretization, the temporal discretization, and the domain decomposition error components. We design an adaptive space-time domain decomposition algorithm, wherein the iterations are stopped when the domain decomposition error does not affect significantly the global error. Thus, a guaranteed bound on the overall error is obtained on each iteration of the space-time domain decomposition algorithm, and simultaneously important savings in terms of the number of domain decomposition iterations can be achieved. Numerical results for two-dimensional problems with strong heterogeneities and local time stepping are presented to illustrate the performance of our adaptive domain decomposition algorithm.","PeriodicalId":50536,"journal":{"name":"Electronic Transactions on Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42301645","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}
引用次数: 10
Asymptotic consistent exponential-type integrators for Klein-Gordon-Schrödinger systems from relativistic to non-relativistic regimes 从相对论到非相对论的Klein-Gordon-Schrödinger系统的渐近一致指数型积分器
IF 1.3 4区 数学 Q2 Mathematics Pub Date : 2018-03-26 DOI: 10.1553/ETNA_VOL48S63
Simon Baumstark, Georgia Kokkala, Katharina Schratz
In this paper we propose asymptotic consistent exponential-type integrators for the Klein-Gordon-Schrodinger system. This novel class of integrators allows us to solve the system from slowly varying relativistic up to challenging highly oscillatory non-relativistic regimes without any step size restriction. In particular, our first- and second-order exponential-type integrators are asymptotically consistent in the sense of asymptotically converging to the corresponding decoupled free Schrodinger limit system.
本文提出了Klein-Gordon薛定谔系统的渐近一致指数型积分器。这类新型积分器使我们能够在没有任何步长限制的情况下求解从缓慢变化的相对论性到具有挑战性的高振荡非相对论性系统。特别地,我们的一阶和二阶指数型积分器在渐近收敛到相应的解耦自由薛定谔极限系统的意义上是渐近一致的。
{"title":"Asymptotic consistent exponential-type integrators for Klein-Gordon-Schrödinger systems from relativistic to non-relativistic regimes","authors":"Simon Baumstark, Georgia Kokkala, Katharina Schratz","doi":"10.1553/ETNA_VOL48S63","DOIUrl":"https://doi.org/10.1553/ETNA_VOL48S63","url":null,"abstract":"In this paper we propose asymptotic consistent exponential-type integrators for the Klein-Gordon-Schrodinger system. This novel class of integrators allows us to solve the system from slowly varying relativistic up to challenging highly oscillatory non-relativistic regimes without any step size restriction. In particular, our first- and second-order exponential-type integrators are asymptotically consistent in the sense of asymptotically converging to the corresponding decoupled free Schrodinger limit system.","PeriodicalId":50536,"journal":{"name":"Electronic Transactions on Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44363121","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}
引用次数: 5
Varying the s in Your s-step GMRES 改变你的s步GMRES中的s
IF 1.3 4区 数学 Q2 Mathematics Pub Date : 2018-01-01 DOI: 10.1553/ETNA_VOL47S206
David Imberti, J. Erhel
Krylov subspace methods are commonly used iterative methods for solving large sparse linear systems, however they suffer from communication bottlenecks on parallel computers. Therefore, $s$-step methods have been developed where the Krylov subspace is built block by block, so that $s$ matrix-vector multiplications can be done before orthonormalizing the block. Then Communication-Avoiding algorithms can be used for both kernels. This paper introduces a new variation on $s$-step GMRES in order to reduce the number of iterations necessary to ensure convergence, with a small overhead in the number of communications. Namely, we develop a $s$-step GMRES algorithm, where the block size is variable and increases gradually. Our numerical experiments show a good agreement with our analysis of condition numbers and demonstrate the efficiency of our variable $s$-step approach.
Krylov子空间方法是求解大型稀疏线性系统的常用迭代方法,但在并行计算机上存在通信瓶颈。因此,已经开发了$s$步方法,其中逐个块构建Krylov子空间,因此可以在对块进行正交规格化之前进行$s$矩阵向量乘法。然后,两个核都可以使用通信避免算法。为了减少保证收敛所需的迭代次数,在通信数量上的开销很小,本文引入了$ 5 $步GMRES的一种新变体。也就是说,我们开发了一个$s$步的GMRES算法,其中块大小是可变的,并且逐渐增加。我们的数值实验结果与我们的条件数分析结果很好地吻合,并证明了我们的变量$s$步方法的有效性。
{"title":"Varying the s in Your s-step GMRES","authors":"David Imberti, J. Erhel","doi":"10.1553/ETNA_VOL47S206","DOIUrl":"https://doi.org/10.1553/ETNA_VOL47S206","url":null,"abstract":"Krylov subspace methods are commonly used iterative methods for solving large sparse linear systems, however they suffer from communication bottlenecks on parallel computers. Therefore, $s$-step methods have been developed where the Krylov subspace is built block by block, so that $s$ matrix-vector multiplications can be done before orthonormalizing the block. Then Communication-Avoiding algorithms can be used for both kernels. This paper introduces a new variation on $s$-step GMRES in order to reduce the number of iterations necessary to ensure convergence, with a small overhead in the number of communications. Namely, we develop a $s$-step GMRES algorithm, where the block size is variable and increases gradually. Our numerical experiments show a good agreement with our analysis of condition numbers and demonstrate the efficiency of our variable $s$-step approach.","PeriodicalId":50536,"journal":{"name":"Electronic Transactions on Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"67352006","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}
引用次数: 0
On generalized iterated Tikhonov regularization with operator-dependent seminorms 算子相关半模的广义迭代Tikhonov正则化
IF 1.3 4区 数学 Q2 Mathematics Pub Date : 2018-01-01 DOI: 10.1553/ETNA_VOL47S73
D. Bianchi, M. Donatelli
We investigate the recently introduced Tikhonov regularization filters with penalty terms having seminorms that depend on the operator itself. Exploiting the singular value decomposition of the operator, we provide optimal order conditions, smoothing properties, and a general condition (with a minor condition of the seminorm) for the saturation level. Moreover, we introduce and analyze both stationary and nonstationary iterative counterparts of the generalized Tikhonov method with operator-dependent seminorms. We establish their convergence rate under conditions affecting only the iteration parameters, proving that they overcome the saturation result. Finally, some selected numerical results confirm the effectiveness of the proposed regularization filters.
我们研究了最近引入的吉洪诺夫正则化滤波器,其惩罚项具有依赖于算子本身的半正规。利用算子的奇异值分解,我们提供了饱和水平的最优阶条件、平滑性质和一般条件(含半正规的次要条件)。此外,我们还引入并分析了具有算子相关半模的广义Tikhonov方法的平稳迭代对应物和非平稳迭代对应物。在只影响迭代参数的条件下,建立了它们的收敛速度,证明了它们克服了饱和结果。最后,选定的一些数值结果验证了所提正则化滤波器的有效性。
{"title":"On generalized iterated Tikhonov regularization with operator-dependent seminorms","authors":"D. Bianchi, M. Donatelli","doi":"10.1553/ETNA_VOL47S73","DOIUrl":"https://doi.org/10.1553/ETNA_VOL47S73","url":null,"abstract":"We investigate the recently introduced Tikhonov regularization filters with penalty terms having seminorms that depend on the operator itself. Exploiting the singular value decomposition of the operator, we provide optimal order conditions, smoothing properties, and a general condition (with a minor condition of the seminorm) for the saturation level. Moreover, we introduce and analyze both stationary and nonstationary iterative counterparts of the generalized Tikhonov method with operator-dependent seminorms. We establish their convergence rate under conditions affecting only the iteration parameters, proving that they overcome the saturation result. Finally, some selected numerical results confirm the effectiveness of the proposed regularization filters.","PeriodicalId":50536,"journal":{"name":"Electronic Transactions on Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"67352120","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}
引用次数: 14
Low-Rank Solvers for Fractional Differential Equations 分数阶微分方程的低秩解
IF 1.3 4区 数学 Q2 Mathematics Pub Date : 2016-01-01 DOI: 10.17617/2.2270973
T. Breiten, V. Simoncini, M. Stoll
Many problems in science and technology can be cast using differential equations with both fractional time and spatial derivatives. To accurately simulate natural phenomena using this technology fine spatial and temporal discretizations are required, leading to large-scale linear systems or matrix equations, especially whenever more than one space dimension is considered. The discretization of fractional differential equations typically involves dense matrices with a Toeplitz or in the variable coefficient case Toeplitz-like structure . We combine the fast evaluation of Toeplitz matrices and their circulant preconditioners with state-of-the-art linear matrix equation methods to efficiently solve these problems, both in terms of CPU time and memory requirements. Additionally, we illustrate how these techniqes can be adapted when variable coefficients are present. Numerical experiments on typical differential problems with fractional derivatives in both space and time showing the effectiveness of the approaches are reported.
许多科学技术问题都可以用具有分数阶时间导数和空间导数的微分方程来解决。为了使用该技术精确地模拟自然现象,需要精细的空间和时间离散化,从而导致大规模的线性系统或矩阵方程,特别是当考虑多个空间维度时。分数阶微分方程的离散化通常涉及具有Toeplitz结构或在变系数情况下具有Toeplitz结构的密集矩阵。我们将Toeplitz矩阵及其循环前置条件的快速评估与最先进的线性矩阵方程方法相结合,以有效地解决这些问题,无论是在CPU时间和内存需求方面。此外,我们还说明了在存在可变系数时如何调整这些技术。对典型的含分数阶导数的微分问题在空间和时间上的数值实验表明了该方法的有效性。
{"title":"Low-Rank Solvers for Fractional Differential Equations","authors":"T. Breiten, V. Simoncini, M. Stoll","doi":"10.17617/2.2270973","DOIUrl":"https://doi.org/10.17617/2.2270973","url":null,"abstract":"Many problems in science and technology can be cast using differential equations with both fractional time and spatial derivatives. To accurately simulate natural phenomena using this technology fine spatial and temporal discretizations are required, leading to large-scale linear systems or matrix equations, especially whenever more than one space dimension is considered. The discretization of fractional differential equations typically involves dense matrices with a Toeplitz or in the variable coefficient case Toeplitz-like structure . We combine the fast evaluation of Toeplitz matrices and their circulant preconditioners with state-of-the-art linear matrix equation methods to efficiently solve these problems, both in terms of CPU time and memory requirements. Additionally, we illustrate how these techniqes can be adapted when variable coefficients are present. Numerical experiments on typical differential problems with fractional derivatives in both space and time showing the effectiveness of the approaches are reported.","PeriodicalId":50536,"journal":{"name":"Electronic Transactions on Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2016-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"67654664","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}
引用次数: 47
A MATRIX-FREE LEGENDRE SPECTRAL METHOD FOR INITIAL-BOUNDARY VALUE PROBLEMS 初边值问题的无矩阵Legendre谱法
IF 1.3 4区 数学 Q2 Mathematics Pub Date : 2016-01-01 DOI: 10.15496/PUBLIKATION-27275
Bernd Brumm, Emil Kieri
We study numerical methods for time-dependent partial differential equations describing wave propagation, primarily applied to problems in quantum dynamics governed by the time-dependent Schrodinger equation (TDSE). We consider both methods for spatial approximation and for time stepping. In most settings, numerical solution of the TDSE is more challenging than solving a hyperbolic wave equation. This is mainly because the dispersion relation of the TDSE makes it very sensitive to dispersion error, and infers a stringent time step restriction for standard explicit time stepping schemes. The TDSE is also often posed in high dimensions, where standard methods are intractable.The sensitivity to dispersion error makes spectral methods advantageous for the TDSE. We use spectral or pseudospectral methods in all except one of the included papers. In Paper III we improve and analyse the accuracy of the Fourier pseudospectral method applied to a problem with limited regularity, and in Paper V we construct a matrix-free spectral method for problems with non-trivial boundary conditions. Due to its stiffness, the TDSE is most often solved using exponential time integration. In this thesis we use exponential operator splitting and Krylov subspace methods. We rigorously prove convergence for force-gradient operator splitting methods in Paper IV. One way of making high-dimensional problems computationally tractable is low-rank approximation. In Paper VI we prove that a splitting method for dynamical low-rank approximation is robust to singular values in the approximation approaching zero, a situation which is difficult to handle since it implies strong curvature of the approximation space.
我们研究了描述波传播的时变偏微分方程的数值方法,主要应用于由时变薛定谔方程(TDSE)控制的量子动力学问题。我们考虑了空间逼近和时间步进两种方法。在大多数情况下,TDSE的数值解比求解双曲波动方程更具挑战性。这主要是因为TDSE的色散关系使其对色散误差非常敏感,并且对标准显式时间步进方案推导出严格的时间步长限制。TDSE也经常是高维的,在这里标准方法是难以处理的。对色散误差的敏感性使光谱方法对TDSE具有优势。除了一篇纳入的论文外,我们在所有论文中都使用了光谱或伪光谱方法。在第三篇论文中,我们改进并分析了傅里叶伪谱方法用于有限正则性问题的准确性,在第五篇论文中,我们构造了一个无矩阵谱方法用于具有非平凡边界条件的问题。由于其刚度,TDSE通常使用指数时间积分来求解。在本文中,我们使用指数算子分裂和Krylov子空间方法。本文严格证明了力梯度算子分裂方法的收敛性。使高维问题在计算上易于处理的一种方法是低秩近似。在第六篇论文中,我们证明了一种动态低秩逼近的分裂方法对于逼近于零的奇异值具有鲁棒性,这种情况由于意味着逼近空间的强曲率而难以处理。
{"title":"A MATRIX-FREE LEGENDRE SPECTRAL METHOD FOR INITIAL-BOUNDARY VALUE PROBLEMS","authors":"Bernd Brumm, Emil Kieri","doi":"10.15496/PUBLIKATION-27275","DOIUrl":"https://doi.org/10.15496/PUBLIKATION-27275","url":null,"abstract":"We study numerical methods for time-dependent partial differential equations describing wave propagation, primarily applied to problems in quantum dynamics governed by the time-dependent Schrodinger equation (TDSE). We consider both methods for spatial approximation and for time stepping. In most settings, numerical solution of the TDSE is more challenging than solving a hyperbolic wave equation. This is mainly because the dispersion relation of the TDSE makes it very sensitive to dispersion error, and infers a stringent time step restriction for standard explicit time stepping schemes. The TDSE is also often posed in high dimensions, where standard methods are intractable.The sensitivity to dispersion error makes spectral methods advantageous for the TDSE. We use spectral or pseudospectral methods in all except one of the included papers. In Paper III we improve and analyse the accuracy of the Fourier pseudospectral method applied to a problem with limited regularity, and in Paper V we construct a matrix-free spectral method for problems with non-trivial boundary conditions. Due to its stiffness, the TDSE is most often solved using exponential time integration. In this thesis we use exponential operator splitting and Krylov subspace methods. We rigorously prove convergence for force-gradient operator splitting methods in Paper IV. One way of making high-dimensional problems computationally tractable is low-rank approximation. In Paper VI we prove that a splitting method for dynamical low-rank approximation is robust to singular values in the approximation approaching zero, a situation which is difficult to handle since it implies strong curvature of the approximation space.","PeriodicalId":50536,"journal":{"name":"Electronic Transactions on Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2016-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"67154485","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}
引用次数: 1
Self-Generating and Efficient Shift Parameters in ADI Methods for Large Lyapunov and Sylvester Equations 大型Lyapunov和Sylvester方程的自生成和有效移位参数的ADI方法
IF 1.3 4区 数学 Q2 Mathematics Pub Date : 2014-01-01 DOI: 10.17617/2.2071065
P. Benner, P. Kürschner, J. Saak
Low-rank versions of the alternating direction implicit (ADI) iteration are popular and well estab- lished methods for the numerical solution of large-scale Sylvester and Lyapunov equations. Probably the biggest disadvantage of these methods is their dependence on a set of shift parameters that are crucial for fast convergence. Here we firstly review existing shift generation strategies that compute a number of shifts before the actual itera- tion. These approaches come with several disadvantages such as, e.g., expensive numerical computations and the difficulty to obtain necessary spectral information or data n eeded to initially setup their generation. Secondly, we propose two novel shift selection strategies with the motivation to resolve these issues at least partially. Both ap- proaches generate shifts automatically in the course of the ADI iterations. Extensive numerical tests show that one of these new approaches, based on a Galerkin projection onto the space spanned by the current ADI data, is superior to other approaches in the majority of cases both in terms of convergence speed and required execution time.
交替方向隐式迭代(ADI)的低秩迭代是求解大规模Sylvester和Lyapunov方程数值解的常用方法。这些方法最大的缺点可能是它们依赖于一组对快速收敛至关重要的移位参数。在这里,我们首先回顾了现有的位移生成策略,这些策略在实际迭代之前计算了许多位移。这些方法有几个缺点,例如,昂贵的数值计算和难以获得必要的光谱信息或数据,需要初始设置它们的生成。其次,我们提出了两种新颖的转移选择策略,至少部分地解决了这些问题。这两种方法都在ADI迭代过程中自动产生移位。大量的数值测试表明,其中一种基于Galerkin投影到当前ADI数据所跨越的空间的新方法,在大多数情况下,在收敛速度和所需的执行时间方面都优于其他方法。
{"title":"Self-Generating and Efficient Shift Parameters in ADI Methods for Large Lyapunov and Sylvester Equations","authors":"P. Benner, P. Kürschner, J. Saak","doi":"10.17617/2.2071065","DOIUrl":"https://doi.org/10.17617/2.2071065","url":null,"abstract":"Low-rank versions of the alternating direction implicit (ADI) iteration are popular and well estab- lished methods for the numerical solution of large-scale Sylvester and Lyapunov equations. Probably the biggest disadvantage of these methods is their dependence on a set of shift parameters that are crucial for fast convergence. Here we firstly review existing shift generation strategies that compute a number of shifts before the actual itera- tion. These approaches come with several disadvantages such as, e.g., expensive numerical computations and the difficulty to obtain necessary spectral information or data n eeded to initially setup their generation. Secondly, we propose two novel shift selection strategies with the motivation to resolve these issues at least partially. Both ap- proaches generate shifts automatically in the course of the ADI iterations. Extensive numerical tests show that one of these new approaches, based on a Galerkin projection onto the space spanned by the current ADI data, is superior to other approaches in the majority of cases both in terms of convergence speed and required execution time.","PeriodicalId":50536,"journal":{"name":"Electronic Transactions on Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2014-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"67654632","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}
引用次数: 73
Rational Interpolation Methods for Symmetric Sylvester Equations 对称Sylvester方程的有理插值方法
IF 1.3 4区 数学 Q2 Mathematics Pub Date : 2014-01-01 DOI: 10.17617/2.2060858
P. Benner, T. Breiten
We discuss low-rank approximation methods for large-scale symmetric Sylvester equations. Follow- ing similar discussions for the Lyapunov case, we introduce an energy norm by the symmetric Sylvester operator. Given a rank nr, we derive necessary conditions for an approximation being optimal with respect to this norm. We show that the norm minimization problem is related to an objective function based on the H2-inner product for sym- metric state space systems. This objective function leads to first-order optimality conditions that are equivalent to the ones for the norm minimization problem. We further propose an iterative procedure and demonstrate its efficiency by means of some numerical examples.
讨论了大规模对称Sylvester方程的低秩逼近方法。在对Lyapunov情况进行类似讨论之后,我们引入了对称Sylvester算子的能量范数。给定一个秩n,我们得到了关于这个范数的最优逼近的必要条件。我们证明了范数最小化问题与一个基于h2 -内积的对称状态空间系统的目标函数有关。这个目标函数导致一阶最优性条件,等同于范数最小化问题的条件。我们进一步提出了一种迭代方法,并通过一些数值算例证明了它的有效性。
{"title":"Rational Interpolation Methods for Symmetric Sylvester Equations","authors":"P. Benner, T. Breiten","doi":"10.17617/2.2060858","DOIUrl":"https://doi.org/10.17617/2.2060858","url":null,"abstract":"We discuss low-rank approximation methods for large-scale symmetric Sylvester equations. Follow- ing similar discussions for the Lyapunov case, we introduce an energy norm by the symmetric Sylvester operator. Given a rank nr, we derive necessary conditions for an approximation being optimal with respect to this norm. We show that the norm minimization problem is related to an objective function based on the H2-inner product for sym- metric state space systems. This objective function leads to first-order optimality conditions that are equivalent to the ones for the norm minimization problem. We further propose an iterative procedure and demonstrate its efficiency by means of some numerical examples.","PeriodicalId":50536,"journal":{"name":"Electronic Transactions on Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2014-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"67655019","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}
引用次数: 3
期刊
Electronic Transactions on Numerical Analysis
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1