Results 1  10
of
39
On the Complexity of Isolating Real Roots and Computing with Certainty the Topological Degree
, 2002
"... In this contribution the isolation of real roots and the computation of the topological degree... ..."
Abstract

Cited by 42 (20 self)
 Add to MetaCart
In this contribution the isolation of real roots and the computation of the topological degree...
Linear systems with large uncertainties, with applications to truss structures
, 2006
"... Linear systems whose coefficients have large uncertainties arise routinely in finite element calculations for structures with uncertain geometry, material properties, or loads. However, a true worst case analysis of the influence of such uncertainties was previously possible only for very small sys ..."
Abstract

Cited by 20 (6 self)
 Add to MetaCart
Linear systems whose coefficients have large uncertainties arise routinely in finite element calculations for structures with uncertain geometry, material properties, or loads. However, a true worst case analysis of the influence of such uncertainties was previously possible only for very small systems and uncertainties, or in special cases where the coefficients do not exhibit dependence. This paper presents a method for computing rigorous bounds on the solution of such systems, with a computable overestimation factor that is frequently quite small. The merits of the new approach are demonstrated by computing realistic bounds for some large, uncertain truss structures, some leading to linear systems with over 5000 variables and over 10000 interval parameters, with excellent bounds for up to about 10 % input uncertainty. Also discussed are some counterexamples for the performance of traditional approximate methods for worst case uncertainty analysis.
Decomposition of Arithmetic Expressions to Improve the Behavior of Interval Iteration for Nonlinear Systems
, 1991
"... Interval iteration can be used, in conjunction with other techniques, for rigorously bounding all solutions to a nonlinear system of equations within a given region, or for verifying approximate solutions. However, because of overestimation which occurs when the interval Jacobian matrix is accumul ..."
Abstract

Cited by 20 (9 self)
 Add to MetaCart
Interval iteration can be used, in conjunction with other techniques, for rigorously bounding all solutions to a nonlinear system of equations within a given region, or for verifying approximate solutions. However, because of overestimation which occurs when the interval Jacobian matrix is accumulated and applied, straightforward linearization of the original nonlinear system sometimes leads to nonconvergent iteration. In this paper, we examine interval iterations based on an expanded system obtained from the intermediate quantities in the original system. In this system, there is no overestimation in entries of the interval Jacobi matrix, and nonlinearities can be taken into account to obtain sharp bounds. We present an example in detail, algorithms, and detailed experimental results obtained from applying our algorithms to the example.
Empirical Evaluation Of Innovations In Interval Branch And Bound Algorithms For Nonlinear Systems
 SIAM J. Sci. Comput
, 1994
"... . Interval branch and bound algorithms for finding all roots use a combination of a computational existence / uniqueness procedure and a tesselation process (generalized bisection). Such algorithms identify, with mathematical rigor, a set of boxes that contains unique roots and a second set within w ..."
Abstract

Cited by 20 (10 self)
 Add to MetaCart
(Show Context)
. Interval branch and bound algorithms for finding all roots use a combination of a computational existence / uniqueness procedure and a tesselation process (generalized bisection). Such algorithms identify, with mathematical rigor, a set of boxes that contains unique roots and a second set within which all remaining roots must lie. Though each root is contained in a box in one of the sets, the second set may have several boxes in clusters near a single root. Thus, the output is of higher quality if there are relatively more boxes in the first set. In contrast to previously implemented similar techniques, a box expansion technique in this paper, based on using an approximate root finder, fflinflation and exact set complementation, decreases the size of the second set, increases the size of the first set, and never loses roots. In addition to the expansion technique, use of secondorder extensions to eliminate small boxes that do not contain roots, and interval slopes versus interval d...
Rigorous sensitivity analysis for systems of linear and nonlinear equations
 Mathematics of Computation 54
, 1990
"... Methods are presented for performing a rigorous sensitivity analysis of numerical problems with independent, noncorrelated data for general systems of linear and nonlinear equations. The methods may serve for the following two purposes. First, to bound the dependency of the solution on changes in t ..."
Abstract

Cited by 17 (2 self)
 Add to MetaCart
(Show Context)
Methods are presented for performing a rigorous sensitivity analysis of numerical problems with independent, noncorrelated data for general systems of linear and nonlinear equations. The methods may serve for the following two purposes. First, to bound the dependency of the solution on changes in the input data. In contrast to condition numbers a componentwise sensitivity analysis of the solution vector is performed. Second, to estimate the true solution set for problems the input data of which are afflicted with tolerances. The methods presented are very effective with the additional property that, due to an automatic error control mechanism, every computed result is guaranteed to be correct. Examples are given for linear systems demonstrating that the computed bounds are in general very sharp. Interesting comparisons to traditional condition numbers are given.
Automatic Computation of a Linear Interval Enclosure
 Reliable Computing
, 2001
"... Abstract. Recently, an alternative interval approximation F ( X) for enclosing a factorable function f(x) in a given box X has been suggested. The enclosure is in the form of an affine interval function n i=1 F ( X) = a X + B where only the additive term B is an interval, the coefficients ai being r ..."
Abstract

Cited by 14 (6 self)
 Add to MetaCart
