21 January 2020 11:07:20 AM

TOMS917_TEST
  C version
  Test the TOMS917 library.

DRIVER:
  Demonstrate simple and extended Wright Omega evaluators.

  Calling:
    w = wrightomega(z);
  returns:
    w = omega(                       0                       +0 I)=       0.567143290409784                       +0 I

  Calling:
    wrightomega_ext ( z, &w, &e, &r, &condest );
  returns:
  w = omega(0.000000 +0.000000 I)=0.567143 +0.000000 I
  e = last update step: -4.92506e-05 +0 I
  r = penultimate residual: -7.71854e-05 +0 I
  condest = condition number estimate: 0 +0 I

  ultimate residual: -1.17257e-17 +0 I

DRIVER:
  Demonstrate simple and extended Wright Omega evaluators.

  Calling:
    w = wrightomega(z);
  returns:
    w = omega(                       1                       +0 I)=                       1                       +0 I

  Calling:
    wrightomega_ext ( z, &w, &e, &r, &condest );
  returns:
  w = omega(1.000000 +0.000000 I)=1.000000 +0.000000 I
  e = last update step: 0 +0 I
  r = penultimate residual: 0 +0 I
  condest = condition number estimate: 0.5 +0 I

  ultimate residual: -0 +0 I

DRIVER:
  Demonstrate simple and extended Wright Omega evaluators.

  Calling:
    w = wrightomega(z);
  returns:
    w = omega(       3.718281828459045                       +0 I)=       2.718281828459045                       +0 I

  Calling:
    wrightomega_ext ( z, &w, &e, &r, &condest );
  returns:
  w = omega(3.718282 +0.000000 I)=2.718282 +0.000000 I
  e = last update step: -7.24546e-13 +0 I
  r = penultimate residual: -2.69407e-12 +0 I
  condest = condition number estimate: 1 +0 I

  ultimate residual: -1.58832e-49 +0 I

DRIVER:
  Demonstrate simple and extended Wright Omega evaluators.

  Calling:
    w = wrightomega(z);
  returns:
    w = omega(                      -1       +3.141592653589793 I)=                      -1                       +0 I

  Calling:
    wrightomega_ext ( z, &w, &e, &r, &condest );
  returns:
  w = omega(-1.000000 +3.141593 I)=-1.000000 +0.000000 I
  e = last update step: 0 +0 I
  r = penultimate residual: 0 +0 I
  condest = condition number estimate: 1 +0 I

  ultimate residual: -nan -nan I

DRIVER:
  Demonstrate simple and extended Wright Omega evaluators.

  Calling:
    w = wrightomega(z);
  returns:
    w = omega(                      -1       -3.141592653589793 I)=                      -1                       +0 I

  Calling:
    wrightomega_ext ( z, &w, &e, &r, &condest );
  returns:
  w = omega(-1.000000 -3.141593 I)=-1.000000 +0.000000 I
  e = last update step: 0 +0 I
  r = penultimate residual: 0 +0 I
  condest = condition number estimate: 1 +0 I

  ultimate residual: -nan -nan I

DRIVER:
  Demonstrate simple and extended Wright Omega evaluators.

  Calling:
    w = wrightomega(z);
  returns:
    w = omega(      -1.306852819440055       +3.141592653589793 I)=                      -2                       -0 I

  Calling:
    wrightomega_ext ( z, &w, &e, &r, &condest );
  returns:
  w = omega(-1.306853 +3.141593 I)=-2.000000 -0.000000 I
  e = last update step: 4.81044e-06 +0 I
  r = penultimate residual: -4.8104e-06 +0 I
  condest = condition number estimate: 1.30685 -0 I

  ultimate residual: 1.23156e-20 -9.04933e-36 I

DRIVER:
  Demonstrate simple and extended Wright Omega evaluators.

  Calling:
    w = wrightomega(z);
  returns:
    w = omega(      -1.306852819440055       -3.141592653589793 I)=       -0.40637573995996                       -0 I

  Calling:
    wrightomega_ext ( z, &w, &e, &r, &condest );
  returns:
  w = omega(-1.306853 -3.141593 I)=-0.406376 -0.000000 I
  e = last update step: 4.65417e-05 +0 I
  r = penultimate residual: 2.76281e-05 +0 I
  condest = condition number estimate: -2.20148 -0 I

  ultimate residual: 3.43693e-17 +0 I

DRIVER:
  Demonstrate simple and extended Wright Omega evaluators.

  Calling:
    w = wrightomega(z);
  returns:
    w = omega(                       0       +2.570796326794897 I)=   4.618701254790458e-17                       +1 I

  Calling:
    wrightomega_ext ( z, &w, &e, &r, &condest );
  returns:
  w = omega(0.000000 +2.570796 I)=0.000000 +1.000000 I
  e = last update step: -2.04459e-15 -2.43787e-14 I
  r = penultimate residual: 2.23341e-14 -2.64233e-14 I
  condest = condition number estimate: 1.2854 +1.2854 I

  ultimate residual: -1.17683e-54 +9.78138e-55 I

DRIVER:
  Demonstrate simple and extended Wright Omega evaluators.

  Calling:
    w = wrightomega(z);
  returns:
    w = omega(                       0       +3.141592653589793 I)=      -0.318131505204764       +1.337235701430689 I

  Calling:
    wrightomega_ext ( z, &w, &e, &r, &condest );
  returns:
  w = omega(0.000000 +3.141593 I)=-0.318132 +1.337236 I
  e = last update step: -7.62322e-16 +1.92451e-16 I
  r = penultimate residual: -7.77156e-16 -8.88178e-16 I
  condest = condition number estimate: 1.86453 +0.95074 I

  ultimate residual: 4.06691e-61 +2.08648e-60 I

DRIVER:
  Demonstrate simple and extended Wright Omega evaluators.

  Calling:
    w = wrightomega(z);
  returns:
    w = omega(                       1                       +1 I)=      0.9372082083733697      +0.5054213160131512 I

  Calling:
    wrightomega_ext ( z, &w, &e, &r, &condest );
  returns:
  w = omega(1.000000 +1.000000 I)=0.937208 +0.505421 I
  e = last update step: 2.85027e-05 -1.38839e-05 I
  r = penultimate residual: 6.22316e-05 -1.24894e-05 I
  condest = condition number estimate: 0.609404 +0.357212 I

  ultimate residual: 8.54308e-19 +1.70422e-18 I

TEST_BOUNDARY:
  Test wrightomega_ext() near approximation region boundaries.
  Store results in a file for comparison with benchmark data.

TEST_BOUNDARY:
  Results saved in file 'results.txt'

TOMS917_TEST
  Normal end of execution.

21 January 2020 11:07:20 AM