##############ECHO OF PROBLEM################# ##############temp/mtest6_revpostode.ode################# diff (x2,t,2) = 3.0 * diff(x2,t,1) - 2.0 * x2 - diff(x1,t,2) - diff (x1,t,1) + x1; diff (x1,t,1) = 4.0 * x2 - 2.0 * diff (x2,t ,1) - 2.0 * x1; ! #BEGIN FIRST INPUT BLOCK # Digits:=32; ELIMINATED in preodein.rb max_terms=30 ! #END FIRST INPUT BLOCK #BEGIN SECOND INPUT BLOCK ## ## was complicatedrev.ode ## t_start=0.5 t_end=5.0 $array_x1_init[0 + 1] = exact_soln_x1(t_start) $array_x1_init[1 + 1] = exact_soln_x1p(t_start) $array_x2_init[0 + 1] = exact_soln_x2(t_start) $array_x2_init[1 + 1] = exact_soln_x2p(t_start) $glob_h=0.00001 $glob_look_poles=true $glob_max_iter=10 #END SECOND INPUT BLOCK #BEGIN OVERRIDE BLOCK $glob_desired_digits_correct=10 $glob_display_interval=0.001 $glob_look_poles=true $glob_max_iter=10000000 $glob_max_minutes=3 #END OVERRIDE BLOCK ! #BEGIN USER DEF BLOCK def exact_soln_x1 (t) c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return(1.0 + 2.0 * c1 + 6.0 * c3 * exp(-t)) end def exact_soln_x1p (t) c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return( - 6.0 * c3 * exp(-t)) end def exact_soln_x2 (t) c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return(1.0 + c1 + c2 * exp(2.0 * t) + c3 * exp(-t)) end def exact_soln_x2p (t) c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return( 2.0 * c2 * exp(2.0 * t) - c3 * exp(-t)) end #END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Optimize min_size 0.0 min_size 1.0 opt_iter 1 $glob_desired_digits_correct 10 desired_abs_gbl_error 1.0e-10 range 4.5 estimated_steps 4500.0 step_error 2.2222222222222224e-14 est_needed_step_err 2.2222222222222224e-14 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 value3 8.112591994867362e+48 value3 1.2266152631992113e+52 max_value3 1.2266152631992113e+52 value3 1.2266152631992113e+52 best_h 0.0 opt_iter 2 $glob_desired_digits_correct 10 desired_abs_gbl_error 1.0e-10 range 4.5 estimated_steps 9000.0 step_error 1.1111111111111112e-14 est_needed_step_err 1.1111111111111112e-14 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 value3 3.0223946222800794e+40 value3 4.569804843345226e+43 max_value3 4.569804843345226e+43 value3 4.569804843345226e+43 best_h 0.0 opt_iter 3 $glob_desired_digits_correct 10 desired_abs_gbl_error 1.0e-10 range 4.5 estimated_steps 18000.0 step_error 5.555555555555556e-15 est_needed_step_err 5.555555555555556e-15 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 value3 1.1260926465501067e+32 value3 1.7026136812708273e+35 max_value3 1.7026136812708273e+35 value3 1.7026136812708273e+35 best_h 0.0 opt_iter 4 $glob_desired_digits_correct 10 desired_abs_gbl_error 1.0e-10 range 4.5 estimated_steps 36000.0 step_error 2.777777777777778e-15 est_needed_step_err 2.777777777777778e-15 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 value3 419623642279730300000000.0 value3 634443512637305100000000000.0 max_value3 634443512637305100000000000.0 value3 634443512637305100000000000.0 best_h 0.0 opt_iter 5 $glob_desired_digits_correct 10 desired_abs_gbl_error 1.0e-10 range 4.5 estimated_steps 72000.0 step_error 1.388888888888889e-15 est_needed_step_err 1.388888888888889e-15 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 value3 1564125639234969.2 value3 2364757743377717000.0 max_value3 2364757743377717000.0 value3 2364757743377717000.0 best_h 0.0 opt_iter 6 $glob_desired_digits_correct 10 desired_abs_gbl_error 1.0e-10 range 4.5 estimated_steps 144000.0 step_error 6.944444444444445e-16 est_needed_step_err 6.944444444444445e-16 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 value3 5833584.610357152 value3 8818900726.11272 max_value3 8818900726.11272 value3 8818900726.11272 best_h 0.0 opt_iter 7 $glob_desired_digits_correct 10 desired_abs_gbl_error 1.0e-10 range 4.5 estimated_steps 288000.0 step_error 3.4722222222222225e-16 est_needed_step_err 3.4722222222222225e-16 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 value3 0.021782391645970806 value3 32.923957420665886 max_value3 32.923957420665886 value3 32.923957420665886 best_h 0.0 opt_iter 8 $glob_desired_digits_correct 10 desired_abs_gbl_error 1.0e-10 range 4.5 estimated_steps 576000.0 step_error 1.7361111111111112e-16 est_needed_step_err 1.7361111111111112e-16 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 value3 8.15258367778097e-11 value3 1.2318431186418186e-07 max_value3 1.2318431186418186e-07 value3 1.2318431186418186e-07 best_h 0.0 opt_iter 9 $glob_desired_digits_correct 10 desired_abs_gbl_error 1.0e-10 range 4.5 estimated_steps 1152000.0 step_error 8.680555555555556e-17 est_needed_step_err 8.680555555555556e-17 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 value3 3.065861785112164e-19 value3 4.629293774222638e-16 max_value3 4.629293774222638e-16 value3 4.629293774222638e-16 best_h 0.0 opt_iter 10 $glob_desired_digits_correct 10 desired_abs_gbl_error 1.0e-10 range 4.5 estimated_steps 2304000.0 step_error 4.340277777777778e-17 est_needed_step_err 4.340277777777778e-17 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 value3 1.1642679562470447e-27 value3 1.7554934852987215e-24 max_value3 1.7554934852987215e-24 value3 1.7554934852987215e-24 best_h 1.953125e-06 START of Soultion t[1] 0.5 x2[1] (analytic) 2.000725615563606 x2[1] (numeric) 2.000725615563606 absolute error 0.0 relative error 0.0% Correct digits 16 h 1.953125e-06 x1[1] (analytic) 3.001091755187483 x1[1] (numeric) 3.001091755187483 absolute error 0.0 relative error 0.0% Correct digits 16 h 1.953125e-06 t[1] 0.5 x2[1] (analytic) 2.000725615563606 x2[1] (numeric) 2.000725615563606 absolute error 0.0 relative error 0.0% Correct digits 16 h 1.953125e-06 x1[1] (analytic) 3.001091755187483 x1[1] (numeric) 3.001091755187483 absolute error 0.0 relative error 0.0% Correct digits 16 h 1.953125e-06 TOP MAIN SOLVE Loop NO POLE NO POLE t[1] 0.5010019531249763 x2[1] (analytic) 2.000726523869006 x2[1] (numeric) 2.0007260221262637 absolute error 5.017427424647281e-07 relative error 2.5078027230551113e-05% Correct digits 6 h 1.953125e-06 x1[1] (analytic) 3.00109066184779 x1[1] (numeric) 3.001091163768893 absolute error 5.019211029022586e-07 relative error 1.672462312728742e-05% Correct digits 6 h 1.953125e-06 TOP MAIN SOLVE Loop NO POLE NO POLE t[1] 0.5020019531249527 x2[1] (analytic) 2.000727432766899 x2[1] (numeric) 2.000725428369249 absolute error 2.004397650079426e-06 relative error 0.00010018344414398572% Correct digits 5 h 1.953125e-06 x1[1] (analytic) 3.0010895717310913 x1[1] (numeric) 3.001091572836983 absolute error 2.0011058916580282e-06 relative error 6.667931242397901e-05% Correct digits 6 h 1.953125e-06 TOP MAIN SOLVE Loop NO POLE NO POLE t[1] 0.503001953124929 x2[1] (analytic) 2.0007283440297376 x2[1] (numeric) 2.000723834063017 absolute error 4.50996672052284e-06 relative error 0.00022541624573774752% Correct digits 5 h 1.953125e-06 x1[1] (analytic) 3.0010884827039646 x1[1] (numeric) 3.0010929762663126 absolute error 4.493562347995805e-06 relative error 0.0001497310850344249% Correct digits 5 h 1.953125e-06 TOP MAIN SOLVE Loop NO POLE NO POLE t[1] 0.5040019531249054 x2[1] (analytic) 2.000729257661712 x2[1] (numeric) 2.00072123721366 absolute error 8.020448052370455e-06 relative error 0.0004008762315868812% Correct digits 5 h 1.953125e-06 x1[1] (analytic) 3.0010873947653205 x1[1] (numeric) 3.001095370059875 absolute error 7.975294554540824e-06 relative error 0.00026574682791483573% Correct digits 5 h 1.953125e-06 TOP MAIN SOLVE Loop NO POLE NO POLE t[1] 0.5050019531248817 x2[1] (analytic) 2.00073017366702 x2[1] (numeric) 2.0007176358292735 absolute error 1.2537837746684488e-05 relative error 0.0006266631008870439% Correct digits 5 h 1.953125e-06 x1[1] (analytic) 3.0010863079140715 x1[1] (numeric) 3.0010987502187 absolute error 1.2442304628379475e-05 relative error 0.00041459336226246676% Correct digits 5 h 1.953125e-06 TOP MAIN SOLVE Loop NO POLE NO POLE t[1] 0.5060019531248581 x2[1] (analytic) 2.00073109204987 x2[1] (numeric) 2.0007130279199643 absolute error 1.8064129905681625e-05 relative error 0.0009028764523859042% Correct digits 5 h 1.953125e-06 x1[1] (analytic) 3.0010852221491304 x1[1] (numeric) 3.001103112741822 absolute error 1.789059269174942e-05 relative error 0.0005961374425394574% Correct digits 5 h 1.953125e-06 TOP MAIN SOLVE Loop NO POLE NO POLE t[1] 0.5070019531248344 x2[1] (analytic) 2.000732012814477 x2[1] (numeric) 2.0007074114978454 absolute error 2.460131663184484e-05 relative error 0.001229615784336733% Correct digits 4 h 1.953125e-06 x1[1] (analytic) 3.0010841374694115 x1[1] (numeric) 3.001108453626325 absolute error 2.431615691333988e-05 relative error 0.0008102457578494906% Correct digits 5 h 1.953125e-06 TOP MAIN SOLVE Loop NO POLE NO POLE Finished! Maximum Time Reached before Solution Completed! diff (x2,t,2) = 3.0 * diff(x2,t,1) - 2.0 * x2 - diff(x1,t,2) - diff (x1,t,1) + x1; diff (x1,t,1) = 4.0 * x2 - 2.0 * diff (x2,t ,1) - 2.0 * x1; Iterations 3738 Total Elapsed Time 3 Minutes 0 Seconds Elapsed Time(since restart) 2 Minutes 59 Seconds Expected Time Remaining 1 Days 6 Hours 45 Minutes 37.35758486409031 Seconds Optimized Time Remaining 1 Days 6 Hours 35 Minutes 22.150042725959793 Seconds Expected Total Time 1 Days 6 Hours 38 Minutes 22.150042725959793 Seconds Time to Timeout 0.0 Seconds Percent Done 0.16228298610727362%