4,423 research outputs found
Harmonic Exponential Families on Manifolds
In a range of fields including the geosciences, molecular biology, robotics
and computer vision, one encounters problems that involve random variables on
manifolds. Currently, there is a lack of flexible probabilistic models on
manifolds that are fast and easy to train. We define an extremely flexible
class of exponential family distributions on manifolds such as the torus,
sphere, and rotation groups, and show that for these distributions the gradient
of the log-likelihood can be computed efficiently using a non-commutative
generalization of the Fast Fourier Transform (FFT). We discuss applications to
Bayesian camera motion estimation (where harmonic exponential families serve as
conjugate priors), and modelling of the spatial distribution of earthquakes on
the surface of the earth. Our experimental results show that harmonic densities
yield a significantly higher likelihood than the best competing method, while
being orders of magnitude faster to train.Comment: fixed typ
Methods for suspensions of passive and active filaments
Flexible filaments and fibres are essential components of important complex
fluids that appear in many biological and industrial settings. Direct
simulations of these systems that capture the motion and deformation of many
immersed filaments in suspension remain a formidable computational challenge
due to the complex, coupled fluid--structure interactions of all filaments, the
numerical stiffness associated with filament bending, and the various
constraints that must be maintained as the filaments deform. In this paper, we
address these challenges by describing filament kinematics using quaternions to
resolve both bending and twisting, applying implicit time-integration to
alleviate numerical stiffness, and using quasi-Newton methods to obtain
solutions to the resulting system of nonlinear equations. In particular, we
employ geometric time integration to ensure that the quaternions remain unit as
the filaments move. We also show that our framework can be used with a variety
of models and methods, including matrix-free fast methods, that resolve low
Reynolds number hydrodynamic interactions. We provide a series of tests and
example simulations to demonstrate the performance and possible applications of
our method. Finally, we provide a link to a MATLAB/Octave implementation of our
framework that can be used to learn more about our approach and as a tool for
filament simulation
- …