Loading in 5 sec....

Testing COSY’s Taylor Model Arithmetic PowerPoint Presentation

Testing COSY’s Taylor Model Arithmetic

- By
**zea** - Follow User

- 123 Views
- Uploaded on

Download Presentation
## PowerPoint Slideshow about 'Testing COSY’s Taylor Model Arithmetic' - zea

**An Image/Link below is provided (as is) to download presentation**

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript

### Testing COSY’s Taylor Model Arithmetic

### What is “correct”?

### Goal

### Verification Process

### Verification process 1

### Verification process 2

### Testing Scope

### Error Report

### Error 1 (cos):

### Error 2 (cos):

### Error 3 (cos):

### Error 4 (sin):

### Error 5 (sin):

### Error 6 (sin):

### Error 7 (asin):

### Error 8 (asin):

### Error 9 (asin):

### Error 10 (asin):

### Error 11 (asin):

Jun Yu and George F. CorlissMarquette University, Dec. 16, 2002

Supported by Martin Berz, MSU

Domain: X

Expression: f

COSY’s Taylor model expression of f: TM_F

COSY’s interval expression of f: INV_F

For all x X, f (x) TM_F ([x]).

For all [x] X, INV_F ([x]) intersects TM_F([x]) .

Find test cases when x X, f (x) TM_F ([x])

Or find test cases when [x] X, INV_F ([x]) and TM_F ([x]) are disjoint.

Process 1: Use COSY’s interval arithmetic to verify COSY’s Taylor model arithmetic.

Process 2: Use Maple as a referee.

Evaluate the function f in Taylor model arithmetic over a domain.E.g., f = COS (3.14 + 1.57 * X), domain [-1, 1].

Construct the Taylor model expression of f (TM_EXPR) in COSYTM_EXPR := COS(-3.14 * TM_ONE + (1.57 * TM_ONE) * TM_INDEP)

Construct the interval expression of f (INL_EXPR) in COSY.IVL_EXPR := COS(INTV(-3.14, -3.14) + INTV(1.57, 1.57) * VAR1)

Pick a point , say A, inside the domain and convert it to a tight interval [A].

Evaluate the polynomial part of TM_EXPR on [A] and add the remainder bound.

Evaluate IVL_EXPR on [A].

Compare the results of 5 and 6. If the intervals are disjoint, there is an error.

Follow the first four steps of verification process 1.

Output the information about Taylor model expression in binary form.

In Maple, we code the function f.

Read the binary information of the Taylor model expression.

Evaluate f at the challenge point. If the point result of f is not in the interval result of Taylor model expression, we have a violation of containment.

Test cases are designed to evaluate the COSY operations of addition, subtraction, multiplication, sin, cos, tan, asin, acos, atan, sinh, cosh, tanh, log, exp, sqrt, sqr, isqrt, unary +, unary -.

We test Taylor models at

dimension 1: order 1, 7, 15, 17, and 20

dimension 2: order 1, 7, 15, 17, and 18

dimension 7: order 1, 2, 3, and 4

Using process 1 (all in COSY)

Observations: (For all the test cases we ran)

1. Errors are found in 3 functions: COS, SIN, and ASIN.

2. The COS and the SIN errors only occur at order 17 of dimension 1 or dimension 2.

3. The ACOS errors only occur at order 1 of dimension 1 and dimension 2, and order 1 or order 3 of dimension 7.

IMD: 0

Dimension: 1

Order: 17

TM_EXPR := COS(-3.141592653590006 * TM_ONE

+ (1.570796326794687 * TM_ONE) * TM_INDEP (1));

Point value: [-0.9999990000000004, -0.9999989999999995]

ERROR: Lower bd of expression > upper bd of Taylor model

Lower bound of expression: -0.1570796328686757E-05

Upper bound of Taylor model: -0.1570796816860323E-05

Difference: 0.4881735668679223E-12

IMD: 1

Dimension: 2

Order: 17

TM_EXPR := COS(-3.141592653590006 * TM_ONE

+ (1.570796326794687 * TM_ONE) * TM_INDEP (1));

Point value: [0.9999899999999996, 0.9999900000000005]

Point value: [0.9999899999999996, 0.9999900000000005]

ERROR: Lower bd of expression > upper bd of Taylor model

Lower bound of expression: -0.1570796369223327E-04

Upper bound of Taylor model: -0.1570796418440198E-04

Difference: 0.4921687142306775E-12

IMD: 1

Dimension: 2

Order: 17

TM_EXPR := COS(-3.141592653589796 * TM_ONE + (0.7853981633974484 * TM_ONE) * TM_INDEP(1) + (0.7853981633974484 * TM_ONE) * TM_INDEP(2));

Point value: [-0.9999900000000005, -0.9999899999999996]

Point value: [-0.9999900000000005, -0.9999899999999996]

ERROR: Lower bd of expression > upper bd of Taylor model

Lower bound of expression: -0.1570796327002688E-04

Upper bound of Taylor model: -0.1570796373588637E-04

Difference: 0.4658594880300668E-12

