updated default dee_rel_tol

doc/shared/sundomeigest/SUNDomEigEst_Power.rst

@@ -85,7 +85,7 @@ routines:
                      need too tight tolerances and consider a safety factor,
                      and (2) an early (less costly) termination will be a good
                      indicator whether the power iteration is compatible.
-   :param rel_tol: relative tolerance for convergence checks (default 0.01). A
+   :param rel_tol: relative tolerance for convergence checks (default 0.005). A
                    value :math:`\leq 0` will result in the default value. The
                    default has been found to small enough for many internal
                    applications.
@@ -149,7 +149,7 @@ information:
 
 * ``num_ATimes`` - number of calls to the ``ATimes`` function,
 
-* ``rel_tol`` - relative tolerance for the convergence criteria (default is 0.01),
+* ``rel_tol`` - relative tolerance for the convergence criteria (default is 0.005),
 
 * ``res`` - the residual from the last :c:func:`SUNDomEigEstimator_Estimate`
   call.

examples/arkode/CXX_lapack/ark_heat2D_lsrk_domeigest.cpp

@@ -602,7 +602,7 @@ static int InitUserData(UserData* udata)
   udata->dee_num_succ_wups = 5;
   udata->dee_max_iters     = 100;
   udata->dee_krylov_dim    = 3;
-  udata->dee_reltol        = 0.01;
+  udata->dee_reltol        = 0.005;
 
   // Timing variables
   udata->timing     = false;

examples/arkode/CXX_lapack/ark_heat2D_lsrk_domeigest.out

@@ -29,7 +29,7 @@
  dee num_succ_wups = 5
  dee_max_iters     = 100
  dee_krylov_dim    = 3
- dee_reltol        = 0.01
+ dee_reltol        = 0.005
  ------------------------------------ 
 
           t                     ||u||_rms                max error      

examples/arkode/C_serial/ark_analytic_lsrk_domeigest.c

@@ -125,7 +125,7 @@ int main(int argc, char* argv[])
   SUNDomEigEstimator DEE = NULL; /* domeig estimator object */
   sunindextype max_iters = 100;  /* max number of power iterations (PI)*/
   sunindextype numwarmup = 10;   /* number of preprocessing warmups */
-  sunrealtype rel_tol    = SUN_RCONST(1.0e-2); /* relative error for PI*/
+  sunrealtype rel_tol    = SUN_RCONST(5.0e-3); /* relative error for PI*/
   N_Vector q             = NULL;               /* random initial eigenvector */
 
   /* Create the SUNDIALS context object for this simulation */

src/sundomeigest/power/sundomeigest_power.c

@@ -31,7 +31,7 @@
 #define DEE_NUM_OF_WARMUPS_PI_DEFAULT 100
 
 /* Default Power Iteration parameters */
-#define DEE_TOL_DEFAULT      SUN_RCONST(0.01)
+#define DEE_TOL_DEFAULT      SUN_RCONST(0.005)
 #define DEE_MAX_ITER_DEFAULT 100
 
 /*