test_simeq_newton3.c running first A generated A[ 0][ 0]=0.000000 A[ 0][ 1]=0.000985 A[ 0][ 2]=0.041631 A[ 0][ 3]=0.176643 A[ 0][ 4]=0.364602 A[ 0][ 5]=0.091331 A[ 0][ 6]=0.092298 A[ 0][ 7]=0.487217 A[ 0][ 8]=0.526750 A[ 0][ 9]=0.454433 A[ 1][ 0]=0.233178 A[ 1][ 1]=0.831292 A[ 1][ 2]=0.931731 A[ 1][ 3]=0.568060 A[ 1][ 4]=0.556094 A[ 1][ 5]=0.050832 A[ 1][ 6]=0.767051 A[ 1][ 7]=0.018915 A[ 1][ 8]=0.252360 A[ 1][ 9]=0.298197 A[ 2][ 0]=0.875981 A[ 2][ 1]=0.531557 A[ 2][ 2]=0.920261 A[ 2][ 3]=0.515431 A[ 2][ 4]=0.810429 A[ 2][ 5]=0.188420 A[ 2][ 6]=0.886314 A[ 2][ 7]=0.570614 A[ 2][ 8]=0.076775 A[ 2][ 9]=0.815274 A[ 3][ 0]=0.984891 A[ 3][ 1]=0.118352 A[ 3][ 2]=0.893906 A[ 3][ 3]=0.784484 A[ 3][ 4]=0.100916 A[ 3][ 5]=0.253311 A[ 3][ 6]=0.019842 A[ 3][ 7]=0.378377 A[ 3][ 8]=0.678877 A[ 3][ 9]=0.680923 solve system of equations A X = Y for X = X0, X1, X2, X3, X0 * X1 X0 * X2 X0 * X3 X1 * X2 X1 * X3 X2 * X3 X_soln[ 0]=1.100000 X_soln[ 1]=1.200000 X_soln[ 2]=1.400000 X_soln[ 3]=1.500000 X_soln[ 4]=1.320000 X_soln[ 5]=1.540000 X_soln[ 6]=1.650000 X_soln[ 7]=1.680000 X_soln[ 8]=1.800000 X_soln[ 9]=2.100000 Y[ 0]=3.819631 Y[ 1]=6.600759 Y[ 2]=9.294213 Y[ 3]=7.497234 simeq_newton3 running 4 equations, 4 variables, with 6 nonlinear terms X[ 0]=1.000000 X[ 1]=1.000000 X[ 2]=1.000000 X[ 3]=1.000000 X[ 4]=1.000000 X[ 5]=1.000000 X[ 6]=1.000000 X[ 7]=1.000000 X[ 8]=1.000000 X[ 9]=1.000000 residual X_resid[ 0]=-1.58374 residual X_resid[ 1]=-2.09305 residual X_resid[ 2]=-3.10316 residual X_resid[ 3]=-2.60336 simeq_newton3 itr=0, prev=1e+12, residual=9.3833 Ja computed Ja[ 0][ 0]=0.548231 Ja[ 0][ 1]=1.379555 Ja[ 0][ 2]=1.074612 Ja[ 0][ 3]=1.250124 Ja[ 1][ 0]=1.607156 Ja[ 1][ 1]=1.658661 Ja[ 1][ 2]=1.299675 Ja[ 1][ 3]=1.885668 Ja[ 2][ 0]=2.761145 Ja[ 2][ 1]=1.989375 Ja[ 2][ 2]=2.494569 Ja[ 2][ 3]=2.293794 Ja[ 3][ 0]=1.358960 Ja[ 3][ 1]=1.276523 Ja[ 3][ 2]=2.206517 Ja[ 3][ 3]=2.164126 Ja inverted invJa[ 0][ 0]=-0.810424 invJa[ 0][ 1]=0.382608 invJa[ 0][ 2]=0.489295 invJa[ 0][ 3]=-0.383842 invJa[ 1][ 0]=1.495549 invJa[ 1][ 1]=-0.382441 invJa[ 1][ 2]=0.390619 invJa[ 1][ 3]=-0.944707 invJa[ 2][ 0]=0.761398 invJa[ 2][ 1]=-1.680925 invJa[ 2][ 2]=0.674920 invJa[ 2][ 3]=0.309454 invJa[ 3][ 0]=-1.149567 invJa[ 3][ 1]=1.699178 invJa[ 3][ 2]=-1.225801 invJa[ 3][ 3]=0.944840 change X_tmp2[ 0]=-0.036397 change X_tmp2[ 1]=-0.320837 change X_tmp2[ 2]=-0.587600 change X_tmp2[ 3]=-0.391746 X_next[ 0]=1.036397 X_next[ 1]=1.320837 X_next[ 2]=1.587600 X_next[ 3]=1.391746 X_next[ 4]=1.368912 X_next[ 5]=1.645384 X_next[ 6]=1.442401 X_next[ 7]=2.096962 X_next[ 8]=1.838270 X_next[ 9]=2.209536 residual X_resid[ 0]=0.270191 residual X_resid[ 1]=0.122444 residual X_resid[ 2]=0.331023 residual X_resid[ 3]=0.32028 simeq_newton3 itr=1, prev=9.3833, residual=1.04394 Ja computed Ja[ 0][ 0]=0.755032 Ja[ 0][ 1]=1.885467 Ja[ 0][ 2]=1.412276 Ja[ 0][ 3]=1.689510 Ja[ 1][ 0]=2.115930 Ja[ 1][ 1]=1.788876 Ja[ 1][ 2]=1.424412 Ja[ 1][ 3]=2.169773 Ja[ 2][ 0]=3.479087 Ja[ 2][ 1]=2.384241 Ja[ 2][ 2]=3.003882 Ja[ 2][ 3]=2.829740 Ja[ 3][ 0]=1.547956 Ja[ 3][ 1]=1.768478 Ja[ 3][ 2]=2.603883 Ja[ 3][ 3]=2.782768 Ja inverted invJa[ 0][ 0]=-0.400658 invJa[ 0][ 1]=0.310176 invJa[ 0][ 2]=0.338071 invJa[ 0][ 3]=-0.342375 invJa[ 1][ 0]=1.181951 invJa[ 1][ 1]=-0.204299 invJa[ 1][ 2]=0.212128 invJa[ 1][ 3]=-0.774014 invJa[ 2][ 0]=0.198673 invJa[ 2][ 1]=-1.323387 invJa[ 2][ 2]=0.650718 invJa[ 2][ 3]=0.249545 invJa[ 3][ 0]=-0.714172 invJa[ 3][ 1]=1.195610 invJa[ 3][ 2]=-0.931755 invJa[ 3][ 3]=0.808196 change X_tmp2[ 0]=-0.068021 change X_tmp2[ 1]=0.116655 change X_tmp2[ 2]=0.186966 change X_tmp2[ 3]=-0.096151 X_next[ 0]=1.104418 X_next[ 1]=1.204182 X_next[ 2]=1.400635 X_next[ 3]=1.487897 X_next[ 4]=1.329921 X_next[ 5]=1.546886 X_next[ 6]=1.643260 X_next[ 7]=1.686619 X_next[ 8]=1.791699 X_next[ 9]=2.083999 residual X_resid[ 0]=-0.0069021 residual X_resid[ 1]=-0.00782104 residual X_resid[ 2]=-0.00610213 residual X_resid[ 3]=-0.0154954 simeq_newton3 itr=2, prev=1.04394, residual=0.0363206 Ja computed Ja[ 0][ 0]=0.704298 Ja[ 0][ 1]=1.869822 Ja[ 0][ 2]=1.405347 Ja[ 0][ 3]=1.549376 Ja[ 1][ 0]=2.115307 Ja[ 1][ 1]=1.847430 Ja[ 1][ 2]=1.454335 Ja[ 1][ 3]=2.136757 Ja[ 2][ 0]=3.434538 Ja[ 2][ 1]=2.340064 Ja[ 2][ 2]=3.028522 Ja[ 2][ 3]=2.728644 Ja[ 3][ 0]=1.490731 Ja[ 3][ 1]=1.769873 Ja[ 3][ 2]=2.642445 Ja[ 3][ 3]=2.577614 Ja inverted invJa[ 0][ 0]=-0.333337 invJa[ 0][ 1]=0.209999 invJa[ 0][ 2]=0.404735 invJa[ 0][ 3]=-0.402167 invJa[ 1][ 0]=1.213306 invJa[ 1][ 1]=-0.228268 invJa[ 1][ 2]=0.233477 invJa[ 1][ 3]=-0.787235 invJa[ 2][ 0]=0.228471 invJa[ 2][ 1]=-1.225402 invJa[ 2][ 2]=0.604162 invJa[ 2][ 3]=0.238924 invJa[ 3][ 0]=-0.874532 invJa[ 3][ 1]=1.291509 invJa[ 3][ 2]=-1.013744 invJa[ 3][ 3]=0.916152 change X_tmp2[ 0]=0.004420 change X_tmp2[ 1]=0.004185 change X_tmp2[ 2]=0.000618 change X_tmp2[ 3]=-0.012075 X_next[ 0]=1.099998 X_next[ 1]=1.199998 X_next[ 2]=1.400016 X_next[ 3]=1.499972 X_next[ 4]=1.319995 X_next[ 5]=1.540015 X_next[ 6]=1.649966 X_next[ 7]=1.680016 X_next[ 8]=1.799962 X_next[ 9]=2.099985 residual X_resid[ 0]=-2.66809e-05 residual X_resid[ 1]=-4.54444e-05 residual X_resid[ 2]=-4.02892e-05 residual X_resid[ 3]=-3.69075e-05 simeq_newton3 itr=3, prev=0.0363206, residual=0.000149322 Ja computed Ja[ 0][ 0]=0.703830 Ja[ 0][ 1]=1.874270 Ja[ 0][ 2]=1.408391 Ja[ 0][ 3]=1.546483 Ja[ 1][ 0]=2.122211 Ja[ 1][ 1]=1.848008 Ja[ 1][ 2]=1.457632 Ja[ 1][ 3]=2.132126 Ja[ 2][ 0]=3.441732 Ja[ 2][ 1]=2.337056 Ja[ 2][ 2]=3.035146 Ja[ 2][ 3]=2.723901 Ja[ 3][ 0]=1.490391 Ja[ 3][ 1]=1.777391 Ja[ 3][ 2]=2.647964 Ja[ 3][ 3]=2.574264 Ja inverted invJa[ 0][ 0]=-0.325565 invJa[ 0][ 1]=0.208677 invJa[ 0][ 2]=0.403487 invJa[ 0][ 3]=-0.404194 invJa[ 1][ 0]=1.214039 invJa[ 1][ 1]=-0.227164 invJa[ 1][ 2]=0.232845 invJa[ 1][ 3]=-0.787563 invJa[ 2][ 0]=0.227464 invJa[ 2][ 1]=-1.223766 invJa[ 2][ 2]=0.606005 invJa[ 2][ 3]=0.235701 invJa[ 3][ 0]=-0.883716 invJa[ 3][ 1]=1.294831 invJa[ 3][ 2]=-1.017723 invJa[ 3][ 3]=0.923793 change X_tmp2[ 0]=-0.000002 change X_tmp2[ 1]=-0.000002 change X_tmp2[ 2]=0.000016 change X_tmp2[ 3]=-0.000028 X_next[ 0]=1.100000 X_next[ 1]=1.200000 X_next[ 2]=1.400000 X_next[ 3]=1.500000 X_next[ 4]=1.320000 X_next[ 5]=1.540000 X_next[ 6]=1.650000 X_next[ 7]=1.680000 X_next[ 8]=1.800000 X_next[ 9]=2.100000 residual X_resid[ 0]=-1.9096e-10 residual X_resid[ 1]=-7.51292e-11 residual X_resid[ 2]=-3.45795e-10 residual X_resid[ 3]=-2.93348e-10 simeq_newton3 itr=4, prev=0.000149322, residual=9.05232e-10 simeq_newton3 finished Returned solution X[ 0]=1.100000, err=2.55367e-11 X[ 1]=1.200000, err=-6.42522e-11 X[ 2]=1.400000, err=-2.30189e-10 X[ 3]=1.500000, err=1.52406e-10 first test finished second A generated A[ 0][ 0]=0.752707 A[ 0][ 1]=0.006229 A[ 0][ 2]=0.624407 A[ 0][ 3]=0.126462 A[ 0][ 4]=0.617550 A[ 0][ 5]=0.770928 A[ 0][ 6]=0.186536 A[ 0][ 7]=0.496722 A[ 0][ 8]=0.509500 A[ 0][ 9]=0.315939 A[ 1][ 0]=0.701703 A[ 1][ 1]=0.367247 A[ 1][ 2]=0.876312 A[ 1][ 3]=0.370995 A[ 1][ 4]=0.525252 A[ 1][ 5]=0.427835 A[ 1][ 6]=0.542743 A[ 1][ 7]=0.667979 A[ 1][ 8]=0.171142 A[ 1][ 9]=0.830779 A[ 2][ 0]=0.861905 A[ 2][ 1]=0.769221 A[ 2][ 2]=0.052579 A[ 2][ 3]=0.966730 A[ 2][ 4]=0.334284 A[ 2][ 5]=0.082097 A[ 2][ 6]=0.270315 A[ 2][ 7]=0.074262 A[ 2][ 8]=0.676226 A[ 2][ 9]=0.810159 A[ 3][ 0]=0.464211 A[ 3][ 1]=0.569185 A[ 3][ 2]=0.512203 A[ 3][ 3]=0.586510 A[ 3][ 4]=0.738759 A[ 3][ 5]=0.871689 A[ 3][ 6]=0.577967 A[ 3][ 7]=0.964704 A[ 3][ 8]=0.333529 A[ 3][ 9]=0.392828 solve system of equations A X = Y for X = X0, X1, X2, X3, X0 * X0 X1 * X1 X2 * X2 X3 * X3 X0 * X1 X2 * X3 X_soln[ 0]=1.100000 X_soln[ 1]=1.200000 X_soln[ 2]=1.400000 X_soln[ 3]=1.500000 X_soln[ 4]=1.210000 X_soln[ 5]=1.440000 X_soln[ 6]=1.960000 X_soln[ 7]=2.250000 X_soln[ 8]=1.320000 X_soln[ 9]=2.100000 Y[ 0]=6.575934 Y[ 1]=8.784807 Y[ 2]=7.208428 Y[ 3]=9.508230 simeq_newton3 running 4 equations, 4 variables, with 6 nonlinear terms X[ 0]=1.000000 X[ 1]=1.000000 X[ 2]=1.000000 X[ 3]=1.000000 X[ 4]=1.000000 X[ 5]=1.000000 X[ 6]=1.000000 X[ 7]=1.000000 X[ 8]=1.000000 X[ 9]=1.000000 residual X_resid[ 0]=-2.16895 residual X_resid[ 1]=-3.30282 residual X_resid[ 2]=-2.31065 residual X_resid[ 3]=-3.49665 simeq_newton3 itr=0, prev=1e+12, residual=11.2791 Ja computed Ja[ 0][ 0]=2.497307 Ja[ 0][ 1]=2.057586 Ja[ 0][ 2]=1.313417 Ja[ 0][ 3]=1.435845 Ja[ 1][ 0]=1.923349 Ja[ 1][ 1]=1.394058 Ja[ 1][ 2]=2.792578 Ja[ 1][ 3]=2.537731 Ja[ 2][ 0]=2.206699 Ja[ 2][ 1]=1.609640 Ja[ 2][ 2]=1.403367 Ja[ 2][ 3]=1.925413 Ja[ 3][ 0]=2.275258 Ja[ 3][ 1]=2.646092 Ja[ 3][ 2]=2.060966 Ja[ 3][ 3]=2.908745 Ja inverted invJa[ 0][ 0]=0.122150 invJa[ 0][ 1]=-0.082815 invJa[ 0][ 2]=1.180439 invJa[ 0][ 3]=-0.769425 invJa[ 1][ 0]=0.810727 invJa[ 1][ 1]=-0.178166 invJa[ 1][ 2]=-1.605828 invJa[ 1][ 3]=0.818202 invJa[ 2][ 0]=0.753112 invJa[ 2][ 1]=0.836562 invJa[ 2][ 2]=-1.403467 invJa[ 2][ 3]=-0.172607 invJa[ 3][ 0]=-1.366679 invJa[ 3][ 1]=-0.365882 invJa[ 3][ 2]=1.531884 invJa[ 3][ 3]=0.323624 change X_tmp2[ 0]=-0.028593 change X_tmp2[ 1]=-0.320434 change X_tmp2[ 2]=-0.550016 change X_tmp2[ 3]=-0.498540 X_next[ 0]=1.028593 X_next[ 1]=1.320434 X_next[ 2]=1.550016 X_next[ 3]=1.498540 X_next[ 4]=1.058004 X_next[ 5]=1.743547 X_next[ 6]=2.402550 X_next[ 7]=2.245621 X_next[ 8]=1.358190 X_next[ 9]=2.322761 residual X_resid[ 0]=0.350849 residual X_resid[ 1]=0.60394 residual X_resid[ 2]=0.33728 residual X_resid[ 3]=0.615493 simeq_newton3 itr=1, prev=11.2791, residual=1.90756 Ja computed Ja[ 0][ 0]=2.695884 Ja[ 0][ 1]=2.566218 Ja[ 0][ 2]=1.676121 Ja[ 0][ 3]=2.104887 Ja[ 1][ 0]=2.008226 Ja[ 1][ 1]=1.673137 Ja[ 1][ 2]=3.803789 Ja[ 1][ 3]=3.660701 Ja[ 2][ 0]=2.442501 Ja[ 2][ 1]=1.681589 Ja[ 2][ 2]=2.104618 Ja[ 2][ 3]=2.445059 Ja[ 3][ 0]=2.424379 Ja[ 3][ 1]=3.214267 Ja[ 3][ 2]=2.892589 Ja[ 3][ 3]=4.086693 Ja inverted invJa[ 0][ 0]=-0.137256 invJa[ 0][ 1]=-0.334546 invJa[ 0][ 2]=1.153181 invJa[ 0][ 3]=-0.319577 invJa[ 1][ 0]=1.096860 invJa[ 1][ 1]=0.309285 invJa[ 1][ 2]=-1.549202 invJa[ 1][ 3]=0.084892 invJa[ 2][ 0]=1.072353 invJa[ 2][ 1]=1.087096 invJa[ 2][ 2]=-1.385302 invJa[ 2][ 3]=-0.697280 invJa[ 3][ 0]=-1.540295 invJa[ 3][ 1]=-0.814247 invJa[ 3][ 2]=1.514895 invJa[ 3][ 3]=0.861052 change X_tmp2[ 0]=-0.057954 change X_tmp2[ 1]=0.101357 change X_tmp2[ 2]=0.136369 change X_tmp2[ 3]=0.008748 X_next[ 0]=1.086548 X_next[ 1]=1.219078 X_next[ 2]=1.413647 X_next[ 3]=1.489792 X_next[ 4]=1.180586 X_next[ 5]=1.486150 X_next[ 6]=1.998398 X_next[ 7]=2.219480 X_next[ 8]=1.324586 X_next[ 9]=2.106040 residual X_resid[ 0]=0.010885 residual X_resid[ 1]=0.0162894 residual X_resid[ 2]=0.00399304 residual X_resid[ 3]=0.0207677 simeq_newton3 itr=2, prev=1.90756, residual=0.0519352 Ja computed Ja[ 0][ 0]=2.715822 Ja[ 0][ 1]=2.439469 Ja[ 0][ 2]=1.622482 Ja[ 0][ 3]=2.053112 Ja[ 1][ 0]=2.051761 Ja[ 1][ 1]=1.596328 Ja[ 1][ 2]=3.648495 Ja[ 1][ 3]=3.535721 Ja[ 2][ 0]=2.412708 Ja[ 2][ 1]=1.704137 Ja[ 2][ 2]=2.023806 Ja[ 2][ 3]=2.333279 Ja[ 3][ 0]=2.476202 Ja[ 3][ 1]=3.056893 Ja[ 3][ 2]=2.731518 Ja[ 3][ 3]=4.016245 Ja inverted invJa[ 0][ 0]=-0.492087 invJa[ 0][ 1]=-0.505078 invJa[ 0][ 2]=1.692560 invJa[ 0][ 3]=-0.287106 invJa[ 1][ 0]=1.603413 invJa[ 1][ 1]=0.522691 invJa[ 1][ 2]=-2.317881 invJa[ 1][ 3]=0.066775 invJa[ 2][ 0]=1.360665 invJa[ 2][ 1]=1.212035 invJa[ 2][ 2]=-1.865851 invJa[ 2][ 3]=-0.678610 invJa[ 3][ 0]=-1.842426 invJa[ 3][ 1]=-0.910758 invJa[ 3][ 2]=1.989669 invJa[ 3][ 3]=0.836713 change X_tmp2[ 0]=-0.012788 change X_tmp2[ 1]=0.018099 change X_tmp2[ 2]=0.013011 change X_tmp2[ 3]=-0.009569 X_next[ 0]=1.099336 X_next[ 1]=1.200979 X_next[ 2]=1.400636 X_next[ 3]=1.499361 X_next[ 4]=1.208539 X_next[ 5]=1.442350 X_next[ 6]=1.961782 X_next[ 7]=2.248083 X_next[ 8]=1.320279 X_next[ 9]=2.100059 residual X_resid[ 0]=0.000273324 residual X_resid[ 1]=0.000236037 residual X_resid[ 2]=-0.000123259 residual X_resid[ 3]=0.000466419 simeq_newton3 itr=3, prev=0.0519352, residual=0.00109904 Ja computed Ja[ 0][ 0]=2.722395 Ja[ 0][ 1]=2.418078 Ja[ 0][ 2]=1.620651 Ja[ 0][ 3]=2.058508 Ja[ 1][ 0]=2.062098 Ja[ 1][ 1]=1.583029 Ja[ 1][ 2]=3.642322 Ja[ 1][ 3]=3.537696 Ja[ 2][ 0]=2.409019 Ja[ 2][ 1]=1.709813 Ja[ 2][ 2]=2.024525 Ja[ 2][ 3]=2.324159 Ja[ 3][ 0]=2.489060 Ja[ 3][ 1]=3.029605 Ja[ 3][ 2]=2.720238 Ja[ 3][ 3]=4.029597 Ja inverted invJa[ 0][ 0]=-0.623364 invJa[ 0][ 1]=-0.569679 invJa[ 0][ 2]=1.884409 invJa[ 0][ 3]=-0.268294 invJa[ 1][ 0]=1.795124 invJa[ 1][ 1]=0.613083 invJa[ 1][ 2]=-2.598117 invJa[ 1][ 3]=0.043244 invJa[ 2][ 0]=1.480004 invJa[ 2][ 1]=1.267790 invJa[ 2][ 2]=-2.045617 invJa[ 2][ 3]=-0.689230 invJa[ 3][ 0]=-1.963693 invJa[ 3][ 1]=-0.964892 invJa[ 3][ 2]=2.170298 invJa[ 3][ 3]=0.846649 change X_tmp2[ 0]=-0.000662 change X_tmp2[ 1]=0.000976 change X_tmp2[ 2]=0.000634 change X_tmp2[ 3]=-0.000637 X_next[ 0]=1.099998 X_next[ 1]=1.200003 X_next[ 2]=1.400002 X_next[ 3]=1.499998 X_next[ 4]=1.209995 X_next[ 5]=1.440007 X_next[ 6]=1.960005 X_next[ 7]=2.249994 X_next[ 8]=1.320001 X_next[ 9]=2.100000 residual X_resid[ 0]=8.24623e-07 residual X_resid[ 1]=6.80915e-07 residual X_resid[ 2]=-4.00722e-07 residual X_resid[ 3]=1.40386e-06 simeq_newton3 itr=4, prev=0.00109904, residual=3.31012e-06 Ja computed Ja[ 0][ 0]=2.722716 Ja[ 0][ 1]=2.416911 Ja[ 0][ 2]=1.620616 Ja[ 0][ 3]=2.058941 Ja[ 1][ 0]=2.062626 Ja[ 1][ 1]=1.582308 Ja[ 1][ 2]=3.642162 Ja[ 1][ 3]=3.538020 Ja[ 2][ 0]=2.408802 Ja[ 2][ 1]=1.710100 Ja[ 2][ 2]=2.024698 Ja[ 2][ 3]=2.323740 Ja[ 3][ 0]=2.489713 Ja[ 3][ 1]=3.028125 Ja[ 3][ 2]=2.719755 Ja[ 3][ 3]=4.030577 Ja inverted invJa[ 0][ 0]=-0.631623 invJa[ 0][ 1]=-0.573770 invJa[ 0][ 2]=1.896360 invJa[ 0][ 3]=-0.267000 invJa[ 1][ 0]=1.807213 invJa[ 1][ 1]=0.618856 invJa[ 1][ 2]=-2.615603 invJa[ 1][ 3]=0.041560 invJa[ 2][ 0]=1.487559 invJa[ 2][ 1]=1.271356 invJa[ 2][ 2]=-2.056815 invJa[ 2][ 3]=-0.690069 invJa[ 3][ 0]=-1.971356 invJa[ 3][ 1]=-0.968404 invJa[ 3][ 2]=2.181576 invJa[ 3][ 3]=0.847453 change X_tmp2[ 0]=-0.000002 change X_tmp2[ 1]=0.000003 change X_tmp2[ 2]=0.000002 change X_tmp2[ 3]=-0.000002 X_next[ 0]=1.100000 X_next[ 1]=1.200000 X_next[ 2]=1.400000 X_next[ 3]=1.500000 X_next[ 4]=1.210000 X_next[ 5]=1.440000 X_next[ 6]=1.960000 X_next[ 7]=2.250000 X_next[ 8]=1.320000 X_next[ 9]=2.100000 residual X_resid[ 0]=7.88436e-12 residual X_resid[ 1]=6.50147e-12 residual X_resid[ 2]=-3.8245e-12 residual X_resid[ 3]=1.34026e-11 simeq_newton3 itr=5, prev=3.31012e-06, residual=3.16129e-11 simeq_newton3 finished Returned solution X[ 0]=1.100000, err=-1.95419e-11 X[ 1]=1.200000, err=2.8834e-11 X[ 2]=1.400000, err=1.86133e-11 X[ 3]=1.500000, err=-1.88267e-11 second test finished third A generated A[ 0][ 0]=1.000000 A[ 0][ 1]=1.000000 A[ 0][ 2]=1.000000 A[ 0][ 3]=1.000000 A[ 0][ 4]=1.000000 A[ 0][ 5]=0.000000 A[ 0][ 6]=0.000000 A[ 0][ 7]=0.000000 A[ 0][ 8]=0.000000 A[ 0][ 9]=0.000000 A[ 1][ 0]=1.000000 A[ 1][ 1]=1.000000 A[ 1][ 2]=1.000000 A[ 1][ 3]=1.000000 A[ 1][ 4]=0.000000 A[ 1][ 5]=1.000000 A[ 1][ 6]=0.000000 A[ 1][ 7]=0.000000 A[ 1][ 8]=0.000000 A[ 1][ 9]=0.000000 A[ 2][ 0]=1.000000 A[ 2][ 1]=1.000000 A[ 2][ 2]=1.000000 A[ 2][ 3]=1.000000 A[ 2][ 4]=0.000000 A[ 2][ 5]=0.000000 A[ 2][ 6]=1.000000 A[ 2][ 7]=0.000000 A[ 2][ 8]=0.000000 A[ 2][ 9]=0.000000 A[ 3][ 0]=1.000000 A[ 3][ 1]=1.000000 A[ 3][ 2]=1.000000 A[ 3][ 3]=1.000000 A[ 3][ 4]=0.000000 A[ 3][ 5]=0.000000 A[ 3][ 6]=0.000000 A[ 3][ 7]=1.000000 A[ 3][ 8]=0.000000 A[ 3][ 9]=0.000000 solve system of equations A X = Y for X = X0, X1, X2, X3, X0 * X0* X0 X1 * X1* X1 X2 * X2* X2 X3 * X3* X3 X0 * X1* X2 X1 * X2* X3 X_soln[ 0]=1.100000 X_soln[ 1]=1.200000 X_soln[ 2]=1.400000 X_soln[ 3]=1.500000 X_soln[ 4]=1.331000 X_soln[ 5]=1.728000 X_soln[ 6]=2.744000 X_soln[ 7]=3.375000 X_soln[ 8]=1.848000 X_soln[ 9]=2.520000 Y[ 0]=6.531000 Y[ 1]=6.928000 Y[ 2]=7.944000 Y[ 3]=8.575000 simeq_newton3 running 4 equations, 4 variables, with 6 nonlinear terms X[ 0]=0.500000 X[ 1]=0.500000 X[ 2]=0.500000 X[ 3]=0.500000 X[ 4]=0.125000 X[ 5]=0.125000 X[ 6]=0.125000 X[ 7]=0.125000 X[ 8]=0.125000 X[ 9]=0.125000 residual X_resid[ 0]=-4.406 residual X_resid[ 1]=-4.803 residual X_resid[ 2]=-5.819 residual X_resid[ 3]=-6.45 simeq_newton3 itr=0, prev=1e+12, residual=21.478 Ja computed Ja[ 0][ 0]=1.750000 Ja[ 0][ 1]=1.000000 Ja[ 0][ 2]=1.000000 Ja[ 0][ 3]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=1.750000 Ja[ 1][ 2]=1.000000 Ja[ 1][ 3]=1.000000 Ja[ 2][ 0]=1.000000 Ja[ 2][ 1]=1.000000 Ja[ 2][ 2]=1.750000 Ja[ 2][ 3]=1.000000 Ja[ 3][ 0]=1.000000 Ja[ 3][ 1]=1.000000 Ja[ 3][ 2]=1.000000 Ja[ 3][ 3]=1.750000 Ja inverted invJa[ 0][ 0]=1.052632 invJa[ 0][ 1]=-0.280702 invJa[ 0][ 2]=-0.280702 invJa[ 0][ 3]=-0.280702 invJa[ 1][ 0]=-0.280702 invJa[ 1][ 1]=1.052632 invJa[ 1][ 2]=-0.280702 invJa[ 1][ 3]=-0.280702 invJa[ 2][ 0]=-0.280702 invJa[ 2][ 1]=-0.280702 invJa[ 2][ 2]=1.052632 invJa[ 2][ 3]=-0.280702 invJa[ 3][ 0]=-0.280702 invJa[ 3][ 1]=-0.280702 invJa[ 3][ 2]=-0.280702 invJa[ 3][ 3]=1.052632 change X_tmp2[ 0]=0.154246 change X_tmp2[ 1]=-0.375088 change X_tmp2[ 2]=-1.729754 change X_tmp2[ 3]=-2.571088 X_next[ 0]=0.345754 X_next[ 1]=0.875088 X_next[ 2]=2.229754 X_next[ 3]=3.071088 X_next[ 4]=0.041334 X_next[ 5]=0.670123 X_next[ 6]=11.085903 X_next[ 7]=28.965209 X_next[ 8]=0.674647 X_next[ 9]=5.992401 residual X_resid[ 0]=0.0320178 residual X_resid[ 1]=0.263808 residual X_resid[ 2]=9.66359 residual X_resid[ 3]=26.9119 simeq_newton3 itr=1, prev=21.478, residual=36.8713 Ja computed Ja[ 0][ 0]=1.358638 Ja[ 0][ 1]=1.000000 Ja[ 0][ 2]=1.000000 Ja[ 0][ 3]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=3.297336 Ja[ 1][ 2]=1.000000 Ja[ 1][ 3]=1.000000 Ja[ 2][ 0]=1.000000 Ja[ 2][ 1]=1.000000 Ja[ 2][ 2]=15.915414 Ja[ 2][ 3]=1.000000 Ja[ 3][ 0]=1.000000 Ja[ 3][ 1]=1.000000 Ja[ 3][ 2]=1.000000 Ja[ 3][ 3]=29.294739 Ja inverted invJa[ 0][ 0]=0.991109 invJa[ 0][ 1]=-0.280564 invJa[ 0][ 2]=-0.043214 invJa[ 0][ 3]=-0.022780 invJa[ 1][ 0]=-0.280564 invJa[ 1][ 1]=0.391488 invJa[ 1][ 2]=-0.006746 invJa[ 1][ 3]=-0.003556 invJa[ 2][ 0]=-0.043214 invJa[ 2][ 1]=-0.006746 invJa[ 2][ 2]=0.066006 invJa[ 2][ 3]=-0.000548 invJa[ 3][ 0]=-0.022780 invJa[ 3][ 1]=-0.003556 invJa[ 3][ 2]=-0.000548 invJa[ 3][ 3]=0.035054 change X_tmp2[ 0]=-1.072931 change X_tmp2[ 1]=-0.066601 change X_tmp2[ 2]=0.619948 change X_tmp2[ 3]=0.936396 X_next[ 0]=1.418685 X_next[ 1]=0.941688 X_next[ 2]=1.609807 X_next[ 3]=2.134692 X_next[ 4]=2.855343 X_next[ 5]=0.835068 X_next[ 6]=4.171778 X_next[ 7]=9.727595 X_next[ 8]=2.150637 X_next[ 9]=3.236057 residual X_resid[ 0]=2.42922 residual X_resid[ 1]=0.0119402 residual X_resid[ 2]=2.33265 residual X_resid[ 3]=7.25747 simeq_newton3 itr=2, prev=36.8713, residual=12.0313 Ja computed Ja[ 0][ 0]=7.038004 Ja[ 0][ 1]=1.000000 Ja[ 0][ 2]=1.000000 Ja[ 0][ 3]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=3.660332 Ja[ 1][ 2]=1.000000 Ja[ 1][ 3]=1.000000 Ja[ 2][ 0]=1.000000 Ja[ 2][ 1]=1.000000 Ja[ 2][ 2]=8.774433 Ja[ 2][ 3]=1.000000 Ja[ 3][ 0]=1.000000 Ja[ 3][ 1]=1.000000 Ja[ 3][ 2]=1.000000 Ja[ 3][ 3]=14.670726 Ja inverted invJa[ 0][ 0]=0.149883 invJa[ 0][ 1]=-0.035711 invJa[ 0][ 2]=-0.012220 invJa[ 0][ 3]=-0.006949 invJa[ 1][ 0]=-0.035711 invJa[ 1][ 1]=0.294842 invJa[ 1][ 2]=-0.027735 invJa[ 1][ 3]=-0.015773 invJa[ 2][ 0]=-0.012220 invJa[ 2][ 1]=-0.027735 invJa[ 2][ 2]=0.119136 invJa[ 2][ 3]=-0.005397 invJa[ 3][ 0]=-0.006949 invJa[ 3][ 1]=-0.015773 invJa[ 3][ 2]=-0.005397 invJa[ 3][ 3]=0.070080 change X_tmp2[ 0]=0.284733 change X_tmp2[ 1]=-0.262394 change X_tmp2[ 2]=0.208717 change X_tmp2[ 3]=0.478941 X_next[ 0]=1.133952 X_next[ 1]=1.204083 X_next[ 2]=1.401090 X_next[ 3]=1.655751 X_next[ 4]=1.458090 X_next[ 5]=1.745698 X_next[ 6]=2.750414 X_next[ 7]=4.539258 X_next[ 8]=1.913010 X_next[ 9]=2.793299 residual X_resid[ 0]=0.321966 residual X_resid[ 1]=0.212574 residual X_resid[ 2]=0.20129 residual X_resid[ 3]=1.35913 simeq_newton3 itr=3, prev=12.0313, residual=2.09496 Ja computed Ja[ 0][ 0]=4.857544 Ja[ 0][ 1]=1.000000 Ja[ 0][ 2]=1.000000 Ja[ 0][ 3]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=5.349447 Ja[ 1][ 2]=1.000000 Ja[ 1][ 3]=1.000000 Ja[ 2][ 0]=1.000000 Ja[ 2][ 1]=1.000000 Ja[ 2][ 2]=6.889160 Ja[ 2][ 3]=1.000000 Ja[ 3][ 0]=1.000000 Ja[ 3][ 1]=1.000000 Ja[ 3][ 2]=1.000000 Ja[ 3][ 3]=9.224531 Ja inverted invJa[ 0][ 0]=0.221490 invJa[ 0][ 1]=-0.033474 invJa[ 0][ 2]=-0.024722 invJa[ 0][ 3]=-0.017702 invJa[ 1][ 0]=-0.033474 invJa[ 1][ 1]=0.200226 invJa[ 1][ 2]=-0.021926 invJa[ 1][ 3]=-0.015700 invJa[ 2][ 0]=-0.024722 invJa[ 2][ 1]=-0.021926 invJa[ 2][ 2]=0.153610 invJa[ 2][ 3]=-0.011595 invJa[ 3][ 0]=-0.017702 invJa[ 3][ 1]=-0.015700 invJa[ 3][ 2]=-0.011595 invJa[ 3][ 3]=0.113285 change X_tmp2[ 0]=0.035161 change X_tmp2[ 1]=0.006033 change X_tmp2[ 2]=0.002540 change X_tmp2[ 3]=0.142598 X_next[ 0]=1.098792 X_next[ 1]=1.198049 X_next[ 2]=1.398550 X_next[ 3]=1.513153 X_next[ 4]=1.326618 X_next[ 5]=1.719587 X_next[ 6]=2.735483 X_next[ 7]=3.464562 X_next[ 8]=1.841061 X_next[ 9]=2.535336 residual X_resid[ 0]=0.00416219 residual X_resid[ 1]=0.000131274 residual X_resid[ 2]=2.71001e-05 residual X_resid[ 3]=0.0981054 simeq_newton3 itr=4, prev=2.09496, residual=0.102426 Ja computed Ja[ 0][ 0]=4.622029 Ja[ 0][ 1]=1.000000 Ja[ 0][ 2]=1.000000 Ja[ 0][ 3]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=5.305968 Ja[ 1][ 2]=1.000000 Ja[ 1][ 3]=1.000000 Ja[ 2][ 0]=1.000000 Ja[ 2][ 1]=1.000000 Ja[ 2][ 2]=6.867827 Ja[ 2][ 3]=1.000000 Ja[ 3][ 0]=1.000000 Ja[ 3][ 1]=1.000000 Ja[ 3][ 2]=1.000000 Ja[ 3][ 3]=7.868893 Ja inverted invJa[ 0][ 0]=0.234306 invJa[ 0][ 1]=-0.035146 invJa[ 0][ 2]=-0.025791 invJa[ 0][ 3]=-0.022032 invJa[ 1][ 0]=-0.035146 invJa[ 1][ 1]=0.202672 invJa[ 1][ 2]=-0.021694 invJa[ 1][ 3]=-0.018533 invJa[ 2][ 0]=-0.025791 invJa[ 2][ 1]=-0.021694 invJa[ 2][ 2]=0.154501 invJa[ 2][ 3]=-0.013600 invJa[ 3][ 0]=-0.022032 invJa[ 3][ 1]=-0.018533 invJa[ 3][ 2]=-0.013600 invJa[ 3][ 3]=0.133966 change X_tmp2[ 0]=-0.001192 change X_tmp2[ 1]=-0.001938 change X_tmp2[ 2]=-0.001440 change X_tmp2[ 3]=0.013048 X_next[ 0]=1.099983 X_next[ 1]=1.199988 X_next[ 2]=1.399990 X_next[ 3]=1.500104 X_next[ 4]=1.330939 X_next[ 5]=1.727948 X_next[ 6]=2.743943 X_next[ 7]=3.375705 X_next[ 8]=1.847940 X_next[ 9]=2.520133 residual X_resid[ 0]=4.68199e-06 residual X_resid[ 1]=1.35123e-05 residual X_resid[ 2]=8.70578e-06 residual X_resid[ 3]=0.000770657 simeq_newton3 itr=5, prev=0.102426, residual=0.000797558 Ja computed Ja[ 0][ 0]=4.629889 Ja[ 0][ 1]=1.000000 Ja[ 0][ 2]=1.000000 Ja[ 0][ 3]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=5.319913 Ja[ 1][ 2]=1.000000 Ja[ 1][ 3]=1.000000 Ja[ 2][ 0]=1.000000 Ja[ 2][ 1]=1.000000 Ja[ 2][ 2]=6.879918 Ja[ 2][ 3]=1.000000 Ja[ 3][ 0]=1.000000 Ja[ 3][ 1]=1.000000 Ja[ 3][ 2]=1.000000 Ja[ 3][ 3]=7.750940 Ja inverted invJa[ 0][ 0]=0.233908 invJa[ 0][ 1]=-0.034940 invJa[ 0][ 2]=-0.025670 invJa[ 0][ 3]=-0.022358 invJa[ 1][ 0]=-0.034940 invJa[ 1][ 1]=0.202127 invJa[ 1][ 2]=-0.021570 invJa[ 1][ 3]=-0.018787 invJa[ 2][ 0]=-0.025670 invJa[ 2][ 1]=-0.021570 invJa[ 2][ 2]=0.154223 invJa[ 2][ 3]=-0.013803 invJa[ 3][ 0]=-0.022358 invJa[ 3][ 1]=-0.018787 invJa[ 3][ 2]=-0.013803 invJa[ 3][ 3]=0.136106 change X_tmp2[ 0]=-0.000017 change X_tmp2[ 1]=-0.000012 change X_tmp2[ 2]=-0.000010 change X_tmp2[ 3]=0.000104 X_next[ 0]=1.100000 X_next[ 1]=1.200000 X_next[ 2]=1.400000 X_next[ 3]=1.500000 X_next[ 4]=1.331000 X_next[ 5]=1.728000 X_next[ 6]=2.744000 X_next[ 7]=3.375000 X_next[ 8]=1.848000 X_next[ 9]=2.520000 simeq_newton3 finished Returned solution X[ 0]=1.100000, err=-9.06942e-10 X[ 1]=1.200000, err=-8.56501e-10 X[ 2]=1.400000, err=-6.51591e-10 X[ 3]=1.500000, err=6.64206e-09 third test finished Fourth test case fourth A generated solve system of equations A X = Y for X = X0, X1, X_soln[ 0]=1.000000 X_soln[ 1]=2.000000 X_soln[ 2]=1.000000 X_soln[ 3]=4.000000 Y[ 0]=4.000000 Y[ 1]=7.000000 simeq_newton3 running 2 equations, 2 variables, with 2 nonlinear terms X[ 0]=0.500000 X[ 1]=0.500000 X[ 2]=0.250000 X[ 3]=0.250000 residual X_resid[ 0]=-2.75 residual X_resid[ 1]=-5.75 simeq_newton3 itr=0, prev=1e+12, residual=8.5 Ja computed Ja[ 0][ 0]=2.000000 Ja[ 0][ 1]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=2.000000 Ja inverted invJa[ 0][ 0]=0.666667 invJa[ 0][ 1]=-0.333333 invJa[ 1][ 0]=-0.333333 invJa[ 1][ 1]=0.666667 change X_tmp2[ 0]=0.083333 change X_tmp2[ 1]=-2.916667 X_next[ 0]=0.416667 X_next[ 1]=3.416667 X_next[ 2]=0.173611 X_next[ 3]=11.673611 residual X_resid[ 0]=0.00694444 residual X_resid[ 1]=8.50694 simeq_newton3 itr=1, prev=8.5, residual=8.51389 Ja computed Ja[ 0][ 0]=1.833333 Ja[ 0][ 1]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=7.833333 Ja inverted invJa[ 0][ 0]=0.586279 invJa[ 0][ 1]=-0.074844 invJa[ 1][ 0]=-0.074844 invJa[ 1][ 1]=0.137214 change X_tmp2[ 0]=-0.632623 change X_tmp2[ 1]=1.166753 X_next[ 0]=1.049290 X_next[ 1]=2.249913 X_next[ 2]=1.101009 X_next[ 3]=5.062110 residual X_resid[ 0]=0.400212 residual X_resid[ 1]=1.36131 simeq_newton3 itr=2, prev=8.51389, residual=1.76153 Ja computed Ja[ 0][ 0]=3.098579 Ja[ 0][ 1]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=5.499827 Ja inverted invJa[ 0][ 0]=0.342847 invJa[ 0][ 1]=-0.062338 invJa[ 1][ 0]=-0.062338 invJa[ 1][ 1]=0.193158 change X_tmp2[ 0]=0.052350 change X_tmp2[ 1]=0.238001 X_next[ 0]=0.996940 X_next[ 1]=2.011913 X_next[ 2]=0.993888 X_next[ 3]=4.047792 residual X_resid[ 0]=0.00274054 residual X_resid[ 1]=0.0566444 simeq_newton3 itr=3, prev=1.76153, residual=0.0593849 Ja computed Ja[ 0][ 0]=2.993879 Ja[ 0][ 1]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=5.023825 Ja inverted invJa[ 0][ 0]=0.357804 invJa[ 0][ 1]=-0.071221 invJa[ 1][ 0]=-0.071221 invJa[ 1][ 1]=0.213228 change X_tmp2[ 0]=-0.003054 change X_tmp2[ 1]=0.011883 X_next[ 0]=0.999993 X_next[ 1]=2.000030 X_next[ 2]=0.999986 X_next[ 3]=4.000118 residual X_resid[ 0]=9.32518e-06 residual X_resid[ 1]=0.000141206 simeq_newton3 itr=4, prev=0.0593849, residual=0.000150531 Ja computed Ja[ 0][ 0]=2.999986 Ja[ 0][ 1]=1.000000 Ja[ 1][ 0]=1.000000 Ja[ 1][ 1]=5.000059 Ja inverted invJa[ 0][ 0]=0.357144 invJa[ 0][ 1]=-0.071428 invJa[ 1][ 0]=-0.071428 invJa[ 1][ 1]=0.214283 change X_tmp2[ 0]=-0.000007 change X_tmp2[ 1]=0.000030 X_next[ 0]=1.000000 X_next[ 1]=2.000000 X_next[ 2]=1.000000 X_next[ 3]=4.000000 residual X_resid[ 0]=4.56382e-11 residual X_resid[ 1]=8.7568e-10 simeq_newton3 itr=5, prev=0.000150531, residual=9.21318e-10 simeq_newton3 finished Returned solution X[ 0]=1.000000, err=-4.6249e-11 X[ 1]=2.000000, err=1.84386e-10 fourth test finished test_simeq_newton3.c finished