##############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; 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 double exact_soln_x1 (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return(1.0 + 2.0 * c1 + 6.0 * c3 * exp(-t)); } double exact_soln_x1p (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return( - 6.0 * c3 * exp(-t)); } double exact_soln_x2 (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return(1.0 + c1 + c2 * exp(2.0 * t) + c3 * exp(-t)); } double exact_soln_x2p (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return( 2.0 * c2 * exp(2.0 * t) - c3 * exp(-t)); } // END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Optimize min_size = 0 min_size = 1 opt_iter = 1 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 4.5 estimated_steps = 4500 step_error = 2.222222222222222e-14 est_needed_step_err = 2.222222222222222e-14 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 8.11259199486736e+48 value3 = 1.226615263199211e+52 max_value3 = 1.226615263199211e+52 value3 = 1.226615263199211e+52 best_h = 0 opt_iter = 2 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 4.5 estimated_steps = 9000 step_error = 1.111111111111111e-14 est_needed_step_err = 1.111111111111111e-14 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 3.022394622280079e+40 value3 = 4.569804843345226e+43 max_value3 = 4.569804843345226e+43 value3 = 4.569804843345226e+43 best_h = 0 opt_iter = 3 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 4.5 estimated_steps = 18000 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.126092646550107e+32 value3 = 1.702613681270827e+35 max_value3 = 1.702613681270827e+35 value3 = 1.702613681270827e+35 best_h = 0 opt_iter = 4 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 4.5 estimated_steps = 36000 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 = 4.196236422797302e+23 value3 = 6.344435126373051e+26 max_value3 = 6.344435126373051e+26 value3 = 6.344435126373051e+26 best_h = 0 opt_iter = 5 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 4.5 estimated_steps = 72000 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 value3 = 2.364757743377717e+18 max_value3 = 2.364757743377717e+18 value3 = 2.364757743377717e+18 best_h = 0 opt_iter = 6 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 4.5 estimated_steps = 144000 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.610357151 value3 = 8818900726.11272 max_value3 = 8818900726.11272 value3 = 8818900726.11272 best_h = 0 opt_iter = 7 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 4.5 estimated_steps = 288000 step_error = 3.472222222222222e-16 est_needed_step_err = 3.472222222222222e-16 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 0.0217823916459708 value3 = 32.92395742066589 max_value3 = 32.92395742066589 value3 = 32.92395742066589 best_h = 0 opt_iter = 8 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 4.5 estimated_steps = 576000 step_error = 1.736111111111111e-16 est_needed_step_err = 1.736111111111111e-16 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 8.152583677780969e-11 value3 = 1.231843118641819e-07 max_value3 = 1.231843118641819e-07 value3 = 1.231843118641819e-07 best_h = 0 opt_iter = 9 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 4.5 estimated_steps = 1152000 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.065861785112163e-19 value3 = 4.629293774222638e-16 max_value3 = 4.629293774222638e-16 value3 = 4.629293774222638e-16 best_h = 0 opt_iter = 10 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 4.5 estimated_steps = 2304000 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.164267956247045e-27 value3 = 1.755493485298721e-24 max_value3 = 1.755493485298721e-24 value3 = 1.755493485298721e-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 relative error = 0 % Correct digits = 16 h = 1.953125e-06 x1[1] (analytic) = 3.001091755187483 x1[1] (numeric) = 3.001091755187483 absolute error = 0 relative error = 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 relative error = 0 % Correct digits = 16 h = 1.953125e-06 x1[1] (analytic) = 3.001091755187483 x1[1] (numeric) = 3.001091755187483 absolute error = 0 relative error = 0 % Correct digits = 16 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 = 123 Total Elapsed Time = 3 Minutes 2 Seconds Elapsed Time(since restart) = 2 Minutes 47 Seconds Expected Time Remaining = 39 Days 3 Hours 18 Minutes 15 Seconds Optimized Time Remaining = 35 Days 21 Hours 53 Minutes 20 Seconds Expected Total Time = 35 Days 21 Hours 56 Minutes 22 Seconds Time to Timeout Unknown Percent Done = 0.005382 %