pith. sign in

arxiv: 2406.19799 · v4 · submitted 2024-06-28 · 🧮 math.NA · cs.NA

Spectral approximation of a new class of stochastic fractional evolution equations

Pith reviewed 2026-05-24 00:00 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords spectral approximationstochastic evolution equationsfractional parabolic equationserror boundsnumerical approximationSPDE methodspatial statistics
0
0 comments X

The pith

A spectral truncation method with quadrature yields strong error bounds for a new class of fractional parabolic stochastic evolution equations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a numerical method that discretizes a recently proposed class of fractional parabolic stochastic evolution equations by truncating their spectral basis expansion in space and applying quadrature to the resulting temporal evolution of each coefficient. Strong error bounds are established for both the spatial spectral approximation and the temporal quadrature step. A sympathetic reader would care because these equations arise as space-time extensions of the SPDE-method used in spatial statistics, so a convergent approximation scheme with explicit error control makes their solutions computable. The method is tested on numerical examples that verify the theoretical rates.

Core claim

The paper introduces and analyses a method for numerical approximation of a new class of fractional parabolic stochastic evolution equations proposed as a space-time extension of the SPDE-method. Space is discretised by truncating the spectral basis function expansion, after which quadrature approximates the temporal evolution of each basis coefficient. Strong error bounds are proved for both the spectral truncation and the temporal quadrature, and the approach is verified by numerical experiments.

What carries the argument

Truncation of the spectral basis function expansion combined with quadrature on the temporal evolution of the coefficients.

If this is right

  • The spatial spectral truncation error decreases strongly as more basis functions are retained.
  • The temporal quadrature error can be bounded independently of the spatial truncation.
  • The combined scheme produces a convergent approximation to solutions of the target class of equations.
  • Numerical experiments confirm that observed errors match the proved strong bounds.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same truncation-plus-quadrature structure could be tested on related fractional SPDEs whose well-posedness is already known.
  • Explicit error bounds may allow construction of adaptive algorithms that choose truncation level and quadrature points to meet a target accuracy.
  • The method supplies a practical route to simulate the space-time models that extend the SPDE-method, enabling direct comparison with existing spatial-statistics techniques.

Load-bearing premise

The new class of fractional parabolic stochastic evolution equations is well-posed and admits a spectral expansion whose coefficients evolve in a manner amenable to quadrature.

What would settle it

A specific instance of the equations where the strong error between the true solution and the spectral-plus-quadrature approximation fails to decrease at the predicted rate as the truncation level or quadrature nodes increase.

Figures

Figures reproduced from arXiv: 2406.19799 by S. Knutsen Furset.

