Nonlinear Reduction using the Extended Group Finite Element Method

In this paper, we develop a nonlinear reduction framework based on our recently introduced extended group finite element method. By interpolating nonlinearities onto approximation spaces defined with the help of finite elements, the extended group finite element formulation achieves a noticeable reduction in the computational overhead associated with nonlinear finite element problems. However, the problem's size still leads to long solution times in most applications. Aiming to make real-time and/or many-query applications viable, we apply model order reduction and complexity reduction techniques in order to reduce the problem size and efficiently handle the reduced nonlinear terms, respectively. For this work, we focus on the proper orthogonal decomposition and discrete empirical interpolation methods. While similar approaches based on the group finite element method only focus on semilinear problems, our proposed framework is also compatible with quasilinear problems. Compared to existing methods, our reduced models prove to be superior in many different aspects as demonstrated in three numerical benchmark problems.


Introduction
Many nonlinear (finite element) problems are numerically expensive to solve with respect to computing time and storage requirements.In scenarios where many solutions need to be computed, such as for various parameter choices or initial conditions, the computational overhead quickly becomes a limiting factor.Therefore, efficient approaches that decrease the costs of solving the problem are vital.The group finite element method (GFEM) interpolates nonlinearities onto the finite element space used for the solution, which reduces overhead by allowing the time-consuming numerical integration to be precomputed and evaluating the nonlinear functions only at the nodes, see, for example, [10,12].Recently, we extended the GFEM in [18] to allow for better approximations and handle nonlinearities that cannot be directly interpolated onto the solution space.In particular, the extended group finite element method (EGFEM) allows for discontinuous nonlinearities, e.g., functions that depend on derivatives of the solution.It also provides more control in the approximation of the nonlinear functions, since general finite element spaces can be used instead of reusing the solution space.Additionally, we developed a tensor structure for quasilinear problems using EGFEM so that such problems can also be solved more efficiently than with the standard Galerkin approach (SGA).However, while the GFEM and EGFEM reduce the computational overhead, the size of the discrete problems is still large and impacts the solution time.Thus, projection-based model reduction techniques are crucial in order to further speed up the computations with the help of reduced order Date: June 7, 2021 † Corresponding author, email: tolle@uni-trier.de,phone: +49 651 201 3472 1 Universität Trier, FB IV -Mathematik, Lehrstuhl Modellierung und Numerik, D-54286 Trier, Germany.models (ROM), see, for example, [5,6,16,7] for more details on a variety of reduction approaches.After interpolating nonlinear functions onto finite element spaces via EGFEM, we employ the proper orthogonal decomposition (POD) method to reduce the problem's dimension with the help of a reduced basis for the solution.Finally, we apply complexity reduction through the discrete empirical interpolation method (DEIM) in order to further reduce the computational costs associated with evaluating the nonlinear terms.Our proposed approach takes advantage of the benefits of the EGFEM in comparison to the GFEM.These advantages deliver a framework that is easy to implement, directly applicable to quasilinear problems and more versatile than existing techniques.For example, the group proper orthogonal decomposition (GPOD) method presented in [11] and the finite element method with interpolated coefficients (FEIC) for ROMs obtained from POD from [20] are based on the GFEM.As such, these methods only focus on semilinear problems.Our approach -in contrast -is able to replicate the methods from [11] and [20], while also being able to handle problems outside their scope.The ROMs and cROMs attained using our approach prove to be powerful approximations of the full order models (FOM), while requiring only a fraction of the computational overhead.
The structure of this paper is as follows.In Section 2 we briefly present the problem setting and introduce the notation and conventions used in Section 3. Our three-step approach, which combines the EGFEM, POD and DEIM in order to attain efficient, reduced models, is described in Section 3. Finally, Section 4 illustrates the advantages of our method in comparison to established reduction techniques with the help of three benchmark problems.

