An Application of a Modiﬁed Gappy Proper Orthogonal Decomposition on Complexity Reduction of Allen-Cahn Equation

: This work considers model reduction techniques that can substantially decrease computational cost in simulating parmetrized Allen–Cahn equation. We first employ the proper orthogonal decomposition (POD) approach to reduce the number of unknowns in the full-order discretized system. Since POD cannot reduce the computational complexity of nonlinearity in Allen–Cahn equation, we also apply discrete empirical interpolation method (DEIM) to approximate the nonlinear term for a substantial reduction in overall simulation time. However, in general, the POD-DEIM approach is less accurate than the POD approach, since it further approximates the nonlinear term. To increase the accuracy of the POD-DEIM approach, this work introduces an extension of the DEIM approximation based on the concept of Gappy POD (GPOD), which is optimal in the least-squares sense. The POD-GPOD approach is tested and compared with the POD and POD-DEIM approaches on Allen–Cahn equation for both cases of fixed parameter value and varying parameter values. The modified GPOD approximation introduced in this work is demonstrated to improve accuracy of DEIM without sacrificing too much efficiency on the computational speedup, e.g., in one of our numerical tests, the POD-GPOD approach provides an approximate solution to the parmetrized Allen–Cahn equation 200 times faster than the full-order system with average error of order O ( 10 − 4 ) . The POD-GPOD approach is therefore shown to be a promising technique that compromises between the accuracy of POD approach and the efficiency of POD-DEIM approach.


Introduction
In the mathematical and computational field, model order reduction (MOR) is a technique for constructing efficient lower dimensional models for reducing the computational cost of large scale dynamical systems from partial differential equations (PDEs) in numerical simulation.Many modern mathematical models are useful for describing physical phenomena in real life.However, the process for solving the numerical solution of these mathematical models often has high computational complexity.To overcome this problem, A reduced order model can be used to approximate the original system while preserving the main properties and containing the important features of the original model.This paper considers nonlinear partial differential equation called Allen-Cahn equation, which was first introduced by Cahn and Allen to describe the motion of anti-phase boundaries in metallic alloys [1].This equation is used in mathematical physics, to describe the process of phase separation.The Allen-Cahn equation has been widely used in many applications, such as image analysis [2,3], crystal growth [4], and motion by mean curvature flow [5].In particular, it has become a basic model equation for the diffuse interface approach developed to study phase transitions and interfacial dynamics in material science [6].
In a general model of Allen-Cahn equation used for describing phase transition in materials science, there is an important parameter called an "interaction length" [5].The change of this parameter can significantly affect the numerical solutions and the behavior of the phase transition.In particular, when this parameter becomes very small (approaching zero), there will be a sharp interface limit (the mean curvature flow).This work aims to construct low-dimensional models for Allen-Cahn equation with different values of this interaction length parameter by using model reduction approach with common projection bases that can be generated only once and reused for all various parameter values.In the mathematical and computational field, model order reduction (MOR) is a technique for constructing a low dimensional model that preserve the necessary information from the full discretized system, while saving computational time and storage.The advantage of a reduced model often lies in the decrease of simulation time and the provision accurate results when compared to the high dimensional system.There are many MOR techniques such as balance truncation [7], Krylov subspace MOR methods [8], projection based MOR [9], and the Moment-matching method [10].This work uses proper orthogonal decomposition (POD) as a starting point for constructing reduced-order model for Allen-Cahn equation.POD approach generally gives high accurate reduced models with much smaller dimensions by generating a set of optimal basis that represents a given data set with minimum error.However, it cannot truly reduce the computational complexity of nonlinear term.We therefore employ the discrete empirical interpolation method (DEIM), as well as a modified Gappy POD (GPOD) for efficiently constructing reduced order systems.
POD is a popular method for model reduction, which was first proposed by Lumley [11].It is also known as Karhunen-Loéve decomposition [12] or principal component analysis [13].It provides a technique for analyzing multidimensional data.This method constructs an orthonormal basis for representing the given data in a certain least squares optimal sense.The basic properties of the POD method have been studied in [14] as it is applied to data compression and used as a model reduction technique for finite dimensional linear systems.Prajna [15] purposed sufficient conditions for preserving stability in POD model reduction.POD has been widely applied in many application such as data analysis, data compression and model reduction in various fields of engineering and science.Applications of POD include image processing [16], data compression, signal analysis [17], turbulence modeling [18], control of fluids [19], electrical power grids [20], and modeling and control of chemical reaction systems [21,22].Noor et al. [23] presented a reduced-basis technique and a computational algorithm within the context for analyzing nonlinear structure.Peterson [24] demonstrated that the reduced basis method based on POD can efficiently approximate the solution of incompressible viscous flow.Ravindran [25] demonstrated a reduced-order modeling approach using POD for simulation and control of incompressible flows.Gunzburger, Peterson, and Shadid [26] focused on reduced order modeling of time-dependent PDEs with multiple parameters on POD approach to the reduced order model.Carlberg and Ferhat [27] used a compact POD to compute a reduced basis for an optimization-oriented reduced order model and applied this to compute a reduced basis for model reduction of static systems [28].
POD also has been used for various nonlinear dynamical systems.In [29], POD has been used with global optimum search framework to construct a temporally-local reduced model for nonlinear parabolic PDEs.An adaptive framework for POD was introduced in [30], when new snapshots are included for feedback control of dissipative nonlinear PDEs.In [31], a modification of data ensemble revision approach was developed and used with adaptive POD framework to construct reduced order models for output feedback control of distributed processes.The notion of sparse POD has been introduced and successfully used for constructing reduced order models of nonlinear parabolic PDEs with moving boundaries [32].
DEIM is a nonlinear technique for approximating nonlinear terms of a dynamical system to further decrease computational complexity of POD reduced systems.The DEIM approach was introduced in [33], and approximates a nonlinear function by combining projection with interpolation.This approach is a discrete variant of the empirical interpolation method (EIM) [34], which was originally described in a function space setting.The DEIM approximation was extended to nonlinear reduced models by localization (LDEIM) [35].There are other nonlinear techniques for approximating nonlinear terms, such as the Gauss-Newton with approximated tensors (GNAT) method [36], the best points interpolation method [37] and missing point estimation [38].There are many works combining a reduced model by POD with DEIM such as the photovoltaic module [39], the shallow water equations model [40], the Navier-Stokes equations [41] and drift-diffusion equations [42].
Gappy POD (GPOD) was first purposed by Everson and Sirovich [43].It has been used as an approach to restore incomplete data.Thanh, Damodaran, and Willcox [44] presented this in the context of fluid dynamic application.Bos et al., [45] presented this in the context of solving the numerical simulation of nonlinear system.Willcox [46] extented GPOD to handle unsteady flow reconstruction problem. Lee and Mavris [47] applied this method for various problems in aerospace engineering.Murray and Ukeiley [48] used GPOD in the context of particle image velocimetry (PIV) data for subsonic cavity flow.

