PCPRB7: PITCON sample program Number of equations is 25 Number of variables is 26 The materially nonlinear problem. Piecewise linears 1 continuity condition at breakpoints. 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 CHECKW - Machine epsilon = 2.22044604925031308E-16 START - Checking the initial point. Fixing variable number 5 Step FX X DX 0 0.141421E-02 0.100000E-02 19 5 43.5309 26 1 0.849584E-02 43.5309 7 26 CORRECTOR - Warning! The residual FX is not decreasing. START - Warning! The starting point needs to be corrected. The first correction attempt failed. Correction will be retried once. Step FX X DX 0 0.141421E-02 0.100000E-02 19 5 43.5309 26 1 0.849584E-02 43.5309 7 26 7.90641 26 2 0.190847E-02 35.6245 11 26 7.54099 26 3 0.695215E-03 43.1655 11 26 1.54283 26 4 0.309331E-04 41.6227 11 26 0.758943E-01 26 5 0.737807E-07 41.5468 11 26 0.181523E-03 26 6 0.423899E-12 41.5466 11 26 1 corrected 41.5466 TANPAR - Tangent vector has limit component = 1.4644191972016257E-3 TRYSTP - Predictor using stepsize 0.25 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.247182E-02 41.5478 3 26 0.274926E-02 26 1 0.655695E-06 41.5506 7 26 TRYSTP - DEBUG - about to return. 2 continuation 41.5497 TANPAR - Tangent vector has limit component = 2.34209226939633662E-2 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.984907E-02 41.5685 3 26 0.106918E-01 26 1 0.104051E-04 41.5792 7 26 0.252490E-04 26 2 0.188437E-10 41.5792 11 26 TRYSTP - DEBUG - about to return. 3 continuation 41.5721 TANPAR - Tangent vector has limit component = 6.55981398514218772E-2 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.965804E-02 41.6242 3 26 0.961276E-02 26 1 0.101033E-04 41.6338 7 26 0.172126E-04 26 2 0.190994E-10 41.6338 11 26 TRYSTP - DEBUG - about to return. 4 continuation 41.6146 TANPAR - Tangent vector has limit component = 0.10271971187222523 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.932711E-02 41.7031 3 26 0.781534E-02 26 1 0.960825E-05 41.7109 7 26 0.131942E-04 26 2 0.167646E-10 41.7110 11 26 TRYSTP - DEBUG - about to return. 5 continuation 41.6737 Current rod coordinates: X, U(X), U'(X) 0.625000E-01 0.437878E-01 0.700604 0.187500 0.103514 0.255017 0.312500 0.966687E-01 -0.364545 0.437500 0.277755E-01 -0.737745 0.562500 -0.594968E-01 -0.658612 0.687500 -0.114651 -0.223857 0.812500 -0.108487 0.322487 0.937500 -0.441656E-01 0.706650 TANPAR - Tangent vector has limit component = 0.13201279011273664 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.888363E-02 41.8007 3 26 0.540755E-02 26 1 0.897048E-05 41.8061 7 26 0.118950E-04 25 2 0.135127E-10 41.8061 11 26 TRYSTP - DEBUG - about to return. 6 continuation 41.7452 TANPAR - Tangent vector has limit component = 0.15128762260777448 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.835695E-02 41.9110 3 26 0.500641E-02 19 1 0.823837E-05 41.9135 7 26 0.107697E-04 25 2 0.101606E-10 41.9135 11 26 TRYSTP - DEBUG - about to return. 7 continuation 41.8233 TANPAR - Tangent vector has limit component = 0.15887547748674813 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.777115E-02 42.0278 3 26 0.476251E-02 19 1 0.744447E-05 42.0270 7 26 0.957465E-05 25 2 0.715487E-11 42.0270 11 26 TRYSTP - DEBUG - about to return. 8 continuation 41.9020 TANPAR - Tangent vector has limit component = 0.15346393184935536 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.714179E-02 42.1441 3 26 0.442411E-02 19 1 0.659976E-05 42.1398 7 26 0.833611E-05 25 2 0.470582E-11 42.1398 9 26 TRYSTP - DEBUG - about to return. 9 continuation 41.9744 TANPAR - Tangent vector has limit component = 0.13397141278568234 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.647708E-02 42.2528 3 26 0.804572E-02 26 1 0.569608E-05 42.2448 7 26 0.705227E-05 25 2 0.288288E-11 42.2448 9 26 TRYSTP - DEBUG - about to return. 10 continuation 42.0334 Current rod coordinates: X, U(X), U'(X) 0.625000E-01 0.104194 1.66710 0.187500 0.238404 0.480250 0.312500 0.199790 -1.09808 0.437500 0.173595E-01 -1.82081 0.562500 -0.185404 -1.42342 0.687500 -0.296851 -0.359736 0.812500 -0.266018 0.853068 0.937500 -0.106351 1.70161 TANPAR - Tangent vector has limit component = 9.95828289664027894E-2 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 20 Step FX X DX 0 0.578308E-02 42.3464 3 26 0.119488E-01 26 1 0.524932E-05 42.3342 5 26 0.137785E-04 26 2 0.333593E-11 42.3342 5 26 TRYSTP - DEBUG - about to return. 11 continuation 42.0712 TANPAR - Tangent vector has limit component = 5.03817876487735183E-2 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 20 Step FX X DX 0 0.507796E-02 42.4171 3 26 0.153552E-01 26 1 0.393540E-05 42.4017 5 26 0.733829E-05 26 2 0.155157E-11 42.4016 5 26 TRYSTP - DEBUG - about to return. 12 continuation 42.0810 TANPAR - Tangent vector has limit component = -1.31035416309155175E-2 LIMIT - Attempt correction of approximate limit point. LIMIT - For S = 0.E+0 Tan(X(S))(Lim)= 5.03817876487735183E-2 LIMIT - For S = 1. Tan(X(S))(Lim)= -1.31035416309155175E-2 Step FX X DX 0 0.785168E-03 42.3873 3 26 0.264401E-02 26 1 0.958062E-07 42.3899 5 26 LIMIT - For S = 0.79359732745203304 Tan(X(S))(Lim)= 1.0219936508729583E-3 Step FX X DX 0 0.129668E-04 42.3908 3 26 LIMIT - For S = 0.80853072407447413 Tan(X(S))(Lim)= 3.03977063087538449E-5 13 limit point 42.0816 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 20 Step FX X DX 0 0.439162E-02 42.4588 3 26 0.184176E-01 26 1 0.263974E-05 42.4405 5 26 0.239989E-05 26 2 0.538560E-12 42.4405 5 26 TRYSTP - DEBUG - about to return. 14 continuation 42.0561 TANPAR - Tangent vector has limit component = -8.82969313457015875E-2 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 20 Step FX X DX 0 0.379593E-02 42.4659 3 26 0.207586E-01 26 1 0.143165E-05 42.4454 5 26 0.132630E-05 26 2 0.868827E-13 42.4454 7 26 TRYSTP - DEBUG - about to return. 15 continuation 41.9912 Current rod coordinates: X, U(X), U'(X) 0.625000E-01 0.151954 2.43126 0.187500 0.339558 0.570406 0.312500 0.262600 -1.80175 0.437500 -0.172616E-01 -2.67603 0.562500 -0.306740 -1.95563 0.687500 -0.453680 -0.395402 0.812500 -0.395673 1.32351 0.937500 -0.156477 2.50363 TANPAR - Tangent vector has limit component = -0.17044876744628307 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 20 Step FX X DX 0 0.339321E-02 42.4354 11 26 0.219981E-01 26 1 0.558401E-06 42.4138 11 26 0.393152E-05 26 2 0.513709E-13 42.4138 3 26 TRYSTP - DEBUG - about to return. 16 continuation 41.8839 TANPAR - Tangent vector has limit component = -0.25320859022097564 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 20 Step FX X DX 0 0.326297E-02 42.3674 11 26 0.219059E-01 26 1 0.714175E-06 42.3461 3 26 0.544770E-05 26 2 0.283184E-13 42.3461 11 26 TRYSTP - DEBUG - about to return. 17 continuation 41.7354 TANPAR - Tangent vector has limit component = -0.33017333487516726 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 20 Step FX X DX 0 0.339479E-02 42.2657 11 26 0.205067E-01 26 1 0.108512E-05 42.2460 5 26 0.597819E-05 26 2 0.100604E-12 42.2460 5 26 TRYSTP - DEBUG - about to return. 18 continuation 41.5498 TANPAR - Tangent vector has limit component = -0.39652585631857784 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 20 Step FX X DX 0 0.371100E-02 42.1364 11 26 0.180425E-01 26 1 0.123100E-05 42.1193 3 26 0.573431E-05 26 2 0.135404E-12 42.1193 5 26 TRYSTP - DEBUG - about to return. 19 continuation 41.3335 TANPAR - Tangent vector has limit component = -0.44971528059323929 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 26 Step FX X DX 0 0.414944E-02 41.9870 11 26 0.113276E-01 20 1 0.148412E-04 41.9823 11 26 0.635542E-04 20 2 0.496237E-09 41.9822 11 26 TRYSTP - DEBUG - about to return. 20 continuation 41.1087 Current rod coordinates: X, U(X), U'(X) 0.625000E-01 0.207647 3.32236 0.187500 0.456017 0.651558 0.312500 0.332523 -2.62747 0.437500 -0.603170E-01 -3.65796 0.562500 -0.451599 -2.60256 0.687500 -0.642135 -0.446018 0.812500 -0.550624 1.91020 0.937500 -0.215618 3.44989 TANPAR - Tangent vector has limit component = -0.48705910958338661 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 26 Step FX X DX 0 0.466246E-02 41.8349 11 26 0.801142E-02 20 1 0.784694E-05 41.8316 11 26 0.278170E-04 20 2 0.102864E-09 41.8316 11 26 TRYSTP - DEBUG - about to return. 21 continuation 40.8652 TANPAR - Tangent vector has limit component = -0.51255732572531798 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 26 Step FX X DX 0 0.531103E-02 41.6756 3 26 0.500615E-02 20 1 0.331100E-05 41.6737 11 26 0.103780E-04 20 2 0.151290E-10 41.6737 11 26 TRYSTP - DEBUG - about to return. 22 continuation 40.6089 TANPAR - Tangent vector has limit component = -0.52566639305080909 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 26 Step FX X DX 0 0.615713E-02 41.5157 3 26 0.430050E-02 19 1 0.217723E-05 41.5154 7 26 0.331884E-05 17 2 0.133005E-11 41.5154 3 26 TRYSTP - DEBUG - about to return. 23 continuation 40.3460 TANPAR - Tangent vector has limit component = -0.52516274712012556 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.730077E-02 41.3631 3 26 0.482999E-02 19 1 0.333327E-05 41.3676 7 26 0.564118E-05 26 2 0.945473E-12 41.3676 11 26 TRYSTP - DEBUG - about to return. 24 continuation 40.0881 TANPAR - Tangent vector has limit component = -0.50867154784990987 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.886860E-02 41.2299 3 26 0.109250E-01 26 1 0.512281E-05 41.2403 7 26 0.819509E-05 26 2 0.227857E-11 41.2402 11 26 TRYSTP - DEBUG - about to return. 25 continuation 39.8447 Current rod coordinates: X, U(X), U'(X) 0.625000E-01 0.258191 4.13105 0.187500 0.565573 0.787067 0.312500 0.411620 -3.25033 0.437500 -0.760949E-01 -4.55311 0.562500 -0.566772 -3.29773 0.687500 -0.808756 -0.574012 0.812500 -0.691742 2.44624 0.937500 -0.269426 4.31081 TANPAR - Tangent vector has limit component = -0.47187968365066074 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.110606E-01 41.1283 3 26 0.185919E-01 26 1 0.797610E-05 41.1459 7 26 0.114148E-04 26 2 0.509869E-11 41.1459 13 26 TRYSTP - DEBUG - about to return. 26 continuation 39.6274 TANPAR - Tangent vector has limit component = -0.40675927321220612 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 17 Step FX X DX 0 0.141603E-01 41.0748 3 26 0.283682E-01 26 1 0.127765E-04 41.1017 7 26 0.154628E-04 26 2 0.111433E-10 41.1017 13 26 TRYSTP - DEBUG - about to return. 27 continuation 39.4523 TANPAR - Tangent vector has limit component = -0.29990120164068318 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 21 Step FX X DX 0 0.184599E-01 41.0924 3 26 0.418188E-01 26 1 0.183930E-04 41.1299 7 26 0.189522E-04 17 2 0.134541E-09 41.1299 3 26 TRYSTP - DEBUG - about to return. 28 continuation 39.3442 TANPAR - Tangent vector has limit component = -0.13664214246113671 SETSTP - Next stepsize HTAN = 0.5 TRYSTP - Predictor using stepsize 0.5 TRYSTP - Corrector is fixing index 21 Step FX X DX 0 0.237452E-01 41.2082 3 26 0.546764E-01 26 1 0.287896E-04 41.2572 7 26 0.271138E-04 17 2 0.343103E-09 41.2572 3 26 TRYSTP - DEBUG - about to return. 29 continuation 39.3306 TANPAR - Tangent vector has limit component = 8.90843383117656124E-2 LIMIT - Attempt correction of approximate limit point. LIMIT - For S = 0.E+0 Tan(X(S))(Lim)= -0.13664214246113671 LIMIT - For S = 1. Tan(X(S))(Lim)= 8.90843383117656124E-2 Step FX X DX 0 0.572224E-02 41.2063 3 26 0.133754E-01 26 1 0.166847E-05 41.1943 7 26 0.162726E-05 17 2 0.113148E-11 41.1943 3 26 LIMIT - For S = 0.60534387455678673 Tan(X(S))(Lim)= -4.77980513518408943E-3 Step FX X DX 0 0.183826E-03 41.1975 3 26 0.431020E-03 26 1 0.170602E-08 41.1971 7 26 LIMIT - For S = 0.62544078596248742 Tan(X(S))(Lim)= -1.17750226395618685E-4 Step FX X DX 0 0.866058E-03 41.2271 3 26 0.202213E-02 26 1 0.376063E-07 41.2252 7 26 LIMIT - For S = 0.81272039298124366 Tan(X(S))(Lim)= 4.39949689398772825E-2 Step FX X DX 0 0.161419E-02 41.2141 3 26 0.376010E-02 26 1 0.132610E-06 41.2107 7 26 LIMIT - For S = 0.71908058947186548 Tan(X(S))(Lim)= 2.17958169206427278E-2 Step FX X DX 0 0.750401E-03 41.2053 3 26 0.174854E-02 26 1 0.287476E-07 41.2038 7 26 LIMIT - For S = 0.6722606877171764 Tan(X(S))(Lim)= 1.08008042019544451E-2 Step FX X DX 0 0.688377E-03 41.1986 3 26 0.160495E-02 26 1 0.242423E-07 41.1972 7 26 LIMIT - For S = 0.62594571123145171 Tan(X(S))(Lim)= -4.00258671556225105E-7 30 limit point 39.3225 Current rod coordinates: X, U(X), U'(X) 0.625000E-01 0.297735 4.76376 0.187500 0.656115 0.970321 0.312500 0.495043 -3.54747 0.437500 -0.559411E-01 -5.26828 0.562500 -0.632143 -3.95095 0.687500 -0.927856 -0.780450 0.812500 -0.799982 2.82643 0.937500 -0.311665 4.98664 Jacobians 104 Factorizations 104 Solves 104 Functions 105