diff --git a/ode/impl/KokkosODE_Newton_impl.hpp b/ode/impl/KokkosODE_Newton_impl.hpp
index 507a1b191f..6ebcf20278 100644
--- a/ode/impl/KokkosODE_Newton_impl.hpp
+++ b/ode/impl/KokkosODE_Newton_impl.hpp
@@ -32,7 +32,7 @@ namespace Impl {
 
 template <class system_type, class mat_type, class ini_vec_type, class rhs_vec_type, class update_type,
           class scale_type>
-KOKKOS_FUNCTION KokkosODE::Experimental::newton_solver_status NewtonSolve(
+[[nodiscard]] KOKKOS_FUNCTION KokkosODE::Experimental::newton_solver_status NewtonSolve(
     system_type& sys, const KokkosODE::Experimental::Newton_params& params, mat_type& J, mat_type& tmp,
     ini_vec_type& y0, rhs_vec_type& rhs, update_type& update, const scale_type& scale) {
   using newton_solver_status = KokkosODE::Experimental::newton_solver_status;