PCPRB5 PITCON test problem Two point boundary value problem. Number of equations is 21 Number of variables is 22 This program will be run three times. This is run number 1 Using the full Newton method. The jacobian is updated on every iteration. Step Type of point Lambda 0 Start point 0.00000 PITCON 7.0 University of Pittsburgh continuation code Modified: 12 November 1999 Linear algebra: LAPACK Precision: double 1 Corrected 0.00000 TRYSTP - DEBUG - about to return. 2 Continuation 0.264435E-01 TRYSTP - DEBUG - about to return. 3 Continuation 0.104600 TRYSTP - DEBUG - about to return. 4 Continuation 0.219655 TRYSTP - DEBUG - about to return. 5 Continuation 0.327553 TRYSTP - DEBUG - about to return. 6 Continuation 0.425474 TRYSTP - DEBUG - about to return. 7 Continuation 0.513404 TRYSTP - DEBUG - about to return. 8 Continuation 0.591468 TRYSTP - DEBUG - about to return. 9 Continuation 0.659916 TRYSTP - DEBUG - about to return. 10 Continuation 0.719109 TRYSTP - DEBUG - about to return. 11 Continuation 0.769501 TRYSTP - DEBUG - about to return. 12 Continuation 0.811615 13 Target point 0.800000 Complete value for target point 1 0.00000 0.602686E-01 0.118413 0.174305 0.227816 0.278815 0.327169 0.372748 0.415422 0.455064 0.491551 0.524765 0.554597 0.580942 0.603707 0.622811 0.638182 0.649762 0.657507 0.661387 0.661387 0.800000 TRYSTP - DEBUG - about to return. 14 Continuation 0.846018 TRYSTP - DEBUG - about to return. 15 Continuation 0.873302 TRYSTP - DEBUG - about to return. 16 Continuation 0.894067 TRYSTP - DEBUG - about to return. 17 Continuation 0.908904 TRYSTP - DEBUG - about to return. 18 Continuation 0.918385 TRYSTP - DEBUG - about to return. 19 Continuation 0.923057 TRYSTP - DEBUG - about to return. 20 Continuation 0.923436 21 Limit point 0.923283 TRYSTP - DEBUG - about to return. 22 Continuation 0.920003 TRYSTP - DEBUG - about to return. 23 Continuation 0.913206 TRYSTP - DEBUG - about to return. 24 Continuation 0.904857 TRYSTP - DEBUG - about to return. 25 Continuation 0.891133 TRYSTP - DEBUG - about to return. 26 Continuation 0.877626 TRYSTP - DEBUG - about to return. 27 Continuation 0.860121 TRYSTP - DEBUG - about to return. 28 Continuation 0.842782 TRYSTP - DEBUG - about to return. 29 Continuation 0.822585 TRYSTP - DEBUG - about to return. 30 Continuation 0.802513 TRYSTP - DEBUG - about to return. 31 Continuation 0.780519 32 Target point 0.800000 Complete value for target point 2 0.00000 0.152806 0.303281 0.451047 0.595672 0.736669 0.873486 1.00551 1.13207 1.25242 1.36578 1.47130 1.56811 1.65533 1.73208 1.79754 1.85093 1.89160 1.91902 1.93283 1.93283 0.800000 Jacobians 78 Factorizations 78 Solves 78 Functions 81 This is run number 2 Using the modified Newton method. The jacobian is held fixed while correcting a point. Step Type of point Lambda 0 Start point 0.00000 PITCON 7.0 University of Pittsburgh continuation code Modified: 12 November 1999 Linear algebra: LAPACK Precision: double 1 Corrected 0.00000 TRYSTP - DEBUG - about to return. 2 Continuation 0.264435E-01 TRYSTP - DEBUG - about to return. 3 Continuation 0.104600 TRYSTP - DEBUG - about to return. 4 Continuation 0.219655 TRYSTP - DEBUG - about to return. 5 Continuation 0.327553 TRYSTP - DEBUG - about to return. 6 Continuation 0.425474 TRYSTP - DEBUG - about to return. 7 Continuation 0.513404 TRYSTP - DEBUG - about to return. 8 Continuation 0.591468 TRYSTP - DEBUG - about to return. 9 Continuation 0.659916 TRYSTP - DEBUG - about to return. 10 Continuation 0.719109 TRYSTP - DEBUG - about to return. 11 Continuation 0.769501 TRYSTP - DEBUG - about to return. 12 Continuation 0.811615 13 Target point 0.800000 Complete value for target point 1 0.00000 0.602686E-01 0.118413 0.174305 0.227816 0.278815 0.327169 0.372748 0.415422 0.455064 0.491551 0.524765 0.554597 0.580942 0.603707 0.622811 0.638182 0.649762 0.657507 0.661387 0.661387 0.800000 TRYSTP - DEBUG - about to return. 14 Continuation 0.846018 TRYSTP - DEBUG - about to return. 15 Continuation 0.873302 TRYSTP - DEBUG - about to return. 16 Continuation 0.894067 TRYSTP - DEBUG - about to return. 17 Continuation 0.908904 TRYSTP - DEBUG - about to return. 18 Continuation 0.918385 TRYSTP - DEBUG - about to return. 19 Continuation 0.923057 TRYSTP - DEBUG - about to return. 20 Continuation 0.923436 21 Limit point 0.923283 TRYSTP - DEBUG - about to return. 22 Continuation 0.920003 TRYSTP - DEBUG - about to return. 23 Continuation 0.913206 TRYSTP - DEBUG - about to return. 24 Continuation 0.904857 TRYSTP - DEBUG - about to return. 25 Continuation 0.891133 TRYSTP - DEBUG - about to return. 26 Continuation 0.877626 TRYSTP - DEBUG - about to return. 27 Continuation 0.860121 TRYSTP - DEBUG - about to return. 28 Continuation 0.842782 TRYSTP - DEBUG - about to return. 29 Continuation 0.822585 TRYSTP - DEBUG - about to return. 30 Continuation 0.802513 TRYSTP - DEBUG - about to return. 31 Continuation 0.780519 32 Target point 0.800000 Complete value for target point 2 0.00000 0.152806 0.303281 0.451047 0.595672 0.736669 0.873486 1.00551 1.13207 1.25242 1.36578 1.47130 1.56811 1.65533 1.73208 1.79754 1.85093 1.89160 1.91902 1.93283 1.93283 0.800000 Jacobians 56 Factorizations 56 Solves 78 Functions 81 This is run number 3 Using the "cheap" Newton method. The jacobian is held fixed as long as possible, perhaps over multiple points, until convergence fails. Step Type of point Lambda 0 Start point 0.00000 PITCON 7.0 University of Pittsburgh continuation code Modified: 12 November 1999 Linear algebra: LAPACK Precision: double 1 Corrected 0.00000 TRYSTP - DEBUG - about to return. 2 Continuation 0.264435E-01 TRYSTP - DEBUG - about to return. 3 Continuation 0.105774 TRYSTP - DEBUG - about to return. 4 Continuation 0.237992 TRYSTP - DEBUG - about to return. 5 Continuation 0.370209 TRYSTP - DEBUG - about to return. 6 Continuation 0.502427 TRYSTP - DEBUG - about to return. 7 Continuation 0.634644 TRYSTP - DEBUG - about to return. 8 Continuation 0.766862 TRYSTP - DEBUG - about to return. 9 Continuation 0.899079 10 Target point 0.800000 Complete value for target point 1 0.00000 0.603078E-01 0.118491 0.174422 0.227971 0.279007 0.327397 0.373010 0.415717 0.455390 0.491906 0.525147 0.555003 0.581370 0.604155 0.623274 0.638658 0.650248 0.658000 0.661883 0.661883 0.800000 TRYSTP - DEBUG - about to return. UPDATE: Predictor stepsize reductions: 2 11 Continuation 0.913770 TRYSTP - Retrying step with IWORK(4) = 1 TRYSTP - DEBUG - about to return. UPDATE: Predictor stepsize reductions: 2 12 Continuation 0.912943 TRYSTP - DEBUG - about to return. 13 Continuation 0.913845 TRYSTP - DEBUG - about to return. 14 Continuation 0.916548 TRYSTP - DEBUG - about to return. 15 Continuation 0.921919 TRYSTP - DEBUG - about to return. 16 Continuation 0.923950 TRYSTP - DEBUG - about to return. 17 Continuation 0.922335 TRYSTP - DEBUG - about to return. 18 Continuation 0.915652 TRYSTP - DEBUG - about to return. 19 Continuation 0.908083 TRYSTP - DEBUG - about to return. 20 Continuation 0.895757 TRYSTP - DEBUG - about to return. 21 Continuation 0.882348 TRYSTP - DEBUG - about to return. 22 Continuation 0.868383 TRYSTP - DEBUG - about to return. 23 Continuation 0.849582 TRYSTP - DEBUG - about to return. UPDATE: Predictor stepsize reductions: 1 24 Continuation 0.844720 TRYSTP - DEBUG - about to return. 25 Continuation 0.839566 TRYSTP - DEBUG - about to return. UPDATE: Predictor stepsize reductions: 2 26 Continuation 0.835549 TRYSTP - DEBUG - about to return. 27 Continuation 0.835708 TRYSTP - DEBUG - about to return. UPDATE: Predictor stepsize reductions: 1 28 Continuation 0.834379 TRYSTP - DEBUG - about to return. 29 Continuation 0.832980 TRYSTP - DEBUG - about to return. 30 Continuation 0.828594 TRYSTP - Retrying step with IWORK(4) = 1 TRYSTP - DEBUG - about to return. UPDATE: Predictor stepsize reductions: 1 31 Continuation 0.827412 LIMIT - Warning! Iteration did not reach limit point after taking 25 steps. PITCON - Warning! An error occurred during the limit point computation. The computed limit point does not satisfy the accuracy requirements. However, the code can continue. 32 Limit point 0.828594 PITCON returned an error code: IERROR = 8 The computation failed.