271,630 research outputs found
Adaptive Consensus: A network pruning approach for decentralized optimization
We consider network-based decentralized optimization problems, where each
node in the network possesses a local function and the objective is to
collectively attain a consensus solution that minimizes the sum of all the
local functions. A major challenge in decentralized optimization is the
reliance on communication which remains a considerable bottleneck in many
applications. To address this challenge, we propose an adaptive randomized
communication-efficient algorithmic framework that reduces the volume of
communication by periodically tracking the disagreement error and judiciously
selecting the most influential and effective edges at each node for
communication. Within this framework, we present two algorithms: Adaptive
Consensus (AC) to solve the consensus problem and Adaptive Consensus based
Gradient Tracking (AC-GT) to solve smooth strongly convex decentralized
optimization problems. We establish strong theoretical convergence guarantees
for the proposed algorithms and quantify their performance in terms of various
algorithmic parameters under standard assumptions. Finally, numerical
experiments showcase the effectiveness of the framework in significantly
reducing the information exchange required to achieve a consensus solution.Comment: 35 pages, 3 figure
New numerical approaches for modeling thermochemical convection in a compositionally stratified fluid
Seismic imaging of the mantle has revealed large and small scale
heterogeneities in the lower mantle; specifically structures known as large low
shear velocity provinces (LLSVP) below Africa and the South Pacific. Most
interpretations propose that the heterogeneities are compositional in nature,
differing in composition from the overlying mantle, an interpretation that
would be consistent with chemical geodynamic models. Numerical modeling of
persistent compositional interfaces presents challenges, even to
state-of-the-art numerical methodology. For example, some numerical algorithms
for advecting the compositional interface cannot maintain a sharp compositional
boundary as the fluid migrates and distorts with time dependent fingering due
to the numerical diffusion that has been added in order to maintain the upper
and lower bounds on the composition variable and the stability of the advection
method. In this work we present two new algorithms for maintaining a sharper
computational boundary than the advection methods that are currently openly
available to the computational mantle convection community; namely, a
Discontinuous Galerkin method with a Bound Preserving limiter and a
Volume-of-Fluid interface tracking algorithm. We compare these two new methods
with two approaches commonly used for modeling the advection of two distinct,
thermally driven, compositional fields in mantle convection problems; namely,
an approach based on a high-order accurate finite element method advection
algorithm that employs an artificial viscosity technique to maintain the upper
and lower bounds on the composition variable as well as the stability of the
advection algorithm and the advection of particles that carry a scalar quantity
representing the location of each compositional field. All four of these
algorithms are implemented in the open source FEM code ASPECT
On Meshfree GFDM Solvers for the Incompressible Navier-Stokes Equations
Meshfree solution schemes for the incompressible Navier--Stokes equations are
usually based on algorithms commonly used in finite volume methods, such as
projection methods, SIMPLE and PISO algorithms. However, drawbacks of these
algorithms that are specific to meshfree methods have often been overlooked. In
this paper, we study the drawbacks of conventionally used meshfree Generalized
Finite Difference Method~(GFDM) schemes for Lagrangian incompressible
Navier-Stokes equations, both operator splitting schemes and monolithic
schemes. The major drawback of most of these schemes is inaccurate local
approximations to the mass conservation condition. Further, we propose a new
modification of a commonly used monolithic scheme that overcomes these problems
and shows a better approximation for the velocity divergence condition. We then
perform a numerical comparison which shows the new monolithic scheme to be more
accurate than existing schemes
Discrete Imaging Models for Three-Dimensional Optoacoustic Tomography using Radially Symmetric Expansion Functions
Optoacoustic tomography (OAT), also known as photoacoustic tomography, is an
emerging computed biomedical imaging modality that exploits optical contrast
and ultrasonic detection principles. Iterative image reconstruction algorithms
that are based on discrete imaging models are actively being developed for OAT
due to their ability to improve image quality by incorporating accurate models
of the imaging physics, instrument response, and measurement noise. In this
work, we investigate the use of discrete imaging models based on Kaiser-Bessel
window functions for iterative image reconstruction in OAT. A closed-form
expression for the pressure produced by a Kaiser-Bessel function is calculated,
which facilitates accurate computation of the system matrix.
Computer-simulation and experimental studies are employed to demonstrate the
potential advantages of Kaiser-Bessel function-based iterative image
reconstruction in OAT
- …