11,838 research outputs found
Analysis of Iterative Methods for the Steady and Unsteady Stokes Problem: Application to Spectral Element Discretizations
A new and detailed analysis of the basic Uzawa algorithm for decoupling of the pressure and the velocity in the steady and unsteady Stokes operator is presented. The paper focuses on the following new aspects: explicit construction of the Uzawa pressure-operator spectrum for a semiperiodic model problem; general relationship of the convergence rate of the Uzawa procedure to classical inf-sup discretization analysis; and application of the method to high-order variational discretization
Cluster Variation Method in Statistical Physics and Probabilistic Graphical Models
The cluster variation method (CVM) is a hierarchy of approximate variational
techniques for discrete (Ising--like) models in equilibrium statistical
mechanics, improving on the mean--field approximation and the Bethe--Peierls
approximation, which can be regarded as the lowest level of the CVM. In recent
years it has been applied both in statistical physics and to inference and
optimization problems formulated in terms of probabilistic graphical models.
The foundations of the CVM are briefly reviewed, and the relations with
similar techniques are discussed. The main properties of the method are
considered, with emphasis on its exactness for particular models and on its
asymptotic properties.
The problem of the minimization of the variational free energy, which arises
in the CVM, is also addressed, and recent results about both provably
convergent and message-passing algorithms are discussed.Comment: 36 pages, 17 figure
High-order regularized regression in Electrical Impedance Tomography
We present a novel approach for the inverse problem in electrical impedance
tomography based on regularized quadratic regression. Our contribution
introduces a new formulation for the forward model in the form of a nonlinear
integral transform, that maps changes in the electrical properties of a domain
to their respective variations in boundary data. Using perturbation theory the
transform is approximated to yield a high-order misfit unction which is then
used to derive a regularized inverse problem. In particular, we consider the
nonlinear problem to second-order accuracy, hence our approximation method
improves upon the local linearization of the forward mapping. The inverse
problem is approached using Newton's iterative algorithm and results from
simulated experiments are presented. With a moderate increase in computational
complexity, the method yields superior results compared to those of regularized
linear regression and can be implemented to address the nonlinear inverse
problem
- …