Results 1  10
of
1,242
Iterative Methods For Total Variation Denoising
 SIAM J. SCI. COMPUT
"... Total Variation (TV) methods are very effective for recovering "blocky", possibly discontinuous, images from noisy data. A fixed point algorithm for minimizing a TVpenalized least squares functional is presented and compared with existing minimization schemes. A variant of the cellcentered finite ..."
Abstract

Cited by 230 (7 self)
 Add to MetaCart
Total Variation (TV) methods are very effective for recovering "blocky", possibly discontinuous, images from noisy data. A fixed point algorithm for minimizing a TVpenalized least squares functional is presented and compared with existing minimization schemes. A variant of the cellcentered finite difference multigrid method of Ewing and Shen is implemented for solving the (large, sparse) linear subproblems. Numerical results are presented for one and twodimensional examples; in particular, the algorithm is applied to actual data obtained from confocal microscopy.
Unified analysis of discontinuous Galerkin methods for elliptic problems
 SIAM J. Numer. Anal
, 2001
"... Abstract. We provide a framework for the analysis of a large class of discontinuous methods for secondorder elliptic problems. It allows for the understanding and comparison of most of the discontinuous Galerkin methods that have been proposed over the past three decades for the numerical treatment ..."
Abstract

Cited by 220 (19 self)
 Add to MetaCart
Abstract. We provide a framework for the analysis of a large class of discontinuous methods for secondorder elliptic problems. It allows for the understanding and comparison of most of the discontinuous Galerkin methods that have been proposed over the past three decades for the numerical treatment of elliptic problems.
Finite Element Methods for Active Contour Models and Balloons for 2D and 3D Images
 IEEE Transactions on Pattern Analysis and Machine Intelligence
, 1991
"... The use of energyminimizing curves, known as "snakes" to extract features of interest in images has been introduced by Kass, Witkin and Terzopoulos [23]. A balloon model was introduced in [12] as a way to generalize and solve some of the problems encountered with the original method. We present a 3 ..."
Abstract

Cited by 161 (22 self)
 Add to MetaCart
The use of energyminimizing curves, known as "snakes" to extract features of interest in images has been introduced by Kass, Witkin and Terzopoulos [23]. A balloon model was introduced in [12] as a way to generalize and solve some of the problems encountered with the original method. We present a 3D generalization of the balloon model as a 3D deformable surface, which evolves in 3D images. It is deformed under the action of internal and external forces attracting the surface toward detected edgels by means of an attraction potential. We also show properties of energyminimizing surfaces concerning their relationship with 3D edge points. To solve the minimization problem for a surface, two simplified approaches are shown first, defining a 3D surface as a series of 2D planar curves. Then, after comparing Finite Element Method and Finite Difference Method in the 2D problem, we solve the 3D model using the Finite Element Method yielding greater stability and faster convergence. We have a...
The Local Discontinuous Galerkin Method For TimeDependent ConvectionDiffusion Systems
"... In this paper, we study the Local Discontinuous Galerkin methods for nonlinear, timedependent convectiondiffusion systems. These methods are an extension of the RungeKutta Discontinuous Galerkin methods for purely hyperbolic systems to convectiondiffusion systems and share with those methods the ..."
Abstract

Cited by 119 (18 self)
 Add to MetaCart
In this paper, we study the Local Discontinuous Galerkin methods for nonlinear, timedependent convectiondiffusion systems. These methods are an extension of the RungeKutta Discontinuous Galerkin methods for purely hyperbolic systems to convectiondiffusion systems and share with those methods their high parallelizability, their highorder formal accuracy, and their easy handling of complicated geometries, for convection dominated problems. It is proven that for scalar equations, the Local Discontinuous Galerkin methods are L2stable in the nonlinear case. Moreover, in the linear case, it is shown that if polynomials of degree k are used, the methods are kth order accurate for general triangulations; although this order of convergence is suboptimal, it is sharp for the LDG methods. Preliminary numerical examples displaying the performance of the method are shown.
ANALYSIS OF MULTISCALE METHODS
, 2004
"... The heterogeneous multiscale method gives a general framework for the analysis of multiscale methods. In this paper, we demonstrate this by applying this framework to two canonical problems: The elliptic problem with multiscale coefficients and the quasicontinuum method. ..."
Abstract

Cited by 118 (13 self)
 Add to MetaCart
The heterogeneous multiscale method gives a general framework for the analysis of multiscale methods. In this paper, we demonstrate this by applying this framework to two canonical problems: The elliptic problem with multiscale coefficients and the quasicontinuum method.
Schwarz Analysis Of Iterative Substructuring Algorithms For Elliptic Problems In Three Dimensions
 SIAM J. Numer. Anal
, 1993
"... . Domain decomposition methods provide powerful preconditioners for the iterative solution of the large systems of algebraic equations that arise in finite element or finite difference approximations of partial differential equations. The preconditioners are constructed from exact or approximate sol ..."
Abstract

Cited by 110 (26 self)
 Add to MetaCart
