Nonlinear problems with sparse Jacobian: Chemical equilibrium problem N L E Q 1 S ***** V e r s i o n 2 . 3 . 2 *** Newton-Method for the solution of nonlinear systems Real Workspace declared as 248 is used up to 248 (100.0 percent) Integer Workspace declared as 344 is used up to 344 (100.0 percent) N = 6 Maximum number of nonzero Jacobian elements : 18 Prescribed relative precision 0.10D-04 Automatic row scaling of the Jacobian is allowed Rank-1 updates are inhibited Problem is specified as being highly nonlinear Bounded damping strategy is off Maximum permitted number of iteration steps : 50 Internal parameters: Starting value for damping factor FCSTART = 0.10D-01 Minimum allowed damping factor FCMIN = 0.10D-03 Rank-1 updates decision parameter SIGMA = 0.10D+06 ****************************************************************** It Normf Normx Damp.Fct. New 0 0.170D+05 0.820D+00 0 1 0.168D+05 * 0.812D+00 0.01000 Returned from call 1 of NLEQ1S 1 0.168D+05 0.810D+00 0 2 0.420D+04 * 0.112D+00 1.00000 Returned from call 2 of NLEQ1S 2 0.420D+04 0.255D+00 0 3 0.105D+04 * 0.637D-01 1.00000 Returned from call 3 of NLEQ1S 3 0.105D+04 0.239D+00 0 4 0.263D+03 * 0.597D-01 1.00000 Returned from call 4 of NLEQ1S 4 0.263D+03 0.236D+00 0 5 0.657D+02 * 0.591D-01 1.00000 Returned from call 5 of NLEQ1S 5 0.657D+02 0.236D+00 0 6 0.164D+02 * 0.589D-01 1.00000 Returned from call 6 of NLEQ1S 6 0.164D+02 0.236D+00 0 7 0.410D+01 * 0.589D-01 1.00000 Returned from call 7 of NLEQ1S 7 0.410D+01 0.214D+00 0 8 0.103D+01 * 0.535D-01 1.00000 Returned from call 8 of NLEQ1S 8 0.103D+01 0.198D+00 0 9 0.257D+00 * 0.494D-01 1.00000 Returned from call 9 of NLEQ1S 9 0.257D+00 0.194D+00 0 10 0.641D-01 * 0.483D-01 1.00000 Returned from call 10 of NLEQ1S 10 0.641D-01 0.193D+00 0 11 0.160D-01 * 0.481D-01 1.00000 Returned from call 11 of NLEQ1S 11 0.160D-01 0.192D+00 0 12 0.401D-02 * 0.483D-01 1.00000 Returned from call 12 of NLEQ1S 12 0.401D-02 0.193D+00 0 13 0.100D-02 * 0.493D-01 1.00000 Returned from call 13 of NLEQ1S 13 0.100D-02 0.196D+00 0 14 0.251D-03 * 0.532D-01 1.00000 Returned from call 14 of NLEQ1S 14 0.251D-03 0.206D+00 0 15 0.638D-04 * 0.664D-01 1.00000 Returned from call 15 of NLEQ1S 15 0.638D-04 0.242D+00 0 16 0.198D-04 * 0.103D+00 1.00000 Returned from call 16 of NLEQ1S 16 0.198D-04 0.346D+00 0 17 0.115D-04 * 0.189D+00 0.71943 Returned from call 17 of NLEQ1S 17 0.115D-04 0.470D+00 0 18 0.836D-05 * 0.278D+00 0.61195 Returned from call 18 of NLEQ1S 18 0.836D-05 0.427D+00 0 19 0.734D-05 * 0.232D+00 0.72178 Returned from call 19 of NLEQ1S 19 0.734D-05 0.395D+00 0 20 0.652D-05 * 0.201D+00 0.76771 Returned from call 20 of NLEQ1S 20 0.652D-05 0.348D+00 0 21 0.553D-05 * 0.129D+00 1.00000 Returned from call 21 of NLEQ1S 21 0.553D-05 0.243D+00 0 22 0.140D-05 * 0.376D-01 1.00000 Returned from call 22 of NLEQ1S 22 0.140D-05 0.109D+00 0 23 0.344D-07 * 0.204D-02 1.00000 Returned from call 23 of NLEQ1S 23 0.344D-07 0.516D-02 0 24 0.925D-11 * 0.121D-05 1.00000 Solution of nonlinear system of equations obtained within 24 iteration steps Achieved relative accuracy 0.121D-05 ############################################################################# # # # # # Results from time monitor program for: NLEQ1S # # # # Total time: 0.006 Sum of parts: 0.005 # # # # Name Calls Time Av-time % Total % Sum # # NLEQ1S 1 0.001 0.0010 16.46 # # FCN 25 0.000 0.0000 0.43 0.51 # # Jacobi 24 0.000 0.0000 4.42 5.29 # # NIFACT 24 0.000 0.0000 4.71 5.63 # # NISOLV 48 0.000 0.0000 1.63 1.95 # # Output 73 0.005 0.0001 72.35 86.61 # # # ############################################################################# ****** Statistics * NLEQ1S ******* *** Newton iterations : 24 *** *** Corrector steps : 0 *** *** Rejected rk-1 st. : 0 *** *** Jacobian eval. : 24 *** *** Function eval. : 25 *** ************************************* Returned from call 24 of NLEQ1S Time used = 0.008 Sec ******************************************************************