Figure 1
Figure 1. Figure 1: Convergence plots for the left-point method for a selection of the values of γ. The regression line is a linear regression of all the observed errors. The theoretical order line emanates from the last point and has slope equal to the theoretical convergence order. Both axes are in log10-scale. Note that the regression line is almost completely covered by the theoretical order line for the theoretical order… view at source ↗
Figure 2
Figure 2. Figure 2: Plot of estimated convergence orders for the left-points method. The slope of the linear regressions of observed errors for all the tested values of γ can be seen in black. The theoretical convergence order is also plotted [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Convergence plots for the projection method for a selection of the values of γ. The regression line is a linear regression of all the observed errors. The theoretical order line emanates from the last point and has slope equal to the theoretical convergence order. Both axes are in log10-scale. Note that the regression line is almost completely covered by the theoretical order line in some of the plots [PI… view at source ↗
Figure 4
Figure 4. Figure 4: Plot of estimated convergence orders for the projection method. The slope of the linear regressions of observed errors for all the tested values of γ can be seen in black. The theoretical convergence order is also plotted [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Convergence plots for the relative spectral mean square error for a selection of the values of νs. The other parameters have been fixed. The regression line is a linear regression of four of the observed errors. The theoretical order line emanates from the second to last point and has slope equal to the theoretical convergence order. Both axes are in log10-scale. Note that the regression line is almost com… view at source ↗
Figure 6
Figure 6. Figure 6: Plot of estimated spectral convergence orders. The slope of the linear regres￾sions of four of the observed errors for all the tested values of νs can be seen in black. The other parameters have been fixed. The theoretical convergence order is also plotted. less accuracy close to γ = 1.5. Indeed an inspection of the proof of Theorem 3.4 shows that the constant C(CΠ, t, γ) blows up as γ approaches m + 3 2 .… view at source ↗
Figure 7
Figure 7. Figure 7: A simulation of Equation (3) with the operators in (2) on the unit sphere with α = 0.5, β = 1, γ = 1.5, κ = 2.828, r = 10.0, and σ = 10. The plots use a Mollweide projection. γ 1.250 α 0.500 β 0.750 r 0.816 κ 2.000 γ 2.750 α 0.167 β 0.750 r 0.297 κ 2.000 γ 2.000 α 0.417 β 1.250 r 0.813 κ 3.464 γ 3.000 α 0.250 β 1.250 r 0.416 κ 3.464 [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Realisations of temporal and spatial traces for a selection of parameter values. Each row is from the same simulation, containing a list of the chosen parameter values (left), a plot of the temporal trace (middle), and a plot of the spatial trace (right) [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
read the original abstract

A method for numerical approximation of a new class of fractional parabolic stochastic evolution equations is introduced and analysed. This class of equations has recently been proposed as a space-time extension of the SPDE-method in spatial statistics. A truncation of the spectral basis function expansion is used to discretise in space, and then a quadrature is used to approximate the temporal evolution of each basis coefficient. Strong error bounds are proved both for the spectral and temporal approximations. The method is tested and the results are verified by several numerical experiments.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The manuscript introduces and analyzes a numerical method for a new class of fractional parabolic stochastic evolution equations, recently proposed as a space-time extension of the SPDE-method. Spatial discretization proceeds via truncation of a spectral basis expansion, after which quadrature approximates the temporal evolution of the resulting coefficients. Strong error bounds are proved for the spectral truncation and the temporal quadrature steps, and the method is tested via numerical experiments.

Significance. If the well-posedness, mild-solution existence, and coefficient regularity assumptions hold in the precise setting used for the bounds, the strong-error analysis supplies a rigorous justification for the combined spectral-temporal scheme. This would constitute a useful contribution to numerical methods for fractional SPDEs arising in spatial statistics.

major comments (2)
  1. [Introduction] Introduction (and abstract): the well-posedness of the new fractional SPDE class, the existence of a mild solution, and the regularity needed for the spectral expansion and quadrature are taken from the cited prior proposal rather than re-established or verified here. Because the strong error bounds rest directly on these properties, the manuscript should either derive the required trace-class and moment conditions or state explicitly which hypotheses from the reference are invoked and confirm they match the discretization setting.
  2. [Abstract and main error-analysis sections] The abstract claims strong error bounds are proved, yet the derivation of those bounds is not inspectable from the provided summary; if the analysis in the main body invokes additional regularity not guaranteed by the external reference, the claimed rates may not hold uniformly.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major point below and will revise the manuscript to improve clarity regarding the invoked assumptions while preserving the paper's focus on the numerical analysis.

read point-by-point responses
  1. Referee: [Introduction] Introduction (and abstract): the well-posedness of the new fractional SPDE class, the existence of a mild solution, and the regularity needed for the spectral expansion and quadrature are taken from the cited prior proposal rather than re-established or verified here. Because the strong error bounds rest directly on these properties, the manuscript should either derive the required trace-class and moment conditions or state explicitly which hypotheses from the reference are invoked and confirm they match the discretization setting.

    Authors: We agree that the well-posedness, mild-solution existence, and coefficient regularity are taken from the cited prior proposal. Re-deriving these foundational results would duplicate the reference and fall outside the scope of this work on numerical approximation. In the revised manuscript, we will explicitly list the specific hypotheses invoked (trace-class conditions on the noise covariance and moment bounds on the initial data) in a dedicated remark in the introduction. We will also add a short verification that these hypotheses are compatible with the spectral truncation and temporal quadrature settings used for the error analysis. revision: yes

  2. Referee: [Abstract and main error-analysis sections] The abstract claims strong error bounds are proved, yet the derivation of those bounds is not inspectable from the provided summary; if the analysis in the main body invokes additional regularity not guaranteed by the external reference, the claimed rates may not hold uniformly.

    Authors: The strong error bounds for both the spectral truncation and temporal quadrature are derived in full detail in Sections 3 and 4, relying only on the regularity assumptions stated in the reference (no additional regularity is introduced). The abstract will be revised to explicitly note that the bounds hold under the hypotheses of the cited work. We will also insert forward references to the precise assumptions at each step of the proofs to make the dependence transparent. revision: partial

Circularity Check

0 steps flagged

No significant circularity; error analysis is independent of inputs

full rationale

The manuscript proves strong error bounds for spectral truncation and temporal quadrature applied to a class of fractional SPDEs. The class definition and well-posedness are referenced to prior work, but the derivation of the bounds proceeds from standard mild-solution estimates, eigenfunction expansions, and quadrature error analysis on the resulting coefficient SDEs. No equation reduces by construction to a fitted parameter, no ansatz is smuggled via self-citation, and the central claims do not rely on a self-referential uniqueness theorem. The analysis is self-contained against external benchmarks once the SPDE class is granted, which is the normal situation for approximation papers and yields a score of 0.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The analysis rests on the well-posedness of the recently proposed equation class and on standard properties of the spectral basis for the spatial operator; no free parameters or invented entities are mentioned.

axioms (1)
  • domain assumption The fractional parabolic stochastic evolution equations admit a well-posed mild solution whose spatial operator possesses a complete eigenbasis suitable for truncation.
    Required for the spectral method to be applicable; referenced via the 'new class' and 'SPDE-method in spatial statistics' in the abstract.

pith-pipeline@v0.9.0 · 5600 in / 1197 out tokens · 19456 ms · 2026-05-24T00:00:03.823801+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages

  1. [1]

    Andersson and S

    A. Andersson and S. Larsson, Weak convergence for a spatial approximation of the nonlinear stochastic heat equation , Mathematics of Computation 85.299 (2015), pp. 1335–1358, doi: 10 . 1090/mcom/3016

  2. [2]

    J. M. Angulo et al., Spatiotemporal random fields associated with stochastic fractional Helmholtz and heat equations, Stochastic Environmental Research and Risk Assessment 22.S1 (2007), pp. 3–13, doi: 10.1007/s00477-007-0157-7

  3. [3]

    Bakka et al., Spatial modeling with R-INLA: A review , WIREs Computational Statistics 10.6 (2018), doi: 10.1002/wics.1443

    H. Bakka et al., Spatial modeling with R-INLA: A review , WIREs Computational Statistics 10.6 (2018), doi: 10.1002/wics.1443

  4. [4]

    Baˇnas, H

    L. Baˇnas, H. Yang, and R. Zhu, Sharp Interface Limit of Stochastic Cahn-Hilliard Equation with Singular Noise, Potential Analysis 59.2 (2022), pp. 497–518,doi: 10.1007/s11118-021-09976-3

  5. [5]

    Beesack, A general form of the remainder in Taylor’s theorem , The American Mathematical Monthly 73.1 (1966), p

    P.R. Beesack, A general form of the remainder in Taylor’s theorem , The American Mathematical Monthly 73.1 (1966), p. 64, doi: 10.2307/2313928

  6. [6]

    Biswas, M

    A. Biswas, M. De Le ´on-Contreras, and P.R. Stinga, Harnack inequalities and H ¨older estimates for master equations , SIAM Journal on Mathematical Analysis 53.2 (2021), pp. 2319–2348, doi: 10.1137/20m1339520

  7. [7]

    Biswasand P.R.Stinga, Regularity estimates for nonlocal space-time master equations in bounded domains, Journal of Evolution Equations 21.1 (2020), pp

    A. Biswasand P.R.Stinga, Regularity estimates for nonlocal space-time master equations in bounded domains, Journal of Evolution Equations 21.1 (2020), pp. 503–565, doi: 10.1007/s00028-020- 00590-1

  8. [8]

    Bolin and K

    D. Bolin and K. Kirchner, The rational SPDE approach for Gaussian random fields with general smoothness, Journal of Computational and Graphical Statistics 29.2 (2019), pp. 274–285, doi: 10. 1080/10618600.2019.1665537

  9. [9]

    Bolin, K

    D. Bolin, K. Kirchner, and M. Kov´acs, Numerical solution of fractional elliptic stochastic PDEs with spatial white noise , IMA Journal of Numerical Analysis 40.2 (2018), pp. 1051–1073, doi: 10.1093/imanum/dry091

  10. [10]

    D. Bolin, K.Kirchner, and M.Kov´acs, Weak convergence of Galerkin approximations for fractional elliptic stochastic PDEs with spatial white noise, BIT Numerical Mathematics 58.4 (2018), pp. 881– 906, doi: 10.1007/s10543-018-0719-8

  11. [11]

    Bonaccorsi, Fractional stochastic evolution equations with L´evy noise, Differential and Integral Equations 22.11/12 (2009), doi: 10.57262/die/1356019409

    S. Bonaccorsi, Fractional stochastic evolution equations with L´evy noise, Differential and Integral Equations 22.11/12 (2009), doi: 10.57262/die/1356019409

  12. [12]

    Bonito, W

    A. Bonito, W. Lei, and J.E. Pasciak, On sinc quadrature approximations of fractional powers of regularly accretive operators , Journal of Numerical Mathematics 27.2 (2019), pp. 57–68, doi: 10.1515/jnma-2017-0116

  13. [13]

    Br´ehier, Approximation of the Invariant Measure with an Euler Scheme for Stochastic PDEs Driven by Space-Time White Noise, Potential Analysis 40.1 (2013), pp

    C. Br´ehier, Approximation of the Invariant Measure with an Euler Scheme for Stochastic PDEs Driven by Space-Time White Noise, Potential Analysis 40.1 (2013), pp. 1–40,doi: 10.1007/s11118- 013-9338-9

  14. [14]

    Br´ehier, J

    C. Br´ehier, J. Cui, and X. Wang, Weak error estimates of fully-discrete schemes for the stochastic Cahn-Hilliard equation, 2022, doi: 10.48550/ARXIV.2207.09266. 20 REFERENCES

  15. [15]

    Breit and A

    D. Breit and A. Prohl, Weak error analysis for the stochastic Allen–Cahn equation, Stochastics and Partial Differential Equations: Analysis and Computations (2024), doi: 10.1007/s40072- 024- 00326-z

  16. [16]

    Cohen and A

    D. Cohen and A. Lang, Numerical approximation and simulation of the stochastic wave equation on the sphere, Calcolo 59.3 (2022), doi: 10.1007/s10092-022-00472-7

  17. [17]

    Cressie and C.K

    N. Cressie and C.K. Wikle, Statistics for Spatio-Temporal Data, en, Wiley-Blackwell, Chichester, England, 2011

  18. [18]

    Da Prato and J

    G. Da Prato and J. Zabczyk, Stochastic Equations in Infinite Dimensions , Cambridge University Press, 2014, doi: 10.1017/cbo9781107295513

  19. [19]

    Davies, Spectral Theory and Differential Operators , Cambridge University Press, 1995, doi: 10.1017/cbo9780511623721

    E.B. Davies, Spectral Theory and Differential Operators , Cambridge University Press, 1995, doi: 10.1017/cbo9780511623721

  20. [20]

    Magnetic resonance imaging (MRI): considerations and applications in radiotherapy treatment planning

    A. Debussche and J. Printems, Numerical simulation of the stochastic Korteweg–de Vries equation, Physica D: Nonlinear Phenomena 134.2 (1999), pp. 200–226, doi: 10.1016/s0167- 2789(99) 00072-x

  21. [21]

    Debussche and J

    A. Debussche and J. Vovelle, Diffusion limit for a stochastic kinetic problem, Communications on Pure and Applied Analysis 11.6 (2011), pp. 2305–2326, doi: 10.3934/cpaa.2012.11.2305

  22. [22]

    Desch and S

    W. Desch and S. Londen, Semilinear stochastic integral equations in Lp, in: Parabolic Problems, Springer Basel, 2011, pp. 131–166, doi: 10.1007/978-3-0348-0075-4_8

  23. [23]

    Fahim, E

    K. Fahim, E. Hausenblas, and M. Kov´acs, Some approximation results for mild solutions of stochastic fractional order evolution equations driven by Gaussian noise , Stochastics and Partial Differential Equations: Analysis and Computations 11.3 (2022), pp. 1044–1088, doi: 10 . 1007 / s40072-022-00250-0

  24. [24]

    Fuglstad et al., Exploring a New Class of Non-stationary Spatial Gaussian Random Fields with Varying Local Anisotropy, Statistica Sinica (2014), doi: 10.5705/ss.2013.106w

    G. Fuglstad et al., Exploring a New Class of Non-stationary Spatial Gaussian Random Fields with Varying Local Anisotropy, Statistica Sinica (2014), doi: 10.5705/ss.2013.106w

  25. [25]

    Jansson, M

    E. Jansson, M. Kov´acs, and A. Lang, Surface finite element approximation of spherical Whittle– Mat´ern Gaussian random fields, SIAM Journal on Scientific Computing 44.2 (2022), A825–A842, doi: 10.1137/21m1400717

  26. [26]

    Jentzen and P.E

    A. Jentzen and P.E. Kloeden, Taylor Approximations for Stochastic Partial Differential Equations, Society for Industrial and Applied Mathematics, 2011, doi: 10.1137/1.9781611972016

  27. [27]

    Kelbert, N.N

    M.Ya. Kelbert, N.N. Leonenko, and M.D. Ruiz-Medina, Fractional random fields associated with stochastic fractional heat equations, Advances in Applied Probability 37.1 (2005), pp. 108–133,doi: 10.1239/aap/1113402402

  28. [28]

    K Kirchner and J Willems, Regularity theory for a new class of fractional parabolic stochastic evolution equations, 2023, arXiv: 2205.00248 [math.PR]

  29. [29]

    Kirchner, A

    K. Kirchner, A. Lang, and S. Larsson, Covariance structure of parabolic stochastic partial dif- ferential equations with multiplicative L ´evy noise, Journal of Differential Equations 262.12 (2017), pp. 5896–5927, doi: 10.1016/j.jde.2017.02.021

  30. [30]

    Kirchner and J

    K. Kirchner and J. Willems, Multiple and weak Markov properties in Hilbert spaces with applica- tions to fractional stochastic evolution equations, 2023, doi: 10.48550/ARXIV.2310.13536

  31. [31]

    Kov´acs, A

    M. Kov´acs, A. Lang, and A. Petersson, Approximation of SPDE covariance operators by finite elements: a semigroup approach, IMA Journal of Numerical Analysis 43.3 (2022), pp. 1324–1357, doi: 10.1093/imanum/drac020

  32. [32]

    Kov´acs and S

    M. Kov´acs and S. Larsson, Introduction to stochastic partial differential(2010)

  33. [33]

    M. Kov´acs, S.Larsson, and F.Lindgren, Strong convergence of the finite element method with trun- cated noise for semilinear parabolic stochastic equations with additive noise, Numerical Algorithms 53.2–3 (2009), pp. 309–320, doi: 10.1007/s11075-009-9281-4

  34. [34]

    Kov´acs, S

    M. Kov´acs, S. Larsson, and F. Lindgren, Weak convergence of finite element approximations of linear stochastic evolution equations with additive noise, BIT Numerical Mathematics 52.1 (2011), pp. 85–108, doi: 10.1007/s10543-011-0344-2

  35. [35]

    Kov´acs, S

    M. Kov´acs, S. Larsson, and F. Saedpanah, Mittag–Leffler Euler Integrator for a Stochastic Frac- tional Order Equation with Additive Noise, SIAM Journal on Numerical Analysis 58.1 (2020), pp. 66– 85, doi: 10.1137/18m1177895

  36. [36]

    Labl´ee, Spectral Theory in Riemannian Geometry, EMS Press, Zuerich, Switzerland, 2015

    O. Labl´ee, Spectral Theory in Riemannian Geometry, EMS Press, Zuerich, Switzerland, 2015

  37. [37]

    Lang and I

    A. Lang and I. Motschan-Armen, Euler–Maruyama approximations of the stochastic heat equation on the sphere , Journal of Computational Dynamics 11.1 (2024), pp. 23–42, doi: 10.3934/jcd. 2023012

  38. [38]

    Lang and C

    A. Lang and C. Schwab, Isotropic Gaussian random fields on the sphere: Regularity, fast simulation and stochastic partial differential equations , The Annals of Applied Probability 25.6 (2015), doi: 10.1214/14-aap1067. REFERENCES 21

  39. [39]

    Lindgren, D

    F. Lindgren, D. Bolin, and H. Rue, The SPDE approach for Gaussian and non-Gaussian fields: 10 years and still running , Spatial Statistics 50 (2022), p. 100599, doi: 10.1016/j.spasta.2022. 100599

  40. [40]

    A., Hekker, S., Stello, D., Guti ´errez-Soto, J., Handberg, R., Huber, D., et al

    F. Lindgren, H. Rue, and J. Lindstr¨om, An explicit link between Gaussian fields and Gaussian Markov random fields: The stochastic partial differential equation approach , Journal of the Royal Statistical Society Series B: Statistical Methodology 73.4 (2011), pp. 423–498, doi: 10.1111/j. 1467-9868.2011.00777.x

  41. [41]

    Lindgren et al., A diffusion-based spatio-temporal extension of Gaussian Mat ´ern fields, SORT 48.1 (2024), pp

    F. Lindgren et al., A diffusion-based spatio-temporal extension of Gaussian Mat ´ern fields, SORT 48.1 (2024), pp. 3–66, doi: 10.57645/20.8080.02.13

  42. [42]

    Litsg˚ard and K

    M. Litsg˚ard and K. Nystr¨om, On local regularity estimates for fractional powers of parabolic operators with time-dependent measurable coefficients, Journal of Evolution Equations 23.1 (2022), doi: 10.1007/s00028-022-00844-0

  43. [43]

    Lord, C.E

    G.J. Lord, C.E. Powell, and T. Shardlow, An Introduction to Computational Stochastic PDEs , Cambridge University Press, 2014, doi: 10.1017/cbo9781139017329

  44. [44]

    Nystr¨om and O

    K. Nystr¨om and O. Sande, Extension properties and boundary estimates for a fractional heat operator, Nonlinear Analysis 140 (2016), pp. 29–37, doi: 10.1016/j.na.2016.02.027

  45. [45]

    A. Prohl and Y.Wang, Strong error estimates for a space-time discretization of the linear-quadratic control problem with the stochastic heat equation with linear noise , IMA Journal of Numerical Analysis 42.4 (2021), pp. 3386–3429, doi: 10.1093/imanum/drab069

  46. [46]

    Stinga and J.L

    P.R. Stinga and J.L. Torrea, Regularity theory and extension problem for fractional nonlocal parabolic equations and the master equation, SIAM Journal on Mathematical Analysis 49.5 (2017), pp. 3893–3924, doi: 10.1137/16m1104317

  47. [47]

    Strang and G.J

    G. Strang and G.J. Fix, An analysis of the finite element method, Englewood Cliffs: Prentice-Hall, 1969

  48. [48]

    Willems, Dirichlet problems associated to abstract nonlocal space-time differential operators , 2024, doi: 10.48550/ARXIV.2404.11289

    J. Willems, Dirichlet problems associated to abstract nonlocal space-time differential operators , 2024, doi: 10.48550/ARXIV.2404.11289

  49. [49]

    Yan, Galerkin finite element methods for stochastic parabolic partial differential equations, SIAM Journal on Numerical Analysis 43.4 (2005), pp

    Y. Yan, Galerkin finite element methods for stochastic parabolic partial differential equations, SIAM Journal on Numerical Analysis 43.4 (2005), pp. 1363–1384, doi: 10.1137/040605278