27,297 research outputs found
Fourier-Splitting methods for the dynamics of rotating Bose-Einstein condensates
We present a new method to propagate rotating Bose-Einstein condensates
subject to explicitly time-dependent trapping potentials. Using algebraic
techniques, we combine Magnus expansions and splitting methods to yield any
order methods for the multivariate and nonautonomous quadratic part of the
Hamiltonian that can be computed using only Fourier transforms at the cost of
solving a small system of polynomial equations. The resulting scheme solves the
challenging component of the (nonlinear) Hamiltonian and can be combined with
optimized splitting methods to yield efficient algorithms for rotating
Bose-Einstein condensates. The method is particularly efficient for potentials
that can be regarded as perturbed rotating and trapped condensates, e.g., for
small nonlinearities, since it retains the near-integrable structure of the
problem. For large nonlinearities, the method remains highly efficient if
higher order p > 2 is sought. Furthermore, we show how it can adapted to the
presence of dissipation terms. Numerical examples illustrate the performance of
the scheme.Comment: 15 pages, 4 figures, as submitted to journa
Sparse Automatic Differentiation for Large-Scale Computations Using Abstract Elementary Algebra
Most numerical solvers and libraries nowadays are implemented to use
mathematical models created with language-specific built-in data types (e.g.
real in Fortran or double in C) and their respective elementary algebra
implementations. However, built-in elementary algebra typically has limited
functionality and often restricts flexibility of mathematical models and
analysis types that can be applied to those models. To overcome this
limitation, a number of domain-specific languages with more feature-rich
built-in data types have been proposed. In this paper, we argue that if
numerical libraries and solvers are designed to use abstract elementary algebra
rather than language-specific built-in algebra, modern mainstream languages can
be as effective as any domain-specific language. We illustrate our ideas using
the example of sparse Jacobian matrix computation. We implement an automatic
differentiation method that takes advantage of sparse system structures and is
straightforward to parallelize in MPI setting. Furthermore, we show that the
computational cost scales linearly with the size of the system.Comment: Submitted to ACM Transactions on Mathematical Softwar
Matrix products for the synthesis of stationary time series with a priori prescribed joint distributions
Inspired from non-equilibrium statistical physics models, a general framework
enabling the definition and synthesis of stationary time series with a priori
prescribed and controlled joint distributions is constructed. Its central
feature consists of preserving for the joint distribution the simple product
struc- ture it has under independence while enabling to input con- trolled and
prescribed dependencies amongst samples. To that end, it is based on products
of d-dimensional matrices, whose entries consist of valid distributions. The
statistical properties of the thus defined time series are studied in details.
Having been able to recast this framework into that of Hidden Markov Models
enabled us to obtain an efficient synthesis procedure. Pedagogical well-chosen
examples (time series with the same marginal distribution, same covariance
function, but different joint distributions) aim at illustrating the power and
potential of the approach and at showing how targeted statistical prop- erties
can be actually prescribed.Comment: 4 pages, 2 figures, conference publication published in IEEE
International Conference on Acoustics, Speech and Signal Processing (ICASSP),
201
Hybrid approximate message passing
Gaussian and quadratic approximations of message passing algorithms on graphs have attracted considerable recent attention due to their computational simplicity, analytic tractability, and wide applicability in optimization and statistical inference problems. This paper presents a systematic framework for incorporating such approximate message passing (AMP) methods in general graphical models. The key concept is a partition of dependencies of a general graphical model into strong and weak edges, with the weak edges representing interactions through aggregates of small, linearizable couplings of variables. AMP approximations based on the Central Limit Theorem can be readily applied to aggregates of many weak edges and integrated with standard message passing updates on the strong edges. The resulting algorithm, which we call hybrid generalized approximate message passing (HyGAMP), can yield significantly simpler implementations of sum-product and max-sum loopy belief propagation. By varying the partition of strong and weak edges, a performance--complexity trade-off can be achieved. Group sparsity and multinomial logistic regression problems are studied as examples of the proposed methodology.The work of S. Rangan was supported in part by the National Science Foundation under Grants 1116589, 1302336, and 1547332, and in part by the industrial affiliates of NYU WIRELESS. The work of A. K. Fletcher was supported in part by the National Science Foundation under Grants 1254204 and 1738286 and in part by the Office of Naval Research under Grant N00014-15-1-2677. The work of V. K. Goyal was supported in part by the National Science Foundation under Grant 1422034. The work of E. Byrne and P. Schniter was supported in part by the National Science Foundation under Grant CCF-1527162. (1116589 - National Science Foundation; 1302336 - National Science Foundation; 1547332 - National Science Foundation; 1254204 - National Science Foundation; 1738286 - National Science Foundation; 1422034 - National Science Foundation; CCF-1527162 - National Science Foundation; NYU WIRELESS; N00014-15-1-2677 - Office of Naval Research
Finite volume approach for the instationary Cosserat rod model describing the spinning of viscous jets
The spinning of slender viscous jets can be described asymptotically by
one-dimensional models that consist of systems of partial and ordinary
differential equations. Whereas the well-established string models possess only
solutions for certain choices of parameters and set-ups, the more sophisticated
rod model that can be considered as -regularized string is generally
applicable. But containing the slenderness ratio explicitely in the
equations complicates the numerical treatment. In this paper we present the
first instationary simulations of a rod in a rotational spinning process for
arbitrary parameter ranges with free and fixed jet end, for which the hitherto
investigations longed. So we close an existing gap in literature. The numerics
is based on a finite volume approach with mixed central, up- and down-winded
differences, the time integration is performed by stiff accurate Radau methods
- …