Related papers: Fast and accurate approximations to fractional pow…
This paper introduces a very fast method for the computation of the resolvent of fractional powers of operators. The analysis is kept in the continuous setting of (potentially unbounded) self adjoint positive operators in Hilbert spaces.…
We investigate the rational approximation of fractional powers of unbounded positive operators attainable with a specific integral representation of the operator function. We provide accurate error bounds by exploiting classical results in…
In this paper we are interested in the approximation of fractional powers of self-adjoint positive operators. Starting from the integral representation of the operators, we apply the trapezoidal rule combined with a single-exponential and a…
Some mathematical models of applied problems lead to the need of solving boundary value problems with a fractional power of an elliptic operator. In a number of works, approximations of such a nonlocal operator are constructed on the basis…
We consider the finite element approximation of fractional powers of regularly accretive operators via the Dunford-Taylor integral approach. We use a sinc quadrature scheme to approximate the Balakrishnan representation of the negative…
Problems of the numerical solution of the Cauchy problem for a first-order differential-operator equation are discussed. A fundamental feature of the problem under study is that the equation includes a fractional power of the self-adjoint…
We study the numerical approximation of fractional powers of accretive operators in this paper. Namely, if $A$ is the accretive operator associated with an accretive sesquilinear form $A(\cdot,\cdot)$ defined on a Hilbert space $\mathbb V$…
We study a reliable pole selection for the rational approximation of the resolvent of fractional powers of operators in both the finite and infinite dimensional setting. The analysis exploits the representation in terms of hypergeometric…
Boundary value problem for a fractional power of an elliptic operator is considered. An integral representation by means of a standard solution problem for parabolic equations is used to solve such problems. Quadrature generalized…
We present a new formula for umbral operators that yields three main insights. First, it makes explicit a connection between umbral calculus and iteration theory. Second, it leads naturally to a definition of fractional exponents of umbral…
We show how rational function approximations to the logarithm, such as $\log z \approx (z^2 - 1)/(z^2 + 6z + 1)$, can be turned into fast algorithms for approximating the determinant of a very large matrix. We empirically demonstrate that…
In this paper, we introduce the new construction of fractional derivatives and integrals with respect to a function, based on a matrix approach. We believe that this is a powerful tool in both analytical and numerical calculations. We begin…
We present efficient approximation of the error function obtained by Fourier expansion of the exponential function $\exp [{- {(t - 2 \sigma)^2}/4}]$. The error analysis reveals that it is highly accurate and can generate numbers that match…
Classical Laguerre spectral approximations are highly effective on the half-line when the target function is smooth in the usual polynomial scale. However, their accuracy deteriorates for nonsmooth functions. Such behavior appears naturally…
To broaden the range of applicability of variable-order fractional differential models, reliable numerical approaches are needed to solve the model equation. In this paper, we develop Laguerre spectral collocation methods for solving…
We consider the approximation of the inverse square root of regularly accretive operators in Hilbert spaces. The approximation is of rational type and comes from the use of the Gauss-Legendre rule applied to a special integral formulation…
In this paper we explore the theory of fractional powers of non-negative (and not necessarily self-adjoint) operators and its amazing relationship with the Chebyshev polynomials of the second kind to obtain results of existence, regularity…
An efficient algorithm and a Fortran 90 module (LaguerrePol) for computing Laguerre polynomials $L^{(\alpha)}_n(z)$ are presented. The standard three-term recurrence relation satisfied by the polynomials and different types of asymptotic…
Operator convex functions defined on the positive half-line play a prominent role in the theory of quantum information, where they are used to define quantum $f$-divergences. Such functions admit integral representations in terms of…
We introduce an efficient algorithm for computing fractional integrals and derivatives and apply it for solving problems of the calculus of variations of fractional order. The proposed approximations are particularly useful for solving…