(Show Context)
Abstract. Recently, an alternative interval approximation F ( X) for enclosing a factorable function f(x) in a given box X has been suggested. The enclosure is in the form of an affine interval function n i=1 F ( X) = a X + B where only the additive term B is an interval, the coefficients ai being real i i numbers. The approximation is applicable to continuously differentiable, continuous and even discontinuous functions. In this paper, a new algorithm for determining the coefficients ai and the interval B of F(X) is proposed. It is based on the introduction of a specific generalized representation of intervals which permits the computation of the enclosure considered to be fully automated. 1.
Exclusion Regions for Systems of Equations
 SIAM J. NUM. ANALYSIS
, 2003
"... Branch and bound methods for finding all zeros of a nonlinear system of equations in a box frequently have the difficulty that subboxes containing no solution cannot be easily eliminated if there is a nearby zero outside the box. This has the effect that near each zero, many small boxes are created ..."
Abstract

Cited by 11 (7 self)
 Add to MetaCart
(Show Context)
Branch and bound methods for finding all zeros of a nonlinear system of equations in a box frequently have the difficulty that subboxes containing no solution cannot be easily eliminated if there is a nearby zero outside the box. This has the effect that near each zero, many small boxes are created by repeated splitting, whose processing may dominate the total work spent on the global search. This paper discusses the reasons for the occurrence of this socalled cluster effect, and how to reduce the cluster effect by defining exclusion regions around each zero found, that are guaranteed to contain no other zero and hence can safely be discarded. Such exclusion regions are traditionally constructed using uniqueness tests based on the Krawczyk operator or the Kantorovich theorem. These results are reviewed; moreover, refinements are proved that significantly enlarge the size of the exclusion region. Existence and uniqueness tests are also given.
Expansion and estimation of the range of nonlinear functions
 Math. Comp
, 1996
"... Abstract. Many verification algorithms use an expansion f(x) ∈ f(˜x)+S· (x−˜x), f: Rn → Rn for x ∈ X, where the set of matrices S is usually computed as a gradient or by means of slopes. In the following, an expansion scheme is described which frequently yields sharper inclusions for S. This allows ..."
Abstract

Cited by 11 (1 self)
 Add to MetaCart
(Show Context)
Abstract. Many verification algorithms use an expansion f(x) ∈ f(˜x)+S· (x−˜x), f: Rn → Rn for x ∈ X, where the set of matrices S is usually computed as a gradient or by means of slopes. In the following, an expansion scheme is described which frequently yields sharper inclusions for S. This allows also to compute sharper inclusions for the range of f over a domain. Roughly speaking, f has to be given by means of a computer program. The process of expanding f can then be fully automatized. The function f need not be differentiable. For locally convex or concave functions special improvements are described. Moreover, in contrast to other methods, ˜x ∩ X may be empty without implying large overestimations for S. This may be advantageous in practical applications. 0. Notation We denote by IR the set of real intervals X ∈ IR ⇒ X = [inf(X), sup(X)] = { x ∈ R  inf(X) ≤ x ≤ sup(X)}. By PT we denote the power set over a given set T, and we use the canonical embedding IR ⊆ PR. Thesetofndimensional interval vectors is denoted by IRn, i.e., X ∈ IR n ⇒ X = { (xi) ∈ R n  xi ∈ Xi} with Xi ∈ IR, 1 ≤ i ≤ n. Interval vectors are compact. Interval operations and power set operations are defined in the usual way. Details can be found in standard books on interval analysis, among others [10, 2, 11]. If not explicitly noted otherwise, all operations are power set operations. 1. Expansion of nonlinear functions A differentiable function f: D ⊆ Rn → R can be locally expanded by its gradient. For ˜x ∈ D, X ⊆ D, and[g]∈IRn with ∇f(˜x ∪ X) ⊆ [g] there holds ∀ x ∈ X ∃ g ∈ [g] : f(x)−f(˜x)=g T (1.0) ·(x−˜x). Here, ∪ denotes the convex hull, and ∇f(˜x ∪ X) denotes the range of ∇f over ˜x ∪ X. The gradient, for real and for interval arguments, can be computed using automatic differentiation [4, 12]. This process is fully automatized. This approach has three disadvantages: (1) f needs to be differentiable,
Taylor Forms  Use and Limits
 Reliable Computing
, 2002
"... This review is a response to recent discussions on the reliable computing mailing list, and to continuing uncertainties about the properties and merits of Taylor forms, multivariate higher degree generalizations of centered forms. They were invented around 1980 by Lanford, documented in detail in 19 ..."
Abstract

Cited by 8 (2 self)
 Add to MetaCart
(Show Context)
This review is a response to recent discussions on the reliable computing mailing list, and to continuing uncertainties about the properties and merits of Taylor forms, multivariate higher degree generalizations of centered forms. They were invented around 1980 by Lanford, documented in detail in 1984 by Eckmann, Koch and Wittwer, and independently studied and popularized since 1996 by Berz, Makino and Hoefkens. A highlight is their application to the verified integration of asteroid dynamics in the solar system in 2001, although the details given are not sufficient to check the validity of their claims.