Problem Setting
For this work, we develop a reduction scheme for nonlinear finite element problems.Let Ω ⊂ R d , d ∈ N, denote a polygonal domain with boundary Γ.Consider the following nonlinear problem: where µ ∈ P denotes some parameter setting.Since our focus is on the numerical solution of the finite element problem based on (2.1), we assume that the scalar-valued functions a, c and q are sufficiently regular and ensure the existence of a weak solution.Neumannand Robin-type boundary conditions result in extra terms in the weak formulation but otherwise present no additional difficulties.Inhomogeneous Dirichlet conditions require purely technical adaptations to the problem formulation, while the ideas remain the same.Note that time-dependent problems can be handled analogously as seen in one of the examples considered in Section 4. The finite element method introduces a finite-dimensional approximation space V h using a triangulation of Ω. Assume that a nodal basis {φ 1 , . . ., φ Nu } exists with φ i (x j ) = δ ij for all degrees of freedom x j ∈ Ω.We are interested in efficiently computing an approximate solution where the coefficients u = (u 1 , . . ., u Nu ) T solve the weak formulation of (2.1) on V h : (2.3) Throughout this paper, we refer to the direct solution of (2.2) as the standard Galerkin approach (SGA).Many practical applications result in a large number of variables, i.e., N u 1.At the same time, the discrete nonlinear problem (2.2) is solved numerically through iterative methods, which require the repeated integration of the terms in (2.3).This is often time-intensive because of the computational overhead associated with the repeated assembly, i.e., numerical integration in (2.3) for each iteration.In order to mitigate the computational costs, we investigate a combination of reduction techniques along with the EGFEM in order to derive efficient reduced models.
Before describing our reduction approach in detail, we briefly introduce the notation that is used throughout Section 3 and Section 4: Notation 2.1 (Tensor Calculus).Vectors, matrices and third-order tensors are denoted as For the manipulation of the third-order tensors, we use the following rules for the single and double contractions: where the single contraction eliminates one degree of freedom and results in a matrix, while the double contraction reduces the order of the tensor by two and results in a vector.In contrast, the tensor product ⊗ combines two vectors into a matrix through For the reduction, we also use the n-mode product × n for n ∈ {1, 2, 3}: where the result is a third-order tensor with the n-th dimension changed according to the matrix used in the product.The definition of the n-mode product is based on [15], which we refer the reader to for more details.

Reduction Approach
Model order reduction techniques serve a vital role in efficient solution strategies.Many nonlinear reduction approaches rely on data-driven methods.However, it has been shown that model order reduction alone is not effective, since nonlinear functions still depend on the full problem dimension, see, for example, [9] for details.Therefore, complexity reduction in the form of (D)EIM [4,9] or empirical cubature [14] is necessary in order to efficiently evaluate the reduced nonlinear terms.In this paper, we focus on POD for the model reduction and DEIM for the complexity reduction.In combination with the EGFEM, our scheme delivers a powerful tool for the reduction of nonlinear finite element models that is more versatile than the approaches described in [11,20].For example, discontinuous nonlinearities can be handled efficiently with the help of EGFEM.When the DEIM is applied to point-wise schemes, such as finite difference discretizations, the complexity reduction is non-intrusive.This advantage is carried over for the reformulation via GFEM and EGFEM, since the nonlinearities are evaluated point-wise.In contrast, applying the DEIM to the SGA requires modifications to the finite element implementation in order to ensure that the complexity reduction actually improves the evaluation of the reduced nonlinear terms.Furthermore, the tensor structure introduced in [18] allows for a straight-forward reduction of the quasilinear terms.As such, our approach can be seen as an extension of [11,20], which uses the GFEM to approximate nonlinear terms before reducing with POD and optionally DEIM.The reduction approach can be separated into the following three steps, which are explained in more detail in what follows: (1) apply the EGFEM in order to rewrite nonlinear terms, (2) compute projection matrices using POD in order to derive a ROM, (3) use the DEIM to attain a cROM, which efficiently handles reduced nonlinear terms.
3.1.Extended Group Finite Element Method.Most of the computational overhead for nonlinear finite element problems lies in the assembly of nonlinear forms.Because the problems must be solved using iterative methods, such as Newton's methods, these forms have to be assembled in every iteration.Afterwards, linear systems are solved and the approximate solution is updated.In order to reduce the overhead and consequently speed up the computations, the (extended) group finite element method interpolates nonlinearities onto finite element spaces.The major difference between the GFEM and EGFEM is the choice of the approximation space for the nonlinear functions, i.e., GFEM interpolates nonlinear terms onto the solution space V h , while EGFEM introduces an additional finite element space W h that can be tailored to the nonlinearity.As such, EGFEM encapsulates GFEM for W h = V h but also allows for more accurate approximations and the treatment of more general nonlinearities, such as functions that depend on derivatives of the solution.For a more detailed investigation into the EGFEM including guidelines and numerical studies, we refer the reader to [18].
The basic idea is to introduce a finite element approximation space W h with a nodal basis {η f 1 , . . ., η f N f } that is used to interpolate the nonlinear function f (u h , ∇u h ; µ).Instead of evaluating the nonlinear terms in (2.2) in every iteration, we approximate The coefficients f j (u; µ), j = 1, . . ., N f , are determined by evaluating the function at the degrees of freedom x f j , i.e., f j (u; µ) In the case of GFEM, we have the special choice of W h = V h and, therefore, η f i = φ i , i = 1, . . ., N u , with the coefficients evaluated at the same degrees of freedom as the solution.For the nonlinear terms in (2.3), this procedure results in the approximations where the integrals can be precomputed and result in the third-order tensor K a and matrix L c , respectively.By introducing these approximations to (2.2), the EGFEM solves the following problem: , where the coefficient vectors a and c are defined point-wise.This leads to the introduction of the functions as well as the interpolation operators Π f u : R Nu → R N f and Π f ∇u : R Nu → R d×N f , which are used to evaluate the function u h and its gradient ∇u h at the degrees of freedom x f = (x f 1 , . . ., x f N f ) T associated to the approximation space W f h for f ∈ {a, c}.Remark 3.1.Note that the extended group finite element approximation can also be used for nonlinear boundary terms.In this case, the approximation space only has to be defined on the boundary, which results in an additional dimension reduction.Furthermore, it can also be applied to source terms in time-dependent problems, which have to be reassembled in each time step.

