Search

Solving large-scale continuous-time algebraic Riccati equations is a significant challenge in various control theory applications. This work demonstrates that when the matrix coefficients of the equation are quasiseparable, the solution also exhibits numerical quasiseparability. This property enables us to develop two efficient Riccati solvers. The first solver is applicable to the general quasiseparable case, while the second is tailored to the particular case of banded coefficients. Numerical experiments confirm the effectiveness of the proposed algorithms on both synthetic examples and case studies from the control of partial differential equations and agent-based models.

We present an approach for the optimization of irrigation in a Richards' equation framework. We introduce a proper cost functional, aimed at minimizing the amount of water provided by irrigation, at the same time maximizing the root water uptake, which is modeled by a sink term in the continuity equation. The control is acting on the boundary of the dynamics and due to the nature of the mathematical problem we use a State-Dependent Riccati approach which provides suboptimal control in feedback form, applied to the system of ODEs resulting from the Richards' equation semidiscretization in space. The problem is tested with existing hydraulic parameters, also considering proper root water uptake functions. The numerical simulations also consider the presence of noise in the model to further validate the use of a feedback control approach.

Approximation of high dimensional functions is in the focus of machine learning and data-based scientific computing. In many applications, empirical risk minimisation techniques over nonlinear model classes are employed. Neural networks, kernel methods and tensor decomposition techniques are among the most popular model classes. We provide a numerical study comparing the performance of these methods on various high-dimensional functions with focus on optimal control problems, where the collection of the dataset is based on the application of the State-Dependent Riccati Equation.