•
A novel model reduction algorithm that improves the accuracy of the existing standard model reduction technique that combines the orthogonal decomposition (POD) and discrete empirical interpolation method (DEIM) by using a modification of GPOD that is based on least-squares approximation to the nonlinear term.

•
Comparison of computational complexity of the proposed method and the standard methods in both theoretical and numerical aspects.

•
Numerical experiments that demonstrate the accuracy and efficiency of the proposed method for solving the parametrized Allen-Cahn equation for both homogeneous and non-homogeneous boundary conditions with different types of initial conditions.
The remainder of this paper is organized as follow.In Section 2, the general form of Allend-Cahn equation and its corresponding discretization are given.Section 3 provides some background on POD and POD-DEIM approaches, as well as introduces the POD-GPOD approach.The accuracy and efficiency of the POD-GPOD framework are compared with the standard POD approach and POD-DEIM approach in Section 4 for simulating the solution of parametrized Allen-Cahn equation.Finally, conclusion and remarks are discussed in Section 5.

Model Problem of the Allen-Cahn Equation and Finite Difference Approximation
This section considers a general form of Allen-Cahn equation and provides the corresponding difference discretization that will be used to find numerical solution.

Model Problem
Consider the Allen-Cahn equation of the form with F(u) = u − u 3 , and Ω = [a, b] , for some real number a, b with a < b where the initial condition : u(x, 0) = f (x), and the boundary conditions : u(a, t) = g 1 (t), u(b, t) = g 2 (t), t > 0. Here, we have variables t for time and x for position, and F(u) is the nonlinear reaction term.The boundary conditions may depend on the problem for a given domain Ω.In this context, u = u(x, t) denotes the phase function of two kinds of liquid at each position x in space, and time t.The parameter is a small positive number, which is related to the width of the interface of two kinds of liquid or the interaction length and can be used to describe the thickness of phase boundary in the laboratory scale [1].The term ∂ 2 u ∂x 2 represents the diffusion of the liquid, and the term F(u) is the kinetic potential of the liquid.For u(x 1 , t 1 ) = 1, it means that at the time t 1 , the position x 1 is filled only one kind of liquid; for u(x 2 , t 2 ) = −1, it means that at the time t 2 , the position x 2 is filled only the other kind of liquid.For t = 0, we will provide the initial data of u(x, 0).This equation has three constant steady states, u = −1, u = 0, and u = 1.The middle state is unstable, but the states u = 1, and u = −1 are attracting, and solutions tend to exhibit at areas close to these values separated by interfaces that may coalesce or vanish on a long time scale, a phenomenon known as metastability.Thus it is called a bistable reaction-diffusion equation or a reaction-diffusion equation with bistable nonlinearity.

