Related papers: Improved Matlab code for Lyapunov exponents of fra…
We propose a novel numerical algorithm utilizing model reduction for computing solutions to stationary partial differential equations involving the spectral fractional Laplacian. Our approach utilizes a known characterization of the…
We point out that two of Milne's fourth-order integrators are well-suited to bit-reversible simulations. The fourth-order method improves on the accuracy of Levesque and Verlet's algorithm and simplifies the definition of the velocity $v$…
We extend Lyapunov--type Mittag--Leffler stability analysis for fuzzy nonlinear fractional differential equations (Caputo sense) and introduce a family of stronger, broadly applicable stability results. In particular, we develop (i) uniform…
LU and Cholesky matrix factorization algorithms are core subroutines used to solve systems of linear equations (SLEs) encountered while solving an optimization problem. Standard factorization algorithms are highly efficient but remain…
Atangana and Baleanu proposed a new fractional derivative with non-local and no-singular Mittag-Leffler kernel to solve some problems proposed by researchers in the field of fractional calculus. This new derivative is better to describe…
We analyze stochastic partial differential equations (SPDEs) with quadratic nonlinearities close to a change of stability. To this aim we compute finite-time Lyapunov exponents (FTLEs), observing a change of sign based on the interplay…
A higher-order numerical method is presented for scalar valued, coupled forward-backward stochastic differential equations. Unlike most classical references, the forward component is not only discretized by an Euler-Maruyama approximation…
We present a new linearly implicit exponential integrator that preserves the polynomial first integrals or Lyapunov functions for the conservative and dissipative stiff equations, respectively. The method is tested by both oscillated…
We present LinApart, a routine designed for efficiently performing the univariate partial fraction decomposition of large symbolic expressions. Our method is based on an explicit closed formula for the decomposition of rational functions…
Computationally efficient numerical methods for high-order approximations of convolution integrals involving weakly singular kernels find many practical applications including those in the development of fast quadrature methods for…
In this paper, we investigate fast algorithms in the small fraction order regime to approximate the Caputo derivative $^C_0D_t^\alpha u(t)$ when $\alpha$ is small. We focus on two fast algorithms, i.e. FIR and FIDR, both relying on the…
In this paper, we consider some aspects of the numerical analysis of the mathematical model of fractional Duffing with a derivative of variable fractional order of the Riemann-Liouville type. Using numerical methods: an explicit…
We present two new adaptive quadrature routines. Both routines differ from previously published algorithms in many aspects, most significantly in how they represent the integrand, how they treat non-numerical values of the integrand, how…
A very simple and efficient local variational iteration method for solving problems of nonlinear science is proposed in this paper. The analytical iteration formula of this method is derived first using a general form of first order…
The efficient numerical solution of fractional differential equations has been recently tackled through the definition of Fractional HBVMs (FHBVMs), a class of Runge-Kutta type methods. Corresponding Matlab (c) codes have been also made…
This research investigates using a mixed-precision iterative refinement method using posit numbers instead of the standard IEEE floating-point format. The method is applied to solve a general linear system represented by the equation $Ax =…
In this paper, we present a framework for Stability Analysis of Systems of Coupled Linear Partial-Differential Equations. The class of PDE systems considered in this paper includes parabolic, elliptic and hyperbolic systems with Dirichelet,…
Contrary to integer order derivative, the fractional-order derivative of a non-constant periodic function is not a periodic function with the same period, as a consequence of this property the time-invariant fractional order system does not…
We propose two coding schemes for distributed matrix multiplication in the presence of stragglers. These coding schemes are adaptations of LT codes and Raptor codes to distributed matrix multiplication and are termed \emph{factored LT (FLT)…
FTLE (Finite Time Lyapunov Exponent) computation is one of the standard approaches to Lagrangian flow analysis. The main features of interest in FTLE fields are ridges that represent hyperbolic Lagrangian Coherent Structures. FTLE ridges…