Towards an Efficient Finite Element Method for the Integral Fractional Laplacian on Polygonal Domains

We explore the connection between fractional order partial differential equations in two or more spatial dimensions with boundary integral operators to develop techniques that enable one to efficiently tackle the integral fractional Laplacian. In particular, we develop techniques for the treatment of the dense stiffness matrix including the computation of the entries, the efficient assembly and storage of a sparse approximation and the efficient solution of the resulting equations. The main idea consists of generalising proven techniques for the treatment of boundary integral equations to general fractional orders. Importantly, the approximation does not make any strong assumptions on the shape of the underlying domain and does not rely on any special structure of the matrix that could be exploited by fast transforms. We demonstrate the flexibility and performance of this approach in a couple of two-dimensional numerical examples

Statistics of non-linear stochastic dynamical systems under L\'evy noises by a convolution quadrature approach

This paper describes a novel numerical approach to find the statistics of the non-stationary response of scalar non-linear systems excited by L\'evy white noises. The proposed numerical procedure relies on the introduction of an integral transform of Wiener-Hopf type into the equation governing the characteristic function. Once this equation is rewritten as partial integro-differential equation, it is then solved by applying the method of convolution quadrature originally proposed by Lubich, here extended to deal with this particular integral transform. The proposed approach is relevant for two reasons: 1) Statistics of systems with several different drift terms can be handled in an efficient way, independently from the kind of white noise; 2) The particular form of Wiener-Hopf integral transform and its numerical evaluation, both introduced in this study, are generalizations of fractional integro-differential operators of potential type and Gr\"unwald-Letnikov fractional derivatives, respectively.Comment: 20 pages, 5 figure

Trees and asymptotic developments for fractional stochastic differential equations

In this paper we consider a n-dimensional stochastic differential equation driven by a fractional Brownian motion with Hurst parameter H>1/3. After solving this equation in a rather elementary way, following the approach of Gubinelli, we show how to obtain an expansion for E[f(X\_t)] in terms of t, where X denotes the solution to the SDE and f:R^n->R is a regular function. With respect to the work by Baudoin and Coutin, where the same kind of problem is considered, we try an improvement in three different directions: we are able to take a drift into account in the equation, we parametrize our expansion with trees (which makes it easier to use), and we obtain a sharp control of the remainder.Comment: 46 page

Good (and Not So Good) practices in computational methods for fractional calculus

The solution of fractional-order differential problems requires in the majority of cases the use of some computational approach. In general, the numerical treatment of fractional differential equations is much more difficult than in the integer-order case, and very often non-specialist researchers are unaware of the specific difficulties. As a consequence, numerical methods are often applied in an incorrect way or unreliable methods are devised and proposed in the literature. In this paper we try to identify some common pitfalls in the use of numerical methods in fractional calculus, to explain their nature and to list some good practices that should be followed in order to obtain correct results