Discretization
In this subsection, we present the discretization of the Allen-Cahn Equation (1) using the Crank-Nicolson finite difference method for the linear term and the forward Euler method for the nonlinear term.The resulting discretized system is given by where ∆x and ∆t are the step-sizes of space and time discretizations, respectively.The notations i and j are the indices of the discretization in space and time.Equation (2) can be written as − 2(∆x) 2 u(x i+1 , t j+1 ) + In matrix notation, (3) is given by where A, B ∈ R n×n are a finite difference matrices, and F is nonlinear function that maps a vector in R n to R n , for n spatial grid points on the domain.We can also write (4) in the form In general, when the solution of the full discretized system is required to be highly accurate, solving the system may require high computational complexity, or has large discretized dimension n.

Model Order Reduction
This section presents three methods for efficiently reducing original model: proper orthogonal decomposition (POD), the discrete empirical interpolation method (DEIM), and a modified Gappy POD (GPOD) that can improve the accuracy of DEIM approximation.

POD Reduced System
In this section, we will reduce the dimension of the full discretized system for the Allen-Cahn equation using POD-Galerkin method.To construct a POD reduced system, we have to first compute POD basis from the n s snapshots of the solution u(t) from the full system (4): where u j = u(t j ).We will next explain how POD basis can be obtained based on optimal orthogonal approximation.
Let Φ be a matrix whose columns consist of orthonormal vectors of rank k ≤ r, where r = rank(S).Then, recall that the best approximation of u j by using this matrix Φ is given by POD basis is the optimal orthonormal basis that minimizes these approximation errors as shown in the following definition.

Definition 1 ([49]
).The POD basis matrix of dimension k < r, denoted by Consider the singular value decomposition (SVD) of the snapsthot matrix S in the form S = VΣW T , where V = [v 1 , ..., v r ] ∈ R n×r and W = [w 1 , ..., w r ] ∈ R n s ×r are matrices with orthonormal columns, and Σ =diag(σ 1 , ..., σ r ) ∈ R r×r is a diagonal matrix with It can be shown [49] that the POD basis matrix V k = [v 1 , ..., v k ] consists of the first k columns of the left singular vectors in V corresponding to the k largest singular values and this minimum sum of approximation errors can be written in terms of neglected singular values as This error can be used to identify an appropriate dimension k for a given threshold of accuracy.We can construct the POD basis V k from the following Algorithm 1 [50], derived from solution snapshots of the full system in S.

2.
Set a threshold to select the k largest modes from the diagonal matrix Σ based on (6).

3.
Select the columns in matrix V which correspond to modes selected in 2 to generate the POD basis {v 1 , ..., v k } ∈ R n×k and construct the POD projection matrix We can use POD basis matrix V k for constructing the reduced system for the full discretized system, by the following two steps.
Step 1 We consider the representation u ≈ V k ũ, where Step 2 We apply the Galerkin projection along the subspace generated by the columns of V k , which forms an orthonormal basis (V T k V k = I ∈ R k×k ).The resulting system is given by Finally, we obtain the following POD reduced system where C = V T k A −1 BV k ∈ R k×k can be precomputed and reused in each iteration.For the nonlinear term, although we can precompute V T k A −1 , the computational complexity for multiplying it with F(V k ũ(t j )) still depends on the original dimension n.This issue can be overcome by applying DEIM or GPOD approximation as explained next.

