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.

RBF-FD Solution for a Financial Partial-Integro Differential Equation Utilizing the Generalized Multiquadric Function

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

Efficient Scheme for the Economic Heston–Hull–White Problem Using Novel RBF-FD Coefficients Derived from Multiquadric Function Integrals

An efficient radial basis function generated finite difference meshfree scheme to price multi-dimensional PDEs in financial options

Multiquadric based RBF-HFD approximation formulas and convergence properties

Numerical investigation of high-dimensional option pricing PDEs by utilizing a hybrid radial basis function - finite difference procedure

An unfitted RBF-FD method in a least-squares setting for elliptic PDEs on complex geometries

A RBF partition of unity collocation method based on finite difference for initial-boundary value problems

A novel Hermite RBF-based differential quadrature method for solving two-dimensional variable-order time fractional advection-diffusion equation with Neumann boundary condition

Numerical treatment of multi-dimensional time-fractional Benjamin-Bona-Mahony-Burgers equations in arbitrary domains with a novel improvised RBF-based method

An Algorithm for the Solution of Integro-Fractional Differential Equations with a Generalized Symmetric Singular Kernel

Method of Order Reduction for the High-Dimensional Convection-Diffusion-Reaction Equation with Robin Boundary Conditions Based on MQ RBF-FD

Fast Solution Methods for Convex Quadratic Optimization of Fractional Differential Equations

Solving High-Dimensional Partial Integral Differential Equations: The Finite Expression Method

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

A stabilized radial basis-finite difference (RBF-FD) method with hybrid kernels

An Efficient QSC Approximation of Variable-Order Time-Fractional Mobile-Immobile Diffusion Equations with Variably Diffusive Coefficients

Hyperviscosity-Based Stabilization for Radial Basis Function-Finite Difference (RBF-FD) Discretizations of Advection-Diffusion Equations

A New Rbf-Trefftz Meshless Method For Partial Differential Equations

A novel local Hermite radial basis function‐based differential quadrature method for solving two‐dimensional variable‐order time fractional advection–diffusion equation with Neumann boundary conditions

An IRK‐QCD Scheme for the Space Fractional Diffusion Equation and Block‐Lower‐Triangle Preconditioners for the Corresponding Linear Systems