diff --git a/test/unit_tests/logging/CMakeLists.txt b/test/unit_tests/logging/CMakeLists.txt
new file mode 100644
index 0000000000..eb00894b9f
--- /dev/null
+++ b/test/unit_tests/logging/CMakeLists.txt
@@ -0,0 +1,111 @@
+# ------------------------------------------------------------------------------
+# Programmer(s): David J. Gardner @ LLNL
+# ------------------------------------------------------------------------------
+# SUNDIALS Copyright Start
+# Copyright (c) 2002-2024, Lawrence Livermore National Security
+# and Southern Methodist University.
+# All rights reserved.
+#
+# See the top-level LICENSE and NOTICE files for details.
+#
+# SPDX-License-Identifier: BSD-3-Clause
+# SUNDIALS Copyright End
+# ------------------------------------------------------------------------------
+
+# List of test tuples of the form "name\;args"
+set(unit_tests)
+
+# if(BUILD_ARKODE)
+#   list(APPEND unit_tests "test_logging_arkode.cpp\;")
+# endif()
+
+if(BUILD_CVODE)
+  list(APPEND unit_tests "test_logging_cvode.cpp\;1 1") # Newton + Dense
+  list(APPEND unit_tests "test_logging_cvode.cpp\;1 0") # Newton + GMRES
+  list(APPEND unit_tests "test_logging_cvode.cpp\;0")   # Fixed-point
+endif()
+
+# if(BUILD_CVODES)
+#   list(APPEND unit_tests "test_logging_cvodes.cpp\;")
+# endif()
+
+# if(BUILD_IDA)
+#   list(APPEND unit_tests "test_logging_ida.cpp\;")
+# endif()
+
+# if(BUILD_IDAS)
+#   list(APPEND unit_tests "test_logging_idas.cpp\;")
+# endif()
+
+# if(BUILD_KINSOL)
+#   list(APPEND unit_tests "test_logging_kinsol.cpp\;")
+# endif()
+
+# Add the build and install targets for each test
+foreach(test_tuple ${unit_tests})
+
+  # parse the test tuple
+  list(GET test_tuple 0 test_src)
+  list(GET test_tuple 1 test_args)
+
+  # extract the file name without extension
+  get_filename_component(test_target ${test_src} NAME_WE)
+
+  string(REGEX MATCH "arkode|cvode|cvodes|ida|idas|kinsol" _pkg ${test_target})
+
+  # check if this test has already been added, only need to add test source
+  # files once for testing with different inputs
+  if(NOT TARGET ${test_target})
+
+    # test source files
+    add_executable(${test_target} ${test_src})
+
+    set_target_properties(${test_target} PROPERTIES FOLDER "unit_tests")
+
+    # include location of public and private header files
+    target_include_directories(
+      ${test_target} PRIVATE
+      ${CMAKE_SOURCE_DIR}/include
+      ${CMAKE_SOURCE_DIR}/test/unit_tests)
+
+    # libraries to link against
+    target_link_libraries(${test_target} sundials_${_pkg} sundials_nvecserial
+      ${EXE_EXTRA_LINK_LIBS})
+
+  endif()
+
+  # Set the test name
+  if(${SUNDIALS_LOGGING_LEVEL} GREATER 2)
+    set(test_name "${test_target}_lvl${SUNDIALS_LOGGING_LEVEL}")
+  else()
+    set(test_name "${test_target}")
+  endif()
+
+  if("${test_args}" STREQUAL "")
+    set(test_name ${test_name})
+  else()
+    string(REPLACE " " "_" test_name "${test_name}_${test_args}")
+    string(REPLACE " " ";" test_args "${test_args}")
+  endif()
+
+  # For now, only diff results with double precision
+  # TODO(DJG): Update Jenkins diff handling for testing single/extended
+  if(SUNDIALS_PRECISION MATCHES "DOUBLE")
+    set(diff_output "")
+  else()
+    set(diff_output "NODIFF")
+  endif()
+
+  # add test to regression tests
+  sundials_add_test(
+    ${test_name} ${test_target}
+    TEST_ARGS ${test_args}
+    ANSWER_DIR ${CMAKE_CURRENT_SOURCE_DIR}
+    ANSWER_FILE ${test_name}.out
+    EXAMPLE_TYPE "develop"
+    LABELS "UnitTest" "Logging"
+    ${diff_output})
+
+endforeach()
+
+message(STATUS "Added logging units tests")
diff --git a/test/unit_tests/logging/test_logging_cvode.cpp b/test/unit_tests/logging/test_logging_cvode.cpp
new file mode 100644
index 0000000000..7a36829068
--- /dev/null
+++ b/test/unit_tests/logging/test_logging_cvode.cpp
@@ -0,0 +1,198 @@
+/* -----------------------------------------------------------------------------
+ * Programmer(s): David J. Gardner @ LLNL
+ * -----------------------------------------------------------------------------
+ * SUNDIALS Copyright Start
+ * Copyright (c) 2002-2024, Lawrence Livermore National Security
+ * and Southern Methodist University.
+ * All rights reserved.
+ *
+ * See the top-level LICENSE and NOTICE files for details.
+ *
+ * SPDX-License-Identifier: BSD-3-Clause
+ * SUNDIALS Copyright End
+ * -----------------------------------------------------------------------------
+ * Test logging output
+ * ---------------------------------------------------------------------------*/
+
+#include <cmath>
+#include <cstdio>
+#include <iomanip>
+#include <iostream>
+#include <limits>
+
+// Include desired integrators, vectors, linear solvers, and nonlinear solvers
+#include "cvode/cvode.h"
+#include "nvector/nvector_serial.h"
+#include "sundials/sundials_context.hpp"
+#include "sundials/sundials_iterative.h"
+#include "sundials/sundials_logger.h"
+#include "sundials/sundials_matrix.h"
+#include "sundials/sundials_nonlinearsolver.h"
+#include "sunlinsol/sunlinsol_dense.h"
+#include "sunmatrix/sunmatrix_dense.h"
+#include "sunlinsol/sunlinsol_spgmr.h"
+#include "sunnonlinsol/sunnonlinsol_fixedpoint.h"
+
+#include "problems/kpr.hpp"
+#include "utilities/check_return.hpp"
+
+using namespace std;
+
+int main(int argc, char* argv[])
+{
+  cout << "Start CVODE Logging test" << endl;
+
+  // SUNDIALS context object for this simulation
+  sundials::Context sunctx;
+
+  // Use Newton (1) otherwise use fixed-point
+  bool newton = true;
+  if (argc > 1) { newton = stoi(argv[1]); }
+
+  // Use direct dense solver (1) otherwise use GMRES
+  bool direct = true;
+  if (argc > 2) { direct = stoi(argv[2]); }
+
+  // Ensure logging output goes to stdout
+  SUNLogger logger;
+  int flag = SUNContext_GetLogger(sunctx, &logger);
+  if (check_flag(flag, "SUNContext_GetLogger")) { return 1; }
+
+  SUNLogger_SetErrorFilename(logger, "stdout");
+  SUNLogger_SetWarningFilename(logger, "stdout");
+  SUNLogger_SetInfoFilename(logger, "stdout");
+  SUNLogger_SetDebugFilename(logger, "stdout");
+
+  // Create initial condition
+  N_Vector y = N_VNew_Serial(2, sunctx);
+  if (check_ptr(y, "N_VNew_Serial")) { return 1; }
+
+  sunrealtype utrue, vtrue;
+  flag = kpr_true_sol(ZERO, &utrue, &vtrue);
+  if (check_flag(flag, "true_sol")) { return 1; }
+
+  sunrealtype* ydata = N_VGetArrayPointer(y);
+  ydata[0]           = utrue;
+  ydata[1]           = vtrue;
+
+  // Create CVODE memory structure
+  void* cvode_mem = CVodeCreate(CV_BDF, sunctx);
+  if (check_ptr(cvode_mem, "CVodeCreate")) { return 1; }
+
+  flag = CVodeInit(cvode_mem, kpr_rhs, ZERO, y);
+  if (check_flag(flag, "CVodeInit")) { return 1; }
+
+  flag = CVodeSetUserData(cvode_mem, &kpr_udata);
+  if (check_flag(flag, "CVodeSetUserData")) { return 1; }
+
+  // Relative and absolute tolerances
+  const sunrealtype rtol = SUN_RCONST(1.0e-6);
+  const sunrealtype atol = SUN_RCONST(1.0e-10);
+
+  flag = CVodeSStolerances(cvode_mem, rtol, atol);
+  if (check_flag(flag, "CVodeSStolerances")) { return 1; }
+
+  SUNNonlinearSolver NLS = nullptr;
+
+  if (!newton)
+  {
+    cout << "Using fixed-point nonlinear solver" << endl;
+
+    NLS = SUNNonlinSol_FixedPoint(y, 0, sunctx);
+    if (check_ptr(NLS, "SUNNonlinSol_FixedPoint")) { return 1; }
+
+    flag = CVodeSetNonlinearSolver(cvode_mem, NLS);
+    if (check_flag(flag, "CVodeSetLinearSolver")) { return 1; }
+  }
+  else
+  {
+    cout << "Using Newton nonlinear solver" << endl;
+  }
+
+  SUNMatrix A        = nullptr;
+  SUNLinearSolver LS = nullptr;
+
+  if (newton && direct)
+  {
+    cout << "Using dense direct linear solver" << endl;
+
+    A = SUNDenseMatrix(2, 2, sunctx);
+    if (check_ptr(A, "SUNDenseMatrix")) { return 1; }
+
+    LS = SUNLinSol_Dense(y, A, sunctx);
+    if (check_ptr(LS, "SUNLinSol_Dense")) { return 1; }
+
+    flag = CVodeSetLinearSolver(cvode_mem, LS, A);
+    if (check_flag(flag, "CVodeSetLinearSolver")) { return 1; }
+
+    flag = CVodeSetJacFn(cvode_mem, kpr_rhs_jac);
+    if (check_flag(flag, "CVodeSetJacFn")) { return 1; }
+  }
+  else if (newton)
+  {
+    cout << "Using GMRES iterative linear solver" << endl;
+
+    LS = SUNLinSol_SPGMR(y, SUN_PREC_NONE, 0, sunctx);
+    if (check_ptr(LS, "SUNLinSol_SPGMR")) { return 1; }
+
+    flag = CVodeSetLinearSolver(cvode_mem, LS, A);
+    if (check_flag(flag, "CVodeSetLinearSolver")) { return 1; }
+  }
+
+  // Initial time and fist output time
+  const sunrealtype dtout = ONE; // output interval
+  const int nout          = 3;   // number of outputs
+  sunrealtype tret        = ZERO;
+  sunrealtype tout        = tret + dtout;
+
+  // Output initial contion
+  cout << scientific;
+  cout << setprecision(numeric_limits<sunrealtype>::digits10);
+  cout << "           t              ";
+  cout << "          u              ";
+  cout << "          v              ";
+  cout << "        u err            ";
+  cout << "        v err            " << endl;
+  for (int i = 0; i < 9; i++) { cout << "--------------"; }
+  cout << endl;
+
+  cout << setw(22) << tret << setw(25) << ydata[0] << setw(25) << ydata[1]
+       << setw(25) << abs(ydata[0] - utrue) << setw(25) << abs(ydata[1] - vtrue)
+       << endl;
+
+  // Advance in time
+  for (int i = 0; i < nout; i++)
+  {
+    flag = CVode(cvode_mem, tout, y, &tret, CV_ONE_STEP);
+    if (check_flag(flag, "CVode")) { return 1; }
+
+    flag = kpr_true_sol(tret, &utrue, &vtrue);
+    if (check_flag(flag, "true_sol")) { return 1; }
+
+    cout << setw(22) << tret << setw(25) << ydata[0] << setw(25) << ydata[1]
+         << setw(25) << abs(ydata[0] - utrue) << setw(25)
+         << abs(ydata[1] - vtrue) << endl;
+
+    // update output time
+    tout += dtout;
+  }
+  for (int i = 0; i < 9; i++) { cout << "--------------"; }
+  cout << endl;
+
+  // Print some final statistics
+  flag = CVodePrintAllStats(cvode_mem, stdout, SUN_OUTPUTFORMAT_TABLE);
+  if (check_flag(flag, "CVodePrintAllStats")) { return 1; }
+
+  // Clean up and return with successful completion
+  N_VDestroy(y);
+  SUNMatDestroy(A);
+  SUNLinSolFree(LS);
+  SUNNonlinSolFree(NLS);
+  CVodeFree(&cvode_mem);
+
+  cout << "End CVODE Logging test" << endl;
+
+  return 0;
+}
+
+/*---- end of file ----*/
diff --git a/test/unit_tests/logging/test_logging_cvode_lvl3_0.out b/test/unit_tests/logging/test_logging_cvode_lvl3_0.out
new file mode 100644
index 0000000000..8f5875afd4
--- /dev/null
+++ b/test/unit_tests/logging/test_logging_cvode_lvl3_0.out
@@ -0,0 +1,54 @@
+Start CVODE Logging test
+Using fixed-point nonlinear solver
+           t                        u                        v                      u err                    v err            
+------------------------------------------------------------------------------------------------------------------------------
+ 0.000000000000000e+00    1.224744871391589e+00    1.732050807568877e+00    0.000000000000000e+00    0.000000000000000e+00
+[INFO][rank 0][cvStep][begin-step-attempt] step = 1, t_n = 0, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499621899552651
+ 1.029860256095084e-04    1.224744869195305e+00    1.732049582942475e+00    1.113801051388918e-09    6.122818017040288e-07
+[INFO][rank 0][cvStep][begin-step-attempt] step = 2, t_n = 0.0001029860256095084, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.249962018079078
+ 2.059720512190168e-04    1.224744864802767e+00    1.732047133691379e+00    2.258890852147033e-09    1.224500395968775e-06
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.001029860256095084, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4064516129032258
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = failed error test, dsm = 8.198472955140451, eflag = 5
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.0002810714755410129, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4217683315129505
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.5885921138418445
+ 4.870435267600296e-04    1.224744841922327e+00    1.732034367982330e+00    5.258975033228808e-09    2.744234303353466e-06
+------------------------------------------------------------------------------------------------------------------------------
+Current time                 = 0.0004870435267600296
+Steps                        = 3
+Error test fails             = 1
+NLS step fails               = 0
+Initial step size            = 0.0001029860256095084
+Last step size               = 0.0002810714755410129
+Current step size            = 0.0002810714755410129
+Last method order            = 2
+Current method order         = 2
+Stab. lim. order reductions  = 0
+RHS fn evals                 = 11
+NLS iters                    = 8
+NLS fails                    = 0
+NLS iters per step           = 2.666666666666667
+LS setups                    = 0
+Root fn evals                = 0
+End CVODE Logging test
diff --git a/test/unit_tests/logging/test_logging_cvode_lvl3_1_0.out b/test/unit_tests/logging/test_logging_cvode_lvl3_1_0.out
new file mode 100644
index 0000000000..d816801faa
--- /dev/null
+++ b/test/unit_tests/logging/test_logging_cvode_lvl3_1_0.out
@@ -0,0 +1,122 @@
+Start CVODE Logging test
+Using Newton nonlinear solver
+Using GMRES iterative linear solver
+           t                        u                        v                      u err                    v err            
+------------------------------------------------------------------------------------------------------------------------------
+ 0.000000000000000e+00    1.224744871391589e+00    1.732050807568877e+00    0.000000000000000e+00    0.000000000000000e+00
+[INFO][rank 0][cvStep][begin-step-attempt] step = 1, t_n = 0, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.499975381317856, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 0.7070719651123717, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 5.130293423419247e-05
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.4999243811465705
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 3.627663976220471e-05, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499621905732853
+ 1.029860256095084e-04    1.224744869258150e+00    1.732049582942250e+00    1.050956877080012e-09    6.122820268572582e-07
+[INFO][rank 0][cvStep][begin-step-attempt] step = 2, t_n = 0.0001029860256095084, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.4999746801363614, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 0.7070709734919923, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 5.129626861394919e-05
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.4999240374084951
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 3.627192645756874e-05, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499620187042475
+ 2.059720512190168e-04    1.224744864928434e+00    1.732047133690926e+00    2.133223819811292e-09    1.224500849161814e-06
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.001029860256095084, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4064516129032258
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 33.34377627547665, b-tol = 0.02032258064516129, res-tol = 0.02874046917080807
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 47.15522062951332, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 0.02280194342842776
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 33.32282742612932
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.01612377930633392, b-tol = 0.02032258064516129, res-tol = 0.02874046917080807
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = failed error test, dsm = 8.19847341436515, eflag = 5
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.0002810714702930874, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4217683318751525
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 2.482943969856862, b-tol = 0.02108841659375763, res-tol = 0.02982352475586586
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 3.511413036784067, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 0.0004635068923188323
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 2.482495053589632
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.0003277492125626657, b-tol = 0.02108841659375763, res-tol = 0.02982352475586586
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.5885920933306285
+ 4.870435215121043e-04    1.224744842672905e+00    1.732034367979995e+00    4.508398321334539e-09    2.744236933471811e-06
+------------------------------------------------------------------------------------------------------------------------------
+Current time                 = 0.0004870435215121043
+Steps                        = 3
+Error test fails             = 1
+NLS step fails               = 0
+Initial step size            = 0.0001029860256095084
+Last step size               = 0.0002810714702930874
+Current step size            = 0.0002810714702930874
+Last method order            = 2
+Current method order         = 2
+Stab. lim. order reductions  = 0
+RHS fn evals                 = 11
+NLS iters                    = 8
+NLS fails                    = 0
+NLS iters per step           = 2.666666666666667
+LS setups                    = 0
+Jac fn evals                 = 0
+LS RHS fn evals              = 4
+Prec setup evals             = 0
+Prec solves                  = 0
+LS iters                     = 4
+LS fails                     = 0
+Jac-times setups             = 0
+Jac-times evals              = 4
+LS iters per NLS iter        = 0.5
+Jac evals per NLS iter       = 0
+Prec evals per NLS iter      = 0
+Root fn evals                = 0
+End CVODE Logging test
diff --git a/test/unit_tests/logging/test_logging_cvode_lvl3_1_1.out b/test/unit_tests/logging/test_logging_cvode_lvl3_1_1.out
new file mode 100644
index 0000000000..6c8b51edf4
--- /dev/null
+++ b/test/unit_tests/logging/test_logging_cvode_lvl3_1_1.out
@@ -0,0 +1,93 @@
+Start CVODE Logging test
+Using Newton nonlinear solver
+Using dense direct linear solver
+           t                        u                        v                      u err                    v err            
+------------------------------------------------------------------------------------------------------------------------------
+ 0.000000000000000e+00    1.224744871391589e+00    1.732050807568877e+00    0.000000000000000e+00    0.000000000000000e+00
+[INFO][rank 0][cvStep][begin-step-attempt] step = 1, t_n = 0, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.4999243864411202
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 2.207987840210415e-11
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499621932116032
+ 1.029860256095084e-04    1.224744869195325e+00    1.732049582942459e+00    1.113781733508290e-09    6.122818176912403e-07
+[INFO][rank 0][cvStep][begin-step-attempt] step = 2, t_n = 0.0001029860256095084, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.499923682675401
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 1
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499618413377005
+ 2.059720512190168e-04    1.224744864802806e+00    1.732047133692213e+00    2.258852216385776e-09    1.224499561969239e-06
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.001029860256095084, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4064516129032258
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 33.32111802218621
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0.001726792321093216
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = failed error test, dsm = 8.198477691401965, eflag = 5
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.0002810714214160092, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4217683352485581
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 2.482483008202045
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 1.190779226370788e-05
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.5885920559914323
+ 4.870434726350260e-04    1.224744841922706e+00    1.732034367987314e+00    5.258602220337139e-09    2.744232363793842e-06
+------------------------------------------------------------------------------------------------------------------------------
+Current time                 = 0.000487043472635026
+Steps                        = 3
+Error test fails             = 1
+NLS step fails               = 0
+Initial step size            = 0.0001029860256095084
+Last step size               = 0.0002810714214160092
+Current step size            = 0.0002810714214160092
+Last method order            = 2
+Current method order         = 2
+Stab. lim. order reductions  = 0
+RHS fn evals                 = 10
+NLS iters                    = 7
+NLS fails                    = 0
+NLS iters per step           = 2.333333333333333
+LS setups                    = 3
+Jac fn evals                 = 1
+LS RHS fn evals              = 0
+Prec setup evals             = 0
+Prec solves                  = 0
+LS iters                     = 0
+LS fails                     = 0
+Jac-times setups             = 0
+Jac-times evals              = 0
+LS iters per NLS iter        = 0
+Jac evals per NLS iter       = 0.1428571428571428
+Prec evals per NLS iter      = 0
+Root fn evals                = 0
+End CVODE Logging test
diff --git a/test/unit_tests/logging/test_logging_cvode_lvl4_0.out b/test/unit_tests/logging/test_logging_cvode_lvl4_0.out
new file mode 100644
index 0000000000..85729f3727
--- /dev/null
+++ b/test/unit_tests/logging/test_logging_cvode_lvl4_0.out
@@ -0,0 +1,81 @@
+Start CVODE Logging test
+Using fixed-point nonlinear solver
+           t                        u                        v                      u err                    v err            
+------------------------------------------------------------------------------------------------------------------------------
+ 0.000000000000000e+00    1.224744871391589e+00    1.732050807568877e+00    0.000000000000000e+00    0.000000000000000e+00
+[INFO][rank 0][cvStep][begin-step-attempt] step = 1, t_n = 0, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 1, update-norm = 0.499975381317856
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 2, update-norm = 6.259026232096693e-05
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 0, h = 0.0001029860256095084, dsm = 0.2499621899552651
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499621899552651
+[DEBUG][rank 0][cvCompleteStep][return] nst = 1, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0001029860256095084, qprime = 1, qwait = 1
+ 1.029860256095084e-04    1.224744869195305e+00    1.732049582942475e+00    1.113801051388918e-09    6.122818017040288e-07
+[INFO][rank 0][cvStep][begin-step-attempt] step = 2, t_n = 0.0001029860256095084, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 1, update-norm = 0.4999746839964076
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 2, update-norm = 6.230194580154422e-05
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 1, h = 0.0001029860256095084, dsm = 0.249962018079078
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.249962018079078
+[DEBUG][rank 0][cvCompleteStep][return] nst = 2, nscon = 2
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 10, hprime = 0.001029860256095084, qprime = 2, qwait = 2
+ 2.059720512190168e-04    1.224744864802767e+00    1.732047133691379e+00    2.258890852147033e-09    1.224500395968775e-06
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.001029860256095084, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4064516129032258
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 1, update-norm = 33.34377648386529
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 2, update-norm = 0.02644643520889491
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.001029860256095084, dsm = 8.198472955140451
+[DEBUG][rank 0][cvDoErrorTest][new-step-eta] eta = 0.2729219560396963
+[INFO][rank 0][cvStep][end-step-attempt] status = failed error test, dsm = 8.198472955140451, eflag = 5
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.0002810714755410129, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4217683315129505
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 1, update-norm = 2.482944084048132
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 2, update-norm = 0.0005559126276672459
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.0002810714755410129, dsm = 0.5885921138418445
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.5885921138418445
+[DEBUG][rank 0][cvCompleteStep][return] nst = 3, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0002810714755410129, qprime = 2, qwait = 2
+ 4.870435267600296e-04    1.224744841922327e+00    1.732034367982330e+00    5.258975033228808e-09    2.744234303353466e-06
+------------------------------------------------------------------------------------------------------------------------------
+Current time                 = 0.0004870435267600296
+Steps                        = 3
+Error test fails             = 1
+NLS step fails               = 0
+Initial step size            = 0.0001029860256095084
+Last step size               = 0.0002810714755410129
+Current step size            = 0.0002810714755410129
+Last method order            = 2
+Current method order         = 2
+Stab. lim. order reductions  = 0
+RHS fn evals                 = 11
+NLS iters                    = 8
+NLS fails                    = 0
+NLS iters per step           = 2.666666666666667
+LS setups                    = 0
+Root fn evals                = 0
+End CVODE Logging test
diff --git a/test/unit_tests/logging/test_logging_cvode_lvl4_1_0.out b/test/unit_tests/logging/test_logging_cvode_lvl4_1_0.out
new file mode 100644
index 0000000000..59b12844fd
--- /dev/null
+++ b/test/unit_tests/logging/test_logging_cvode_lvl4_1_0.out
@@ -0,0 +1,133 @@
+Start CVODE Logging test
+Using Newton nonlinear solver
+Using GMRES iterative linear solver
+           t                        u                        v                      u err                    v err            
+------------------------------------------------------------------------------------------------------------------------------
+ 0.000000000000000e+00    1.224744871391589e+00    1.732050807568877e+00    0.000000000000000e+00    0.000000000000000e+00
+[INFO][rank 0][cvStep][begin-step-attempt] step = 1, t_n = 0, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.499975381317856, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 0.7070719651123717, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 5.130293423419247e-05
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.4999243811465705
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 3.627663976220471e-05, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 0, h = 0.0001029860256095084, dsm = 0.2499621905732853
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499621905732853
+[DEBUG][rank 0][cvCompleteStep][return] nst = 1, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0001029860256095084, qprime = 1, qwait = 1
+ 1.029860256095084e-04    1.224744869258150e+00    1.732049582942250e+00    1.050956877080012e-09    6.122820268572582e-07
+[INFO][rank 0][cvStep][begin-step-attempt] step = 2, t_n = 0.0001029860256095084, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.4999746801363614, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 0.7070709734919923, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 5.129626861394919e-05
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.4999240374084951
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 3.627192645756874e-05, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 1, h = 0.0001029860256095084, dsm = 0.2499620187042475
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499620187042475
+[DEBUG][rank 0][cvCompleteStep][return] nst = 2, nscon = 2
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 10, hprime = 0.001029860256095084, qprime = 2, qwait = 2
+ 2.059720512190168e-04    1.224744864928434e+00    1.732047133690926e+00    2.133223819811292e-09    1.224500849161814e-06
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.001029860256095084, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4064516129032258
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 33.34377627547665, b-tol = 0.02032258064516129, res-tol = 0.02874046917080807
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 47.15522062951332, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 0.02280194342842776
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 33.32282742612932
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.01612377930633392, b-tol = 0.02032258064516129, res-tol = 0.02874046917080807
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.001029860256095084, dsm = 8.19847341436515
+[DEBUG][rank 0][cvDoErrorTest][new-step-eta] eta = 0.2729219509439317
+[INFO][rank 0][cvStep][end-step-attempt] status = failed error test, dsm = 8.19847341436515, eflag = 5
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.0002810714702930874, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4217683318751525
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 2.482943969856862, b-tol = 0.02108841659375763, res-tol = 0.02982352475586586
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 3.511413036784067, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 0.0004635068923188323
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 2.482495053589632
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.0003277492125626657, b-tol = 0.02108841659375763, res-tol = 0.02982352475586586
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.0002810714702930874, dsm = 0.5885920933306285
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.5885920933306285
+[DEBUG][rank 0][cvCompleteStep][return] nst = 3, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0002810714702930874, qprime = 2, qwait = 2
+ 4.870435215121043e-04    1.224744842672905e+00    1.732034367979995e+00    4.508398321334539e-09    2.744236933471811e-06
+------------------------------------------------------------------------------------------------------------------------------
+Current time                 = 0.0004870435215121043
+Steps                        = 3
+Error test fails             = 1
+NLS step fails               = 0
+Initial step size            = 0.0001029860256095084
+Last step size               = 0.0002810714702930874
+Current step size            = 0.0002810714702930874
+Last method order            = 2
+Current method order         = 2
+Stab. lim. order reductions  = 0
+RHS fn evals                 = 11
+NLS iters                    = 8
+NLS fails                    = 0
+NLS iters per step           = 2.666666666666667
+LS setups                    = 0
+Jac fn evals                 = 0
+LS RHS fn evals              = 4
+Prec setup evals             = 0
+Prec solves                  = 0
+LS iters                     = 4
+LS fails                     = 0
+Jac-times setups             = 0
+Jac-times evals              = 4
+LS iters per NLS iter        = 0.5
+Jac evals per NLS iter       = 0
+Prec evals per NLS iter      = 0
+Root fn evals                = 0
+End CVODE Logging test
diff --git a/test/unit_tests/logging/test_logging_cvode_lvl4_1_1.out b/test/unit_tests/logging/test_logging_cvode_lvl4_1_1.out
new file mode 100644
index 0000000000..46e9c3743e
--- /dev/null
+++ b/test/unit_tests/logging/test_logging_cvode_lvl4_1_1.out
@@ -0,0 +1,104 @@
+Start CVODE Logging test
+Using Newton nonlinear solver
+Using dense direct linear solver
+           t                        u                        v                      u err                    v err            
+------------------------------------------------------------------------------------------------------------------------------
+ 0.000000000000000e+00    1.224744871391589e+00    1.732050807568877e+00    0.000000000000000e+00    0.000000000000000e+00
+[INFO][rank 0][cvStep][begin-step-attempt] step = 1, t_n = 0, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.4999243864411202
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 2.207987840210415e-11
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 0, h = 0.0001029860256095084, dsm = 0.2499621932116032
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499621932116032
+[DEBUG][rank 0][cvCompleteStep][return] nst = 1, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0001029860256095084, qprime = 1, qwait = 1
+ 1.029860256095084e-04    1.224744869195325e+00    1.732049582942459e+00    1.113781733508290e-09    6.122818176912403e-07
+[INFO][rank 0][cvStep][begin-step-attempt] step = 2, t_n = 0.0001029860256095084, h = 0.0001029860256095084, q = 1
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.499923682675401
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 1
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 1, h = 0.0001029860256095084, dsm = 0.2499618413377005
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499618413377005
+[DEBUG][rank 0][cvCompleteStep][return] nst = 2, nscon = 2
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 10, hprime = 0.001029860256095084, qprime = 2, qwait = 2
+ 2.059720512190168e-04    1.224744864802806e+00    1.732047133692213e+00    2.258852216385776e-09    1.224499561969239e-06
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.001029860256095084, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4064516129032258
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 33.32111802218621
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0.001726792321093216
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.001029860256095084, dsm = 8.198477691401965
+[DEBUG][rank 0][cvDoErrorTest][new-step-eta] eta = 0.2729219034840186
+[INFO][rank 0][cvStep][end-step-attempt] status = failed error test, dsm = 8.198477691401965, eflag = 5
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.0002810714214160092, q = 2
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4217683352485581
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 2.482483008202045
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 1.190779226370788e-05
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.0002810714214160092, dsm = 0.5885920559914323
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.5885920559914323
+[DEBUG][rank 0][cvCompleteStep][return] nst = 3, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0002810714214160092, qprime = 2, qwait = 2
+ 4.870434726350260e-04    1.224744841922706e+00    1.732034367987314e+00    5.258602220337139e-09    2.744232363793842e-06
+------------------------------------------------------------------------------------------------------------------------------
+Current time                 = 0.000487043472635026
+Steps                        = 3
+Error test fails             = 1
+NLS step fails               = 0
+Initial step size            = 0.0001029860256095084
+Last step size               = 0.0002810714214160092
+Current step size            = 0.0002810714214160092
+Last method order            = 2
+Current method order         = 2
+Stab. lim. order reductions  = 0
+RHS fn evals                 = 10
+NLS iters                    = 7
+NLS fails                    = 0
+NLS iters per step           = 2.333333333333333
+LS setups                    = 3
+Jac fn evals                 = 1
+LS RHS fn evals              = 0
+Prec setup evals             = 0
+Prec solves                  = 0
+LS iters                     = 0
+LS fails                     = 0
+Jac-times setups             = 0
+Jac-times evals              = 0
+LS iters per NLS iter        = 0
+Jac evals per NLS iter       = 0.1428571428571428
+Prec evals per NLS iter      = 0
+Root fn evals                = 0
+End CVODE Logging test
diff --git a/test/unit_tests/logging/test_logging_cvode_lvl5_0.out b/test/unit_tests/logging/test_logging_cvode_lvl5_0.out
new file mode 100644
index 0000000000..aa827efba8
--- /dev/null
+++ b/test/unit_tests/logging/test_logging_cvode_lvl5_0.out
@@ -0,0 +1,97 @@
+Start CVODE Logging test
+Using fixed-point nonlinear solver
+           t                        u                        v                      u err                    v err            
+------------------------------------------------------------------------------------------------------------------------------
+ 0.000000000000000e+00    1.224744871391589e+00    1.732050807568877e+00    0.000000000000000e+00    0.000000000000000e+00
+[INFO][rank 0][cvStep][begin-step-attempt] step = 1, t_n = 0, h = 0.0001029860256095084, q = 1
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448713915889e+00
+1.7320508075688772e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 1, update-norm = 0.499975381317856
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 2, update-norm = 6.259026232096693e-05
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 0, h = 0.0001029860256095084, dsm = 0.2499621899552651
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499621899552651
+[DEBUG][rank 0][cvCompleteStep][return] nst = 1, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0001029860256095084, qprime = 1, qwait = 1
+ 1.029860256095084e-04    1.224744869195305e+00    1.732049582942475e+00    1.113801051388918e-09    6.122818017040288e-07
+[INFO][rank 0][cvStep][begin-step-attempt] step = 2, t_n = 0.0001029860256095084, h = 0.0001029860256095084, q = 1
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448669990217e+00
+1.7320483583160731e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 1, update-norm = 0.4999746839964076
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 2, update-norm = 6.230194580154422e-05
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 1, h = 0.0001029860256095084, dsm = 0.249962018079078
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.249962018079078
+[DEBUG][rank 0][cvCompleteStep][return] nst = 2, nscon = 2
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 10, hprime = 0.001029860256095084, qprime = 2, qwait = 2
+ 2.059720512190168e-04    1.224744864802767e+00    1.732047133691379e+00    2.258890852147033e-09    1.224500395968775e-06
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.001029860256095084, q = 2
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448208773850e+00
+1.7320226411804154e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4064516129032258
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 1, update-norm = 33.34377648386529
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 2, update-norm = 0.02644643520889491
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.001029860256095084, dsm = 8.198472955140451
+[DEBUG][rank 0][cvDoErrorTest][new-step-eta] eta = 0.2729219560396963
+[INFO][rank 0][cvStep][end-step-attempt] status = failed error test, dsm = 8.198472955140451, eflag = 5
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.0002810714755410129, q = 2
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448528145659e+00
+1.7320404491473784e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4217683315129505
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-solver] solver = Fixed-Point
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 1, update-norm = 2.482944084048132
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = continue
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][begin-iteration] 
+[DEBUG][rank 0][SUNNonlinSolSolve_FixedPoint][nonlinear-iterate] cur-iter = 2, update-norm = 0.0005559126276672459
+[INFO][rank 0][SUNNonlinSolSolve_FixedPoint][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.0002810714755410129, dsm = 0.5885921138418445
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.5885921138418445
+[DEBUG][rank 0][cvCompleteStep][return] nst = 3, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0002810714755410129, qprime = 2, qwait = 2
+ 4.870435267600296e-04    1.224744841922327e+00    1.732034367982330e+00    5.258975033228808e-09    2.744234303353466e-06
+------------------------------------------------------------------------------------------------------------------------------
+Current time                 = 0.0004870435267600296
+Steps                        = 3
+Error test fails             = 1
+NLS step fails               = 0
+Initial step size            = 0.0001029860256095084
+Last step size               = 0.0002810714755410129
+Current step size            = 0.0002810714755410129
+Last method order            = 2
+Current method order         = 2
+Stab. lim. order reductions  = 0
+RHS fn evals                 = 11
+NLS iters                    = 8
+NLS fails                    = 0
+NLS iters per step           = 2.666666666666667
+LS setups                    = 0
+Root fn evals                = 0
+End CVODE Logging test
diff --git a/test/unit_tests/logging/test_logging_cvode_lvl5_1_0.out b/test/unit_tests/logging/test_logging_cvode_lvl5_1_0.out
new file mode 100644
index 0000000000..cd864d6015
--- /dev/null
+++ b/test/unit_tests/logging/test_logging_cvode_lvl5_1_0.out
@@ -0,0 +1,149 @@
+Start CVODE Logging test
+Using Newton nonlinear solver
+Using GMRES iterative linear solver
+           t                        u                        v                      u err                    v err            
+------------------------------------------------------------------------------------------------------------------------------
+ 0.000000000000000e+00    1.224744871391589e+00    1.732050807568877e+00    0.000000000000000e+00    0.000000000000000e+00
+[INFO][rank 0][cvStep][begin-step-attempt] step = 1, t_n = 0, h = 0.0001029860256095084, q = 1
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448713915889e+00
+1.7320508075688772e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.499975381317856, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 0.7070719651123717, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 5.130293423419247e-05
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.4999243811465705
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 3.627663976220471e-05, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 0, h = 0.0001029860256095084, dsm = 0.2499621905732853
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499621905732853
+[DEBUG][rank 0][cvCompleteStep][return] nst = 1, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0001029860256095084, qprime = 1, qwait = 1
+ 1.029860256095084e-04    1.224744869258150e+00    1.732049582942250e+00    1.050956877080012e-09    6.122820268572582e-07
+[INFO][rank 0][cvStep][begin-step-attempt] step = 2, t_n = 0.0001029860256095084, h = 0.0001029860256095084, q = 1
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448671247101e+00
+1.7320483583156228e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.4999746801363614, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 0.7070709734919923, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 5.129626861394919e-05
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.4999240374084951
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 3.627192645756874e-05, b-tol = 0.01, res-tol = 0.01414213562373095
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 1, h = 0.0001029860256095084, dsm = 0.2499620187042475
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499620187042475
+[DEBUG][rank 0][cvCompleteStep][return] nst = 2, nscon = 2
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 10, hprime = 0.001029860256095084, qprime = 2, qwait = 2
+ 2.059720512190168e-04    1.224744864928434e+00    1.732047133690926e+00    2.133223819811292e-09    1.224500849161814e-06
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.001029860256095084, q = 2
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448216312835e+00
+1.7320226411776802e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4064516129032258
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 33.34377627547665, b-tol = 0.02032258064516129, res-tol = 0.02874046917080807
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 47.15522062951332, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 0.02280194342842776
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 33.32282742612932
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.01612377930633392, b-tol = 0.02032258064516129, res-tol = 0.02874046917080807
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.001029860256095084, dsm = 8.19847341436515
+[DEBUG][rank 0][cvDoErrorTest][new-step-eta] eta = 0.2729219509439317
+[INFO][rank 0][cvStep][end-step-attempt] status = failed error test, dsm = 8.19847341436515, eflag = 5
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.0002810714702930874, q = 2
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448531116913e+00
+1.7320404491464272e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4217683318751525
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 2.482943969856862, b-tol = 0.02108841659375763, res-tol = 0.02982352475586586
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-solver] solver = spgmr
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] cur-iter = 0, total-iters = 0, res-norm = 3.511413036784067, status = continue
+[INFO][rank 0][SUNLinSolSolve_SPGMR][begin-linear-iterate] 
+[INFO][rank 0][SUNLinSolSolve_SPGMR][linear-iterate] cur-iter = 1, total-iters = 1, res-norm = 0.0004635068923188323
+[INFO][rank 0][SUNLinSolSolve_SPGMR][end-linear-iterate] status = success
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 1, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 2.482495053589632
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 1, b-norm = 0.0003277492125626657, b-tol = 0.02108841659375763, res-tol = 0.02982352475586586
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success small rhs
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.0002810714702930874, dsm = 0.5885920933306285
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.5885920933306285
+[DEBUG][rank 0][cvCompleteStep][return] nst = 3, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0002810714702930874, qprime = 2, qwait = 2
+ 4.870435215121043e-04    1.224744842672905e+00    1.732034367979995e+00    4.508398321334539e-09    2.744236933471811e-06
+------------------------------------------------------------------------------------------------------------------------------
+Current time                 = 0.0004870435215121043
+Steps                        = 3
+Error test fails             = 1
+NLS step fails               = 0
+Initial step size            = 0.0001029860256095084
+Last step size               = 0.0002810714702930874
+Current step size            = 0.0002810714702930874
+Last method order            = 2
+Current method order         = 2
+Stab. lim. order reductions  = 0
+RHS fn evals                 = 11
+NLS iters                    = 8
+NLS fails                    = 0
+NLS iters per step           = 2.666666666666667
+LS setups                    = 0
+Jac fn evals                 = 0
+LS RHS fn evals              = 4
+Prec setup evals             = 0
+Prec solves                  = 0
+LS iters                     = 4
+LS fails                     = 0
+Jac-times setups             = 0
+Jac-times evals              = 4
+LS iters per NLS iter        = 0.5
+Jac evals per NLS iter       = 0
+Prec evals per NLS iter      = 0
+Root fn evals                = 0
+End CVODE Logging test
diff --git a/test/unit_tests/logging/test_logging_cvode_lvl5_1_1.out b/test/unit_tests/logging/test_logging_cvode_lvl5_1_1.out
new file mode 100644
index 0000000000..dd4a6329f8
--- /dev/null
+++ b/test/unit_tests/logging/test_logging_cvode_lvl5_1_1.out
@@ -0,0 +1,120 @@
+Start CVODE Logging test
+Using Newton nonlinear solver
+Using dense direct linear solver
+           t                        u                        v                      u err                    v err            
+------------------------------------------------------------------------------------------------------------------------------
+ 0.000000000000000e+00    1.224744871391589e+00    1.732050807568877e+00    0.000000000000000e+00    0.000000000000000e+00
+[INFO][rank 0][cvStep][begin-step-attempt] step = 1, t_n = 0, h = 0.0001029860256095084, q = 1
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448713915889e+00
+1.7320508075688772e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.4999243864411202
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 2.207987840210415e-11
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 0, h = 0.0001029860256095084, dsm = 0.2499621932116032
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499621932116032
+[DEBUG][rank 0][cvCompleteStep][return] nst = 1, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0001029860256095084, qprime = 1, qwait = 1
+ 1.029860256095084e-04    1.224744869195325e+00    1.732049582942459e+00    1.113781733508290e-09    6.122818176912403e-07
+[INFO][rank 0][cvStep][begin-step-attempt] step = 2, t_n = 0.0001029860256095084, h = 0.0001029860256095084, q = 1
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448669990604e+00
+1.7320483583160411e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.2
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 0.499923682675401
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 1
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 1, h = 0.0001029860256095084, dsm = 0.2499618413377005
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.2499618413377005
+[DEBUG][rank 0][cvCompleteStep][return] nst = 2, nscon = 2
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 10, hprime = 0.001029860256095084, qprime = 2, qwait = 2
+ 2.059720512190168e-04    1.224744864802806e+00    1.732047133692213e+00    2.258852216385776e-09    1.224499561969239e-06
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.001029860256095084, q = 2
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448208776175e+00
+1.7320226411897481e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4064516129032258
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 33.32111802218621
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 0.001726792321093216
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.001029860256095084, dsm = 8.198477691401965
+[DEBUG][rank 0][cvDoErrorTest][new-step-eta] eta = 0.2729219034840186
+[INFO][rank 0][cvStep][end-step-attempt] status = failed error test, dsm = 8.198477691401965, eflag = 5
+[INFO][rank 0][cvStep][begin-step-attempt] step = 3, t_n = 0.0002059720512190168, h = 0.0002810714214160092, q = 2
+[DEBUG][rank 0][cvPredict][return] zn_0(:) =
+1.2247448528146598e+00
+1.7320404491518191e+00
+
+[INFO][rank 0][cvNls][begin-nonlinear-solve] tol = 0.4217683352485581
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-solver] solver = Newton
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 1, total-iters = 1, update-norm = 2.482483008202045
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = continue
+[INFO][rank 0][SUNNonlinSolSolve_Newton][begin-nonlinear-iterate] 
+[INFO][rank 0][cvLsSolve][begin-linear-solve] iterative = 0
+[INFO][rank 0][cvLsSolve][end-linear-solve] status = success, iters = 0, p-solves = 0
+[INFO][rank 0][SUNNonlinSolSolve_Newton][nonlinear-iterate] cur-iter = 2, total-iters = 2, update-norm = 1.190779226370788e-05
+[INFO][rank 0][SUNNonlinSolSolve_Newton][end-nonlinear-iterate] status = success
+[INFO][rank 0][cvNls][end-nonlinear-solve] status = success, iters = 2
+[DEBUG][rank 0][cvDoErrorTest][error-test] step = 2, h = 0.0002810714214160092, dsm = 0.5885920559914323
+[INFO][rank 0][cvStep][end-step-attempt] status = success, dsm = 0.5885920559914323
+[DEBUG][rank 0][cvCompleteStep][return] nst = 3, nscon = 1
+[DEBUG][rank 0][cvPrepareNextStep][return] eta = 1, hprime = 0.0002810714214160092, qprime = 2, qwait = 2
+ 4.870434726350260e-04    1.224744841922706e+00    1.732034367987314e+00    5.258602220337139e-09    2.744232363793842e-06
+------------------------------------------------------------------------------------------------------------------------------
+Current time                 = 0.000487043472635026
+Steps                        = 3
+Error test fails             = 1
+NLS step fails               = 0
+Initial step size            = 0.0001029860256095084
+Last step size               = 0.0002810714214160092
+Current step size            = 0.0002810714214160092
+Last method order            = 2
+Current method order         = 2
+Stab. lim. order reductions  = 0
+RHS fn evals                 = 10
+NLS iters                    = 7
+NLS fails                    = 0
+NLS iters per step           = 2.333333333333333
+LS setups                    = 3
+Jac fn evals                 = 1
+LS RHS fn evals              = 0
+Prec setup evals             = 0
+Prec solves                  = 0
+LS iters                     = 0
+LS fails                     = 0
+Jac-times setups             = 0
+Jac-times evals              = 0
+LS iters per NLS iter        = 0
+Jac evals per NLS iter       = 0.1428571428571428
+Prec evals per NLS iter      = 0
+Root fn evals                = 0
+End CVODE Logging test