POD-DEIM Reduced System
Since the nonlinear term of the POD reduced system in (9) may still have a large computational cost due to the dependence on n, we will use DEIM to approximate this nonlinear term.DEIM approximation consists of two main steps.The first step is to project the nonlinear term on a basis set of dimension m, which is generally much less than n.This basis is obtained from POD basis of the nonlinear snapshot matrix F = {F(u(t 1 )), ..., F(u(t n s ))} by using Algorithm 1. Suppose the resulting POD basis matrix of the nonlinear snapshots is denoted by Z = [z 1 , . . ., z m ] ∈ R n×m .The second step of DEIM is to perform interpolation of the projected nonlinear term on its certain components, which are chosen by a greedy selection procedure as shown in Algorithm 2. This procedure is called DEIM algorithm in this work.
Solve (P T Z)c = P T z i for c 5. r = z i − Zc 6.

7.
Z ← Z z i , P ← P e i , ← i

end for
Note that the input to Algorithm 2 can be any set of linearly independent vectors {z i } m i=1 .In the standard POD-DEIM approach, this set of input vectors is obtained from POD basis of nonlinear snapshots.The output is the vector that contains the indices for unity components in matrix P with P = [e 1 , ..., e m ] ∈ R n×m , where e i = [0 ... 0 1 0 ... 0] T ∈ R n is the i th column of the identity matrix I ∈ R n×n for i = 1, ..., m.From Algorithm 2, in line 1, the index 1 is determined by the largest absolute value of elements in z 1 .In line 2, we define Z = z 1 ∈ R n×1 and P = e 1 ∈ R n×1 .The selected index is stored in 1 ∈ R. In line 3, we start the iterative process for i = 2, ..., m and solve P T Z c = P T z i to find the constant c in line 4.In line 5, we consider the residual r = z i − Zc and determine the next index by finding the largest component in |r|.In line 6, the index i is determined by the largest absolute value of element in r.In line 7, this index is used to update the matrix P by appending a new vector e i to the existing one.From the iterative process, we obtain Z, P and .Note that the matrix P T Z in each iteration of Algorithm 2 is a nonsingular matrix as proved in [51].
After obtaining Z = [z 1 , ..., z m ] ∈ R n×m and P = [e 1 , ..., e m ] ∈ R n×m , we are ready to apply the DEIM approximation for the nonlinear term in (9): Suppose that f(t j ) = F(V k ũ(t j )) and define the approximation from projecting on the span of Z as where c(t j ) can be determined by selecting m rows from (11), i.e., c(t j ) can be computed by solving the system where multiplying by P T is equivalent to selecting m rows in this system (13).Since P T Z is nonsingular [51], the approximation of ( 11) is From f(t j ) = F(V k ũ(t j )), (13) can be written as Since F has to be evaluated for the original dimension before being interpolated by the matrix P, we use the fact that F is a componentwise function, which implies P T F(V k ũ(t j )) = F(P T V k ũ(t j )) and Therefore (10) can now be approximated by From the POD reduced system (9), we obtain the following POD-DEIM reduced system where M = V T k A −1 Z(P T Z) −1 and M ∈ R k×m can be precomputed and used in each iteration.The dimension n of the nonlinear term in (9) decreases to m in (17) where m n, which can reduce computational cost for obtaining the approximate numerical solution.

POD-GPOD Reduced System
This subsection provides an extension of DEIM based on Gappy POD to improve the accuracy of the POD-DEIM reduced system while still reducing the computational cost of nonlinear term in a POD system.This method uses a dimension of the nonlinear basis less than the number of the indices used for selecting row in the approximation.Recall the approximation in (11): Instead of using m components of the equation above to determine c(t j ), we can use q components where q > m.These q components can be determined from Algorithm 2 with the input basis vectors z 1 , ..., z m , z m+1 , ..., z q obtained from the left singular vectors of the nonlinear snapshot matrix F = {F(u(t 1 )), ..., F(u(t n s ))} corresponding to the q largest singular values.In this case the maximum value of q is the rank of the matrix F. Suppose P ∈ R n×q be the output from Algorithm 2 from using {z j } q j=1 .Then, the goal is to obtain c(t j ) from such that min c(t j ) The optimal solution is given by where Since F is a componentwise function, c(t j ) = ( P T Z) + F( P T V k Z(t j )), and we obtain the following POD-GPOD reduced system ũ(t j+1 ) = C precomputed:k×k where N = V T k A −1 Z( P T Z) + ∈ R k×q , can be precomputed and used in every iteration.Note that, similar to the DEIM approximation, pre-multiplying P T is equivalent to selecting the rows 1 , . . ., q and, therefore, there is no need to perform the matrix multiplication directly.Notice that the computational complexity for each time step depends on only the dimensions k and q, In general, when k, q n, the POD-GPOD reduced system can significantly speed up the simulation time for solving the Allen-Cahn equation.The steps for constructing the POD-GPOD reduced system are summarized in Algorithm 3.
Table 1 provides the comparison of computational complexity for solving the original full-order Allen-Cahn Equation ( 5), the POD reduced system (9), the POD-DEIM reduced system (17), and the POD-GPOD reduced system (20) in each time step.
From Table 1, the computational cost of the FD full system in each time step is dominated by O(n 3 ) where n is the dimension of spatial domain, which is generally large.The cost of computing the POD reduced system is O(k 3 + nk 2 ) where k is the dimension of POD basis and k n.We reduced the complexity for computing the nonlinear term of the POD reduced system by using DEIM, and GPOD.The resulting complexity of the POD-DEIM and POD-GPOD reduced systems, respectively, are given by O(k 3 + mk 2 ) and O(k 3 + qk 2 ), where m is the dimension of nonlinear POD basis and q is the dimension of selected indices used in GPOD.Notice that when m n and q n, the POD-DEIM and POD-GPOD reduced systems can substantially reduce the computational time of the original full-order system.Note that, although POD-GPOD reduced system has higher complexity than the POD-DEIM reduced system when q > m, it can provide more accurate solution as shown next in the numerical experiments.

