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.

Efficient and Accurate Spectral Method Using Generalized Jacobi Functions for Solving Riesz Fractional Differential Equations

Efficient Chebyshev Spectral Method for Solving Linear Elliptic Pdes Using Quasi-Inverse Technique

Generalized Jacobi Functions and Their Applications to Fractional Differential Equations.

Generalized Finite Difference Method for a Class of Multidimensional Space-Fractional Diffusion Equations

Efficient and Accurate Spectral Method for Solving Fractional Differential Equations on the Half Line Using Orthogonal Generalized Rational Jacobi Functions

Exponentially accurate spectral and spectral element methods for fractional ODEs

An Efficient and Accurate Numerical Method for the Spectral Fractional Laplacian Equation

Efficient spectral–Galerkin methods for fractional partial differential equations with variable coefficients

An efficient spectral method for the fractional Schrödinger equation on the real line

An Efficient Numerical Method for Solving a Class of Nonlinear Fractional Differential Equations and Error Estimates

The Spectral Collocation Method for Efficiently Solving PDEs with Fractional Laplacian.

Adaptive finite element method for fractional differential equations using hierarchical matrices

Efficient Numerical Algorithms for Three-Dimensional Fractional Partial Differential Equations

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

Two Efficient Generalized Laguerre Spectral Algorithms for Fractional Initial Value Problems

Sine-transform-based fast solvers for Riesz fractional nonlinear Schrödinger equations with attractive nonlinearities

Fast Solution Methods for Convex Quadratic Optimization of Fractional Differential Equations

A spectrally accurate step-by-step method for the numerical solution of fractional differential equations

A search for a spectral technique to solve nonlinear fractional differential equations

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

A priori and a posteriori error estimates of a space–time Petrov–Galerkin spectral method for time-fractional diffusion equation