Monte Carlo convergence rates for kth moments in Banach spaces
We formulate standard and multilevel Monte Carlo methods for the kth moment š^k_Īµ[Ī¾] of a Banach space valued random variable Ī¾Ī©ā E, interpreted as an element of the k-fold injective tensor product space ā^k_Īµ E. For the standard Monte Carlo estimator of š^k_Īµ[Ī¾], we prove the k-independent convergence rate 1-1/p in the L_q(Ī©;ā^k_Īµ E)-norm, provided that (i) Ī¾ā L_kq(Ī©;E) and (ii) qā[p,ā), where pā[1,2] is the Rademacher type of E. We moreover derive corresponding results for multilevel Monte Carlo methods, including a rigorous error estimate in the L_q(Ī©;ā^k_Īµ E)-norm and the optimization of the computational cost for a given accuracy. Whenever the type of E is p=2, our findings coincide with known results for Hilbert space valued random variables. We illustrate the abstract results by three model problems: second-order elliptic PDEs with random forcing or random coefficient, and stochastic evolution equations. In these cases, the solution processes naturally take values in non-Hilbertian Banach spaces. Further applications, where physical modeling constraints impose a setting in Banach spaces of type p<2, are indicated.
READ FULL TEXT