Algorithm 3:
Steps for constructing the POD-GPOD reduced system • INPUTS: -Reduced dimensions k, m, q -Coefficient matrices: A and B from the discretized Allen-Cahn Equation ( 4).

2.
Construct the POD basis matrix V k of dimension k from S by using Algorithm 1.

3.
Construct the POD basis matrices Z and Z of dimensions m and q, respectively, from F by using Algorithm 1.

4.
Use a set of vectors in Z as an input to Algorithm 2 to obtain q selected indices (or equivalently the matrix P).

5.
Precompute matrices Use C, N, V k , and P to form the POD-GPOD reduced system (20).
Table 1.Computational complexity for computing numerical solutions for Allen-Cahn equation.

System Complexity (One Iteration Step)
The full discretized system O(n 3 ) The POD reduced system O(k 3 + nk 2 ) The POD-DEIM reduced system O(k 3 + mk 2 ) The POD-GPOD reduced system O(k 3 + qk 2 )

Numerical Results
This section presents the numerical tests that demonstrate the accuracy and efficiency of the resulting POD-GPOD approach on approximating the solution of the parametrized Allen-Cahn equation when compared with the standard POD and POD-DEIM approaches.It mainly considers 3 numerical tests.In Section 4.1, we test these model reduction frameworks on Allen-Cahn equation with non-homogenious boundary conditions and square block initial condition.In Sections 4.2 and 4.3, we test these model reduction approaches with various parameter values for two different initial conditions.In particular, we use only one basis for POD and one basis for DEIM or GPOD to construct several reduced systems for different parameter values that are not previously used in the snapshot sets.In all of these numerical experiments, we use the MATLAB program for solving the numerical solutions.The approximation error of the solution from each of these reduced systems is defined as where u and ũ, respectively, are snapshot matrices from solving the FD full system and reduced system, which could be POD, POD-DEIM, or POD-GPOD reduced system.