. Domain decomposition methods provide powerful preconditioners for the iterative solution of the large systems of algebraic equations that arise in finite element or finite difference approximations of partial differential equations. The preconditioners are constructed from exact or approximate solvers for the same partial differential equation restricted to a set of subregions into which the given region has been divided. In addition, the preconditioner is often augmented by a coarse, second level approximation, that provides additional, global exchange of information, and which can enhance the rate of convergence considerably. The iterative substructuring methods, based on decompositions of the region into nonoverlapping subregions, form one of the main families of such algorithms. Many domain decomposition algorithms can conveniently be described and analyzed as Schwarz methods. These algorithms are fully defined in terms of a set of subspaces and auxiliary bilinear forms. A gener...
Domain Decomposition Algorithms With Small Overlap
, 1994
"... Numerical experiments have shown that twolevel Schwarz methods often perform very well even if the overlap between neighboring subregions is quite small. This is true to an even greater extent for a related algorithm, due to Barry Smith, where a Schwarz algorithm is applied to the reduced linear ..."
Abstract

Cited by 82 (11 self)
 Add to MetaCart
Numerical experiments have shown that twolevel Schwarz methods often perform very well even if the overlap between neighboring subregions is quite small. This is true to an even greater extent for a related algorithm, due to Barry Smith, where a Schwarz algorithm is applied to the reduced linear system of equations that remains after that the variables interior to the subregions have been eliminated. In this paper, a supporting theory is developed.
Sliver Exudation
 ANNUAL SYMPOSIUM ON COMPUTATIONAL GEOMETRY
, 1999
"... A sliver is a tetrahedron whose four vertices lie close to a plane and whose orthogonal projection to that plane is a convex quadrilateral with no short edge. Slivers are notoriously common in 3dimensional Delaunay triangulations even for wellspaced point sets. We show that if the Delaunay triangu ..."
Abstract

Cited by 81 (11 self)
 Add to MetaCart
A sliver is a tetrahedron whose four vertices lie close to a plane and whose orthogonal projection to that plane is a convex quadrilateral with no short edge. Slivers are notoriously common in 3dimensional Delaunay triangulations even for wellspaced point sets. We show that if the Delaunay triangulation has the ratio property introduced in [15] then there is an assignment of weights so the weighted Delaunay triangulation contains no slivers. We also give an algorithm to compute such a weight assignment.
GALERKIN FINITE ELEMENT APPROXIMATIONS OF STOCHASTIC ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS
, 2004
"... We describe and analyze two numerical methods for a linear elliptic problem with stochastic coefficients and homogeneous Dirichlet boundary conditions. Here the aim of the computations is to approximate statistical moments of the solution, and, in particular, we give a priori error estimates for the ..."
Abstract

Cited by 77 (5 self)
 Add to MetaCart
We describe and analyze two numerical methods for a linear elliptic problem with stochastic coefficients and homogeneous Dirichlet boundary conditions. Here the aim of the computations is to approximate statistical moments of the solution, and, in particular, we give a priori error estimates for the computation of the expected value of the solution. The first method generates independent identically distributed approximations of the solution by sampling the coefficients of the equation and using a standard Galerkin finite element variational formulation. The Monte Carlo method then uses these approximations to compute corresponding sample averages. The second method is based on a finite dimensional approximation of the stochastic coefficients, turning the original stochastic problem into a deterministic parametric elliptic problem. A Galerkin finite element method, of either the h or pversion, then approximates the corresponding deterministic solution, yielding approximations of the desired statistics. We present a priori error estimates and include a comparison of the computational work required by each numerical approximation to achieve a given accuracy. This comparison suggests intuitive conditions for an optimal selection of the numerical approximation.
Subdivision Surfaces: A New Paradigm For ThinShell FiniteElement Analysis
 INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING
, 2000
"... We develop a new paradigm for thinshell finiteelement analysis based on the use of subdivision surfaces for: i) describing the geometry of the shell in its undeformed configuration, and ii) generating smooth interpolated displacement fields possessing bounded energy within the strict framework ..."
Abstract

Cited by 74 (19 self)
 Add to MetaCart
We develop a new paradigm for thinshell finiteelement analysis based on the use of subdivision surfaces for: i) describing the geometry of the shell in its undeformed configuration, and ii) generating smooth interpolated displacement fields possessing bounded energy within the strict framework of the KirchhoffLove theory of thin shells. The particular subdivision strategy adopted here is Loop's scheme, with extensions such as required to account for creases and displacement boundary conditions. The displacement fields obtained by subdivision are H 2 and, consequently, have a finite KirchhoffLove energy. The resulting finite elements contain three nodes and element integrals are computed by a onepoint quadrature. The displacement field of the shell is interpolated from nodal displacements only. In particular, no nodal rotations are used in the interpolation. The interpolation scheme induced by subdivision is nonlocal, i. e., the displacement field over one element depend on the nodal displacements of the element nodes and all nodes of immediately neighboring elements. However, the use of subdivision surfaces ensures that all the local displacement fields thus constructed combine conformingly to define one single limit surface.