3.2.
Proper Orthogonal Decomposition.The next major reduction in complexity is achieved by applying the proper orthogonal decomposition technique on the spatial dependence of the system.This means that the size of the discrete problem is reduced with the help of orthogonal projections that select the most important modes of information with the help of a set of snapshots, i.e. a collection of solutions at different times and/or different parameter settings, see, for example, [16].The snapshots are collected in the matrix where n s denotes the number of snapshots and u k represents the solution in a specific configuration, e.g., for a fixed parameter value and/or at a specific time.The singular value decomposition (SVD) of Y delivers with U ∈ R Nu×Nu and W ∈ R ns×ns .Then, the projection matrix V u ∈ R Nu×nu consists of the first n u N u columns from the left projection U.For our problem, the snapshot matrix consists of the coefficient vectors u at discrete parameter settings, i.e., u i := u(µ i ) for 1 ≤ i ≤ n s .After computing the SVD of Y and determining the projection matrix V u , the reduced-order variable u r ∈ R nu with n u N u approximates the full-order variable u through The projection-based reduction of u in the SGA (2.2) and EGFEM (3.4) results in the following respective ROMs: with q r (µ) = V T u q(µ).Irregardless of which reduced model is considered, a major hindrance is the evaluation of the nonlinear functions.This is because these functions require the prolongation of the reduced variables to the full order dimension, evaluation of the function and finally the restriction back to the reduced order dimension.Therefore, the evaluation of the nonlinearities still depends on the size of the full order model.In light of this, we use the DEIM in order to remove the dependence on the full order dimension.