Numerical Test 1 (Fixed Parameter): Non-Homogenious Boundary Conditions with Square Block Initial Data
In this subsection, we used a square block initial condition to test model reduction approaches for Allen-Cahn Equation (1).The initial condition is defined by a piecewise constant function with value a on wave crest and value b on wave trough.We used the following inputs: internal point (n = 100) in [−1, 1], time steps (n t = 600) on [0, 12] and = 0.01.We used a space step-size ∆x = 2 n−1 , and time step-size ∆t = 12 n t −1 .First, we selected wave crest a = 1 and wave trough b = −1.We solved the numerical solution of the FD full system.When the time increases, the evolution of the Allen-Cahn equation still keeps the phase separation after long time calculation in the left plot of Figure 1.From the right plot of Figure 1, the decreasing of singular values of the snapshot matrix and then beginning to stabilize implies that solution information lies within a subspace whose dimension is significantly lower than the full dimension used in the discretization.That is, we can further reduce the dimensions of the full discretized system and still obtain accurate approximate solution.Next, we solved the numerical solution of the reduced system and used dimensions of POD basis k = 50 and DEIM basis m = 60.We obtained the approximate evolution of phase function of Allen-Cahn equation from POD, and POD-DEIM approaches for time t = 0 to t = 12 as shown in the first two plots of Figure 2. Finally, GPOD is used for reducing the complexity of the nonlinear term from the POD reduced system using a dimension of nonlinear basis m that is less than the number of selected rows q.The dimension of POD basis, nonlinear basis, and the number of selected rows in this numerical test are k = 50, m = 60, and q = 90, respectively.We obtained the evolution of the phase function of Allen-Cahn equation by using POD-GPOD method for time t = 0 to t = 12 as shown in the last plot of Figure 2. Notice that the approximations from both POD-DEIM and POD-GPOD seem to be visually indistinguishable from the original solution.More details on errors are given in Table 2 and Figure 3.
The runtime and error from ( 21) are computed for POD-reduced system as shown in Table 2 and for POD-DEIM and POD-GPOD reduced systems as shown in Table 3 for different dimension m of nonlinear basis, and number of selected rows q.It can be observed in Figure 3 that the error decreases as the dimension of POD increases as also shown in Table 2.For a fixed POD dimension of k = 50, the error of POD-DEIM reduced system can be reduced by increasing the dimension of DEIM basis as shown in Table 3 and Figure 3. Similarly for POD-GPOD reduced system, when POD dimension is fixed to be k = 50, the error can be decreased by increasing both the dimension of nonlinear basis and the number of selected rows as shown in Table 3. From Figure 3 and Table 3, for each m, the error of POD-GPOD approach is clearly smaller than the error of POD-DEIM approach.
Note that, from Figure 3, although the accuracy of the POD approach is higher (less error) than the POD-DEIM and POD-GPOD approaches, the reduction of the simulation time seems to be far less than these last two approaches, as shown in Tables 2 and 3. To compare the speedup of the reduced order approaches, the runtimes given in Tables 2 and 3 are all scaled with the simulation time for the full-order system, which is considered to be 100%.From Table 2, the runtime of POD approach is shown to be approximately 20-30% of the original full order system, i.e., around 3.3-5 times faster than the full order system, when the POD dimensions k = 10, 20, ...., 70 are used.From Table 3, the POD-DEIM approach used only 2.51-4.88% of the full order system's runtime, i.e., roughly 20-40 time faster, when POD dimension k = 50 with DEIM dimension m = 10, ..., 70 are used.For POD-GPOD approach, as shown in the second part of Table 3, the runtime is reduced to approximately 5.49-6.57% of the original system, i.e., roughly 15-19 times faster than the original system.This shows that DEIM and GPOD techniques can significantly speed-up the simulation time of POD approach with small trade-off on accuracy.Notice that POD-DEIM approach provides the fastest speedup, but it has slightly less accuracy than the others.As a result, the POD-GPOD approach is shown to balance between the accuracy and the computational speedup.
In the next two numerical experiments, we will consider Allen-Cahn equation with parameter variation for different initial conditions and boundary conditions.

