A subspace-adaptive weights cubature method with application to the local hyperreduction of parameterized finite element models

This paper is concerned with quadrature/cubature rules able to deal with multiple subspaces of functions, in such a way that the integration points are common for all the subspaces, yet the weights are tailored to each specific subspace. These subspace-adaptive weights cubature rules can be used to accelerate computational mechanics applications requiring efficiently evaluating spatial integrals whose integrand function dynamically switches between multiple pre-computed subspaces. One of such applications is local hyperreduced-order modeling (HROM), in which the solution manifold is approximately represented as a collection of basis matrices, each basis matrix corresponding to a different region in parameter space. We pose the optimization problem and propose two algorithms for its resolution. The first one is a greedy strategy based on an enhanced version of the Empirical Cubature Method (ECM) developed by the authors elsewhere (we call it the Subspace-Adaptive Weights ECM, SAW-ECM for short), while the second method is based on a convexification of the cubature problem so that it can be addressed by linear programming algorithms. We show in a toy problem involving integration of polynomial functions that the SAW-ECM clearly outperforms the other method both in terms of computational cost and optimality. We illustrate the performance of the SAW-ECM in the construction of a local HROMs in a highly nonlinear equilibrium problem (large strains regime). We demonstrate that, provided that the subspace-transition errors are negligible, the error associated to hyperreduction using adaptive weights can be controlled by the truncation tolerances of the SVDs used for determining the basis matrices. The Python source codes of the proposed SAW-ECM are openly accessible in the public repository https://github.com/Rbravo555/localECM.