next | previous | forward | backward | up | top | index | toc | Macaulay2 website
NumericalCertification :: krawczykMethod

krawczykMethod -- certify the interval box for square polynomial system

Synopsis

Description

For given interval and polynomial system, this function computes the Krawczyk operator and check that the operator is contained in the input interval.

i1 : R = QQ[x1,x2,y1,y2];
i2 : f = polySystem {3*y1 + 2*y2 -1, 3*x1 + 2*x2 -7/2,x1^2 + y1^2 -1, x2^2 + y2^2 - 1};
i3 : (I1, I2, I3, I4) = (interval(.94,.96), interval(.31,.33), interval(-.31,-.29), interval(.94,.96));

Intervals for certification should be given in advance using other system solvers. After that we set the relationships between variables and intervals.

i4 : o = intervalOptionList {("x1" => "I1"), ("x2" => "I2"), ("y1" => "I3"), ("y2" => "I4")};

If the Krawczyk operator is contained in the input interval, then the function returns the result that the input interval (or the Krawczyk operator) contains a unique root of the system.

i5 : krawczykMethod(f,o)
given interval contains a unique solution

o5 = true

By default, these functions uses the inverse of the matrix obtained by evaluating the system at the midpoint of the input interval.

However, users can choose an invertible matrix by themselves. (See InvertibleMatrix.)

i6 : Y = matrix {{.095, .032, .476, -.100},{-.143, .452, -.714, .150},{.301,.101,-.160, -.317},{.048,-.152,.240,.476}};

                4          4
o6 : Matrix RR    <--- RR
              53         53
i7 : krawczykMethod(f,o,InvertibleMatrix => Y)

o7 = true

Ways to use krawczykMethod :

For the programmer

The object krawczykMethod is a method function with options.