The proper generalized decomposition (PGD) is an iterative numerical method for solving boundary value problems (BVPs), that is, partial differential equations constrained by a set of boundary conditions, such as the Poisson's equation or the Laplace's equation.

The PGD algorithm computes an approximation of the solution of the BVP by successive enrichment. This means that, in each iteration, a new component (or mode) is computed and added to the approximation. In principle, the more modes obtained, the closer the approximation is to its theoretical solution. Unlike POD principal components, PGD modes are not necessarily orthogonal to each other.

By selecting only the most relevant PGD modes, a reduced order model of the solution is obtained. Because of this, PGD is considered a dimensionality reduction algorithm.


The proper generalized decomposition is a method characterized by

  1. a variational formulation of the problem,
  2. a discretization of the domain in the style of the finite element method,
  3. the assumption that the solution can be approximated as a separate representation and
  4. a numerical greedy algorithm to find the solution.[1][2]

Variational formulation

The most implemented variational formulation in PGD is the Bubnov-Galerkin method,[3][4] although other implementations exist.[5][3]

Domain discretization

The discretization of the domain is a well defined set of procedures that cover (a) the creation of finite element meshes, (b) the definition of basis function on reference elements (also called shape functions) and (c) the mapping of reference elements onto the elements of the mesh.

Separate representation

PGD assumes that the solution u of a (multidimensional) problem can be approximated as a separate representation of the form

where the number of addends N and the functional products X1(x1), X2(x2), ..., Xd(xd), each depending on a variable (or variables), are unknown beforehand.

Greedy algorithm

The solution is sought by applying a greedy algorithm, usually the fixed point algorithm, to the weak formulation of the problem. For each iteration i of the algorithm, a mode of the solution is computed. Each mode consists of a set of numerical values of the functional products X1(x1), ..., Xd(xd), which enrich the approximation of the solution. Due to the greedy nature of the algorithm, the term 'enrich' is used rather than 'improve', since some modes may actually worsen the approach. The number of computed modes required to obtain an approximation of the solution below a certain error threshold depends on the stop criterium of the iterative algorithm.


PGD is suitable for solving high-dimensional problems, since it overcomes the limitations of classical approaches. In particular, PGD avoids the curse of dimensionality, as solving decoupled problems is computationally much less expensive than solving multidimensional problems.

Therefore, PGD enables to re-adapt parametric problems into a multidimensional framework by setting the parameters of the problem as extra coordinates:

where a series of functional products K1(k1), K2(k2), ..., Kp(kp), each depending on a parameter (or parameters), has been incorporated to the equation.

In this case, the obtained approximation of the solution is called computational vademecum: a general meta-model containing all the particular solutions for every possible value of the involved parameters.[6]

Sparse Subspace Learning

This section may be confusing or unclear to readers. Please help clarify the section. There might be a discussion about this on the talk page. (April 2021) (Learn how and when to remove this template message)

The Sparse Subspace Learning (SSL) method leverages the use of hierarchical collocation to approximate the numerical solution of parametric models. With respect to traditional projection-based reduced order modeling, the use of a collocation enables non-intrusive approach based on sparse adaptive sampling of the parametric space. This allows to recover the lowdimensional structure of the parametric solution subspace while also learning the functional dependency from the parameters in explicit form. A sparse low-rank approximate tensor representation of the parametric solution can be built through an incremental strategy that only needs to have access to the output of a deterministic solver. Non-intrusiveness makes this approach straightforwardly applicable to challenging problems characterized by nonlinearity or non affine weak forms.[7]


  1. ^ Amine Ammar; Béchir Mokdad; Francisco Chinesta; Roland Keunings (2006). "A New Family of Solvers for Some Classes of Multidimensional Partial Differential Equations Encountered in Kinetic Theory Modeling of Complex Fluids". Journal of Non-Newtonian Fluid Mechanics.
  2. ^ Amine Ammar; Béchir Mokdad; Francisco Chinesta; Roland Keunings (2007). "A new family of solvers for some classes of multidimensional partial differential equations encountered in kinetic theory modelling of complex fluids. Part II: Transient simulation using space-time separated representations". Journal of Non-Newtonian Fluid Mechanics.
  3. ^ a b Croft, Thomas Lloyd David (2015-04-09). Proper generalised decompositions: theory and applications (phd thesis). Cardiff University.
  4. ^ Chinesta, Francisco; Keunings, Roland; Leygue, Adrien (2014). The Proper Generalized Decomposition for Advanced Numerical Simulations: A Primer. SpringerBriefs in Applied Sciences and Technology. Springer International Publishing. ISBN 978-3-319-02864-4.
  5. ^ Aguado, José Vicente (18 Nov 2018). "Advanced strategies for the separated formulation of problems in the Proper Generalized Decomposition framework".
  6. ^ Francisco Chinesta, Adrien Leygue, Felipe Bordeu, Elías Cueto, David Gonzalez, Amine Ammar, Antonio Huerta (2013). "PGD-Based Computational Vademecum for Efficient Design, Optimization and Control". Archives of Computational Methods in Engineering.((cite journal)): CS1 maint: multiple names: authors list (link)
  7. ^ Borzacchiello, Domenico; Aguado, José V.; Chinesta, Francisco (April 2019). "Non-intrusive Sparse Subspace Learning for Parametrized Problems". Archives of Computational Methods in Engineering. 26 (2): 303–326. doi:10.1007/s11831-017-9241-4. ISSN 1134-3060. S2CID 126121268.