Numerical Test 2 (Parameter Variation): Homogeneous Boundary Conditions
The goal of this section is to construct reduced systems that have different values of a given parameter by using one basis set for each POD, DEIM, and GPOD.We consider Allen-Cahn Equation (1) with Ω = 2π n , 2π , n ∈ Z + , n > 1, the initial condition : u(x, 0) = 0.25 sin(x), and the homogeneous boundary conditions : u(0, t) = 0, u(2π, t) = 0, t > 0. In our numerical tests, the number of internal point is n = 600 in 2π, 2π n , the number of time steps is n t = 700 on [0, 5] and = 0.01.We used a space step-size is ∆x = 2π n−1 , time step-size is ∆t = 5 n t −1 and the basis sets used in POD and DEIM approximations were constructed from n s = 700 snapshots.
We used snapshots from the full systems with = 0.01 and = 0.99, as shown in Figure 4, respectively, to construct snapshot matrices S 1 , F 1 and S 2 , F 2 .These snapshot matrices were used to construct basis sets for constructing reduced models with different parameter values , i.e., = 0.2, 0.5, 0.8 as shown in Figure 5.The singular values corresponding to the combined snapshot matrix [S 1 , S 2 ], which is used for computing POD basis, is shown in the last plot of Figure 4.The combined nonlinear snapshot matrix [F 1 , F 1 ] is used for constructing the basis for DEIM and GPOD approximations.Notice from the plots in Figures 4 and 5 that the behaviors of the solutions are clearly different as the value of parameter changes.In Tables 4 and 5, the average of errors from ( 21) is given for POD, POD-DEIM, and POD-GPOD approaches with parameter variation = 0.2, 0.5, and 0.8.The corresponding average runtimes given in these tables are scaled with the runtime of the full order system, whose average runtime is considered to be 100%.Note that in Table 5, the dimension k is fixed to be 50 with different dimension m of nonlinear basis and different dimension q of selected rows.It can be observed that the decrease of the average error in Figure 6 corresponds to the increase of dimension of POD basis, as also shown in Table 4.At a fixed POD dimension of k = 50, the decrease of the average error corresponds to the increase of dimension of DEIM basis, as also shown in Table 5.Similarly, at a fixed POD dimension of k = 50, the decrease of the average error corresponds to the increase of both dimension of nonlinear basis and the number of selected rows in POD-GPOD approach.In Figure 6, the average error for the POD-GPOD reduced system is decreased more than the average error from the POD-DEIM reduced system.The average runtime of the reduced systems by POD approach was ranging between 39.37% to 46.91% of the average runtime of the full-order system, which is roughly 2.1-2.5 times speedup as shown in Table 4.The average runtime in Table 4 shows that POD approach can speedup the simulation of the full-order model to roughly 2.1-2.5 times faster (using 39.37% to 46.91% of the average runtime for the full-order system).In Table 5, the POD-DEIM approach and POD-GPOD approach are shown to give around 300-500 times and 240-255 times, respectively, for the speedups.Although GPOD seems to use more simulation time than DEIM, it has already reduced the runtime of the full-order system (more than 200 time speedup) and, more importantly, it can provide approximation with less error than DEIM around one order of magnitude when the dimension m ≤ 50.Table 6 summarizes the accuracy and the scaled average computational time of the reduced systems when compared with the full-order system (the runtime of full-order system is normalized to be one), for a special case fixed dimensions k = 50, m = 40, and q = 70.

Table 6. [Test 2]
The ratio of the average runtime of these systems to the runtime of the full discretized system with parameter values = 0.2, 0.5, and 0.8, using k = 50, m = 40 and q = 70.

