Adaptivity in Bayesian Inverse Finite Element Problems: Learning and Simultaneous Control of Discretisation and Sampling Errors

The local size of computational grids used in partial differential equation (PDE)-based probabilistic inverse problems can have a tremendous impact on the numerical results. As a consequence, numerical model identification procedures used in structural or material engineering may yield erroneous, mesh-dependent result. In this work, we attempt to connect the field of adaptive methods for deterministic and forward probabilistic finite-element (FE) simulations and the field of FE-based Bayesian inference. In particular, our target setting is that of exact inference, whereby complex posterior distributions are to be sampled using advanced Markov Chain Monte Carlo (MCMC) algorithms. Our proposal is for the mesh refinement to be performed in a goal-oriented manner. We assume that we are interested in a finite subset of quantities of interest (QoI) such as a combination of latent uncertain parameters and/or quantities to be drawn from the posterior predictive distribution. Next, we evaluate the quality of an approximate inversion with respect to these quantities. This is done by running two chains in parallel: (i) the approximate chain and (ii) an enhanced chain whereby the approximate likelihood function is corrected using an efficient deterministic error estimate of the error introduced by the spatial discretisation of the PDE of interest. One particularly interesting feature of the proposed approach is that no user-defined tolerance is required for the quality of the QoIs, as opposed to the deterministic error estimation setting. This is because our trust in the model, and therefore a good measure for our requirement in terms of accuracy, is fully encoded in the prior. We merely need to ensure that the finite element approximation does not impact the posterior distributions of QoIs by a prohibitively large amount. We will also propose a technique to control the error introduced by the MCMC sampler, and demonstrate the validity of the combined mesh and algorithmic quality control strategy.


Introduction
The Bayesian statistical framework has been used extensively in the problem of system identification [1] or model updating based on experimental test data [2,3].The objective of using inverse problems to learn or calibrate model or latent parameters, including model error terms, lies in the fact that underlying parametrised computational models are uncertain, or erroneous [4].The experimental data when assimilated into the model, using the Bayesian inference framework, is expected to provide joint estimates of the model parameters conditional on the data and compensate for uncertainty/bias in the model predictions.However, operating only in the field of model parameters is problematic for a number of reasons.Firstly, the underlying computational models for any practical application is expensive.As a result working with a very high resolution numerical model at every stage is prohibitively expensive.Secondly, when using surrogate models in the parameter space in conjunction with a simulator, adaptive enrichment of the response surface does not take adaptive mesh refinement within its purview.This is not optimal since ensuring that the model error in the energy norm is bounded in the parameter space requires a uniformly high resolution and increases the associated computational overhead.Lastly the advanced adaptive mesh refinement techniques are rarely used in conjunction with parametric learning using Bayesian inference.There is significant room for improvement in this regard since a simultaneous control of both statistical and discretisation error would lead to substantially improved predictive numerical models both in terms of accuracy and computational efficiency.
At the root or our methodology is the estimation of errors due to the finite element approximation of the partial differential equation (PDE) of interest.Classically, spatial resolution of finite element models can be adaptively refined (also known as local h-refinement) based on a posteriori error estimation techniques [5][6][7], combined with re-meshing strategies.Within this, methods focussing on errors estimated in terms of specific quantities of interest (QoI), rather than the classical energy norm, constitute the goal-oriented adaptivity scheme [8][9][10][11][12], which is of particular interest in the present study.Numerical studies have shown that these give better convergence in the local features of the solution compared to traditional approaches.
Integrating model reduction techniques with finite element model updating techniques has received some attention in recent years [13], where the motivation is to use Bayesian model updating framework with an adaptive scheme of enriching the surrogate response surface.Multistage Bayesian inverse problems are quite important in this respect [14,15] and have important applications for system identification of vibrating systems.The prediction error is an important parameter to be calibrated in such cases, as the authors point out.But the improvement in model predictions, if solely focussed on obtaining the posterior probabilistic parameter estimates or for adaptive enrichment of the response surface in the parameter space, without considering simultaneous enhancement of the resolution of the numerical simulator would be unsatisfactory both from the standpoints of computational accuracy and efficiency.
The forward problem of uncertainty propagation has been investigated extensively for the solution of stochastically parametrised partial differential equations.These range from efficient stochastic Galerkin methods using polynomial chaos basis functions [16][17][18][19], stochastic collocation techniques [20,21], Monte-Carlo sampling based methods (and its various improvements) [22][23][24] and other deterministic sampling methods [25][26][27].The main challenge is to obtain a good approximation of the lower order statistical moments of the state vector or specific quantities of interest.On the other hand, advances in the resolution of stochastic inverse problems has become a very active topic in engineering and mathematical research (see e.g., [28]).Scalable Bayesian inversion algorithms for large-scale problems have been investigated [29] and as well as Bayesian inversion to probabilistic robust optimization under uncertainty [30].The use of adaptive sparse-grid surrogates unified with Bayesian inversion for posterior density estimates of the model or design has also been studied parameters [3,31,32].This research focuses mostly on the definition of surrogates of the numerical model and having adaptive methods to control the statistical sampling error in the definition of the response surface.Lately there has been some research [33] that defines adaptivity as the local enrichment of the surrogate model and uses error estimation to bound this error, with no mention of spatial discretisation.However, coupling engineering uncertainty quantification (UQ) with an adaptive scheme for goal-oriented finite element model refinement remains a sparsely studied domain and presents significant challenges.It is so both from the problem formulation perspective, owing to the choice of appropriate candidate estimates based on which adaptive model enrichment can be performed, as well as incorporating it into the general formulation of Bayesian inversion.
The main focus of the paper is to develop a robust methodology for the simultaneous control of errors from multiple sources-the goal-oriented finite element error and the uncertainty-driven statistical error-in a Bayesian identification framework for identification system parameters conditional on data (experimental or otherwise).The novel algorithmic approach proposed in this paper consists in running two Markov Chain Monte Carlo algorithms simultaneously to sample the posterior densities of the quantities of interest (component-wise MCMC).The first chain utilises the current finite element model, whilst the second chain runs with a corrected likelihood function that takes into account the discretisation error.The latter quantity may be constructed by making use of a posteriori finite element error estimates available in the literature.At any time during the sampling process, two empirical densities are available and may be compared to evaluate the effect of the discretisation error onto posterior densities.The second building block of our algorithmic methodology allows us to determine whether enough samples have been drawn by the component-wise MCMC.Using multiple parallel chains combined with bootstrap-based estimates of sampling errors, we automatically stop the MCMC algorithms when either (i) the required level of accuracy is achieved or (ii) we have generated enough statistical confidence in the fact that the discretisation error is too large for our purpose, and therefore that mesh refinement is necessary.While any available error estimate of the discretisation error may be used within the general strategy outlined above (a goal-oriented residual or recovery-based error estimate, for instance [8,9,11,34]), we choose instead to construct this estimate via a dedicated machine learning approach.This interesting feature, inspired by previous work in data-driven error modelling [35][36][37][38], will be briefly outlined in the paper.
The paper is organised as follows.Section 2 introduces the Bayesian inverse problem based on a finite element model of a parametrised vibrating structure.This section also includes discussions on discretisation error.Some numerical examples are presented in Section 3, which aims to demonstrate the convergence of the joint posterior distributions on the model parameters through successive stages of mesh refinement.Section 4 discusses the total error as a combination of statistical and discretisation error that results from the MCMC algorithm used for sampling from posterior distributions.Section 5 gives the methodology for robust, simultaneous control of all error sources within the adaptive inverse problem solver using a component-wise MCMC algorithm in conjunction with bootstrap-aggregated regression model for model parameters.Numerical examples are presented and discussed in Section 6 to demonstrate the capabilities of the proposed methodology.

Finite Element Bayesian Inverse Problems
Although the methodology proposed in this paper is general, we will apply it to the verification of popular inverse finite element procedures used to monitor the integrity of structures during service.We will assume that the structure can be modelled as an elastic body, and that potential structural damage can be modelled by the evolution of a field of elastic constants characterised by a finite number of uncertain parameters.Assuming that the resonance frequencies of the structure have been measured physically, we attempt to identify this set of mathematical parameters, through the solution of an inverse problem.Significant deviations from the parameter values corresponding to an undamaged structural state may reveal structural failure.This is a typical condition monitoring method used for non-destructive structural testing.
In order to clarify the proposed study, we will assume that we have measured the n d first dynamic of the structural dynamics model, here the log of the elastic modulii of the subdomains represented in Figure 1.
We will also aim to predict the remaining n p − n d frequencies p = ln(ω n d +1 ) ... ln(ω n d +n p ) T .
The relationships µ → d(µ) and µ → p(µ), i.e., the computational model, are defined implicitly through the evaluation of a standard finite element model of the steady-state, undamped structural vibrations.

Applied & Computational
Mechanics Group