Discrete Empirical Interpolation Method.
The difficulty with nonlinear reduction is that the evaluation of the nonlinear functions depends on the full-order dimension.In [9], the discrete empirical interpolation method solves this problem via an interpolation approach, which only uses a few components of the nonlinear function as an approximation.
In the scope of this work, we consider f (µ) := f (V u u r ; µ), where f comes from either the nonlinear terms K(a, u; µ) and (c, u; µ) in the case of the SGA or the point-wise evaluations a(u; µ) and c(u; µ) in the EGFEM setting for every µ ∈ P. Note that the matrices K(a, u; µ) ∈ R Nu×Nu are treated as vectors in R N 2 u and reduced as described here.In this case, one speaks of the matrix discrete empirical interpolation method (MDEIM), see [8,21].Similar to the reduction of u, the form of the approximation of f is given through where τ is either a parameter setting µ or a time t, with V f ∈ R N f ×n f and f r (τ ) ∈ R n f .However, the system f (τ ) = V f f r (τ ) is overdetermined.Therefore, in order to get a welldefined approximation, which is efficient to evaluate, n f unique rows must be selected in order to determine f r .This is accomplished with the help of a selection matrix where e i denotes the i-th canonical basis vector in R N f .Then, the coefficient vector f r is uniquely determined through with the DEIM operator D f , where the projection matrix V f and the selection matrix P f still need to be specified.POD is applied to the nonlinear snapshots which results in the projection matrix V f , while the selection matrix P f is chosen inductively using a Greedy-type algorithm from [9].The DEIM-algorithm based on our notation is presented in Algorithm 3.2, where the first two steps compute the projection matrix V f via POD and the rest of the algorithm selects the rows for P f via a Greedy algorithm.
For our model problem, the DEIM is applied to (3.5) and (3.6) in order to efficiently evaluate the nonlinear terms.Beginning with the SGA-ROM in (3.5), DEIM approximates the nonlinear functions K(a, u; µ) and (c, u; µ).This leads to the following SGA-cROM: where the functions T and T −1 denote the transformation of a matrix to a vector (columnwise) and vice versa.For the EGFEM-ROM in (3.6), we attain the EGFEM-cROM: The remarkable difference between (3.7) and (3.8) is the "hidden" overhead contained in (3.7).While the numerical integration in (3.8) is only performed once in the computation of K a and L c , the matrix K(a, V u u r ; µ) and vector (c, V u u r ; µ) require assembly for each evaluation.In light of this, an effective implementation of (3.7) requires a modification of the assembly process in order to take full advantage of the complexity reduction from the DEIM, so that information that is immediately disregarded by the selection operator P T f , f ∈ {K, }, is not generated at all.In contrast, the evaluations of a and c in (3.8), which are point-wise, can take full advantage of the DEIM.Similar to the DEIM for finite difference methods in [9], the selection operator P T f , f ∈ {a, c}, can be applied directly to the arguments of f .This means that the nonlinear functions are only evaluated at the n f points chosen in the DEIM algorithm.