IMD: 0

Dimension: 1

Order: 17

TM_EXPR := SIN(-4.712388980384691 * TM_ONE

+ (1.570796326794690 * TM_ONE) * TM_INDEP(1));

Point value: [-0.9999990000000004, -0.9999989999999995]

ERROR: Upper bd of expression < lower bd of Taylor model

Upper bound of expression: 0.1570796538419324E-05

Lower bound of Taylor model: 0.1570797020611943E-05

Difference: -0.4821926187347235E-12

IMD: 1

Dimension: 2

Order: 17

TM_EXPR := SIN(-4.712388980384691 * TM_ONE

+ (1.570796326794690 * TM_ONE) * TM_INDEP(1));

Point value: [-0.9999900000000005, -0.9999899999999996]

Point value: [-0.5000000000000002, -0.4999999999999998]

ERROR: Upper bd of expression < lower bd of Taylor model

Upper bound of expression: 0.1570796347887127E-04

Lower bound of Taylor model: 0.1570796396102999E-04

Difference: -0.4821587160292514E-12

IMD: 1

Dimension: 2

Order: 17

TM_EXPR := SIN(4.712388980384691 * TM_ONE

+ (0.7853981633974483 * TM_ONE) * TM_INDEP (1)

+ (0.7853981633974483 * TM_ONE)* TM_INDEP (2));

Point value: [0.9999899999999996, 0.9999900000000005]

Point value: [0.9999899999999996, 0.9999900000000005]

ERROR: Lower bd of expression > upper bd of Taylor model

Lower bound of expression: -0.1570796327370206E-04

Upper bound of Taylor model: -0.1570796373932624E-04

Difference: 0.4656241856654513E-12

IMD: 0

Dimension: 1

Order: 1

TM_EXPR := ASIN((0.9999999999999983E-03 * TM_ONE)

* TM_INDEP(1));

Point value: [-0.9999990000000004, -0.9999989999999995]

ERROR: Upper bd of expression < lower bd of Taylor model

Upper bound of expression: -0.9999991666662369E-03

Lower bound of Taylor model: -0.9999990000000064E-03

Difference: -0.1666662305128269E-09

IMD: 0

Dimension: 1

Order: 1

TM_EXPR := ASIN(-0.4935 * TM_ONE + (0.3499999999999997E-02

* TM_ONE) * TM_INDEP(1));

Point value: [-0.9999990000000004, -0.9999989999999995]

ERROR: Upper bd of expression < lower bd of Taylor model

Upper bound of expression: -0.5201381202432389

Lower bound of Taylor model: -0.5201381145553354

Difference: -0.5687903481543799E-08

IMD: 0

Dimension: 2

Order: 1

TM_EXPR := ASIN((0.4999999999999991E-03 * TM_ONE)

* TM_INDEP(1) + (0.4999999999999991E-03

* TM_ONE) * TM_INDEP (2));}

Point value: [-0.9999900000000005, -0.9999899999999996]

Point value: [-0.9999900000000005, -0.9999899999999996]

ERROR: Upper bd of expression < lower bd of Taylor model

Upper bound of expression: -0.9999901666617370E-03

Lower bound of Taylor model: -0.9999900000000064E-03

Difference: -0.1666617306401302E-09

IMD: 1

Dimension: 7

Order: 1

TM_EXPR := ASIN((-0.4935 * TM_ONE) * TM_INDEP(1));

Point value: [0.4999999999999998, 0.5000000000000002]

Point value: [-0.5000000000000002, -0.4999999999999998]

Point value: [0.4999999999999998, 0.5000000000000002]

Point value: [0.4999999999999998, 0.5000000000000002]

Point value: [-0.5000000000000002, -0.4999999999999998]

Point value: [0.4999999999999998, 0.5000000000000002]

Point value: [0.4999999999999998, 0.5000000000000002]

ERROR: Upper bd of expression < lower bd of Taylor model

Upper bound of expression: -0.2493251165458725

Lower bound of Taylor model: -0.2467500000000038

Difference: -0.2575116545868672E-02

IMD: 1

Dimension: 7

Order: 3

TM_EXPR := ASIN((-0.4935 * TM_ONE) * TM_INDEP(1));

Point value: [0.4999999999999998, 0.5000000000000002]

Point value: [0.4999999999999998, 0.5000000000000002]

Point value: [0.4999999999999998, 0.5000000000000002]

Point value: [0.4999999999999998, 0.5000000000000002]

Point value: [0.4999999999999998, 0.5000000000000002]

Point value: [0.4999999999999998, 0.5000000000000002]

Point value: [0.4999999999999998, 0.5000000000000002]

ERROR: Upper bd of expression < lower bd of Taylor model

Upper bound of expression: -0.2493251165458725

Lower bound of Taylor model: -0.2492539187578167

Difference: -0.7119778805578236E-04

Using verification process 2 (Maple as referee)

Observation: (For all the test cases we ran)

Verification process 2 detects the exactly the same errors we have detected in verification process 1.

Download Presentation

Connecting to Server..