Interval contractor

From HandWiki
Short description: Mathematical construct

In mathematics, an interval contractor (or contractor for short)[1] associated to a set [math]\displaystyle{ X }[/math] is an operator [math]\displaystyle{ C }[/math] which associates to a hyperrectangle [math]\displaystyle{ [x] }[/math] in [math]\displaystyle{ \bold{R}^n }[/math] another box [math]\displaystyle{ C([x]) }[/math] of [math]\displaystyle{ \bold{R}^n }[/math] such that the two following properties are always satisfied:

  • [math]\displaystyle{ C([x]) \subset [x] }[/math] (contractance property)
  • [math]\displaystyle{ C([x]) \cap X = [x] \cap X }[/math] (completeness property)

A contractor associated to a constraint (such as an equation or an inequality) is a contractor associated to the set [math]\displaystyle{ X }[/math] of all [math]\displaystyle{ x }[/math] which satisfy the constraint. Contractors make it possible to improve the efficiency of branch-and-bound algorithms classically used in interval analysis.

Properties of contractors

A contractor C is monotonic if we have [math]\displaystyle{ [x] \subset [y] \Rightarrow C([x])\subset C([y]) }[/math].

It is minimal if for all boxes [x], we have [math]\displaystyle{ C([x]) = [[x]\cap X] }[/math], where [A] is the interval hull of the set A, i.e., the smallest box enclosing A.

The contractor C is thin if for all points x, [math]\displaystyle{ C(\{x\}) = \{x\}\cap X }[/math] where {x} denotes the degenerated box enclosing x as a single point.

The contractor C is idempotent if for all boxes [x], we have [math]\displaystyle{ C \circ C([x]) = C([x]). }[/math]

The contractor C is convergent if for all sequences [x](k) of boxes containing x, we have [math]\displaystyle{ [x](k) \rightarrow x \implies C([x](k))\rightarrow \{x\} \cap X. }[/math]

Illustration

Figure 1 represents the set X painted grey and some boxes, some of them degenerated (i.e., they correspond to singletons). Figure 2 represents these boxes after contraction. Note that no point of X has been removed by the contractor. The contractor is minimal for the cyan box but is pessimistic for the green one. All degenerated blue boxes are contracted to the empty box. The magenta box and the red box cannot be contracted.

Figure 1: Boxes before contraction
Figure 2: Boxes after contraction

Contractor algebra

Some operations can be performed on contractors to build more complex contractors. [2] The intersection, the union, the composition and the repetition are defined as follows.

[math]\displaystyle{ (C_1 \cap C_2)([x]) = C_1([x]) \cap C_2([x]) }[/math]
[math]\displaystyle{ (C_1 \cup C_2) ([x]) = [C_1([x]) \cup C_2([x])] }[/math]
[math]\displaystyle{ (C_1 \circ C_2) ([x]) = C_1( C_2([x]) ) }[/math]
[math]\displaystyle{ C^\infty([x]) = C \circ C \circ C \circ \cdots \circ C([x]) }[/math]

Building contractors

There exist different ways to build contractors associated to equations and inequalities, say, f(x) in [y]. Most of them are based on interval arithmetic. One of the most efficient and most simple is the forward/backward contractor (also called as HC4-revise).[3][4]

The principle is to evaluate f(x) using interval arithmetic (this is the forward step). The resulting interval is intersected with [y]. A backward evaluation of f(x) is then performed in order to contract the intervals for the xi (this is the backward step). We now illustrate the principle on a simple example.

Consider the constraint [math]\displaystyle{ (x_1+x_2)\cdot x_3 \in [1,2]. }[/math] We can evaluate the function f(x) by introducing the two intermediate variables a and b, as follows

[math]\displaystyle{ a = x_1+x_2 }[/math]
[math]\displaystyle{ b =a \cdot x_3 }[/math]

The two previous constraints are called forward constraints. We get the backward constraints by taking each forward constraint in the reverse order and isolating each variable on the right hand side. We get

[math]\displaystyle{ x_3 = \frac{b}{a} }[/math]
[math]\displaystyle{ a = \frac{b}{x_3} }[/math]
[math]\displaystyle{ x_1 = a-x_2 }[/math]
[math]\displaystyle{ x_2 = a-x_1 }[/math]

The resulting forward/backward contractor [math]\displaystyle{ C([x_1],[x_2],[x_3]) }[/math] is obtained by evaluating the forward and the backward constraints using interval analysis.

[math]\displaystyle{ [a] = [x_1]+[x_2] }[/math]
[math]\displaystyle{ [b] = [a] \cdot [x_3] }[/math]
[math]\displaystyle{ [b] = [b] \cap [1,2] }[/math]
[math]\displaystyle{ [x_3] = [x_3]\ \cap \ \frac{[b]}{[a]} }[/math]
[math]\displaystyle{ [a] = [a] \cap \frac{[b]}{[x_3]} }[/math]
[math]\displaystyle{ [x_1] = [x_1]\ \cap \ \ [a]-[x_2] }[/math]
[math]\displaystyle{ [x_2] = [x_2]\ \cap \ \ [a]-[x_1] }[/math]

References

  1. Jaulin, Luc; Kieffer, Michel; Didrit, Olivier; Walter, Eric (2001). Applied Interval Analysis. Berlin: Springer. ISBN 1-85233-219-0. 
  2. Chabert, G.; Jaulin, L. (2009). "Contractor programming". Artificial Intelligence 173 (11): 1079–1100. doi:10.1016/j.artint.2009.03.002. http://www.ensta-bretagne.fr/jaulin/paper_chabert_quimper.pdf. 
  3. Messine, F. (1997). Méthode d'optimisation globale basée sur l'analyse d'intervalles pour la résolution de problèmes avec contraintes. Thèse de doctorat, Institut National Polytechnique de Toulouse. http://messine.perso.enseeiht.fr/PUBLIS/these.ps. 
  4. Benhamou, F.; Goualard, F.; Granvilliers, L.; Puget, J.F. (1999). Revising hull and box consistency. In Proceedings of the 1999 international conference on Logic programming. https://frederic.goualard.net/publications/revising-hull-and-box-consistency_benhamou_et-al_iclp1999.pdf.