Interval contractor explained

In mathematics, an interval contractor (or contractor for short)[1] associated to a set

X

is an operator

C

which associates to a hyperrectangle

[x]

in

\bold{R}n

another box

C([x])

of

\bold{R}n

such that the two following properties are always satisfied:

C([x])\subset[x]

(contractance property)

C([x])\capX=[x]\capX

(completeness property)

A contractor associated to a constraint (such as an equation or an inequality) is a contractor associated to the set

X

of all

x

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

[x]\subset[y]C([x])\subsetC([y])

.

It is minimal if for all boxes [''x''], we have

C([x])=[[x]\capX]

,where [''A''] is the interval hull of the set A, i.e., the smallestbox enclosing A.

The contractor C is thin if for all points x,

C(\{x\})=\{x\}\capX

where denotes the degenerated box enclosing x as a single point.

The contractor C is idempotent if for all boxes [''x''], we have

C\circC([x])=C([x]).

The contractor C is convergent if for all sequences [''x''](k) of boxes containing x, we have

[x](k)x\impliesC([x](k))\{x\}\capX.

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 contractoris minimal for the cyan box but is pessimistic for the green one. All degenerated blue boxes are contracted tothe empty box. The magenta box and the red box cannot be contracted.

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.

(C1\capC2)([x])=C1([x])\capC2([x])

(C1\cupC2)([x])=[C1([x])\cupC2([x])]

(C1\circC2)([x])=C1(C2([x]))

Cinfty([x])=C\circC\circC\circ\circC([x])

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

(x1+x2)x3\in[1,2].

We can evaluate the function f(x) by introducing the two intermediatevariables a and b, as follows

a=x1+x2

b=ax3

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

x3=

b
a

a=

b
x3

x1=a-x2

x2=a-x1

The resulting forward/backward contractor

C([x1],[x2],[x3])

is obtained by evaluating the forward and the backward constraints using interval analysis.

[a]=[x1]+[x2]

[b]=[a][x3]

[b]=[b]\cap[1,2]

[x3]=[x3]\cap

[b]
[a]

[a]=[a]\cap

[b]
[x3]

[x1]=[x1]\cap  [a]-[x2]

[x2]=[x2]\cap  [a]-[x1]

Notes and References

  1. Book: Jaulin. Luc. Kieffer. Michel. Didrit. Olivier. Walter. Eric. Applied Interval Analysis. 2001. Springer. Berlin. 1-85233-219-0.
  2. Chabert. G.. Jaulin. L.. Contractor programming. Artificial Intelligence. 2009. 173. 11. 10.1016/j.artint.2009.03.002. 1079–1100. free.
  3. Book: Messine. F.. Méthode d'optimisation globale basée sur l'analyse d'intervalles pour la résolution de problèmes avec contraintes. 1997. Thèse de doctorat, Institut National Polytechnique de Toulouse.
  4. Book: Benhamou. F.. Goualard. F.. Granvilliers. L.. Puget. J.F.. Revising hull and box consistency. 1999. In Proceedings of the 1999 international conference on Logic programming.