Exponential monotonicity of quadratic forms in ODEs and preserving methods
Future Generation Computer Systems - Selected papers on theoretical and computational aspects of structural dynamical systems in linear algebra and control
RCMS: right correction Magnus series approach for oscillatory ODEs
Journal of Computational and Applied Mathematics
Brief paper: Computing the L2 gain for linear periodic continuous-time systems
Automatica (Journal of IFAC)
Dynamic Dual-Tracer PET Reconstruction
IPMI '09 Proceedings of the 21st International Conference on Information Processing in Medical Imaging
Tracer kinetics guided dynamic PET reconstruction
IPMI'07 Proceedings of the 20th international conference on Information processing in medical imaging
Multimodal admittance method in waveguides and singularity behavior at high frequencies
Journal of Computational and Applied Mathematics
Journal of Computational and Applied Mathematics
Cardiac motion recovery: continuous dynamics, discrete measurements, and optimal estimation
MICCAI'06 Proceedings of the 9th international conference on Medical Image Computing and Computer-Assisted Intervention - Volume Part I
Numerical integration of the differential riccati equation: a high performance computing approach
VECPAR'04 Proceedings of the 6th international conference on High Performance Computing for Computational Science
Time-averaging and exponential integrators for non-homogeneous linear IVPs and BVPs
Applied Numerical Mathematics
Magnus integrators for solving linear-quadratic differential games
Journal of Computational and Applied Mathematics
Hi-index | 0.01 |
This paper introduces a new class of methods, which we call Möbius schemes, for the numerical solution of matrix Riccati differential equations. The approach is based on viewing the Riccati equation in its natural geometric setting, as a flow on the Grassmannian of m-dimensional subspaces of an (n+m)-dimensional vector space. Since the Grassmannians are compact differentiable manifolds, and the coefficients of the equation are assumed continuous, there are no singularities or intrinsic instabilities in the associated flow. The presence of singularities and numerical instabilities is an artifact of the coordinate system, but since Möbius schemes are based on the natural geometry, they are able to deal with numerical instability and pass accurately through the singularities. A number of examples are given to demonstrate these properties.