A TENSOR APPROXIMATION METHOD BASED ON IDEAL MINIMAL RESIDUAL FORMULATIONS FOR THE SOLUTION OF HIGH-DIMENSIONAL PROBLEMS

被引:27
作者
Billaud-Friess, M. [1 ]
Nouy, A. [1 ]
Zahm, O. [1 ]
机构
[1] Univ Nantes, Ecole Cent Nantes, CNRS, UMR 6183,GeM, F-44321 Nantes 3, France
来源
ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS-MODELISATION MATHEMATIQUE ET ANALYSE NUMERIQUE | 2014年 / 48卷 / 06期
关键词
High-dimensional problems; nonlinear approximation; low-rank approximation; proper generalized decomposition; minimal residual; stochastic partial differential equation; DECOMPOSITION; CONVERGENCE;
D O I
10.1051/m2an/2014019
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
In this paper, we propose a method for the approximation of the solution of high-dimensional weakly coercive problems formulated in tensor spaces using low-rank approximation formats. The method can be seen as a perturbation of a minimal residual method with a measure of the residual corresponding to the error in a specified solution norm. The residual norm can be designed such that the resulting low-rank approximations are optimal with respect to particular norms of interest, thus allowing to take into account a particular objective in the definition of reduced order approximations of high-dimensional problems. We introduce and analyze an iterative algorithm that is able to provide an approximation of the optimal approximation of the solution in a given low-rank subset, without any a priori information on this solution. We also introduce a weak greedy algorithm which uses this perturbed minimal residual method for the computation of successive greedy corrections in small tensor subsets. We prove its convergence under some conditions on the parameters of the algorithm. The proposed numerical method is applied to the solution of a stochastic partial differential equation which is discretized using standard Galerkin methods in tensor product spaces.
引用
收藏
页码:1777 / 1806
页数:30
相关论文
共 44 条
[1]   A new family of solvers for some, classes of multidimensional partial differential equations encountered in kinetic theory modeling of complex fluids [J].
Ammar, A. ;
Mokdad, B. ;
Chinesta, F. ;
Keunings, R. .
JOURNAL OF NON-NEWTONIAN FLUID MECHANICS, 2006, 139 (03) :153-176
[2]   On the Convergence of a Greedy Rank-One Update Algorithm for a Class of Linear Systems [J].
Ammar, A. ;
Chinesta, F. ;
Falco, A. .
ARCHIVES OF COMPUTATIONAL METHODS IN ENGINEERING, 2010, 17 (04) :473-486
[3]   Adaptive Near-Optimal Rank Tensor Approximation for High-Dimensional Operator Equations [J].
Bachmayr, Markus ;
Dahmen, Wolfgang .
FOUNDATIONS OF COMPUTATIONAL MATHEMATICS, 2015, 15 (04) :839-898
[4]   A projection method to solve linear systems in tensor format [J].
Ballani, Jonas ;
Grasedyck, Lars .
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS, 2013, 20 (01) :27-43
[5]   Algorithms for numerical analysis in high dimensions [J].
Beylkin, G ;
Mohlenkamp, MJ .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2005, 26 (06) :2133-2159
[6]  
Cances E., 2012, ARXIV12106688V1
[7]   CONVERGENCE OF A GREEDY ALGORITHM FOR HIGH-DIMENSIONAL CONVEX NONLINEAR PROBLEMS [J].
Cances, Eric ;
Ehrlacher, Virginie ;
Lelievre, Tony .
MATHEMATICAL MODELS & METHODS IN APPLIED SCIENCES, 2011, 21 (12) :2433-2467
[8]   A Short Review on Model Order Reduction Based on Proper Generalized Decomposition [J].
Chinesta, Francisco ;
Ladeveze, Pierre ;
Cueto, Elias .
ARCHIVES OF COMPUTATIONAL METHODS IN ENGINEERING, 2011, 18 (04) :395-404
[9]   ADAPTIVITY AND VARIATIONAL STABILIZATION FOR CONVECTION-DIFFUSION EQUATIONS [J].
Cohen, Albert ;
Dahmen, Wolfgang ;
Welper, Gerrit .
ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS-MODELISATION MATHEMATIQUE ET ANALYSE NUMERIQUE, 2012, 46 (05) :1247-1273
[10]   ADAPTIVE PETROV-GALERKIN METHODS FOR FIRST ORDER TRANSPORT EQUATIONS [J].
Dahmen, Wolfgang ;
Huang, Chunyan ;
Schwab, Christoph ;
Welper, Gerrit .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2012, 50 (05) :2420-2445