test_simeq_newton3.java running first A generated 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.1 X_soln[1]=1.2 X_soln[2]=1.4 X_soln[3]=1.5 Y[0]=9.21809726972151 Y[1]=8.202810561712287 Y[2]=8.331973328814726 Y[3]=6.572916901000221 simeq_newton3 running n=4, nlin=6 initial guess at solution, can be bad X[0]=1.0 X[1]=1.0 X[2]=1.0 X[3]=1.0 X[4]=1.0 X[5]=1.0 X[6]=1.0 X[7]=1.0 X[8]=1.0 X[9]=1.0 residual X_resid[0]=-3.1928605395392866 residual X_resid[1]=-2.9298350666235926 residual X_resid[2]=-3.3567807621676975 residual X_resid[3]=-2.347267621874928 simeq_newton3 itr 0, prev=1.0E12 residual=11.826743990205504 Ja computed Ja inverted change X_tmp2[0]=-0.07725619886769461 change X_tmp2[1]=-0.2697017769542099 change X_tmp2[2]=-0.4844714276764934 change X_tmp2[3]=-0.5097329549050893 X_next[0]=1.0772561988676945 X_next[1]=1.2697017769542098 X_next[2]=1.4844714276764934 X_next[3]=1.5097329549050893 X_next[4]=1.3677941099372495 X_next[5]=1.599156047506479 X_next[6]=1.626369184306349 X_next[7]=1.8848360095585965 X_next[8]=1.9169106155693219 X_next[9]=2.241155434978209 residual X_resid[0]=0.45357718665657565 residual X_resid[1]=0.34156906043627266 residual X_resid[2]=0.507896152066543 residual X_resid[3]=0.21704980591392875 simeq_newton3 itr 1, prev=11.826743990205504 residual=1.52009220507332 Ja computed Ja inverted change X_tmp2[0]=-0.020744071377467345 change X_tmp2[1]=0.06585912798075273 change X_tmp2[2]=0.0901942253706941 change X_tmp2[3]=0.004992516302355532 X_next[0]=1.098000270245162 X_next[1]=1.203842648973457 X_next[2]=1.3942772023057994 X_next[3]=1.5047404386027337 X_next[4]=1.3218195539055073 X_next[5]=1.530916744928436 X_next[6]=1.652205408234625 X_next[7]=1.6784903606271142 X_next[8]=1.8114707156249965 X_next[9]=2.0980252889314213 residual X_resid[0]=0.0029453695275662994 residual X_resid[1]=-0.0011244865117099323 residual X_resid[2]=0.005458630766581152 residual X_resid[3]=0.0011997114489679106 simeq_newton3 itr 2, prev=1.52009220507332 residual=0.010728198254825294 Ja computed Ja inverted change X_tmp2[0]=-0.0020022112521770937 change X_tmp2[1]=0.0038562771607696346 change X_tmp2[2]=-0.005729653439163246 change X_tmp2[3]=0.004744628070252938 X_next[0]=1.100002481497339 X_next[1]=1.1999863718126875 X_next[2]=1.4000068557449628 X_next[3]=1.4999958105324807 X_next[4]=1.3199879867569448 X_next[5]=1.5400110154327462 X_next[6]=1.649999113821341 X_next[7]=1.6799891473382864 X_next[8]=1.799974530415103 X_next[9]=2.100004418334195 residual X_resid[0]=-3.2354978351634145E-5 residual X_resid[1]=-1.096533857491977E-5 residual X_resid[2]=-2.6421856158975743E-5 residual X_resid[3]=-1.1341734215797317E-5 simeq_newton3 itr 3, prev=0.010728198254825294 residual=8.108390730132697E-5 Ja computed Ja inverted change X_tmp2[0]=2.4814886355210255E-6 change X_tmp2[1]=-1.3628136390542786E-5 change X_tmp2[2]=6.855683970522746E-6 change X_tmp2[3]=-4.189429929215116E-6 X_next[0]=1.1000000000087036 X_next[1]=1.199999999949078 X_next[2]=1.4000000000609922 X_next[3]=1.4999999999624098 X_next[4]=1.31999999995443 X_next[5]=1.5400000000792764 X_next[6]=1.649999999971706 X_next[7]=1.6800000000019 X_next[8]=1.7999999998785088 X_next[9]=2.1000000000388623 residual X_resid[0]=-7.692513293022785E-11 residual X_resid[1]=-2.6147972675971687E-12 residual X_resid[2]=-5.171862937913829E-11 residual X_resid[3]=-3.348077370901592E-11 simeq_newton3 itr 4, prev=8.108390730132697E-5 residual=1.6473933328597923E-10 simeq_newton3.java finished Returned solution X[0]=1.1000000000087036 err=8.703482379246452E-12 X[1]=1.199999999949078 err=-5.092193333666728E-11 X[2]=1.4000000000609922 err=6.099232230383222E-11 X[3]=1.4999999999624098 err=-3.7590153212363475E-11 first test finished second A generated 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.1 X_soln[1]=1.2 X_soln[2]=1.4 X_soln[3]=1.5 X_soln[4]=1.2100000000000002 X_soln[5]=1.44 X_soln[6]=1.9599999999999997 X_soln[7]=2.25 X_soln[8]=1.32 X_soln[9]=2.0999999999999996 Y[0]=8.985761062264652 Y[1]=8.073947421304396 Y[2]=8.068557284299333 Y[3]=7.057936619755721 simeq_newton3 running n=4, nlin=6 initial guess at solution, can be bad X[0]=1.0 X[1]=1.0 X[2]=1.0 X[3]=1.0 X[4]=1.0 X[5]=1.0 X[6]=1.0 X[7]=1.0 X[8]=1.0 X[9]=1.0 residual X_resid[0]=-3.3297334390281748 residual X_resid[1]=-2.66282216827402 residual X_resid[2]=-2.9519104939109404 residual X_resid[3]=-2.6247816047357926 simeq_newton3 itr 0, prev=1.0E12 residual=11.569247705948927 Ja computed Ja inverted change X_tmp2[0]=-0.0882915592236817 change X_tmp2[1]=-0.2007136205607838 change X_tmp2[2]=-0.45730345384182414 change X_tmp2[3]=-0.6390929465805466 X_next[0]=1.0882915592236817 X_next[1]=1.2007136205607838 X_next[2]=1.4573034538418241 X_next[3]=1.6390929465805466 X_next[4]=1.1843785178775124 X_next[5]=1.4417131986001859 X_next[6]=2.1237333565793097 X_next[7]=2.6866256875300984 X_next[8]=1.3067264983012075 X_next[9]=2.388655812219603 residual X_resid[0]=0.6346152615217964 residual X_resid[1]=0.4707274650410618 residual X_resid[2]=0.5878599394550434 residual X_resid[3]=0.5513484542800633 simeq_newton3 itr 1, prev=11.569247705948927 residual=2.244551120297965 Ja computed Ja inverted change X_tmp2[0]=-0.01137167905604182 change X_tmp2[1]=9.420557644631011E-4 change X_tmp2[2]=0.05624582230507605 change X_tmp2[3]=0.1338508774760005 X_next[0]=1.0996632382797236 X_next[1]=1.1997715647963207 X_next[2]=1.401057631536748 X_next[3]=1.505242069104546 X_next[4]=1.209259237623848 X_next[5]=1.439451807693812 X_next[6]=1.9629624868873619 X_next[7]=2.265753686602135 X_next[8]=1.3193446841398533 X_next[9]=2.1089308882290894 residual X_resid[0]=0.01860144472939851 residual X_resid[1]=0.014497322208287144 residual X_resid[2]=0.01651324614037719 residual X_resid[3]=0.017817199629472924 simeq_newton3 itr 2, prev=2.244551120297965 residual=0.06742921270753577 Ja computed Ja inverted change X_tmp2[0]=-3.3647434736902484E-4 change X_tmp2[1]=-2.2997726288035217E-4 change X_tmp2[2]=0.0010567173882211345 change X_tmp2[3]=0.0052369365460719595 X_next[0]=1.0999997126270926 X_next[1]=1.2000015420592012 X_next[2]=1.4000009141485268 X_next[3]=1.500005132558474 X_next[4]=1.2099993677796863 X_next[5]=1.4400037009444608 X_next[6]=1.9600025596167108 X_next[7]=2.250015397701765 X_next[8]=1.3200013514171893 X_next[9]=2.1000085568093456 residual X_resid[0]=2.2998860055523096E-5 residual X_resid[1]=1.857964916851529E-5 residual X_resid[2]=1.8036562773104947E-5 residual X_resid[3]=2.266445078635826E-5 simeq_newton3 itr 3, prev=0.06742921270753577 residual=8.227952278350159E-5 Ja computed Ja inverted change X_tmp2[0]=-2.873724440651923E-7 change X_tmp2[1]=1.5420573772922703E-6 change X_tmp2[2]=9.141479396965531E-7 change X_tmp2[3]=5.132553246119616E-6 X_next[0]=1.0999999999995367 X_next[1]=1.2000000000018238 X_next[2]=1.400000000000587 X_next[3]=1.5000000000052278 X_next[4]=1.2099999999989808 X_next[5]=1.4400000000043771 X_next[6]=1.9600000000016435 X_next[7]=2.2500000000156835 X_next[8]=1.3200000000014502 X_next[9]=2.1000000000081993 residual X_resid[0]=2.2575719071937783E-11 residual X_resid[1]=1.8319568084734783E-11 residual X_resid[2]=1.6971313243629993E-11 residual X_resid[3]=2.3105961588498758E-11 simeq_newton3 itr 4, prev=8.227952278350159E-5 residual=8.097256198880132E-11 simeq_newton3.java finished Returned solution X[0]=1.0999999999995367 err=-4.634070904785403E-13 X[1]=1.2000000000018238 err=1.823874384854207E-12 X[2]=1.400000000000587 err=5.870859354217828E-13 X[3]=1.5000000000052278 err=5.227818178354937E-12 second test finished third A generated 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.1 X_soln[1]=1.2 X_soln[2]=1.4 X_soln[3]=1.5 Y[0]=2.762036571991004 Y[1]=1.832466842414773 Y[2]=2.30950278030278 Y[3]=2.7556783133198763 simeq_newton3 running n=4, nlin=6 initial guess at solution, can be bad X[0]=0.5 X[1]=0.5 X[2]=0.5 X[3]=0.5 X[4]=0.125 X[5]=0.125 X[6]=0.125 X[7]=0.125 X[8]=0.125 X[9]=0.125 residual X_resid[0]=-1.4238369939042248 residual X_resid[1]=-0.7421062342399456 residual X_resid[2]=-0.9546966109209363 residual X_resid[3]=-1.3891239474988688 simeq_newton3 itr 0, prev=1.0E12 residual=4.509763786563976 Ja computed Ja inverted change X_tmp2[0]=-1.1471628069496207 change X_tmp2[1]=-1.246139427494921 change X_tmp2[2]=1.3515949922032149 change X_tmp2[3]=-0.25559761094179256 X_next[0]=1.6471628069496207 X_next[1]=1.746139427494921 X_next[2]=-0.8515949922032149 X_next[3]=0.7555976109417926 X_next[4]=4.468992048761381 X_next[5]=5.323984178676888 X_next[6]=-0.6175886368584559 X_next[7]=0.4313916436641835 X_next[8]=-2.4493370107788666 X_next[9]=-1.1235763616853134 residual X_resid[0]=-0.03222660366706753 residual X_resid[1]=3.204487355506015 residual X_resid[2]=1.4370404518508 residual X_resid[3]=1.2336286399025607 simeq_newton3 itr 1, prev=4.509763786563976 residual=5.907383050926443 Ja computed Ja inverted change X_tmp2[0]=1.0575053915811068 change X_tmp2[1]=0.8056732738825778 change X_tmp2[2]=-2.947905989342368 change X_tmp2[3]=2.8943513929184785 X_next[0]=0.5896574153685139 X_next[1]=0.9404661536123431 X_next[2]=2.0963109971391534 X_next[3]=-2.138753781976686 X_next[4]=0.20502144656381852 X_next[5]=0.8318202928806105 X_next[6]=9.212280177023587 X_next[7]=-9.783232428526945 X_next[8]=1.1625152198809625 X_next[9]=-4.21657348542289 residual X_resid[0]=-5.428298889437429 residual X_resid[1]=-0.8615161467103447 residual X_resid[2]=-9.02418746382526 residual X_resid[3]=-6.106539904205595 simeq_newton3 itr 2, prev=5.907383050926443 residual=21.42054240417863 Ja computed Ja inverted change X_tmp2[0]=1.2495154473403725 change X_tmp2[1]=-0.8554441730881157 change X_tmp2[2]=0.8737572212342676 change X_tmp2[3]=-0.9877840233774071 X_next[0]=-0.6598580319718585 X_next[1]=1.7959103267004588 X_next[2]=1.2225537759048857 X_next[3]=-1.1509697585992789 X_next[4]=-0.287310516084707 X_next[5]=5.792338624436488 X_next[6]=1.827275006508861 X_next[7]=-1.5247257626441366 X_next[8]=-1.4487822831524468 X_next[9]=-2.52706569278223 residual X_resid[0]=-1.9563947359343346 residual X_resid[1]=1.5088988409556248 residual X_resid[2]=-3.164887732704405 residual X_resid[3]=-0.22949743451901927 simeq_newton3 itr 3, prev=21.42054240417863 residual=6.859678744113384 Ja computed Ja inverted change X_tmp2[0]=0.5466386137708379 change X_tmp2[1]=0.09104901052203795 change X_tmp2[2]=1.365134050066624 change X_tmp2[3]=-1.5148528797219574 X_next[0]=-1.2064966457426964 X_next[1]=1.704861316178421 X_next[2]=-0.14258027416173835 X_next[3]=0.3638831211226785 X_next[4]=-1.7562177268697778 X_next[5]=4.955268251365886 X_next[6]=-0.002898533581892476 X_next[7]=0.04818210096466838 X_next[8]=0.2932747146507955 X_next[9]=-0.08845256129808882 residual X_resid[0]=-0.7328005768787618 residual X_resid[1]=1.439661867244481 residual X_resid[2]=-1.2472347029749409 residual X_resid[3]=1.5735709085942382 simeq_newton3 itr 4, prev=6.859678744113384 residual=4.993268055692422 Ja computed Ja inverted change X_tmp2[0]=0.11747875773166139 change X_tmp2[1]=0.5704252261464596 change X_tmp2[2]=-0.3397655408151503 change X_tmp2[3]=-1.2210032132803512 X_next[0]=-1.3239754034743578 X_next[1]=1.1344360900319614 X_next[2]=0.19718526665341196 X_next[3]=1.5848863344030297 X_next[4]=-2.3208108750455994 X_next[5]=1.4599571288360595 X_next[6]=0.007666963332341962 X_next[7]=3.981020025769458 X_next[8]=-0.2961654636811628 X_next[9]=0.35452969509754395 residual X_resid[0]=1.5408609764720218 residual X_resid[1]=2.380618051420166 residual X_resid[2]=1.3331258806256976 residual X_resid[3]=0.954744270314845 simeq_newton3 itr 5, prev=4.993268055692422 residual=6.20934917883273 Ja computed Ja inverted change X_tmp2[0]=0.1280193791357997 change X_tmp2[1]=0.23316910797365112 change X_tmp2[2]=-1.1302488613218846 change X_tmp2[3]=0.42907819852387563 X_next[0]=-1.4519947826101576 X_next[1]=0.9012669820583104 X_next[2]=1.3274341279752966 X_next[3]=1.155808135879154 X_next[4]=-3.061224408586952 X_next[5]=0.7320831025930563 X_next[6]=2.339046938993573 X_next[7]=1.544035358861017 X_next[8]=-1.7371267012409863 X_next[9]=1.3827771272965474 simeq_newton3.java finished Returned solution X[0]=-1.4519947826101576 err=-2.5519947826101577 X[1]=0.9012669820583104 err=-0.2987330179416896 X[2]=1.3274341279752966 err=-0.07256587202470333 X[3]=1.155808135879154 err=-0.3441918641208459 third test finished Fourth test case fourth A generated solve system of equations A X = Y for X = X0, X1, X0*X0, X1*X1, X_soln[0]=1.0 X_soln[1]=2.0 X_soln[2]=1.0 X_soln[3]=4.0 Y[0]=4.0 Y[1]=7.0 simeq_newton3 running n=2, nlin=2 initial guess at solution, can be bad X[0]=0.5 X[1]=0.5 X[2]=0.25 X[3]=0.25 residual X_resid[0]=-2.75 residual X_resid[1]=-5.75 simeq_newton3 itr 0, prev=1.0E12 residual=8.5 Ja computed Ja inverted change X_tmp2[0]=0.08333333333333326 change X_tmp2[1]=-2.9166666666666665 X_next[0]=0.41666666666666674 X_next[1]=3.4166666666666665 X_next[2]=0.17361111111111116 X_next[3]=11.67361111111111 residual X_resid[0]=0.006944444444444642 residual X_resid[1]=8.506944444444443 simeq_newton3 itr 1, prev=8.5 residual=8.513888888888888 Ja computed Ja inverted change X_tmp2[0]=-0.6326230076230074 change X_tmp2[1]=1.1667532917532917 X_next[0]=1.0492896742896742 X_next[1]=2.249913374913375 X_next[2]=1.1010088205709307 X_next[3]=5.0621101946140925 residual X_resid[0]=0.4002118697739796 residual X_resid[1]=1.3613132438171416 simeq_newton3 itr 2, prev=8.513888888888888 residual=1.7615251135911212 Ja computed Ja inverted change X_tmp2[0]=0.05235014632184856 change X_tmp2[1]=0.23800078748599252 X_next[0]=0.9969395279678257 X_next[1]=2.0119125874273824 X_next[2]=0.9938884224247111 X_next[3]=4.0477922594487445 residual X_resid[0]=0.002740537819919453 residual X_resid[1]=0.05664437484395268 simeq_newton3 itr 3, prev=1.7615251135911212 residual=0.059384912663872136 Ja computed Ja inverted change X_tmp2[0]=-0.0030537163884401383 change X_tmp2[1]=0.011882995358037841 X_next[0]=0.9999932443562658 X_next[1]=2.0000295920693447 X_next[2]=0.9999864887581703 X_next[3]=4.0001183691530695 residual X_resid[0]=9.325183780717339E-6 residual X_resid[1]=1.4120557868047E-4 simeq_newton3 itr 4, prev=0.059384912663872136 residual=1.5053076246118735E-4 Ja computed Ja inverted change X_tmp2[0]=-6.755597485054146E-6 change X_tmp2[1]=2.9591884959060132E-5 X_next[0]=0.9999999999537509 X_next[1]=2.000000000184386 X_next[2]=0.9999999999075018 X_next[3]=4.000000000737543 residual X_resid[0]=4.5638159917871235E-11 residual X_resid[1]=8.756799729781051E-10 simeq_newton3 itr 5, prev=1.5053076246118735E-4 residual=9.213181328959763E-10 simeq_newton3.java finished Returned solution X[0]=0.9999999999537509 err=-4.624911564832246E-11 X[1]=2.000000000184386 err=1.84385839929746E-10 fourth test finished test_simeq_newton3 finished