Numerical Results
In this section, we focus on the quality of the (complexity) reduced models in comparison with the FOMs.We consider three examples that focus on different aspects of the considered method(s).The first example is taken from [9] and compares the SGA, GFEM and EGFEM for a semilinear problem.For the second example, we solve the viscous Burgers' equation in two dimensions, which is also considered in [11].This illustrates how our approach can also be applied to time-dependent problems.Finally, the third example highlights the superiority of our approach in comparison to the existing methods based on the GFEM.Note that all examples assume homogeneous Dirichlet boundary conditions.For all the examples, we use piecewise polynomial functions for the approximation space, where P k denotes the piecewise polynomial functions of degree k ∈ N 0 .The solution space V h uses piecewise linear functions, i.e., P 1 .The interpolation of the nonlinear function onto linear finite elements corresponds to the GFEM, since V h = W h .In contrast, we consider the approximation spaces W h using piecewise constant and quadratic functions, which are denoted with P 0 and P 2 , respectively, for the EGFEM.The ROMs are attained using POD, while the cROMs use POD in conjunction with the DEIM.The size of the complexity reduction, i.e., the number of DEIM modes retained, is provided in parenthesis, where applicable.Finally, all computations are performed in Matlab (R2020a) on a server with a Xeon E5-2699V4 processor and 756 GB of memory.The triangulations are generated using the software Gmsh [13], while an external toolbox for tensors [2,3] is used for sparse third-order tensors.The stationary examples are solved using Matlab's fsolve function with all tolerances set to 10 −8 and using the initial guess u 0 = 0.The instationary (second) example is solved using Matlab's ode15s method with backward differentiation formula, where the absolute and relative tolerances are set to 10 −8 and 10 −6 , respectively.In order to ensure optimal performance of Matlab's solvers, we provide the Jacobian matrix analytically.Finally, the runtimes refer to the average time elapsed as measured by Matlab's tic-toc functionality.
4.1.Parameter-dependent Semilinear Equation.The first example is a nonlinear parameter-dependent problem with a homogeneous Dirichlet boundary condition from [9].The equation considered is given by with the parameter setting µ = (µ 1 , µ 2 ) T ∈ P = [0.01,10] 2 .With respect to (2.1), we have where the exponential function acts on each element of its vector-valued argument and diag(•) denotes the matrix with the argument on its diagonal.The reduced models have similar derivatives with the inclusion of the chain rule, which results in additional scaling with the projection matrix V u .Remark 4.1.Note that the constant factor in c could also be shifted to the source term q.However, the current choice matches with [9] and ensures that c ≡ 0 on the boundary.
Following [9], the solution and nonlinearity snapshots are generated by sampling the parameter space P = [0.01,10] 2 uniformly with 144 samples.The singular values belonging to the snapshot matrices are shown in Figure 1.Up to a constant shift, the singular values of the nonlinear snapshot matrices show almost the same decay rate.This shift in the singular values of the nonlinear snapshots between the SGA and EGFEM problems is most likely due to the additional integration step involved in the SGA.We evaluate the ROMs and cROMs on a uniform sampling of the parameter space using 225 samples, which do not coincide with the training parameters.In order to compare the different choices of W h in the EGFEM and because we do not have an analytical solution available, we consider the average absolute error using the Euclidean norm with respect to two different references.The average absolute errors of the reduced solutions with respect to the full order solution are shown in the left column of Figure 2.These errors are comparable to those shown in [9].The right column in Figure 2 shows the error with respect to the SGA-FOM.This choice as a reference is motivated by the fact that the SGA coupled with a fourth-order numerical quadrature rule for the nonlinear term should be the most accurate approximation of the exponential function in the nonlinearity.In contrast, the GFEM and EGFEM use lower-degree polynomial approximations of the nonlinear function c in order to simplify the problem.From the right column of Figure 2, we see that the GFEM (P 1 ) approximation as well as the EGFEM approximation with piece-wise constant functions (P 0 ) poorly describe the reference SGA solution, while the EGFEM with piecewise quadratic functions (P 2 ) captures the behavior of the SGA solution and even performs better in some configurations, e.g., when using 19 DEIM modes.Finally, we see noticeable reductions in the computational overhead for the EGFEM formulations in comparison to the SGA, as shown in Figure 3.The SGA-ROMs require the complexity reduction via DEIM (SGA-cROMs) in order to even show an improvement in the computing times.Our approach with the EGFEM -in contrast to the SGA -turns out to be even faster than most of the SGA-cROMs.Additionally, DEIM reduces the average time per iteration by an additional order of magnitude for the models using the EGFEM.Interestingly, as seen in Figure 3, the average time per iteration for the GFEM-and EGFEM-cROMs appears to be independent of the DEIM dimension.This is most likely due to the fact that Matlab efficiently handles matrix-vector multiplications, and the difference between the DEIM dimensions is negligible at the considered scale.4.2.Burgers' Equation.For the second example, we consider the viscous Burgers' equation in two dimensions with a homogeneous Dirichlet boundary condition: on the domain x = (x 1 , x 2 ) ∈ [0, 1] 2 and time interval t ∈ [0, 10] with u(x, 0) = u 0 (x) and a constant viscosity ν = 1/100.In the spirit of [11], we choose the source q via the method of manufactured solutions and the initial condition u 0 , so that the solution is given through Here, we have a(x) ≡ ν.We modify the weak formulation by partially integrating the terms with derivatives with respect to x 1 and x 2 .This leads to the following new term for the SGA: For the finite element discretization, this results in a term quite similar to (c, u; µ): where d(x, u, ∇u; t) = u(x, t) 2 in this example.However, since d is quadratic, we also consider a SGA using a multi-linear form as discussed in [18].This leads to a simple reduction of the problem without needing complexity reduction.The SGA and multi-linear standard Galerkin approach (ML-SGA) solve the following semi-discrete problems: respectively, where E and K denote the standard mass and stiffness matrices and the thirdorder tensor M is defined through The ML-SGA has the same benefits as the EGFEM, i.e., the cost-intensive numerical integration only has to be performed once.In contrast, the EGFEM solves the following semi-discrete problem:   where d is defined in a similar manner to a and c in (3.4).In addition to the reduced terms in (3.5) and (3.6), we also introduce the reduced mass and stiffness matrices as well as the reduced tensor for the ML-SGA, which read: The advantage of the ML-SGA is that the nonlinearity no longer depends on the full dimension, i.e., M r ∈ R nu×nu×nu , which means that the DEIM is not needed for this formulation.However, the ML-SGA is only possible for polynomial nonlinearities and only practical for low-degree polynomials.
We begin by examining the singular values of the snapshot matrices.The snapshots are generated by solving the FOMs and evaluating the solutions on an equidistant time grid with step size δt = 10 −2 .Figure 4 displays the singular values associated with the respective snapshot matrix.We see that the singular values of the various formulations are qualitatively similar.For this example, we choose the number of DEIM modes to be equal to the number of POD modes.Because the GFEM loses accuracy with respect to the energy norm in comparison to the SGA, see [10], we use the energy norm to compute the errors.Relative errors are shown in Figure 5.The discretization error compares the analytical solution û to the solution of the FOM for increasingly finer triangulations.The reduction error compares the FOMs with N u = 4, 225 to the ROMs.The solid lines without markers show the projection error, i.e., the error arising from restricting and prolongating the fullorder solution.Finally, the total error compares the ROMs to the analytical solution.From the reduction error in Figure 5, we see that the SGA-cROMs perform worse than  the other cROMs as the size of the reduced models increases.However, when considering the total error, which includes the discretization error, we see that the differences become less obvious.We do note that the EGFEM models perform slightly better than the GFEM models, which shows that our method can be better than the existing approaches using GFEM.To summarize, this problem is most easily reduced using the ML-SGA, which is applicable for low-degree polynomial nonlinearities.However, the EGFEM using quadratic (P 2 ) elements is analytically equivalent to the SGA, while also being simpler to reduce and much more efficient in terms of the computational costs.
Remark 4.2 (Instability).When using DEIM, we have found that certain choices of POD and DEIM sizes may lead to unstable approximations.In particular, the cROM may be instable, when the DEIM dimension is smaller than the POD dimension.This phenomenon has also been reported for other models, see, for example, [1] and [17].
4.3.Minimal Surface Equation.Finally, we consider an example for which the GFEM is not applicable.The minimal surface equation with a source term [19] is given through The domain Ω is a polygonal approximation of the unit disk.The nonlinear function a is an example of a nonlinearity, that cannot be approximated using the GFEM.Since V h ⊂ H 1 0 , the gradient ∇u h of u h ∈ V h is not well-defined on the nodes.This means that the interpolation of a(∇u h ) onto V h is not possible and as such any reduction methods based on the GFEM are not applicable.Therefore, we only consider EGFEM with a piecewise constant (P 0 ) approximation of the nonlinearity.Interestingly, this is an exact reformulation of the nonlinearity, since ∇u h is piecewise constant.Similar to the first example, we sample the parameter space P uniformly at 144 points, while the errors are computed on a uniform sampling with 225 points, so that the ROMs and cROMs are evaluated on parameter settings not included in the training set.The singular values for the solution and nonlinear snapshot matrices are shown in Figure 6.From the singular values, we expect good approximations using few modes.The average absolute error over all samples is shown in Figure 7.The errors of both methods are comparable, although the EGFEM shows smaller errors as the size of the cROM increases.We also compare the average time per iteration in Figure 8. Once again, we see that the ROM without complexity reduction takes as long as the FOM for the SGA.The cROMs show a noticeable reduction in the computational overhead associated with evaluating the  nonlinearity.When comparing the two approaches, we see that the EGFEM is noticeably faster than the SGA-counterparts.In summary, the EGFEM leads to faster models that are easier to reduce and prove to be more accurate than the standard approach.

