|
D.8.2.5 applyMatrix
Procedure from library interval.lib (see interval_lib).
- Usage:
- A * b, A ivmat, b box
- Return:
- A*b
Example:
| LIB "interval.lib";
ring R = 0,(x,y,z),lp;
ideal I = xyz3+z2y2+x,x4+y3+2z+3,xyz+1/2;
interval J = bounds2(1/2, 3/2);
box B = list(J,J,J);
ivmat A = evalJacobianAtBox(I, B); A;
==> [17/16, 97/16],[5/16, 189/16],[7/16, 351/16]
==> [1/2, 27/2],[3/4, 27/4],[2, 2]
==> [1/4, 9/4],[1/4, 9/4],[1/4, 9/4]
==>
A*B;
==> [29/32, 1911/32],[13/8, 267/8],[3/8, 81/8]
==>
unitMatrix2(3) * B;
==> [1/2, 3/2],[1/2, 3/2],[1/2, 3/2]
==>
diagMatrix(3, bounds2(0, 1)) * B;
==> [0, 3/2],[0, 3/2],[0, 3/2]
==>
|
|