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.

An Efficient Iteration Method for Toeplitz-Plus-Band Triangular Systems Generated from Fractional Ordinary Differential Equation

A Novel Method for Linear Systems of Fractional Ordinary Differential Equations with Applications to Time-Fractional PDEs

Differential Operator Multiplication Method for Fractional Differential Equations

Diagonal and normal with Toeplitz-block splitting iteration method for space fractional coupled nonlinear Schrödinger equations with repulsive nonlinearities

Numerical algorithms for the fast and reliable solution of periodic tridiagonal Toeplitz linear systems

Multigrid method for symmetric Toeplitz tridiagonal matrix: Convergence analysis and application to fractional Feynman-Kac equation

Efficient numerical algorithms of time fractional telegraph‐type equations involving Hadamard derivatives

A Parallel Algorithm for the Two-Dimensional Time Fractional Diffusion Equation with Implicit Difference Method

Efficiently solving fractional delay differential equations of variable order via an adjusted spectral element approach

Exponentially accurate spectral and spectral element methods for fractional ODEs

A New Operational Matrices-Based Spectral Method For Multi-Order Fractional Problems

Numerical investigation of two fractional operators for time fractional delay differential equation

Fast Solution Methods for Convex Quadratic Optimization of Fractional Differential Equations

Computational Challenge of Fractional Differential Equations and the Potential Solutions: A Survey

Novel operational matrices-based method for solving fractional-order delay differential equations via shifted Gegenbauer polynomials

A fast method for variable-order Caputo fractional derivative with applications to time-fractional diffusion equations

A family of iterative methods to solve nonlinear problems with applications in fractional differential equations

Two Finite Difference Schemes for Multi-Dimensional Fractional Wave Equations with Weakly Singular Solutions

The numerical solution of fractional integral equations via orthogonal polynomials in fractional powers

Superfast Tikhonov Regularization of Toeplitz Systems

Spectral Collocation Method for Fractional Differential/Integral Equations with Generalized Fractional Operator