Abstract:In this paper, we present numerical methods suitable for solving convex quadratic fractional differential equation (FDE) constrained optimization problems, with box constraints on the state and/or control variables. We develop an alternating direction method of multipliers (ADMM) framework, which uses preconditioned Krylov subspace solvers for the resulting subproblems. The latter allows us to tackle a range of partial differential equation (PDE) optimization problems with box constraints, posed on space-time domains, that were previously out of the reach of state-of-the-art preconditioners. In particular, by making use of the powerful generalized locally Toeplitz (GLT) sequences theory, we show that any existing GLT structure present in the problem matrices is preserved by ADMM, and we propose some preconditioning methodologies that could be used within the solver, to demonstrate the generality of the approach. Focusing on convex quadratic programs with time-dependent 2-dimensional FDE constraints, we derive multilevel circulant preconditioners, which may be embedded within Krylov subspace methods, for solving the ADMM subproblems. Discretized versions of FDEs involve large dense linear systems. In order to overcome this difficulty, we design a recursive linear algebra, which is based on the fast Fourier transform (FFT). We manage to keep the storage requirements linear, with respect to the grid size , while ensuring an order computational complexity per iteration of the Krylov solver. We implement the proposed method and demonstrate its scalability, generality, and efficiency through a series of experiments over different setups of the FDE optimization problem.

Fractional BDF Methods for Solving Fractional Differential Matrix Equations

A New Method For Solving Fractional And Classical Differential Equations Based On a New Generalized Fractional Power Series

Operational-matrix-based Algorithm for Differential Equations of Fractional Order with Dirichlet Boundary Conditions

Numerical solution of multidimensional time‐space fractional differential equations of distributed order with Riesz derivative

Approximate Time‐Fractional Differential Equations

Numerical Method for Multi-Dimensional Coupled Forward-Backward Stochastic Differential Equations Based on Fractional Fourier Fast Transform

Higher order numerical methods for fractional delay differential equations

A Convolution Method for Numerical Solution of Backward Stochastic Differential Equations Based on the Fractional FFT

Applications to Riccati and Lorenz Systems

A new numerical method to solve fractional differential equations in terms of Caputo-Fabrizio derivatives

On Failed Methods of Fractional Differential Equations: The Case of Multi-step Generalized Differential Transform Method

Solution of Fractional Order Riccati Differential Equations Using Euler Wavelet Method

A new approach to solving multi-order fractional equations using BEM and Chebyshev matrix

Solving Fractional Order Differential Equations by Using Fractional Radial Basis Function Neural Network

Numerical Study of Caputo Fractional-Order Differential Equations by Developing New Operational Matrices of Vieta–Lucas Polynomials

The BEM and DRBEM schemes for the numerical solution of the two-dimensional time-fractional diffusion-wave equations

Matrix approach to discrete fractional calculus II: partial fractional differential equations

Adaptive finite element method for fractional differential equations using hierarchical matrices

Numerical Approximations for Fractional Differential Equations

Fast Solution Methods for Convex Quadratic Optimization of Fractional Differential Equations

Numerical study of generalized modified Caputo fractional differential equations