Bayesian Inverse Problem
We assume that the quantities d that are measured experimentally are described by the mathematical model, up to an error, which is modelled in a probabilistic manner.Following standard Bayesian procedures, this error is modelled as a white Gaussian noise, with Σ 0 a positive definite covariance matrix.
is the mathematical prediction that corresponds to the experimental observation.Model A d is a deterministic function of n µ uncertain model parameters that we organise in a n µ -dimensional vector µ.Bayesian inversion requires to associate a prior probability density with the model parameters µ.This prior probability density, denoted by π µ (µ) in the following, encodes the knowledge that we possess about µ before making any physical observations.Seen from a different point of view, the probabilistic inversion setting amounts to the definition of a joint probability distribution for µ and d: The formal expression of likelihood function L(µ; It is now possible to formally proceed to the inversion itself by conditioning the joint distribution to actually observed quantities.Applying Bayes' formula, the posterior probability density of the model parameters is where π d (d) is a normalising constant, whose computation requires a usually intractable integration.Bayes formula provides us with an updated knowledge about the uncertain part of our mathematical model.It is now possible to predict unobserved quantities i.e., by formally propagating the posterior uncertainty through model A p : µ → p m := A p (µ).The numerical model is implicitly defined through the solution of a continuum mechanics problem.The components of d m and p are the logarithm of the eigenvalues corresponding to the following parametrised eigenvalue problem: In the previous variational statement, Ω is the domain occupied by the structure of interest, H 1 (Ω) is the space of functions defined over Ω, with values in R 2 , that are zero on part ∂ u Ω of the boundary of the domain, and whose derivatives up to order one are square integrable.∇ s denotes the symmetric part of the gradient operator.D is the fourth order Hooke tensor and ρ is the mass density of the solid material.The problem possesses an infinite number of solutions (ω, u) called free vibration modes.We order the free vibration frequencies in increasing order ω 1 < ω 2 < ... < ω n d +n p .The free vibration modes are functions of uncertain material parameters µ through the definition of the Hooke tensor.Specifically, with the parametrised Lamé constants Solving the continuum mechanics model is equivalent to evaluating the mapping There is, in general, no analytical solution to the continuous vibration problem, and a standard way to obtain approximate solutions is to substitute a finite element space U h (Ω) for infinite dimensional search space H 1 (Ω) ⊃ U h (Ω) [39].Too coarse a finite element discretisation may result in poorly predictive results, while too fine a mesh will lead to numerically intractable results, or, in any case, to a waste of computing resources.

Finite Element Approximation of the Bayesian Inverse Problem
Whilst using the continuum mechanics model exactly would deliver the posterior density π µ|d (µ; d) as solution of the Bayesian inverse problem, we now have an approximate posterior density where the finite element likelihood L is obtained by substituting finite element mapping Ād for A d in Equation (3).Similarly, the approximate posterior density of p is denoted by πp|d (µ; d) and obtained by evaluating finite element mapping Āp instead of A p when propagating the posterior uncertainties forward.

Discretisation Error
The finite element error is the mismatch between πµ|d (µ; d) and π µ|d (µ; d) on the one hand, and πp|d (µ; d) and π p|d (µ; d) on the other hand.Various measures can be used to quantify this mismatch, amongst which the Hellinger distance, defined by the Kullback-Leibler divergence, the total variation distance and the Kolmogorov-Smirnov (KS) distance, defined by where the capital symbols Π .denotes the cumulative probability density corresponding to π ., and Π. is the finite element approximation of Π . .This contribution will make use of the latter measure, in a one-dimensional setting (i.e., it will be applied to control the accuracy of the posterior density of one of the elements of µ or one of the elements of p).The attractiveness of the Kolmogorov-Smirnov distance is its straightforward application in the context of Monte-Carlo procedures, where only empirical densities are available, and its closeness to confidence intervals (CIs), which makes its values relatively easy to interpret within the context of a posteriori error estimation.Notice that the D ks is lower bounded by 0 (identical density functions) and upper bounded by 1 (non-overlapping support for the probability density functions).

Numerical Examples-Part I: Effect of Discretisation Errors Onto Posterior Densities
This section introduces the numerical examples that will be investigated in this paper, and aims to provide a first qualitative understanding of the effect of mesh refinement onto the quality of posterior probability densities.

Forward Stochastic Model
The stochastic field of Young's modulus that will be used to exemplify the error control approach proposed in this paper is defined via a decomposition of domain The domain decomposition is represented in Figure 1.Then, the proposed model is such that Hence, the scalar parameters contained in vector µ ∈ R n µ are the logarithms of the Young's modulii corresponding to each of the subdomains.
The prior density is a multivariate Gaussian, and is given by equation The prior mean µ 0 is the null vector, and the prior variance Σ 0 = σ 2 0 I d , where I d ∈ R n µ × R n µ is the identity matrix, is diagonal and isotropic.The prior density is represented in dimensions (µ 1 , µ 3 ) in Figure 1.
Finally, we model the error as a zero-mean multivariate Gaussian (consistently with what was described in previous section), with independent components and isotropic variance, i.e., Σ n = σ 2 n I d .

Computational Meshes
The evaluation of the likelihood function appearing in solution (11) of the Bayesian inverse problem requires solving the continuum mechanics problem using the finite element method.In this example, we use a sequence of meshes {M i } i∈ n m associated with a monotonically increasing number of degrees of freedom.These meshes are represented in Figure 2.Although the sequence of meshes is not strictly hierarchical, we see that the typically uniform element size is divided by ≈ 2 when moving from mesh M i to mesh M i+1 .The intermediate meshes {M i+ 1 2 } i∈ 1 ñm represented in Figure 2 will be used later on.

Applied & Computational Mechanics Group
Figure 2. Hierarchy of computational meshes used in the study.We aim at selecting the coarsest mesh that delivers the targeted distributions up to a user-defined numerical accuracy.The second line of meshes is only used for error estimation purposes.

Inverse Problems and First Results
Two tests will now be investigated: • Test 1(weakly informative data): only the first eigenvalue is measured, i.e., d is scalar.This can be interpreted as a task of model updating, where new data is used to update an existing knowledge.
• Test 2 (strongly informative data): the first three eigenvalues of the structure are measured.This can be interpreted as an inverse problem, where rich information is used to identify all the unknown of the model, and the probabilistic setting acts as a regulariser.
The two corresponding marginal posterior density πµ 1 ,µ 3 |d obtained when using mesh M 3 are represented in Figure 3.Samples from this density are obtained by a Monte-Carlo sampler, as presented in next section.A Kernel Density Estimate (KDE) is used as a smoother for illustration purposes only.Notice that for the predictive posterior densities, the histograms correspond to the marginal densities of each individual eigenvalue, which explains their overlap.

Applied & Computational
Mechanics Group

Data
Posterior It is interesting to notice that the posterior densities observed in Test 2 are much sharper than that of Test 1, owing to the quality of the data.The symmetry in the results are a consequence of structural and probabilistic symmetries (see Figure 1 and the definition of the prior probability density).
Synthetic data for this problem is generated by computing the average of the spectra delivered by meshes M 2 and M 2+ 1 2 , for a reference parameter vector µ d = µ 0 but situated in the vicinity of the prior mode.The model averaging, together with selecting relatively coarse finite element meshes in sequence {M i } i∈ n m , is meant to circumvent the "inverse crime" problem.

Convergence with Mesh Refinement
The posterior densities corresponding to Test 2 and to various level of mesh refinement are displayed in Figure 4.The difference between the modes of the predicted eigenvalues that are obtained with coarse and fine meshes are very large.The fact that this discrepancy increases with the mode number is to be expected as the spatial wave length of the deformations in the continuum can be shown to increase linearly with the eigenfrequency.As a consequence, a good mesh for the first range of frequencies might be unable to capture the faster spatial variations associated with higher vibration mode.

Applied & Computational
Mechanics Group The difference between the posteriori densities corresponding to meshes M 3 and M 4 is qualitatively small.The solution of the Bayesian inverse problem converges with mesh refinement.Notice that the posterior distribution for the model parameters goes from mono-modal to multi-modal, which may prove a stumbling block when selecting an appropriate Monte-Carlo sampler.

Tempered Metropolis-Hastings Markov Chain Monte Carlo Algorithm
The approximate posterior density described by Equation ( 11) is an arbitrarily complex function of µ, as it depends on the nonlinear computational model (A d (µ) T A p (µ) T ) T .In particular, standard random number generators cannot be used to draw samples from this distribution.Importance Monte-Carlo procedures whereby the prior distribution is used as proposal density, will also fail.This is because the posterior density may be arbitrarily different from the prior density, resulting in unacceptably large variances of importance sampling estimates.Designing better proposal densities a priori is impossible and a successful importance sampling approach would require generating the proposal density using advanced methods such as sequential Monte-Carlo samplers.
One of the simplest and generic generator of samples from posterior densities is the Metropolis-Hastings (MH) Markov Chain Monte Carlo (MCMC) algorithm (see Figure 5).
The algorithm works as follows.Starting from sample µ n , the next sample µ n+1 is drawn from transitional distribution Applied & Computational Mechanics Group This is done by first drawing a random move from µ k using the multivariate Gaussian (any proposal can be used, but the expressions exposed in this section are only valid for symmetric proposals), and then accepting or rejecting the move in order to account for the first term of the transition.Typically, a move ending up in a state of higher posterior density is always accepted, whilst a move ending up in a state of lower density may be accepted or not, depending on the result of a die and the ratio of posterior densities between current and proposed states.
This particular transition is designed such that the following ergodic property holds: As a result, under some assumptions, the Markov chain is guaranteed to have πµ|d as its stationary distribution, which means that µ n ∼ πµ|d as n → ∞.Each sample, taken individually, is distributed according to the posterior density, if the chain is run long enough.Due to the Markov process, the samples are not independent, which makes frequentist error estimation and convergence diagnosis more difficult than in the context of traditional Monte-Carlo algorithm, where tools such as standard error and bootstrap apply without particular difficulty.Practically, a burn-in phase is first observed, whereby the chain "seeks" the regions of high probability densities (see Figure 5).Once found, samples become progressively distributed according to the target posterior density.
Notice that Σ is difficult to choose a priori.Adaptive proposal MCMC have been proposed in [40,41].Alternative solutions to choose good proposal are methods based on particle mechanics (e.g., Langevin diffusion, Hamiltonian dynamics) (see e.g., [42]).In this particular contribution, the proposal densities have been calibrated "by hand".We use a tempered version of the MCMC algorithm, whereby we simultaneously sample multiple tempered replicates of the posterior density, with proposed state exchange between replicates that are subsequently Metropolis corrected [43,44] (see also [45] for an interesting application in the context of structural damage assessment).This relatively classical MCMC algorithm allows us to sample from; densities that are multi-modal or become multi-modal with mesh refinement (at least in the low-dimensional parametric settings that is investigated in this paper, as multi-modality in high-dimensions remains an open problem).
The call to the finite element solver is hidden in the acceptance/rejection test, which requires the evaluation of the finite element likelihood at the proposed state.Therefore, as for a standard MC algorithm, one iteration means one evaluation of the computational model.

Empirical Posterior Densities
Once n samples M = {µ 1 , µ 2 , ... , µ n } have been computed (n larger than 10) by the MH algorithm, we discard the first 25% (burn-in) of these samples, resulting in ñ samples that we hope are located in the ergodic part of the Markov chain.The resulting empirical probability density is given by πµ|d where δ is the Dirac delta function.Accordingly, the empirical predictive posterior density is The cumulative empirical predictive posterior density may be expressed as where I is the indicator function, and the inequality is to be understood in a component-by-component manner.A similar definition holds for the cumulative posterior density of p.

Total Error Measure
It is now clear that both the finite element discretisation and the Monte-Carlo approximate sampling affect the quality of the resulting posterior densities.The total error, for the marginal distribution of a single element µ i of µ, reads as which can be formally decomposed as follows In the last expression, the first term is the pure finite element error, which would occur if we could run the Markov process for an infinite number of iterations, while the second term is a pure statistical error.
We also define the error of an element p i of p as follows:

Simulation of the Discretisation Error
The finite element method introduces an error in computational mapping . Mapping A contains all the scalar quantities that need to be evaluated through calls to the finite element solver, namely the numerical predictions of the physical measurements d, and the posterior predictions p.We define the error in the simulated data as and the error in the posteriori predictions as For now, we assume that both these quantities can be estimated, at affordable numerical cost and in a reliable manner.Therefore, for any value of parameter µ, a corrected computational model is available, which reads as where symbol .denotes computable estimates.

Component-Wise MCMC
It is now posible to sample the corrected posterior distribution using MCMC.It should be clear that the corrected posterior distribution is simply obtained by replacing the corrected computational model ( 26) into the expression of the likelihood function, Equation (3).Notice that the normalising constant is affected by modifications of the computational model.This is of no practical consequence as MCMC samplers work with unnormalised densities, and the KS distance uses empirical cumulative distributions directly, without the need for smoothing or marginalisation.
We formally define a component-wise MCMC were the uncorrected and corrected computational models are sampled at the same time.This will yield an estimate of the effect of the discretisation error on posteriori densities at any stage of the Markov process, which, in turn, will allow us to develop an early-stopping methodology.The Component-wise MCMC iteration proceeds as follows, given a current sample ( μn , μn ) of the uncorrected/corrected finite element posteriori densities,
In the field of a posteriori finite element error estimation, the error estimate is usually a post-processing operation of the coarse finite element solution.This adapts seamlessly to non-Markovian Monte-Carlo samplers, by post-processing each of the independent samples.Here, unfortunately, the finite element model has to be called twice at every iteration of the Markov process.Whether this can be avoided or not, for instance by making use of elements of sequential Monte-Carlo samplers, is unclear to us at this stage.

Machine Learning-Based Simulation of the Discretisation Error
At this stage, any numerical error estimator can be used, provided that it is goal-oriented.A method of choice could be a residual-based [5,6] or smoothing-based a posteriori error estimator [12,46] in conjunction with the adjoint methodology [8,9,11].In addition, nothing prevents us from using a meta-modelling approach, such as projection-based reduced order modelling [47][48][49][50][51][52] or polynomial chaos expansions [16,17] to approximate the variations of the computed quantities of interest with parameter variation.Error estimates also exist for such two-level approximations.
In this contribution, we develop and use a feature-based method, that finds its roots in data-science methodologies and is much more "black-box" than the previously mentioned strategies.The proposed technique is inspired by the work of [35][36][37].
The exact continuum mechanics model is well approximated by a very refined numerical strategy (e.g., no meta-modelling and a very fine mesh).However, this very refined numerical model cannot be used at every iteration of MCMC as its evaluation is very costly.We propose to train a model that will map parameter µ to the output of the generally intractable very fine model through the combination of (i) a dedicated feature extractor and (ii) a weakly parametric regression model.This combination is defined as where the .denote quantities that are delivered by the overkill (but computable) numerical model and R is a weakly parametrised regression model, here a neural network regression (a Gaussian process could be used as well, but a random forest would probably have been the most efficient choice, given the way we bootstrap the regression model to generate estimates of generalisation errors), parametrised by a set of parameters θ ∈ R n θ (we do not explicitly distinguish parameters and hyper-parameters in our notations).F is a mapping from input µ to a feature space.Its careful design is critical to the success of the machine learning procedure.We choose to construct the following features: where Ãd (µ) are slighlty corrected computational models, here generated by refining the mesh by a moderate factor.In our examples, the typical mesh size is divided by 1.5 (see Figures 2 and 6 where we use mesh M i+ 1 2 to correct mesh M i , whilst the overkill solution is computed using M i+2 ).In this fashion, the generation of features will remain of the order of the computation of the finite element solution itself.
where e m j denotes the jth canonical vector of R m .
Each of the regression R l,i is a single-hidden layer bootstrap-aggregated neural network model with n n neurons and n nbs bootstrap replicates.

Training
We sample artificial "data" by running the overkill computational model n ml times, after sampling the training set parameters {µ ml 1 , µ ml 2 , ... , µ ml n ml } from prior π µ .The number of neurons and the cardinality of the training set is chosen automatically by making use of an automatised early-stopping methodology that aims to maximise the predictive coefficient of determination.We will not detail this procedure here.Fitting of the nonlinear regression coefficients is performed by employing standard least-squares method, solved by a gradient descent algorithm with randomised initialisation.Outliers of the set of bootstrap replicates are identified and eliminated to decrease the variance of the boostrap-aggregated regression model.
An example of fitted regression model is represented in Figure 6, where the output is the discretisation error in the first free vibration circular frequency (i.e., regression model R d,1 ).

Bootstrap Confidence Intervals for the MCMC Sampler
At any iteration n of the MCMC sampler, a Monte-Carlo estimate of the discretisation error for the posterior density of the ith component of µ is given by where Π(n) and Π(n) Crucially, D (n) ks is a random variable whose statistics, and in particular its bias and variance, strongly depend on the length of the Markov chain.Unfortunately, evaluating the convergence of any statistics provided by MCMC is difficult, due to the statistical dependency between successively drawn samples.
Following standard diagnostic approaches for MCMC samplers (e.g., the Gelman-Rubin convergence test [53,54]), we will run n c ≥ 10 independent (tempered) MCMC chains in parallel and pool all the resulting samples, after discarding the first 25% of every individual chain as burn-in (see Figure 7 as a visual aid).The pooled samples at iteration n of the multiple-chain MCMC (MC 3 ) algorithm are where ∏ denotes the cartesian product.In the previous expression, the sample set from chain i (at ambient temperature) is D (n)  ks is now the pooled KS distance estimate provided by the MC 3 algorithm (we will keep the same notation for the sake of simplicity).Formally, we simply replace Equation (37) by and perform a similar operation to define the pooled corrected empirical distribution.
The independence of the n c MCMC chains allows us to compute confidence intervals for D (n)  ks by making use of the non-parametric bootstrap.This is done by resampling S(n) and Ŝ(n) with replacement, generating bootstrap replicates of the pooled sample sets { S(n where B k ∈ ( 1 n c ) n bs is such that each element of this set is drawn uniformly over 1 n c and k varies between 1 and a large number n bs , typically set to 1000.For each replicate, statistics D (n) ks , denoted by D (n)  ks,k , can be computed in a straightforward manner by using the bootstrap replicates of the pooled empirical distributions which reads as Finally, the bootstrap confidence intervals are calculated by calculating the Xth and (100 − X)th bootstrap percentiles such that the Xth percentile reads as where Q X is an operator that extracts the Xth and Yth percentile of the set passed as argument.
It is important to understand that the derived bootstrap confidence interval stands for a chain of finite length n, and not for the asymptotic limit.In fact, estimate D (n) ks of D ks is strongly biased (upward for small asymptotic errors D ks and typically downward for large asymptotic errors), which is due to two factors: • the existence of the burn-in phase.For small n, each individual chain will be strongly affected by the initialisation of the chains.For small asymptotic errors, this can be expected to have a strong upward bias effect providing that the initialisation is disperse, which is often the case in practice (e.g., initialisation from the prior, sequential MC approaches with decreasing levels of noise).This can be visualised in Figure 7. Two incomplete chains running on the same probability density may be exploring completely different regions of space, yielding values of D ks that are large, even in the case where the corrected and uncorrected densities are close to one another.

•
the discrete evaluation of the KS distance itself, which generates an additional (upward) bias.
The variance of D (n) ks decreases with the number of chains of the MC 3 algorithm (which is not a free parameter as overall CPU cost increases linearly with n c ).Of course, bias and variance are both expected to decrease with the length n of the run.The results are pooled, and will be boostrapped to derive confidence intervals for the summary statistics of the posterior densities.

Simultaneous Control of All Sources of Errors
We now make use of the CI derived for D (n) ks to construct an adaptive inverse problem solver that jointly controls the quality of the mesh, and that of the statistical evaluation of the posterior densities.The algorithm is as follows.Given a current mesh M i and number of iterations n of the Monte-Carlo solver, do: Convergence of the finite element-based Bayesian inverse problem is achieved when D ks ≤ γ, where γ is a numerical tolerance that will typically be chosen between 0.01 and 0.2.As we only have Monte-Carlo estimates of D ks , we require instead that Criterion C m is an indicator of the combined effect of the discretisation and statistical errors onto the posterior densities.The discretisation error is evaluated through the choice of measure D ks as reliability indicator, whilst the statistical error is taken into account by making use of the upper limit of the bootstrap confidence interval for D (n) ks .The risk of falsely detecting mesh convergence due to a high statistical error is small, due to the fact that D (n) ks is an upwardly biased estimate of D ks for small D ks .Confidence in the result can be increased through a loose Gelman-Rubin convergence test in order to eliminate the risk of stopping the MC 3 algorithm in its non-ergodic phase.However, this has proved to be unnecessary in our numerical tests.This is because C m is a rather strict criterion.

Criterion C s
The second criterion will help us determine whether mesh refinement is actually needed, or whether the statistical error is too large for use to take a robust decision regarding mesh refinement.The ideal criterion is D ks ≥ η, where η < γ.The obvious strategy that would require q (n) X ≤ η does not work.This is due to the previously explained upward bias of estimator D (n) ks , which would eventually lead to the spurious satisfaction C s , and consequently to systematic mesh refinement operations for low n count even when the mesh becomes fine enough for D ks to be well below target γ.In order to derive an appropriate criterion C s , we remark that whilst the ideal criterion D ks ≥ η cannot be statistically evaluated for general values of η, criterion D ks > 0 can be.More precisely, we propose to estimate the statistical error by evaluation whether, at current n, the corrected and uncorrected posterior distributions can be statistically distinguished.
We postulate the following null hypothesis: Hypothesis 0 (H0): the corrected and uncorrected posterior densities are identical.
The rejection of the null hypothesis will indicate that the two densities are significantly different.Now, the criterion for the need for mesh refinement becomes the following: The probability density of D (n) ks under null hypothesis H0 can be approximated by adapting the previously described bootstrap procedure.We will estimate the density of where π(n),1 µ i |d and π(n),2 µ i |d are obtained, respectively, by pooling the result of two independent runs of n c Markov chains of lengths n and corresponding to the uncorrected finite element model only (one can equally choose the corrected one).Notice, to clarify the idea, that D (n)  ks,0 tends to 0 as n tends to infinity: this is a measure of the statistical error only.The desired density can be estimated by resampling S(n) twice, computing the KS distance between the two pooled sets of samples, and repeating the operation n bs times, which generates a sequence of real numbers {D (n)  ks,0,k }.The replicated empirical cumulative distributions are Π(n),j with the sampling sets S(n),j where j = 1 or j = 2 and B 1 k and B 2 k are sets of 1 n c n c constructed such that each element of these two sets is drawn uniformly over 1 n c .We can now extract the (100 corresponding to p-value ξ, and evaluate whether q If so, test C s is true.If not, the statistical sampling error is too large to allow us to decide whether mesh refinement should be performed or not.In this case, and assuming that mesh convergence criterion C m is not satisfied, we need to continue sampling with MCMC. Notice that the use of a larger percentile Z yields a less conservative indicator for the need for mesh refinement, which can be compensated by decreasing p-value ξ.The criterion that is arguably the easiest to interpret is Pr D (n)  ks ≥ q (n) 50 |H0 ≤ ξ, where q (n) 50 is directly the KS distance D (n) ks (i.e., computed without bootstrapping).

Numerical Examples-Part II. Automatised Error Control and Discussion
We now come back to the numerical examples introduced in Section 3 and produce three series of results, illustrating three difference aspects of the proposed approach.

MCMC Iterations Only When Needed
The results reported in Figure 8 correspond to Test 1.For now, we aim to control the quality of the posterior density of the first parameter µ 1 .Each of the graphs shows the evolution of various statistics of the corresponding KS distance as a function of the number of iteration n of the tempered MC 3 algorithm.The mesh used is displayed on the left-hand side of each of the graphs.The solid black line represents the evolution D (n) ks itself.The 5th and and 95th bootstrap percentiles of this quantity are also reported, in solid blue line and solid red line, respectively.We set convergence criterion C m such that γ = 0.2 and Y = 95 for the pure discretisation error.Therefore, convergence with respect to the finite element discretisation is obtained when, for sufficient long chains, the red line (upper limit of the bootstrap CI for D Iteration when needed The data is composed of a unique noisy measurement of the first eigenfrequency.The quantity of interest is here the first elasticity parameter, in log scale.The grey lines are an indicator of the bias of the KS distance estimate, which decreases with the sample size.The bias should be be small for the estimate to be robust, but not too small so as to avoid running unnecessarily long chains using meshes that are too coarse for our needs.Sufficiently accurate results are obtained with the third mesh.Importantly, very few MCMC iterations are needed to determine that the two first meshes are too coarse for our needs.
In grey colours, we report evolution of D (n) ks,0 for the coarse scale (light grey), consistently with the derivations of the previous section, but also, for reference, a similar statistics constructed with the corrected finite element model (dark grey).Both curves have similar evolutions and, of course, tend to 0 as n tends to infinity.We set Z = 5, and η = 0.05.This means that the statistical error is considered small enough once the blue line (lower limit of the bootstrap CI for D (n) ks ) is above the 95th percentile of D (n) ks,0 .Here, "small enough" is to be understood, as small enough for a decision regarding the need to refine the mesh one step further to be taken.
The results are as follows.For the coarsest mesh, we detect a separation of the discretisation and statistical errors at iteration 60 (vertical solid red line).C s is satisfied.As C m is not satisfied, we can stop the MCMC sampler and refine the mesh.A similar behaviour is seen for mesh M 2 .For the last mesh, convergence is achieved after 530 iterations, C s never reaching satisfaction.The discretisation and sampling errors remain entangled, but the sum of them is below the desired target.
We can see here that the proposed algorithm allows us to stop the MCMC sampler early, moving directly to the level of mesh refinement that will yield the desired quality of the posterior densities.

Goal-Oriented Error Control
The second set of results correspond to Test 1, still.However, we now monitor the convergence of the posterior predictive density of the fourth eigenvalue.This is reported in the top two graphs of Figure 9.For M 1 Early-stopping is performed at iteration 40 of the MCMC sampler.For M 2 , global convergence is achieved after 120 iterations, and this error keeps decreasing.The MCMC is allowed to continue generating samples, meaning that the discretisation error is actually a lot smaller than the requested tolerance γ.It is interesting to see that the convergence is much faster than for the first parameter, whose convergence was studied in the previous subsection.This shows that for the same inverse problem, the algorithm may spend more or less resources depending on the engineering quantity of interest.

Applied & Computational
Mechanics Group can be solved appropriately with a coarse mesh, while the data-rich case (c,d) requires a much finer grid to be solved accurately.

Uncertainty-Driven Error Control
Finally, the last set of results concerns Test 2. Here, remember that the three first eigenvalues are used as measurements.We control the convergence of the fourth eigenvalue, which was controlled in the context of Test 1 in the previous subsection.We see here that for mesh M 2 , the posterior density of the fourth eigenvalue is still far from convergence, whilst it was evaluated in a very precise way in Test 1.This shows that the proposed error control algorithm automatically adapts to the level of posterior uncertainty.Qualitatively, a wide posterior density is associated with a high uncertainty concerning the value of QoIs.Consequently, the mesh does not need to be very refined to capture the posterior density correctly.Conversely, for Test 2, the posterior density is sharper, and similar levels of discretisation errors have a much stronger impact on the KS distance.

Concluding Remarks and Discussion
We have presented a methodology to control the various sources of errors arising in finite-element based Bayesian inverse problems.We have focussed on a simple numerical approximation chain consisting of (i) a finite element discretisation of the continuum mechanics problem and (ii) a MCMC solver to draw samples from posterior density distributions.So far we did not consider further error sources such as that engendered by meta-modelling.We have showed that it was possible to drive the mesh refinement process in a goal-oriented manner, by quantifying the impact of the associated error onto posterior density distributions.In order to do so, we run two independent MCMC simultaneously, one of them using a corrected likelihood function that takes the discretisation error.Any a posteriori error estimate available for the PDE under consideration may be used to obtain this correction.Of course, the accuracy (effectivity) of the chosen error estimate will impact the accuracy of the methodology developed in this paper.The study and control of this effect will require further research to be carried out, We have shown that by using multiple replicates of the component-wise MCMC, it was also possible to (i) derive bootstrap-based confidence intervals for the error in posterior densities engendered by the spatial discretisation of the underlying PDE.Importantly, this approach allows us to stop the MCMC iteration as soon as we can be sufficiently confident in the fact that the posterior densities obtained with the current mesh are not accurate enough for our purpose.The approach is goal-oriented: the adaptivity is performed in the sense of the posterior distribution of either some of the latent parameters, or in the sense of the predictive posterior density of engineering QoIs.Finally, we have shown that the approach is uncertainty driven: it will only refine the mesh and/or request additional samples to be drawn by the MCMC algorithm if the effect on the QoIs can be felt when measured using a statistical distance between their posterior densities.In particular, model updating problems with wider priors tend to require less computational effort than parameter estimation problems with very rich observed data and/or narrow prior densities.The proposed methodology establishes a bridge between the field of reliability estimation for the sampling-based algorithms used to solve probabilistic inverse problems, and the field of deterministic, goal-oriented finite element error estimation.
Although the results presented in this paper are encouraging, the proposed approach has shortcomings that need to be addressed in future research work.The error estimation procedure is relatively wasteful for two reasons.Firstly, multiple MCMC algorithms need to run independently for bootstrapping to be possible.These chains all exhibit their own burn-in phase, which increases the amount of discarded samples.Secondly, the finite element error estimation procedure is not merely a post-processing operation any longer; one cannot post-process the finite element results corresponding to the uncorrected chain in order to compute the corrected likelihood.This is because both chains are independent and require running the coarse finite element simulations at different points of the parameter domain.Finally, let us acknowledge that Bayesian finite element inverse problems should not be solved by a MCMC solver without constructing a surrogate model first.The number of finite element computations involved, even if the proposal distribution is well designed, is in the tens of thousands.We are currently investigating the use of Polynomial Chaos surrogates, which adds another layer of numerical approximation that needs to be controlled in a robust and efficient manner.In this context, an elegant approach to separate the sources of errors (i.e., Finite Element discretisation error and Polynomial Chaos error) may be found in [55], which may constitute a solid starting point for the next step of our investigations.
eigenfrequencies d = ln(ω 1 ) ... ln(ω n d ) T of the structure represented in Figure 1.The purpose of the model inversion is to identify the values of n µ parameters µ = µ 1 ... µ n µ T

Figure 1 .
Figure 1.(a) Structure undergoing structural vibrations.The grey area have a known, deterministic Young modulus.This property is piecewise constant in the coloured areas, and Gaussian distributed.(b) The corresponding marginal probability density function for two of the 6 parameters is represented in the top-left corner.(c) Standard Monte-Carlo finite element process is employed to compute the distribution of the 10 first free vibration frequencies and represent the marginal distribution of each of these quantities.

Figure 3 .
Figure 3. (a) Joint posterior probability distribution for two of the 6 unknown elastic modulii and (c) marginal posterior predictive distributions of the first 10 free vibration frequencies.The data set consists of a noisy measurement of the first frequency only, resulting in fat posterior densities.(b) Posterior distribution of the two unknown elastic modulii when the first 3 vibration frequencies are used as dataset, which results as in much thinner posterior densities.(d) Corresponding posterior predictive densities.

3 Figure 4 .
Figure 4. Evolution of the joint posterior density of two of the unknown elasticity parameters as the computational mesh is progressively refined from M 1 to M 4 .

3 Figure 5 .
Figure 5. Three independent MCMC chains sampling a multivariate Gaussian distribution.
t e x i t s h a 1 _ b a s e 6 4 = " l e r O z 5 c y h k D k t Y 2 t S + d r N n L e v H U = " > A A A C I H i c b V D L T g I x F O 3 g C / G F u n T T S E x c k R l j g k t 8 L F x q I m D C w K R T L t D Q e a S 9 Y y S T + R Q 3 / o o b F x q j O / 0 a y 8 B C x Z M 0 P T 3 n 3 v T e 4 8 d S a L T t T 6 u w s L i 0 v F J c L a 2 t b 2 x u l b d 3 m j p K F I c G j 2 S k b n 2 e n b J L 1 h f 3 1 D K p P A = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " l e r O z 5 c y h k D k t Y 2 t S + d r N n L e v H U = " > A A A C I H i c b V D L T g I x F O 3 g C / G F u n T T S E x c k R l j g k t 8 L F x q I m D C w K R T L t D Q e a S 9 Y y S T + R Q 3 / o o b F x q j O / 0 a y 8 B C x Z M 0 P T 3 n 3 v T e 4 8 d S a L T t T 6 u w s L i 0 v F J c L a 2 t b 2 x u l b d 3 m j p K F I c G j 2 S k b n 2 e n b J L 1 h f 3 1 D K p P A = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " l e r O z 5 c y h k D k t Y 2 t S + d r N n L e v H U = " > A A A C I H i c b V D L T g I x F O 3 g C / G F u n T T S E x c k R l j g k t 8 L F x q I m D C w K R T L t D Q e a S 9 Y y S T + R Q 3 / o o b F x q j O / 0 a y 8 B C x Z M 0 P T 3 n 3 v T e 4 8 d S a L T t T 6 u w s L i 0 v F J c L a 2 t b 2 x u l b d 3 m j p K F I c G j 2 S k b n 2 e n b J L 1 h f 3 1 D K p P A = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " l e r O z 5 c y h k D k t Y 2 t S + d r N n L e v H U = " > A A A C I H i c b V D L T g I x F O 3 g C / G F u n T T S E x c k R l j g k t 8 L F x q I m D C w K R T L t D Q e a S 9 Y y S T + R Q 3 / o o b F x q j O / 0 a y 8 B C x Z M 0 P T 3 n 3 v T e 4 8 d S a L T t T 6 u w s L i 0 v F J c L a 2 t b 2 x u l b d 3 m j p K F I c G j 2 S k b n 2 6 D m 5 O B f n k L S o s k b m I / B S 9 i 1 7 H o C s 7 Q U H 7 1 w A 0 S G e Z Q + M 5 V k 7 p 3 I o Q e w 3 z I 6 n 5 k v t w 9 B Y m s G B m P i s J 3 E e 4 x T w u / W r P r 9 r D o O H B K U C N l N f z q m x s m P I s g R i 6 Z 1 h 3 H T t H L m U L B J R Q V N 9 O Q M n 7 D r q F j Y M w i 0 F 4 + P L K g W 4 Y J a T d R 5 s V I h + z P j p x F e r C x c U Y M e / q v N i D / 0 z o Z d v e 9 X M R p h h D z 0 a B u J i k m d J A Y D Y U C j r J v A O N K m F 0 p 7 z H F O J p c K y Y E 5 + / J 4 + B i p + 7 Y d e d 8 t 3 Z 4 X M Y x R z b I J t k m D t k j h + S M N E i L c P J A n s g L e b U e r W f r 3 f o Y W S e s s m e d / C r r 6 x v y j 6 b e < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " j c 3 X B Q t H X 0 s O c N / 4 C c 4 x O b S G J 6 D m 5 O B f n k L S o s k b m I / B S 9 i 1 7 H o C s 7 Q U H 7 1 w A 0 S G e Z Q + M 5 V k 7 p 3 I o Q e w 3 z I 6 n 5 k v t w 9 B Y m s G B m P i s J 3 E e 4 x T w u / W r P r 9 r D o O H B K U C N l N f z q m x s m P I s g R i 6 Z 1 h 3 H T t H L m U L B J R Q V N 9 O Q M n 7 D r q F j Y M w i 0 F 4 + P L K g W 4 Y J a T d R 5 s V I h + z P j p x F e r C x c U Y M e / q v N i D / 0 z o Z d v e 9 X M R p h h D z 0 a B u J i k m d J A Y D Y U C j r J v A O N K m F 0 p 7 z H F O J p c K y Y E 5 + / J 4 + B i p + 7 Y d e d 8 t 3 Z 4 X M Y x R z b I J t k m D t k j h + S M N E i L c P J A n s g L e b U e r W f r 3 f o Y W S e s s m e d / C r r 6 x v y j 6 b e < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " j c 3 X B Q t H X 0 s O c N / 4 C c 4 x O b S G J

1 g
6 D m 5 O B f n k L S o s k b m I / B S 9 i 1 7 H o C s 7 Q U H 7 1 w A 0 S G e Z Q + M 5 V k 7 p 3 I o Q e w 3 z I 6 n 5 k v t w 9 B Y m s G B m P i s J 3 E e 4 x T w u / W r P r 9 r D o O H B K U C N l N f z q m x s m P I s g R i 6 Z 1 h 3 H T t H L m U L B J R Q V N 9 O Q M n 7 D r q F j Y M w i 0 F 4 + P L K g W 4 Y J a T d R 5 s V I h + z P j p x F e r C x c U Y M e / q v N i D / 0 z o Z d v e 9 X M R p h h D z 0 a B u J i k m d J A Y D Y U C j r J v A O N K m F 0 p 7 z H F O J p c K y Y E 5 + / J 4 + B i p + 7 Y d e d 8 t 3 Z 4 X M Y x R z b I J t k m D t k j h + S M N E i L c P J A n s g L e b U e r W f r 3 f o Y W S e s s m e d / C r r 6 x v y j 6 b e < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " j c 3 X B Q t H X 0 s O c N / 4 C c 4 x O b S G J d c = " > A A A C J H i c b V D L S s N A F J 3 4 t r 6 q L t 0 M F s F V S U R Q c O N r 4 b J C q 4 U m h s n k 1 g 5 O H s z c q C X k Y 9 z 4 K 2 5 c + M C F G 7 / F a Z u F W i 8 M c z j n X O 6 9 J 0 i l 0 G j b n 9 b E 5 N T 0 z O z c f G V h c W l 5 p b q 6 d q G T T H F o 8 U Q m q h 0 w D V L E 0 E K B E t q p A h Y F E i 6 D m 5 O B f n k L S o s k b m I / B S 9 i 1 7 H o C s 7 Q U H 7 1 w A 0 S G e Z Q + M 5 V k 7 p 3 I o Q e w 3 z I 6 n 5 k v t w 9 B Y m s G B m P i s J 3 E e 4 x T w u / W r P r 9 r D o O H B K U C N l N f z q m x s m P I s g R i 6 Z 1 h 3 H T t H L m U L B J R Q V N 9 O Q M n 7 D r q F j Y M w i 0 F 4 + P L K g W 4 Y J a T d R 5 s V I h + z P j p x F e r C x c U Y M e / q v N i D / 0 z o Z d v e 9 X M R p h h D z 0 a B u J i k m d J A Y D Y U C j r J v A O N K m F 0 p 7 z H F O J p c K y Y E 5 + / J 4 + B i p + 7 Y d e d 8 t 3 Z 4 X M Y x R z b I J t k m D t k j h + S M N E i L c P J A n s g L e b U e r W f r 3 f o Y W S e s s m e d / C r r 6 x v y j 6 b e < / l a t e x i t > e T A p < l a t e x i t s h a 1 _ b a s e 6 4 = " s r r M n 5 y e t B Z L 7 c P Q G J r B h b D 4 v C d x E e M E 8 L v 9 6 w m 3 Z Z 9 C 9 w J q B B J t X 2 6 0 M 3 T H g W Q Y x c M q 2 7 j p 2 i l z O F g k s o a m 6 m I W W 8 z 2 6 g a 2 D M I t B e X p 5 Z 0 C 3 D h L S X K P N i p C X 7 v S N n k R 5 t b J w R w 1 v 9 W x u R / 2 n d D H v 7 X i 7 i N E O I + X h Q L 5 M U E z r K j I Z C A U c 5 M I B x J c y u l N 8 y x T i a Z G s m B O f 3 y X / B x U 7 T s Z v O 2 W 6 j d T S J o 0 o 2 y C b Z J g 7 Z I y 1 y S t q k Q z h 5 J M 9 k S F 6 t J + v F e r P e x 9 a K N e l Z J z / K + v w C s K e n x w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " s r r M n 5 y e t B Z L 7 c P Q G J r B h b D 4 v C d x E e M E 8 L v 9 6 w m 3 Z Z 9 C 9 w J q B B J t X 2 6 0 M 3 T H g W Q Y x c M q 2 7 j p 2 i l z O F g k s o a m 6 m I W W 8 z 2 6 g a 2 D M I t B e X p 5 Z 0 C 3 D h L S X K P N i p C X 7 v S N n k R 5 t b J w R w 1 v 9 W x u R / 2 n d D H v 7 X i 7 i N E O I + X h Q L 5 M U E z r K j I Z C A U c 5 M I B x J c y u l N 8 y x T i a Z G s m B O f 3 y X / B x U 7 T s Z v O 2 W 6 j d T S J o 0 o 2 y C b Z J g 7 Z I y 1 y S t q k Q z h 5 J M 9 k S F 6 t J + v F e r P e x 9 a K N e l Z J z / K + v w C s K e n x w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " s r r M n 5 y e t B Z L 7 c P Q G J r B h b D 4 v C d x E e M E 8 L v 9 6 w m 3 Z Z 9 C 9 w J q B B J t X 2 6 0 M 3 T H g W Q Y x c M q 2 7 j p 2 i l z O F g k s o a m 6 m I W W 8 z 2 6 g a 2 D M I t B e X p 5 Z 0 C 3 D h L S X K P N i p C X 7 v S N n k R 5 t b J w R w 1 v 9 W x u R / 2 n d D H v 7 X i 7 i N E O I + X h Q L 5 M U E z r K j I Z C A U c 5 M I B x J c y u l N 8 y x T i a Z G s m B O f 3 y X / B x U 7 T s Z v O 2 W 6 j d T S J o 0 o 2 y C b Z J g 7 Z I y 1 y S t q k Q z h 5 J M 9 k S F 6 t J + v F e r P e x 9 a K N e l Z J z / K + v w C s K e n x w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " s r r M n

1 g
5 y e t B Z L 7 c P Q G J r B h b D 4 v C d x E e M E 8 L v 9 6 w m 3 Z Z 9 C 9 w J q B B J t X 2 6 0 M 3 T H g W Q Y x c M q 2 7 j p 2 i l z O F g k s o a m 6 m I W W 8 z 2 6 g a 2 D M I t B e X p 5 Z 0 C 3 D h L S X K P N i p C X 7 v S N n k R 5 t b J w R w 1 v 9 W x u R / 2 n d D H v 7 X i 7 i N E O I + X h Q L 5 M U E z r K j I Z C A U c 5 M I B x J c y u l N 8 y x T i a Z G s m B O f 3 y X / B x U 7 T s Z v O 2 W 6 j d T S J o 0 o 2 y C b Z J g 7 Z I y 1 y S t q k Q z h 5 J M 9 k S F 6 t J + v F e r P e x 9 a K N e l Z J z / K + v w C s K e n x w = = < / l a t e x i t > e T A p < l a t e x i t s h a 1 _ b a s e 6 4 = " s r r M n 5 M I B x J c y u l N 8 y x T i a Z G s m B O f 3 y X / B x U 7 T s Z v O 2 W 6 j d T S J o 0 o 2 y C b Z J g 7 Z I y 1 y S t q k Q z h 5 J M 9 k S F 6 t J + v F e r P e x 9 a K N e l Z J z / K + v w C s K e n x w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " s r r M n 5 M I B x J c y u l N 8 y x T i a Z G s m B O f 3 y X / B x U 7 T s Z v O 2 W 6 j d T S J o 0 o 2 y C b Z J g 7 Z I y 1 y S t q k Q z h 5 J M 9 k S F 6 t J + v F e r P e x 9 a K N e l Z J z / K + v w C s K e n x w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " s r r M n 5 y e t B Z L 7 c P Q G J r B h b D 4 v C d x E e M E 8 L v 9 6 w m 3Z Z 9 C 9 w J q B B J t X 2 6 0 M 3 T H g W Q Y x c M q 2 7 j p 2 i l z O F g k s o a m 6 m I W W 8 z 2 6 g a 2 D M I t B e X p 5 Z 0 C 3 D h L S X K P N i p C X 7 v S N n k R 5 t b J w R w 1 v 9 W x u R / 2 n d D H v 7 X i 7 i N E O I + X h Q L 5 M U E z r K j I Z C A U c 5 M I B x J c y u l N 8 y x T i a Z G s m B O f 3 y X / B x U 7 T s Z v O 2 W 6 j d T S J o 0 o 2 y Cb Z J g 7 Z I y 1 y S t q k Q z h 5 J M 9 k S F 6 t J + v F e r P e x 9 a K N e l Z J z / K + v w C s K e n x w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " s r r M n 5 y e t B Z L 7 c P Q G J r B h b D 4 v C d x E e M E 8 L v 9 6 w m 3Z Z 9 C 9 w J q B B J t X 2 6 0 M 3 T H g W Q Y x c M q 2 7 j p 2 i l z O F g k s o a m 6 m I W W 8 z 2 6 g a 2 D M I t B e X p 5 Z 0 C 3 D h L S X K P N i p C X 7 v S N n k R 5 t b J w R w 1 v 9 W x u R / 2 n d D H v 7 X i 7 i N E O I + X h Q L 5 M U E z r K j I Z C A U c 5 M I B x J c y u l N 8 y x T i a Z G s m B O f 3 y X / B x U 7 T s Z v O 2 W 6 j d T S J o 0 o 2 y C b Z J g 7 Z I y 1 y S t q k Q z h 5 J M9 k S F 6 t J + v F e r P e x 9 a K N e l Z J z / K + v w C s K e n x w = = < / l a t e x i t > e T 1 Āp < l a t e x i t s h a 1 _ b a s e 6 4 = " c y O d D j j W S K E w Z 1 / 4 7 T N Q l s P X D i c c y / 3 3 h O m j C r t O N / W 0 v L K 6 t p 6 Z a O 6 u b W 9 s 2 v v 7 X e U y C Q m b S y Y k L 0 Q K c I o J 2 1 N N S O 9 V B K U h I x 0 w 9 H N x O 8 + E K m o 4 C 0 9 T o m f o A G n M c V I G y m w j 7 x Q s C g n R e D e t 6 A X I p n P l K u i C N L A r j l 1 Z w q 4 S N y S 1 E C J Z m B / e Z H A W U K 4 x g w p 1 X e d V P s 5 k p p i R o q q l y m S I j x C A 9 I 3 l K O E K D + f v l H A E 6 N E M B b S F N d w q v 6 e y F G i 1 D g J T W e C 9 F D N e x P x P 6 + f 6 f j S z y l P M 0 0 4n i 2 K M w a 1 g J N M Y E Q l w Z q N D U F Y U n M r x E M k E d Y m u a o J w Z 1 / e Z F 0 z u q u U 3 f v z m u N 6 z K O C j g E x + A U u O A C N M At a I I 2 w O A R P I N X 8 G Y 9 W S / W u / U x a 1 2 y y p k D 8 A f W 5 w 9 u Q J k X < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " c y O d D j j W S K E w Z 1 / 4 7 T N Q l s P X D i c c y / 3 3 h O m j C r t O N / W 0 v L K 6 t p 6 Z a O 6 u b W 9 s 2 v v 7 X e U y C Q m b S y Y k L 0 Q K c I o J 2 1 N N S O 9 V B K U h I x 0 w 9 H N x O 8 + E K m o 4 C 0 9 T o m f o A G n M c V I G y m w j 7 x Q s C g n R e D e t 6 A X I p n P l K u i C N L A r j l 1 Z w q 4 S N y S 1 E C J Z m B / e Z H A W U K 4 x g w p 1 X e d V P s 5 k p p i R o q q l y m S I j x C A 9 I 3 l K O E K D + f v l H A E 6 N E M B b S F N d w q v 6 e y F G i 1 D g J T W e C 9 F D N e x P x P 6 + f 6 f j S z y l P M 0 0 4n i 2 K M w a 1 g J N M Y E Q l w Z q N D U F Y U n M r x E M k E d Y m u a o J w Z 1 / e Z F 0 z u q u U 3 f v z m u N 6 z K O C j g E x + A U u O A C N M At a I I 2 w O A R P I N X 8 G Y 9 W S / W u / U x a 1 2 y y p k D 8 A f W 5 w 9 u Q J k X < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " c y O d D j j W S K E w Z 1 / 4 7 T N Q l s P X D i c c y / 3 3 h O m j C r t O N / W 0 v L K 6 t p 6 Z a O 6 u b W 9 s 2 v v 7 X e U y C Q m b S y Y k L 0 Q K c I o J 2 1 N N S O 9 V B K U h I x 0 w 9 H N x O 8 + E K m o 4 C 0 9 T o m f o A G n M c V I G y m w j 7 x Q s C g n R e D e t 6 A X I p n P l K u i C N L A r j l 1 Z w q 4 S N y S 1 E C J Z m B / e Z H A W U K 4 x g w p 1 X e d V P s 5 k p p i R o q q l y m S I j x C A 9 I 3 l K O E K D + f v l H A E 6 N E M B b S F N d w q v 6 e y F G i 1 D g J T W e C 9 F D N e x P x P 6 + f 6 f j S z y l P M 0 0 4n i 2 K M w a 1 g J N M Y E Q l w Z q N D U F Y U n M r x E M k E d Y m u a o J w Z 1 / e Z F 0 z u q u U 3 f v z m u N 6 z K O C j g E x + A U u O A C N M At a I I 2 w O A R P I N X 8 G Y 9 W S / W u / U x a 1 2 y y p k D 8 A f W 5 w 9 u Q J k X < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " c y O d D j j W S K E

1 2<
w Z 1 / 4 7 T N Q l s P X D i c c y / 3 3 h O m j C r t O N / W 0 v L K 6 t p 6 Z a O 6 u b W 9 s 2 v v 7 X e U y C Q m b S y Y k L 0 Q K c I o J 2 1 N N S O 9 V B K U h I x 0 w 9 H N x O 8 + E K m o 4 C 0 9 T o m f o A G n M c V I G y m w j 7 x Q s C g n R e D e t 6 A X I p n P l K u i C N L A r j l 1 Z w q 4 S N y S 1 E C J Z m B / e Z H A W U K 4 x g w p 1 X e d V P s 5 k p p i R o q q l y m S I j x C A 9 I 3 l K O E K D + f v l H A E 6 N E M B b S F N d w q v 6 e y F G i 1 D g J T W e C 9 F D N e x P x P 6 + f 6 f j S z y l P M 0 0 4n i 2 K M w a 1 g J N M Y E Q l w Z q N D U F Y U n M r x E M k E d Y m u a o J w Z 1 / e Z F 0 z u q u U 3 f v z m u N 6 z K O C j g E x + A U u O A C N M A t a I I 2 w O A R P I N X 8 G Y 9 W S / W u / U x a 1 2 y y p k D 8 A f W 5 w 9 u Q J k X < / l a t e x i t > M 3+ l a t e x i t s h a 1 _ b a s e 6 4 = " i I 1 B Z w + j 0 A X u X O x 6 8 d F Z i q U W R j 0 = " > A A A C B H i c b V D L S s N A F J 3 U V 6 2 v q M t u B o s g C C W p g i 6 L b t w I F e w D m l A m 0 0 k 7 d D I J M x O h D F m 4 8 V f c u F D E r R / h z r 9 x 0 m a h r Q c u H M 6 5 l 3 v v C R J G p X K c b 6 u 0 s r q 2 v l H e r G x t 7 + z u 2 f s H H R m n A p M 2 j l k s e g G S h F F O 2 o o q R n q J I C g K G O k G k + v c 7 z 4 Q I W n M 7 9 U 0 I X 6 E R p y G F C N l p I F d 9 S K k x h g x f Z s N 9 N m p F w q E t Z v p R p Y N 7 J p T d 2 a A y 8 Q t S A 0 U a A 3 s L 2 8 Y 4 z Q i X G G G p O y 7 T q J 8 j Y S i m J G s 4 q W S J A h P 0 I j 0 D e U o I t L X s y c y e G y U I Q x j Y Y o r O F N / T 2 g U S T m N A t O Z n y w X v V z 8 z + u n K r z 0 N e V J q g j H 8 0 V h y q C K Y Z 4 I H F J B s G J T Q x A W 1 N w K 8 R i Z G J T J r W J C c B d f X i a d Rt 1 1 6 u 7 d e a 1 5 V c R R B l V w B E 6 A C y 5 A E 9 y A F m g D D B 7 B M 3 g F b 9 a T 9 W K 9 W x / z 1 p J V z B y C P 7 A + f w D X O Z g 0 < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " i I 1 B Z w+ j 0 A X u X O x 6 8 d F Z i q U W R j 0 = " > A A A C B H i c b V D L S s N A F J 3 U V 6 2 v q M t u B o s g C C W p g i 6 L b t w I F e w D m l A m 0 0 k 7 d D I J M x O h D F m 4 8 V f c u F D E r R / h z r 9 x 0 m a h r Q c u H M 6 5 l 3 v v C R J G p X K c b 6 u 0 s r q 2 v l H e r G x t 7 + z u 2 f s H H R m n A p M 2 j l k s e g G S h F F O 2 o o q R n q J I C g K G O k G k + v c 7 z 4 Q I W n M 7 9 U 0 I X 6 E R p y G F C N l p I F d 9 S K k x h g x f Z s N 9 N m p F w q E t Z v p R p Y N 7 J p T d 2 a A y 8 Q t S A 0 U a A 3 s L 2 8 Y 4 z Q i X G G G p O y 7 T q J 8 j Y S i m J G s 4 q W S J A h P 0 I j 0 D e U o I t L X s y c y e G y U I Q x j Y Y o r O F N / T 2 g U S T m N A t O Z n y w X v V z 8 z + u n K r z 0 N e V J q g j H 8 0 V h y q C K Y Z 4 I H F J B s G J T Q x A W 1 N w K 8 R i Z G J T J r W J C c B d f X i a d Rt 1 1 6 u 7 d e a 1 5 V c R R B l V w B E 6 A C y 5 A E 9 y A F m g D D B 7 B M 3 g F b 9 a T 9 W K 9 W x / z 1 p J V z B y C P 7 A + f w D X O Z g 0 < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " i I 1 B Z w + j 0 A X u X O x 6 8 d F Z i q U W R j 0 = " > A A A C B H i c b V D L S s N A F J 3 U V 6 2 v q M t u B o s g C C W p g i 6 L b t w I F e w D m l A m 0 0 k 7 d D I J M x O h D F m 4 8 V f c u F D E r R / h z r 9 x 0 m a h r Q c u H M 6 5 l 3 v v C R J G p X K c b 6 u 0 s r q 2 v l H e r G x t 7 + z u 2 f s H H R m n A p M 2 j l k s e g G S h F F O 2 o o q R n q J I C g K G O k G k + v c 7 z 4 Q I W n M 7 9 U 0 I X 6 E R p y G F C N l p I F d 9 S K k x h g x f Z s N 9 N m p F w q E t Z v p R p Y N 7 J p T d 2 a A y 8 Q t S A 0 U a A 3 s L 2 8 Y 4 z Q i X G G G p O y 7 T q J 8 j Y S i m J G s 4 q W S J A h P 0 I j 0 D e U o I t L X s y c y e G y U I t 1 1 6 u 7 d e a 1 5 V c R R B l V w B E 6 A C y 5 A E 9 y A F m g D D B 7 B M 3 g F b 9 a T 9 W K 9 W x / z 1 p J V z B y C P 7 A + f w D X O Z g 0 < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " i I1 B Z w + j 0 A X u X O x 6 8 d F Z i q U W R j 0 = " > A A A C B H i c b V D L S s N A F J 3 U V 6 2 v q M t u B o s g C C W p g i 6 L b t w I F e w D m l A m 0 0 k 7 d D I J M x O h D F m 4 8 V f c u F D E r R / h z r 9 x 0 m a h r Q c u H M 6 5 l 3 v v C R J G p X K c b 6 u 0 s r q 2 v l H e r G x t 7 + z u 2 f s H H R m n A p M 2 j l k s e g G S h F F O 2 o o q R n q J I C g K G O k G k + v c 7 z 4 Q I W n M 7 9 U 0 I X 6 E R p y G F C N l p I F d 9 S K k x h g x f Z s N 9 N m p F w q E t Z v p R p Y N 7 J p T d 2 a A y 8 Q t S A 0 U a A 3 s L 2 8 Y 4 z Q i X G G G pO y 7 T q J 8 j Y S i m J G s 4 q W S J A h P 0 I j 0 D e U o I t L X s y c y e G y U I

3 <
t 1 1 6 u 7 d e a 1 5 V c R R B l V w B E 6 A C y 5 A E 9 y A F m g D D B 7 B M 3 g F b 9 a T 9 W K 9 W x / z 1 p J V z B y C P 7 A + f w D X O Z g 0 < / l a t e x i t > M l a t e x i t s h a 1 _ b a s e 6 4 = " g 9 l g y S u 0 u d w u v s B u S Q 8 r s N L N e D M = " > A A A B + H i c b V D L S s N A F L 3 x W e u j U Z d u B o v g q i Q q 6 L L o x o 1 Q w T 6 g D W E y n b R D J 5 M w M x F q y J e 4 c a G I W z / F n X / j p M 1 C W w 8 M H M 6 5 l 3 v m B A l n S j v O t 7 W y u r a + s V n Z q m 7 v 7 O 7 V 7 P 2 D j o p T S W i b x D y W v Q A r y p m g b c 0 0 p 7 1 E U h w F n H a D y U 3 h d x + p V C w W D 3 q a U C / C I 8 F C R r A 2 k m / X B h H W Y 4 J 5 d p f 7 2 X n u 2 3 W n 4 c y A l o l b k j q U a P n 2 1 2 A Y k z S i Q h O O l e q 7 T q K 9 D E v N C K d 5 d Z A q m m A y w S P a N 1 T g i C o v m w X P 0 Y l R h i i M p X l C o 5 n 6 e y P D k 4 5 W X S O W u 4 T s O 9 v 6 g 3 r 8 s 6 K n A E x 3 A K L l x C E 2 6 h B W 0 g k M I z v M K b 9 W S 9 W O / W x 3 x 0 x S p 3 D u E P r M 8 f 9 7 K T R g = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " g 9 l g y S u 0 u d w u v s B u S Q 8 r s N L N e D M = " > A A A B + H i c b V D L S s N A F L 3 x W e u j U Z d u B o v g q i Q q 6 L L o x o 1 Q w T 6 g D W E y n b R D J 5 M w M x F q y J e 4 c a G I W z / F n X / j p M 1 C W w 8 M H M 6 5 l 3 v m B A l n S j v O t 7 W y u r a + s V n Z q m 7 v 7 O 7 V 7 P 2 D j o p T S W i b x D y W v Q A r y p m g b c 0 0 p 7 1 E U h w F n H a D y U 3 h d x + p V C w W D 3 q a U C / C I 8 F C R r A 2 k m / X B h H W Y 4 J 5 d p f 7 2 X n u 2 3 W n 4 c y A l o l b k j q U a P n 2 1 2 A Y k z S i Q h O O l e q 7 T q K 9 D E v N C K d 5 d Z A q m m A y w S P a N 1 T g i C o v m w X P 0 Y l R h i i M p X l C o 5 n 6 e y P D k 4 5 W X S O W u 4 T s O 9 v 6 g 3 r 8 s 6 K n A E x 3 A K L l x C E 2 6 h B W 0 g k M I z v M K b 9 W S 9 W O / W x 3 x 0 x S p 3 D u E P r M 8 f 9 7 K T R g = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " g 9 l g y S u 0 u d w u v s B u S Q 8 r s N L N e D M = " > A A A B + H i c b V D L S s N A F L 3 x W e u j U Z d u B o v g q i Q q 6 L L o x o 1 Q w T 6 g D W E y n b R D J 5 M w M x F q y J e 4 c a G I W z / F n X / j p M 1 C W w 8 M H M 6 5 l 3 v m B A l n S j v O t 7 W y u r a + s V n Z q m 7 v 7 O 7 V 7 P 2 D j o p T S W i b x D y W v Q A r y p m g b c 0 0 p 7 1 E U h w F n H a D y U 3 h d x + p V C w W D 3 q a U C / C I 8 F C R r A 2 k m / X B h H W Y 4 J 5 d p f 7 2 X n u 2 3 W n 4 c y A l o l b k j q U a P n 2 1 2 A Y k z S i Q h O O l e q 7 T q K 9 D E v N C K d 5 d Z A q m m A y w S P a N 1 T g i C o v m w X P 0 Y l R h i i M p X l C o 5 n 6 e y P D k 4 5 W X S O W u 4 T s O 9 v 6 g 3 r 8 s 6 K n A E x 3 A K L l x C E 2 6 h B W 0 g k M I z v M K b 9 W S 9 W O / W x 3 x 0 x S p 3 D u E P r M 8 f 9 7 K T R g = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " g 9 l g y S u 0 u d w u v s B u S Q 8 r s N L N e D M = " > A A A B + H i c b V D L S s N A F L 3 x W e u j U Z d u B o v g q i Q q 6 L L o x o 1 Q w T 6 g D W E y n b R D J 5 M w M x F q y J e 4 c a G I W z / F n X / j p M 1 C W w 8 M H M 6 5 l 3 v m B A l n S j v O t 7 W y u r a + s V n Z q m 7 v 7 O 7 V 7 P 2 D j o p T S W i b x D y W v Q A r y p m g b c 0 0 p 7 1 E U h w F n H a D y U 3 h d x + p V C w W D 3 q a U C / C I 8 F C R r A 2 k m / X B h H W Y 4 J 5 d p f 7 2 X n u 2 3 W n 4 c y A l o l b k j q U a P n 2 1 2 A Y k z S i Q h O O l e q 7 T q K 9 D E v N C K d 5 d Z A q m m A y w S P a N 1 T g i C o v m w X P 0 Y l R h i i M p X l C o 5 n 6 e y P D k

5 <
4 5 W X S O W u 4 T s O 9 v 6 g 3 r 8 s 6 K n A E x 3 A K L l x C E 2 6 h B W 0 g k M I z v M K b 9 W S 9 W O / W x 3 x 0 x S p 3 D u E P r M 8 f 9 7 K T R g = = < / l a t e x i t > M l a t e x i t s h a 1 _ b a s e 6 4 = " e A 5 A H u 9 B z K j A b i i B k S X 0 D G S G w f E = " > A A A B + H i c b V D L S s N A F L 3 x W e u j U Z d u B o v g q i S i 6 L L o x o 1 Q w T 6 g D W E y n b R D J 5 M w M x F q y J e 4 c a G I W z / F n X / j p M 1 C W w 8 M H M 6 5 l 3 v m B A l n S j v O t 7 W y u r a + s V n Z q m 7 v 7 O 7 V 7 P 2 D j o p T S W i b x D y W v Q A r y p m g b c 0 0 p 7 1 E U h w F n H a D y U 3 h d x + p V C w W D 3 q a U C / C I 8 F C R r A 2 k m / X B h H W Y 4 J 5 d p f 7 2 U X u 2 3 W n 4 c y A l o l b 4 5 W X S O W u 4 T s O 9 P 6 8 3 r 8 s 6 K n A E x 3 A K L l x C E 2 6 h B W 0 g k M I z v M K b 9 W S 9 W O / W x 3 x 0 x S p 3 D u E P r M 8 f + r y T S A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " e A 5 A H u 9 B z K j A b i i B k S X 0 D G S G w f E = " > A A A B + H i c b V D L S s N A F L 3 x W e u j U Z d u B o v g q i S i 6 L L o x o 1 Q w T 6 g D W E y n b R D J 5 M w M x F q y J e 4 c a G I W z / F n X / j p M 1 C W w 8 M H M 6 5 l 3 v m B A l n S j v O t 7 W y u r a + s V n Z q m 7 v 7 O 7 V 7 P 2 D j o p T S W i b x D y W v Q A r y p m g b c 0 0 p 7 1 E U h w F n H a D y U 3 h d x + p V C w W D 3 q a U C / C I 8 F C R r A 2 k m / X B h H W Y 4 J 5 d p f 7 2 U X u 2 3 W n 4 c y A l o l b 4 5 W X S O W u 4 T s O 9 P 6 8 3 r 8 s 6 K n A E x 3 A K L l x C E 2 6 h B W 0 g k M I z v M K b 9 W S 9 W O / W x 3 x 0 x S p 3 D u E P r M 8 f + r y T S A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " e A 5 A H u 9 B z K j A b i i B k S X 0 D G S G w f E = " > A A A B + H i c b V D L S s N A F L 3 x W e u j U Z d u B o v g q i S i 6 L L o x o 1 Q w T 6 g D W E y n b R D J 5 M w M x F q y J e 4 c a G I W z / F n X / j p M 1 C W w 8 M H M 6 5 l 3 v m B A l n S j v O t 7 W y u r a + s V n Z q m 7 v 7 O 7 V 7 P 2 D j o p T S W i b x D y W v Q A r y p m g b c 0 0 p 7 1 E U h w F n H a D y U 3 h d x + p V C w W D 3 q a U C / C I 8 F C R r A 2 k m / X B h H W Y 4 J 5 d p f 7 2 U X u 2 3 W n 4 c y A l o l b 4 5 W X S O W u 4 T s O 9 P 6 8 3 r 8 s 6 K n A E x 3 A K L l x C E 2 6 h B W 0 g k M I z v M K b 9 W S 9 W O / W x 3 x 0 x S p 3 D u E P r M 8 f + r y T S A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " e A 5 A H u 9 B z K j A b

1 d
4 5 W X S O W u 4 T s O 9 P 6 8 3 r 8 s 6 K n A E x 3 A K L l x C E 2 6 h B W 0 g k M I z v M K b 9 W S 9 W O / W x 3 x 0 x S p 3 D u E P r M 8 f + r y T S A = = < / l a t e x i t > e T A p < l a t e x i t s h a 1 _ b a s e 6 4 = " j c 3 s k b m I / B S 9 i 1 7 H o C s 7 Q U H 7 1 w A 0 S G e Z Q + M 5 V k 7 p 3 I o Q e w 3 z I 6 n 5 k v t w 9 B Y m s G B m P i s J 3 E e 4 x T w u / W r P r 9 r D o O HB K U C N l N f z q m x s m P I s g R i 6 Z 1 h 3 H T t H L m U L B J R Q V N 9 O Q M n 7 D r q F j Y M w i 0 F 4 + P L K g W 4 Y J a T d R5 s V I h + z P j p x F e r C x c U Y M e / q v N i D / 0 z o Z d v e 9 X M R p h h D z 0 a B u J i k m d J A Y D Y U C j r J v A O N K m F 0 p 7 z H F O J p c K y Y E 5 + / J 4 + B i p + 7 Y d e d 8 t 3 Z 4 X M Y x R z b I J t k m D t k j h + S M N E i L c P J A n s g L e b U e r W f r 3 f o Y W S e s s m e d / C r r 6 x v y j 6 b e < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " j c 3 s k b m I / B S 9 i 1 7 H o C s 7 Q U H 7 1 w A 0 S G e Z Q + M 5 V k 7 p 3 I o Q e w 3 z I 6 n 5 k v t w 9 B Y m s G B m P i s J 3 E e 4 x T w u / W r P r 9 r D o O HB K U C N l N f z q m x s m P I s g R i 6 Z 1 h 3 H T t H L m U L B J R Q V N 9 O Q M n 7 D r q F j Y M w i 0 F 4 + P L K g W 4 Y J a T d R5 s V I h + z P j p x F e r C x c U Y M e / q v N i D / 0 z o Z d v e 9 X M R p h h D z 0 a B u J i k m d J A Y D Y U C j r J v A O N K m F 0 p 7 z H F O J p c K y Y E 5 + / J 4 + B i p + 7 Y d e d 8 t 3 Z 4 X M Y x R z b I J t k m D t k j h + S M N E i L c P J A n s g L e b U e r W f r 3 f o Y W S e s s m e d / C r r 6 x v y j 6 b e < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " j c 3 s k b m I / B S 9 i 1 7 H o C s 7 Q U H 7 1 w A 0 S G e Z Q + M 5 V k 7 p 3 I o Q e w 3 z I 6 n 5 k v t w 9 B Y m s G B m P i s J 3 E e 4 x T w u / W r P r 9 r D o O HB K U C N l N f z q m x s m P I s g R i 6 Z 1 h 3 H T t H L m U L B J R Q V N 9 O Q M n 7 D r q F j Y M w i 0 F 4 + P L K g W 4 Y J a T d R5 s V I h + z P j p x F e r C x c U Y M e / q v N i D / 0 z o Z d v e 9 X M R p h h D z 0 a B u J i k m d J A Y D Y U C j r J v A O N K m F 0 p 7 z H F O J p c K y Y E 5 + / J 4 + B i p + 7 Y d e d 8 t 3 Z 4 X M Y x R z b I J t k m D t k j h + S M N E i L c P J A n s g L e b U e r W f r 3 f o Y W S e s s m e d / C r r 6 x v y j 6 b e < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " j c 3

1 A
s k b m I / B S 9 i 1 7 H o C s 7 Q U H 7 1 w A 0 S G e Z Q + M 5 V k 7 p 3 I o Q e w 3 z I 6 n 5 k v t w 9 B Y m s G B m P i s J 3 E e 4 x T w u / W r P r 9 r Do O H B K U C N l N f z q m x s m P I s g R i 6 Z 1 h 3 H T t H L m U L B J R Q V N 9 O Q M n 7 D r q F j Y M w i 0 F 4 + P L K g W 4 Y J a T d R 5 s V I h + z P j p x F e r C x c U Y M e / q v N i D / 0 z o Z d v e 9 X M R p h h D z 0 a B u J i k m d J A Y D Y U C j r J v A O N K m F 0 p 7 zH F O J p c K y Y E 5 + / J 4 + B i p + 7 Y d e d 8 t 3 Z 4 X M Y x R z b I J t k m D t k j h + S M N E i L c P J A n s g L e b U e r W f r 3 f o Y W S e s s m e d / C r r 6 x v y j 6 b e < / l a t e x i t > e T ?p < l a t e x i t s h a 1 _ b a s e 6 4 = " l e r O z 5 c y h k D k t Y 2 t S + d r N n L e v H U = " > A A A C I H i c b V D L T g I x F O 3 g C / G F u n T T S E x c k R l j g k t 8 L F x q I m D C w K R T L t D Q e a S 9 Y y S T + R Q 3 / o o b F x q j O / 0 a y 8 B C x Z M 0 P T 3 n 3 v T e 4 8 d S a L T t T 6 u w s L i 0 v F J c L a 2 t b 2 x u l b d 3 m j p K F I c G j 2 S k b n 2 e n b J L 1 h f 3 1 D K p P A = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " l e r O z 5 c y h k D k t Y 2 t S + d r N n L e v H U = " > A A A C I H i c b V D L T g I x F O 3 g C / G F u n T T S E x c k R l j g k t 8 L F x q I m D C w K R T L t D Q e a S 9 Y y S T + R Q 3 / o o b F x q j O / 0 a y 8 B C x Z M 0 P T 3 n 3 v T e 4 8 d S a L T t T 6 u w s L i 0 v F J c L a 2 t b 2 x u l b d 3 m j p K F I c G j 2 S k b n 2 e n b J L 1 h f 3 1 D K p P A = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " l e r O z 5 c y h k D k t Y 2 t S + d r N n L e v H U = " > A A A C I H i c b V D L T g I x F O 3 g C / G F u n T T S E x c k R l j g k t 8 L F x q I m D C w K R T L t D Q e a S 9 Y y S T + R Q 3 / o o b F x q j O / 0 a y 8 B C x Z M 0 P T 3 n 3 v T e 4 8 d S a L T t T 6 u w s L i 0 v F J c L a 2 t b 2 x u l b d 3 m j p K F I c G j 2 S k b n 2

Figure 6 .
Figure6.Machine-learning-based error estimation procedure.The overkill error (difference between the current mesh and much finer mesh as shown in (a)) is postulated as being an unknown function of an error estimate.The inexpensive estimate is the distance between the results obtained when using the current mesh and those obtained when using a slightly refined mesh as shown in.(b) The statistical learning is done offline via a Monte-Carlo sampling of the true error (using the prior as sampling density) and the adaptive fitting of a Neural Network regression.The regression is bootstrapped to provide enhanced stability, and derive confidence interval estimates for the overall error estimation procedure which is shown in (c).Both the size of the dataset and the hyperparameters of the network are found automatically via a Greedy process that is not detailed in this paper.Now, we train n d + n p multivariate neural network regression models

Figure 7 .
Figure 7. Multiple MCMC routines were run in parallel, for both the coarse likelihood function (coarse mesh) as shown in (a) and the enhanced likelihood function obtained after correction by the goal-oriented finite element error estimate, as shown in (b).The results are pooled, and will be boostrapped to derive confidence intervals for the summary statistics of the posterior densities.
perform m 0 iterations of the MC3 algorithm and set n ← n + m 0 .• Evaluate mesh convergence criterion C m .If this criterion is satisfied, exit the adaptation procedure; • If C m is not satisfied, evaluate statistical convergence criterion C s -If C s is satisfied, set i ← i + 1, reinitialise the MC 3 algorithm and set n = 0; -otherwise perform m(n) iterations of the MC 3 algorithm and set n ← n + m(n) (m should be an exponentially increasing function).5.5.1.Criterion C m

Figure 8 .
Figure 8. Convergence of the KS distance as a function of the number of Monte-Carlo samples in each chain (black line, with confidence interval represented by the area between the blue and red lines).The mesh used is represented schematically on the left-hand side of each of the graphs in (a-c).The data is composed of a unique noisy measurement of the first eigenfrequency.The quantity of interest is here the first elasticity parameter, in log scale.The grey lines are an indicator of the bias of the KS distance estimate, which decreases with the sample size.The bias should be be small for the estimate to be robust, but not too small so as to avoid running unnecessarily long chains using meshes that are too coarse for our needs.Sufficiently accurate results are obtained with the third mesh.Importantly, very few MCMC iterations are needed to determine that the two first meshes are too coarse for our needs.

Figure 9 .
Figure 9. Convergence of the KS distance as a function of the number of Monte-Carlo samples in each chain.The (a,b) correspond to the inverse problem where only the first eigenfrequency is observed.(c,d)corresponds to the case where the first three eigenfrequencies are observed.The quantity of interest is the fourth eigenvalue, which is unobserved in both cases and must be inferred.The corresponding posterior density is illustrated above the mesh pictogram.Interestingly, the data-poor problem (a,b) can be solved appropriately with a coarse mesh, while the data-rich case (c,d) requires a much finer grid to be solved accurately.
The posterior predictive probability density of p will be denoted by symbol π p|d (µ; d).