System Average Error Ratio of Average Runtime
The full discretized system -1 The POD reduced system O(10 This section considers parametrized Allen-Cahn Equation ( 1) with non-homogeneous boundary conditions.We used Ω = [−1, 1], with the initial condition: u(x, 0) = 0.53x + 0.47 sin(1.5πx),and the non-homogeneous boundary conditions: u(−1, t) = −1, u(1, t) = 1, t > 0. In this numerical tests, the number of internal points is n = 600 in [−1, 1], the number of time steps is n t = 700 on [0, 60] and = 0.01.We used a space step-size is ∆x = 2 n−1 , time step-size is ∆t = 60 n t −1 .The value of parameter is varying in this numerical test.We used snapshots from the full systems with = 0.011, and = 0.009, from time t = 0 to t = 60, as shown in Figure 7, to construct matrices of solution snapshots and nonlinear snapshots.The decay of singular values of the matrix of solution snapshots is shown in the last plot of Figure 7.These snapshot matrices were used to construct basis sets for constructing reduced models of full systems with different parameter values , i.e., = 0.0095, 0.01, 0.0105, as shown in Figure 8.Notice from Figures 7 and 8 that the behaviors of the solutions are clearly different as the parameter changes, i.e., as the value of parameter gets bigger, the solution reaches the steady state faster.In Tables 7 and 8, we present average runtimes and average errors from the parameter variation using = 0.0105, 0.01, and 0.0095, for POD, POD-DEIM, and POD-GPOD approaches, respectively.We computed the runtime and error from (4.1).
Table 8 shows the average error of the POD-GPOD reduced systems with POD basis of dimension k = 50, where m is the dimension of nonlinear basis, and q is the number of selected rows.The results in this section illustrate the same trends as in the previous numerical tests.From Tables 7 and 8, the POD approach is more accurate, but it is much less efficient in term of computational time than the POD-DEIM and POD-GPOD approaches.From Table 7, when the POD approach uses k = 50, the average error is of order O(10 −9 ) and the average runtime is shown to be 32.957% of the full-order system, which is approximately three times faster than the full-order system.From Table 8, when the POD-DEIM reduced system uses k = 50 and m = 50, the average error is of order O(10 −6 ) and the average runtime is 0.191% of the full-order system, which is approximately 524 times faster than the full-order system.From Table 8, when the POD-GPOD reduced system uses k = 50, m = 50 and q = 80, the the average error is of order O(10 −7 ) and the average runtime is 0.379% of the full-order system, which is approximately 264 times faster than the full-order system.Figure 9 shows that the reduced systems become more accurate as the dimension k for POD and dimension m for DEIM increase.As in the previous numerical test, the second plot in Figure 9 also demonstrates that the POD-GPOD reduced system is more accurate than the POD-DEIM reduced system for each fixed dimension m.
Table 9 summarizes the accuracy and the scaled computational time of the reduced systems when compared with the full-order system.In Table 9, the ratio of the average runtime of the full discretized system is normalized to be 1.It considers a special case fixed dimensions k = 50, m = 60, and q = 90.These results demonstrate that GPOD can further improve the accuracy of the DEIM, while still substantially decrease the computational complexity of the nonlinear term of POD reduced system.Table 9. [Test 3] The ratio of the average runtime of these systems to the runtime of the full discretized system, and average error of these systems with parameter values = 0.0105, 0.01, and 0.0095, using k = 50, m = 60, and q = 90.

System Average Error Ratio of Average Runtime
The

Conclusions
In this paper, we consider model reduction techniques for Allen-Cahn equation.We introduced a modification of DEIM that is based on the concept of GPOD, called POD-GPOD approach, which compromises between the accuracy of POD and computational efficiency of DEIM.We compared the POD-GPOD approach with the standard POD method and POD-DEIM method on the parametrized Allen-Cahn equation.From the numerical results, as expected, the POD reduced system is more accurate than the POD-DEIM reduced system.However, the POD-DEIM reduced system can be used to substantially reduced the computational complexity in the nonlinear terms, and therefore the runtime for solving is significantly less than the POD reduced system.We extend the POD-DEIM reduced system to the POD-GPOD reduced system by increasing the selected rows used in the nonlinear approximation to improve the accuracy.The POD-GPOD reduced system can further reduce the approximation error when compare with the POD-DEIM reduced system with a small increase in simulation time that is still much less than the one used in POD approach.Therefore, POD-GPOD approach is an efficient method that can balance between the accuracy and efficiency for parametrized nonlinear model reduction.This approach can also readily be applied to other nonlinear dynamical

Figure 1 .
Figure 1.Solution of the FD full-order system for Allen-Cahn equation with = 0.01 and the corresponding singular values of the snapshot matrix S from SVD using POD basis.

Figure 3 .
Figure 3. [Test 1] Error plot of the approximation the POD reduced system (left plot) and the POD-DEIM, and POD-GPOD reduced systems (right plot) with dimension of DEIM or nonlinear basis m.

Figure 4 .
Figure 4. [Test 2] Solution of the full discretized system with = 0.01.The last plot shows the singular values of the snapshot matrix [S 1 , S 2 ] from SVD using POD basis.

Figure 6 .
Figure 6.[Test 2] Average error plot of the approximation the POD, POD-DEIM, and POD-GPOD reduced systems with dimension of DEIM or nonlinear basis m.

Figure 7 .
Figure 7. [Test 3] Solution of the full discretized systems with = 0.009 and = 0.011, and the corresponding singular values these snapshots (from left to right).

261 Figure 9 .
Figure 9. [Test 3] Average error plot of the approximation from POD reduced systems with parameter values = 0.0105, 0.01, and 0.0095, using different dimensions of POD basis k (left) and from the POD-DEIM and POD-GPOD reduced systems with different dimensions of nonlinear basis m for a fixed POD dimension k = 50 (right).

Table 2 .
[Test 1] Runtime and error of the POD reduced system.

Table 3 .
[Test 1] Runtime and error of the POD-DEIM reduced system (left table) and POD-GPOD reduced system (right table) with POD basis with dimension k = 50.

Table 4 .
[Test 2] Average runtime and average error of the POD reduced systems with parameter values = 0.2, 0.5 and 0.8.

Table 5 .
[Test 2] Average runtime and average error of the POD-DEIM and POD-GPOD reduced systems with various parameter values and POD basis with dimension k = 50.

Table 7 .
[Test 3] Average runtime and average error of the POD reduced systems with parameter values = 0.0105, 0.01, and 0.0095,

Table 8 .
[Test 3] Average runtime and average error of the POD-DEIM reduced systems with with parameter values = 0.0105, 0.01, and 0.0095.The POD basis has dimension k = 50 for both POD-DEIM and POD-GPOD approaches.