Conclusion
In this work, we combine our recently introduced EGFEM with nonlinear reduction techniques.The presented results use POD and DEIM.However, the choice is illustrative and not fixed.Note that our approach is also applicable with other model order and complexity reduction methods, such as the reduced basis method in place of POD and the empirical interpolation method instead of the DEIM.Due to the nature of the EGFEM, we are able to effectively reduce quasilinear problems as well as handle discontinuous nonlinearities that could not be treated using the methods presented in [11] and [20].At the same time, the same approaches are embedded in our framework, since EGFEM is a generalization of the GFEM, which results for the special choice W h = V h .The superior performance of our method in comparison to existing methods has been numerically demonstrated for different applications in Section 4. Furthermore, it is less intrusive in terms of an efficient implementation of the complexity reduction in comparison to standard reduction approaches.

Figure 1 .
Figure 1.Singular values associated to the snapshot matrices for the solution (left) and the nonlinearity (right).

Figure 3 .
Figure 3. Average time to solve the respective system.

Figure 4 .
Figure 4. Singular values associated to the snapshot matrix containing the full-order solution (left) and nonlinearity (right).

Figure 5 .
Figure 5. Discretization error with respect to mesh refinement (top, left), where the colored marker denotes the mesh used for reduction.The reduction error (top, right) compares the reduced order solutions to the full order solutions, while the total error (bottom, left) compares the reduced order solutions to the analytical solution and includes the discretization error.

Figure 6 .
Figure 6.Singular values associated to the snapshot matrices for the solution (left) and the nonlinearity (right).

Figure 7 .Figure 8 .
Figure 7. Average absolute error over all samples between the reduced order and full order models for the SGA (left) and EGFEM (right).The errors are computed using the Euclidean (point-wise) norm.